Search

We investigate the driving mechanism of Alfv\'en wave-driven stellar winds from red giant stars, Arcturus ($\alpha$ Boo; K1.5 III) and Aldebaran ($\alpha$ Tau; K5 III), with nonideal MHD simulations. Since the atmosphere is not fully ionized, upward propagating Alfv\'enic waves excited by surface convection are affected by ambipolar diffusion. Our fiducial run with the nonideal MHD effect for $\alpha$ Boo gives time-averaged mass-loss rate, $\dot{M}=3.3\times 10^{-11}M_{\odot}/$yr, which is more than one order of magnitude reduced from the result in the ideal MHD run and nicely explains the observational value. Magnetized hot bubbles with $T\gtrsim 10^6$K are occasionally present simultaneously with cool gas with $T\sim$ a few thousands K in the atmosphere because of the thermal instability triggered by radiative cooling; there coexist fully ionized plasma emitting soft X-rays and molecules absorbing/emitting infrared radiations. The inhomogeneity in the atmosphere also causes large temporal variations in the mass-loss rate within an individual magnetic flux tube. We also study the effect of magnetic field strength and metallicity, and find that the wind density, and accordingly the mass-loss rate, positively and sensitively depends on both of them through the ambipolar diffusion of Alfv\'enic waves.The nonideal MHD simulation for $\alpha$ Tau, which is slightly more evolved than $\alpha$ Boo and has weaker magnetic field, results in weaker wind with $\dot{M}=1.5\times 10^{-12}M_{\odot}/$yr with the atmospheric temperature $\lesssim 10^5$K throughout the simulation time. However, given the observations implying the presence of locally strong magnetic fields on the surface of $\alpha$~Tau, we also conduct a simulation with a field strength twice as strong. This results in $\dot{M}=2.0\times 10^{-11}M_{\odot}/$yr -- comparable to the observed value -- with transient magnetized hot bubbles., Comment: 15 pages, 10 figures, submitted to PASJ, comments welcome

Stellar spin is one of the fundamental quantities that characterize a star itself and its planetary system. Nevertheless, stellar spin-down mechanisms in protostellar and pre-main-sequence stellar phases have been a long-standing problem in the star formation theory. To realize the spin-down, previous axisymmetric models based on the conventional magnetospheric paradigm have to assume massive stellar winds or produce highly time-variable magnetospheric ejections. However, this picture has been challenged by both numerical simulations and observations. With a particular focus on the propeller regime for solar-mass stars, we propose a new picture of stellar spin-down based on our recent three-dimensional (3D) magnetohydrodynamic simulation and stellar evolution calculation. We show that failed magnetospheric winds, unique to 3D models, significantly reduce the spin-up accretion torque, which make it easier for the star to spin-down. Additionally, the amplitude of time variability associated with magnetospheric ejections is reduced by 3D effects. Our simulation demonstrates that the star spins down by generating a conical disk wind, driven by a rotating stellar magnetosphere. Our theoretical estimates, inspired by the numerical model, suggest that the conical disk wind is likely to play a crucial role in extracting stellar angular momentum during the protostellar phase. As magnetospheric accretion is expected to occur in other accreting objects such as proto-giant planets, this study will also contribute to the understanding of the angular momentum of such objects., Comment: 21 pages, 14 figures. Accepted for publication in ApJ

We present steady-state solutions for a one-dimensional, magnetically-driven accretion disk wind model based on magnetohydrodynamic equations. We assume a geometrically thin, gas-pressure-dominated accretion disk, incorporating both magnetic braking and turbulent viscosity introduced by an extended alpha-viscosity prescription. Additionally, the vertical stress parameter is assumed to scale with the disk aspect ratio. We confirm that the derived solutions result in standard disk solutions when the wind is absent. We find that the mass accretion rate decreases as the disk mass falls inward, while the mass loss rate increases with radius. The disk spectrum emitted from the magnetically-driven disk wind can be observed without interference from the wind medium because the wind is significantly optically thin. The spectral luminosity is proportional to $\nu^{1/3}$ in the intermediate, multicolor-blackbody wavebands, in the absence of wind, as predicted by standard disk theory. However, in the presence of wind, it follows a different power-law dependence on frequency over the same range. A deviation from the spectral slope of $1/3$, particularly a negative spectral slope, is a clear indicator of the presence of a magnetically driven wind. We also discuss an observational strategy to test our model with multi-wavelength observations., Comment: 28 pages, 5 figures, accepted for publication in Progress of Theoretical and Experimental Physics (PTEP)

We showcase a study on the physical properties of the Circumnuclear Disk surrounding the Sgr A* of the Galactic Center, emphasizing the role of magnetic field (B field) with 0.47 pc spatial resolution, based on the sensitive {\lambda} = 850 {\mu}m polarization data taken with the JCMT. The B field within the CND exhibits a coherent spiral pattern. Applying the model described by Wardle and Ko\ddot{o}nigl 1990 (WK model) to the observed B field pattern, it favors gas-pressure-dominant models without dismissing a gas-and-B field comparable model, leading us to estimate the B-field strength in the ionized cavity around Sgr A* as 0.24 + 0.05 mG. Analysis -0.04 based on the WK model further allows us to derive representative B-field strengths for the radial, azimuthal, and vertical components as (Br,B{\phi},Bz) = (0.4 \pm 0.1,-0.7 \pm 0.2,0.2 \pm 0.05) mG, respectively. A key finding is that the |{\phi}| component is dominant over Br and Bz components, consistent with the spiral morphology, indicating that the CND' s B-field is predominantly toroidal, possibly shaped by accretion dynamics. Considering the turbulent pressure, estimated plasma \{beta} values indicate the effective gas pressure should surpass the magnetic pressure. Assessing the CND of our MWG in the toroidal-and-vertical stability parameter space, we propose that such an "effective" magnetoro-tational instability (MRI) may likely be active. The estimated maximum unstable wavelength, {\lambda}max = 0.1 \pm 0.1 pc, is smaller than the CND' s scale height (0.2 \pm 0.1 pc), which indicates the potential for the effective MRI intermittent cycles of \sim 10^{6} years, which should profoundly affect the CND's evolution, considering the estimated mass accretion rate of 10^{-2}M_{\odot} yr^{-1} to the Sgr A*.

We investigate non-ideal magnetohydrodynamical (MHD) effects in the chromosphere on the solar wind by performing MHD simulations for Alfv\'en-wave driven winds with explicitly including Ohmic and ambipolar diffusion. We find that MHD waves are significantly damped in the chromosphere by ambipolar diffusion so that the Alfv\'enic Poynting flux that reaches the corona is substantially reduced. As a result, the coronal temperature and the mass loss rate of the solar wind are considerably reduced, compared with those obtained from an ideal MHD case, which is indicative of a great importance of the non-ideal MHD effects in the solar atmosphere. However, the temperature and the mass loss rate are recovered by a small increase in the convection-originated velocity perturbation at the photosphere because of the sensitive dependence of the ambipolar diffusion and reflection of Alfv\'en waves on the physical properties of the chromosphere. We also find that density perturbations in the corona are reduced by the ambipolar diffusion of Alfv\'en waves in the chromosphere because the nonlinear generation of compressible perturbations is suppressed., Comment: 17 pages, 11 figures, 1 table, submitted to ApJ

We propose a novel method to constrain the tidal quality factor, $Q'$, from an observed non-synchronized star-planet system consisting of a slowly rotating low-mass star and a close-in Jovian planet, taking into account the co-evolution of stellar spin and planetary orbit due to the tidal interaction and the magnetic braking. On the basis of dynamical system theory, the track of the co-evolution of angular momentum from the fast rotator regime for such a system exhibits the existence of a forbidden region in the $\Omega_\mathrm{orb}$ -- $\Omega_\mathrm{spin}$ plane , where $\Omega_\mathrm{spin}$ and $\Omega_\mathrm{orb}$ denote the angular velocity of the stellar spin and planetary orbit, respectively. The forbidden region is determined primarily by the strength of the tidal interaction. By comparing ($\Omega_\mathrm{orb},\Omega_\mathrm{spin}$) of a single star-planet system to the forbidden region, we can constrain the tidal quality factor regardless of the evolutionary history of the system. The application of this method to the star-planet system, NGTS-10 -- NGTS-10 b, gives $Q' \gtrsim 10^8$, leading to an tight upper bound on the tidal torque. Since this cannot be explained by previous theoretical predictions for non-synchronized star-planet systems, our result requires mechanisms that suppress the tidal interaction in such systems., Comment: 21 pages, 10 figures, 5 tables, accepted to ApJ

We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Els\"asser numbers of the Ohmic resistivity and ambipolar diffusion is divided into two regions: "the transition zone" and "the coherent zone". The coherent zone has the same properties as the conventional dead zone, and extends outside of the transition zone in the radial direction. Between the active and coherent zones, we discover the transition zone, the inner edge of which is identical to that of the conventional dead zone. The transition zone extends out over the regions where thermal ionization determines diffusion coefficients. The transition zone has completely different physical properties than the conventional dead zone, the so-called undead zone, and the zombie zone. The combination of amplification of the radial magnetic field owing to the ambipolar diffusion and a steep radial gradient of the Ohmic diffusivity causes the efficient evacuation of the net vertical magnetic flux from the transition zone within several rotations. Surface gas accretion occurs in the coherent zone but not in the transition zone. The presence of the transition zone prohibits mass and magnetic flux transport from the coherent zone to the active zone. Mass accumulation occurs at both edges of the transition zone as a result of mass supply from the active and coherent zones., Comment: 43 pages, 37 figures, PASJ open access (https://doi.org/10.1093/pasj/psae036)

We present a time-dependent, one-dimensional, magnetically-driven disk wind model based on magnetohydrodynamic (MHD) equations, in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin and gas-pressure dominated, and explicitly accounts for magnetic braking and turbulent viscosity through an extended alpha-viscosity prescription. We find a particular wind solution for a set of basic equations that satisfies the necessary and sufficient conditions for vertically unbound MHD flows. The solution shows that the disk evolves with mass loss due to wind and accretion from the initial Gaussian density distribution. We confirm that the mass accretion rate follows the power law of time $t^{-19/16}$ at late times in the absence of wind, which matches the classical solution of Cannizzo et al. (1990). We find that the mass accretion rate is steeper than the $t^{-19/16}$ curve when the wind is present. Mass accretion is also induced by magnetic braking, known as the wind-driven accretion mechanism, which results in a faster decay with time of both the mass accretion and loss rates. In the disk emission, the ultraviolet (UV) luminosity is the highest among the optical, UV, and X-ray luminosities. While the optical and X-ray emission is observationally insignificant without magnetic braking, the X-ray emission is brighter at late times, especially in the presence of magnetic braking. This provides a possible explanation for observed delayed X-ray flares. Our model predicts that late-time bolometric light curves steeper than $t^{-19/16}$ in UV-bright TDEs are potentially compelling indicators of magnetically driven winds., Comment: 27 pages, 17 figures, Accepted for publication in The Astrophysical Journal

We investigate the role of magnetic field on the gas dynamics in a galactic bulge region by three dimensional simulations with radiative cooling and heating. While high-temperature corona with $T>10^6\ {\rm K}$ is formed in the halo regions, the temperature near the mid-plane is $\lesssim 10^4\ {\rm K}$ following the thermal equilibrium curve determined by the radiative cooling and heating. Although the thermal energy of the interstellar gas is lost by radiative cooling, the saturation level of the magnetic field strength does not significantly depend on the radiative cooling and heating. The magnetic field strength is amplified to $10\ {\rm \mu G}$ on average, and reaches several hundred ${\rm \mu G}$ locally. We find the formation of magnetically dominated regions at mid-latitudes in the case with the radiative cooling and heating, which is not seen in the case without radiative effect. The vertical thickness of the mid-latitude regions is $50-150\ {\rm pc}$ at the radial location of $0.4-0.8\ {\rm kpc}$ from the galactic center, which is comparable to the observed vertical distribution of neutral atomic gas. When we take the average of different components of energy density integrated over the galactic bulge region, the magnetic energy is comparable to the thermal energy. We conclude that the magnetic field plays a substantial role in controlling the dynamical and thermal properties of the galactic bulge region., Comment: Accepted to ApJ; 24 pages, 20 figures 3 tables

By performing ideal magnetohydrodynamical (MHD) simulations with weak vertical magnetic fields in unstratified cylindrical shearing boxes with modified boundary treatment, we investigate MHD turbulence excited by magnetorotational instability. The cylindrical simulation exhibits extremely large temporal variation in the magnetic activity compared to the simulation in a normal Cartesian shearing box, although the time-averaged field strengths are comparable in the cylindrical and Cartesian setups. Detailed analysis of the terms describing magnetic-energy evolution with "triangle diagrams" surprisingly reveals that in the cylindrical simulation the compression of toroidal magnetic field is unexpectedly as important as the winding due to differential rotation in amplifying magnetic fields and triggering intermittent magnetic bursts, which are not seen in the Cartesian simulation. The importance of the compressible amplification is also true for a cylindrical simulation with tiny curvature; the evolution of magnetic fields in the nearly Cartesian shearing box simulation is fundamentally different from that in the exact Cartesian counterpart. The radial gradient of epicycle frequency, $\kappa$, which cannot be considered in the normal Cartesian shearing box model, is the cause of this fundamental difference. An additional consequence of the spatial variation of $\kappa$ is continuous and ubiquitous formation of narrow high(low)-density and weak(strong)-field localized structures; seeds of these ring-gap structures are created by the compressible effect and subsequently amplified and maintained under the marginally unstable condition regarding "viscous-type" instability., Comment: 26 pages, 21 figures, ApJ, in press

Stellar metallicity is a critical factor to characterize the stellar coronae because it directly affects the radiative energy loss from the atmosphere. By extending theoretical relations for solar coronal loops introduced by \cite{Rosner1978}, we analytically derive scaling relations for stellar coronal loops with various metallicities. In order to validate the derived relations, we also perform magnetohydrodyamic simulations for the heating of coronal loops with different metallicities by changing radiative loss functions according to the adopted elemental abundances. The simulation results nicely explain the generalized analytical scaling relations and show a strong dependence of the thermodynamical and radiative properties of the loops on metallicity. Higher density and temperature are obtained in lower-metallicity coronae because of the inefficient radiative cooling, provided that the surface condition is unchanged. Thus, it is estimated that the X-ray radiation from metal-poor coronae is higher because of their denser coronal gas. The generalized scaling laws can also be used as a tool to study the condition of high-energy radiation around magnetically active stars and their impact on planetary environments., Comment: 14 pages, 12 figures, accepted to ApJ

By performing $N$-body simulations, we investigated fundamental processes of collisions between dust aggregates composed of submicron-sized icy dust monomers. We examined the mass distribution of fragments in the collisional outcomes in a wide range of the mass ratio and the collision velocity between colliding dust aggregates. We derived analytic expressions of the mass distribution of large remnants and small fragments by numerical fitting to the simulation results. Our analytic formulae for masses of the large remnants can reproduce the contribution of mass transfer from a large target to a small projectile, which occurs for a mass ratio of $\gtrsim 3$ and is shown in a previous study (Hasegawa et al. 2021). We found that the power-law index of the cumulative mass distribution of the small fragments is independent of the mass ratio and only weakly dependent on the collision velocity. On the other hand, the mass fraction of fragments of individual dust monomers decreases with an increasing total mass of colliding aggregates for a fixed mass ratio. This tendency implies that multiple hierarchical disruptive collisions (i.e., collisions between fragments, collisions between fragments of fragments) are required for producing a large amount of individual dust monomers via collisional fragmentation. Our fragment model suggests that the total geometric cross section integrated over the fragments is estimated to be about the same order of the geometric cross section of the target., Comment: 30 pages, 23 figures, 4 tables, accepted for publication in ApJ

We investigate the role of latitudinal differential rotation (DR) in the spin evolution of solar-type stars. Recent asteroseismic observation detected the strong equator-fast DR in some solar-type stars. Numerical simulations show that the strong equator-fast DR is a typical feature of young fast-rotating stars and that this tendency is gradually reduced with stellar age. Incorporating these properties, we develop a model for the long-term evolution of stellar rotation. The magnetic braking is assumed to be regulated dominantly by the rotation rate in the low-latitude region. Therefore, in our model, stars with the equator-fast DR spin down more efficiently than those with the rigid-body rotation. We calculate the evolution of stellar rotation in ranges of stellar mass, $0.9 \, \mathrm{M}_{\odot} \le M \le 1.2\, \mathrm{M}_{\odot}$, and metallicity, $0.5\, \mathrm{Z}_{\odot} \le Z \le 2\, \mathrm{Z}_{\odot}$, where $\mathrm{M}_{\odot}$ and $\mathrm{Z}_{\odot}$ are the solar mass and metallicity, respectively. Our model, using the observed torque in the present solar wind, nicely explains both the current solar rotation and the average trend of the rotation of solar-type stars, including the dependence on metallicity. In addition, our model naturally reproduces the observed trend of the weakened magnetic braking in old slowly rotating solar-type stars because strong equator-fast DR becomes reduced. Our results indicate that LDR and its transition are essential factors that control the stellar spin down., Comment: 20 pages, 16 figures, 3 table, accepted to MNRAS

We study the role the the p-mode-like vertical oscillation on the photosphere in driving solar winds in the framework of Alfven-wave-driven winds. By performing one-dimensional magnetohydrodynamical numerical simulations from the photosphere to the interplanetary space, we discover that the mass-loss rate is raised up to 4 times as the amplitude of longitudinal perturbations at the photosphere increases. When the longitudinal fluctuation is added, transverse waves are generated by the mode conversion from longitudinal waves in the chromosphere, which increases Alfvenic Poynting flux in the corona. As a result, the coronal heating is enhanced to yield higher coronal density by the chromospheric evaporation, leading to the increase of the mass-loss rate. Our findings clearly show the importance of the p-mode oscillation in the photosphere and the mode conversion in the chromosphere in determining the basic properties of the wind from the sun and solar-type stars., Comment: Appear in Proceedings of IAU 370 "Winds of Stars and Exoplanets" A. A. Vidotto, L. Fossati & J. Vink, eds

We perform three-dimensional magnetohydrodynamic simulations of magnetospheric accretion in a T Tauri star to study the accretion and wind structures in the close vicinity of the star. The gas accreting onto the star consists of the gas from the magnetospheric boundary and the failed disk winds. The accreting gas is commonly found as a multi-column accretion, which is consistent with observations. A significant fraction of the angular momentum of the accreting flows is removed by the magnetic fields of conical disk winds and turbulent failed winds inside and near the magnetosphere. As a result, the accretion torque is significantly reduced compared to the simple estimation based on the mass accretion rate. The stellar spin affects the time variability of the conical disk wind by changing the stability condition of the magnetospheric boundary. However, the time-averaged magnetospheric radius only weakly depends on the stellar spin, which is unlike the prediction of classical theories that the stellar spin controls the magnetospheric radius through the magnetic torque. The ratio of the toroidal to the poloidal field strengths at the magnetospheric boundary, which is a key parameter for the magnetic torque, is also insensitive to the spin; it is rather determined by the disk dynamics. Considering newly found three-dimensional effects, we obtain a scaling relation of the magnetospheric radius very similar to the Ghosh & Lamb relation from the steady angular momentum transport equation., Comment: 38 pages, 21 figures. Accepted for publication in ApJ

Recent observational and numerical studies show a variety of thermal structures in the solar chromosphere. Given that the thermal interplay across the transition region is a key to coronal heating, it is worth investigating how different thermal structures of the chromosphere yield different coronal properties. In this work, by MHD simulations of Alfv\'{e}n-wave heating of coronal loops, we study how the coronal properties are affected by the chromospheric temperature. To this end, instead of solving the radiative transfer equation, we employ a simple radiative loss function so that the chromospheric temperature is easily tuned. When the chromosphere is hotter, because the chromosphere extends to a larger height, the coronal part of the magnetic loop becomes shorter, which enhances the conductive cooling. A larger loop length is therefore required to maintain the high-temperature corona against the thermal conduction. From our numerical simulations we derive a condition for the coronal formation with respect to the half loop length $l_{\rm loop}$ in a simple form: $l_{\rm loop} > a T_{\rm min} + l_{\rm th}$, where $T_{\rm min}$ is the minimum temperature in the atmosphere and parameters $a$ and $l_{\rm th}$ have negative dependencies on the coronal field strength. Our conclusion is that the chromospheric temperature has a non-negligible impact on coronal heating for loops with small length and weak coronal field. In particular, the enhanced chromospheric heating could prevent the formation of the corona., Comment: 18 pages, 19 figures, accepted for publication in ApJ

The TRAPPIST-1 system is an iconic planetary system in various aspects (e.g., habitability, resonant relation, and multiplicity) and hence has attracted considerable attention. The mass distribution of the TRAPPIST-1 planets is characterized by two features: the two inner planets are large, and the masses of the four planets in the outer orbit increase with orbital distance. The origin of these features cannot be explained by previous formation models. We investigate whether the mass distribution of the TRAPPIST-1 system can be reproduced by a planet formation model using N-body simulations. We used a gas disk evolution model around a low-mass star constructed by considering disk winds and followed the growth and orbital migration from planetary embryos with the isolation mass, which increases with orbital distance. As a result, we find that from the initial phase, planets in inner orbits undergo rapid orbital migration, and the coalescence growth near the inner disk edge is enhanced. This allows the inner planets to grow larger. Meanwhile, compared with the inner planets, planets in outer orbits migrate more slowly and do not frequently collide with neighboring planets. Therefore, the trend of increasing mass toward the outer orbit, called reversed mass ranking, is maintained. The final mass distribution approximately agrees with the two features of the mass distribution in the TRAPPIST-1 system. We discover that the mass distribution in the TRAPPIST-1 system can be reproduced when embryos experience rapid migration and become trapped near the disk inner edge, and then more massive embryos undergo slower migration. This migration transition can be achieved naturally in a disk evolution model with disk winds., Comment: 7 pages, 5 figures, 1 table, accepted for publication in A&A, a 5-minute video summary can be found here: https://exoplanet-talks.org/talk/401

We systematically investigated the heating of coronal loops on metal-free stars with various stellar masses and magnetic fields by magnetohydrodynamic simulations. It is found that the coronal property is dependent on the coronal magnetic field strength $B_{\rm c}$ because it affects the difference of the nonlinearity of the Alfv\'{e}nic waves. Weaker $B_{\rm c}$ leads to cooler and less dense coronae because most of the input waves dissipate in the lower atmosphere on account of the larger nonlinearity. Accordingly EUV and X-ray luminosities also correlate with $B_{\rm c}$, while they are emitted in a wide range of the field strength. Finally we extend our results to evaluating the contribution from low-mass Population III coronae to the cosmic reionization. Within the limited range of our parameters on magnetic fields and loop lengths, the EUV and X-ray radiations give a weak impact on the ionization and heating of the gas at high redshifts. However, there still remains a possibility of the contribution to the reionization from energetic flares involving long magnetic loops., Comment: 11 pages, 5 figures. Accepted for publication in MNRAS

We investigated fundamental processes of collisional sticking and fragmentation of dust aggregates by carrying out N-body simulations of submicron-sized icy dust monomers. We examined the condition for collisional growth of two colliding dust aggregates in a wide range of the mass ratio, 1-64. We found that the mass transfer from a larger dust aggregate to a smaller one is a dominant process in collisions with a mass ratio of 2-30 and impact velocity of \approx 30-170 m s^-1. As a result, the critical velocity, v_fra, for fragmentation of the largest body is considerably reduced for such unequal-mass collisions; v_fra of collisions with a mass ratio of 3 is about half of that obtained from equal-mass collisions. The impact velocity is generally higher for collisions between dust aggregates with higher mass ratios because of the difference between the radial drift velocities in the typical condition of protoplanetary disks. Therefore, the reduced v_fra for unequal-mass collisions would delay growth of dust grains in the inner region of protoplanetary disks., Comment: 33 pages, 17 figures, 3 tables, accepted for publication in ApJ

We aim to understand the effect of stellar evolution on the evolution of protoplanetary disks. We focus in particular on the disk evolution around intermediate-mass (IM) stars, which evolve more rapidly than low-mass ones. We numerically solve the long-term evolution of disks around 0.5-5 solar-mass stars considering viscous accretion and photoevaporation (PE) driven by stellar far-ultraviolet (FUV), extreme-ultraviolet (EUV), and X-ray emission. We also take stellar evolution into account and consider the time evolution of the PE rate. We find that the FUV, EUV, and X-ray luminosities of IM stars evolve by orders of magnitude within a few Myr along with the time evolution of stellar structure, stellar effective temperature, or accretion rate. Therefore, the PE rate also evolves with time by orders of magnitude, and we conclude that stellar evolution is crucial for the disk evolution around IM stars., Comment: 21 pages, 15 figures, published in ApJ. Evolutionary models of young 0.5-5 Msun stars are available at https://zenodo.org/record/4595807

Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfv\'en-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $\tau_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( \Omega_\ast / \Omega_\odot \right)^{2.82}$, yielding a spin-down law $\Omega_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfv\'en waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto \Omega_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfv\'en radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfv\'en radius matches previous works which used thermally-driven winds. Our results suggest that the Alfv\'en-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence., Comment: accepted for publication in The Astrophysical Journal

We discovered a new growth mode of dust grains to km-sized bodies in protoplanetary disks that evolve by viscous accretion and magnetically driven disk winds (MDWs). We solved an approximate coagulation equation of dust grains with time-evolving disks that consist of both gas and solid components by a one-dimensional model. With the grain growth, all solid particles initially drift inward toward the central star by the gas drag force. However, the radial profile of gas pressure, $P$, is modified by the MDW that disperses the gas in an inside-out manner. Consequently, a local concentration of solid particles is created by the converging radial flux of drifting dust grains at the location with the convex upward profile of $P$. When the dimensionless stopping time, ${\rm St}$, there exceeds unity, the solid particles spontaneously reach the growth dominated state because of the positive feedback between the suppressed radial drift and the enhanced accumulation of dust particles that drift from the outer part. Once the solid particles are in the drift limited state, the above-mentioned condition of ${\rm St} \gtrsim 1$ for the dust growth is equivalent with \begin{equation} \Sigma_{\rm d}/\Sigma_{\rm g}\gtrsim \eta, \nonumber \end{equation} where $\Sigma_{\rm d}/\Sigma_{\rm g}$ is the dust-to-gas surface-density ratio and $\eta$ is dimensionless radial pressure-gradient force. As a consequence of the successful growth of dust grains, a ring-like structure containing planetesimal-sized bodies is formed at the inner part of the protoplanetary disks. Such a ring-shaped concentration of planetesimals is expected to play a vital role in the subsequent planet formation., Comment: 24 pages, 10 figures, 1 table, accepted for publication in The Astrophysical Journal

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

We develop a framework for magnetohydrodynamical (MHD) simulations in a local cylindrical shearing box by extending the formulation of the Cartesian shearing box. We construct shearing-periodic conditions at the radial boundaries of a simulation box from the conservation relations of the basic MHD equations, taking into account the explicit radial dependence of physical quantities. We demonstrate quasi-steady mass accretion, which cannot be handled by the standard Cartesian shearing box model, with an ideal MHD simulation in a vertically unstratified cylindrical shearing box up to 200 rotations. In this demonstrative run we set up (i) net vertical magnetic flux, (ii) a locally isothermal equation of state, and (iii) a sub-Keplerian equilibrium rotation, whereas the sound velocity and the initial Alfven velocity have the same radial dependence as that of the Keplerian velocity. Inward mass accretion is induced to balance with the outward angular momentum flux of the MHD turbulence triggered by the magnetorotational instability in a self-consistent manner. We discuss detailed physical properties of the saturated magnetic field, in comparison to the results of a Cartesian shearing box simulation., Comment: 27 pages, 10 figures, accepted by PASJ, movie is available at http://ea.c.u-tokyo.ac.jp/astro/Members/stakeru/research/cylshbx/

Recent theoretical studies suggest the existence of low-mass zero-metal stars in the current universe. In order to study the basic properties of the atmosphere of low-mass first stars, we performed one dimensional magnetohydrodynamical simulations for the heating of coronal loops on low-mass stars with various metallicities. While the simulated loops are heated up to >$10^6$ K by the dissipation of Alfvenic waves originating from the convective motion irrespectively of the metallicity, the coronal properties sensitively depend on the metallicity. Lower-metal stars create hotter and denser coronae because the radiative cooling is suppressed. The zero-metal star gives more than ten times higher coronal density than the solar-metallicity counterpart, and as a result, the UV and X-ray fluxes from the loop are (1-5) times higher than those of the solar metallicity star. We also discuss the dependence of the coronal properties on the length of the simulated coronal loops, Comment: 11 pages, 12 figures

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

Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings and gaps form from the same basic mechanism that was identified in the 2D simulations, namely, the redistribution of the poloidal magnetic flux relative to the disk material as a result of the reconnection of a sharply pinched poloidal magnetic field lines. Thus, the less dense gaps are more strongly magnetized with a large poloidal magnetic field compared to the less magnetized (dense) rings. The rings and gaps start out rather smoothly in 3D simulations that have axisymmetric initial conditions. Non-axisymmetric variations arise spontaneously at later times, but they do not grow to such an extent as to disrupt the rings and gaps. These disk substructures persist to the end of the simulations, lasting up to 3000 orbital periods at the inner edge of the simulated disk. The longevity of the perturbed yet still coherent rings make them attractive sites for trapping large grains that would otherwise be lost to rapid radial migration due to gas drag. As the ambipolar diffusivity decreases, both the disk and the wind become increasingly turbulent, driven by the magnetorotational instability, with tightly-wound spiral arms becoming more prominent in the disk., Comment: 19 pages, 13 figures, accepted by MNRAS

The recent discovery of hot dust grains in the vicinity of main-sequence stars has become a hot issue among the scientific community of debris disks. Hot grains must have been enormously accumulated near their sublimation zones, but it is a mystery how such a high concentration of hot grains is sustained. The most difficult conundrum is that the size of hot dust grains is estimated to lie in the submicrometer range, while submicrometer-sized grains are instantly swept away from near-stellar environments by stellar radiation pressure. One and only mechanism proposed for prolonging the residence time of hot grains in the near-stellar environments is trapping of charged nanoparticles by stellar magnetic fields. We revisit the model of magnetic grain trapping around main-sequence stars of various spectral classes by taking into account sublimation and electric charging of the grains. The model of magnetic grain trapping predicts that hot dust grains are present in the vicinity of main-sequence stars with high rotation velocities and intermediate magnetic-field strengths. On the contrary, we find that the detection of hot dust grains has no correlation with the rotation velocities of central stars nor the magnetic field strengths of the stars. Our numerical evaluation of electric grain charging indicates that the surface potential of submicrometer-sized grains in the vicinity of main-sequence stars is typically 4-5 V, which is one order of magnitude smaller than the value assumed by the model of magnetic grain trapping. On the basis of our numerical simulation on sublimation of dust grains in the vicinity of a star, it turns out that their lives end due to sublimation in a timescale much shorter than the period of one revolution at the gyroradius. It is, therefore, infeasible to dynamically extend the dwell time of hot grains inside the sublimation zone by magnetic trapping, ..., Comment: 30 pages, 8 figures, 3 tables, to appear in Planetary and Space Science

We calculate accretion mass of interstellar objects (ISOs) like `Oumuamua onto low-mass population III stars (Pop.~III survivors), and estimate surface pollution of Pop.~III survivors. An ISO number density estimated from the discovery of `Oumuamua is so high ($\sim 0.2$~au$^{-3}$) that Pop.~III survivors have chances at colliding with ISOs $\gtrsim 10^5$ times per $1$~Gyr. `Oumuamua itself would be sublimated near Pop.~III survivors, since it has small size, $\sim 100$~m. However, ISOs with size $\gtrsim 3$~km would reach the Pop.~III survivor surfaces. Supposing an ISO cumulative number density with size larger than $D$ is $n \propto D^{-\alpha}$, Pop.~III survivors can accrete ISO mass $\gtrsim 10^{-16}M_\odot$, or ISO iron mass $\gtrsim 10^{-17}M_\odot$, if $\alpha < 4$. This iron mass is larger than the accretion mass of interstellar medium (ISM) by several orders of magnitude. Taking into account material mixing in a convection zone of Pop.~III survivors, we obtain their surface pollution is typically [Fe/H] $\lesssim -8$ in most cases, however the surface pollution of Pop.~III survivors with $0.8M_\odot$ can be [Fe/H] $\gtrsim -6$ because of the very shallow convective layer. If we apply to Pop.III survivors located at the Galactocentric distance of 8 kpc, the dependence of the metal pollustion is as follows. If $\alpha > 4$, Pop.~III survivors have no chance at colliding with ISOs with $D \gtrsim 3$~km, and keep metal-free. If $3 < \alpha < 4$, Pop.~III survivors would be most polluted by ISOs up to [Fe/H] $\sim -7$. If $\alpha < 3$ up to $D \sim 10$~km, Pop.~III survivors could hide in metal-poor stars so far discovered. Pop.~III survivors would be more polluted with decreasing the Galactocentric distance. Although the metal pollution depends on $\alpha$ and the Galactocentric distance, we first show the importance of ISOs for the metal pollution of Pop.~III survivors., Comment: 6 pages, 3 figures, accepted for publication in PASJ

No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study of disk evolution that includes magnetically driven disk winds showed that the gas disk obtains a positive surface density slope inside ~ 1 au from the central star. In a region with positive midplane pressure gradient, planetesimals undergo outward radial drift. We investigate the radial drift of planetesimals and type I migration of planetary embryos in a disk that viscously evolves with magnetically driven disk winds. We show a case in which no planets remain in the close-in region. Radial drifts of planetesimals are simulated using a recent disk evolution model that includes effects of disk winds. The late stage of planet formation is also examined by performing N-body simulations of planetary embryos. We demonstrate that in the middle stage of disk evolution, planetesimals can undergo convergent radial drift in a magnetorotational instability (MRI)-inactive disk, in which the pressure maximum is created, and accumulate in a narrow ring-like region with an inner edge at ~ 0.7 au from the Sun. We also show that planetary embryos that may grow from the narrow planetesimal ring do not exhibit significant type I migration in the late stage of disk evolution. The origin of the localized configuration of the terrestrial planets of the solar system, in particular the deficit of close-in planets, can be explained by the convergent radial drift of planetesimals in disks with a positive pressure gradient in the close-in region., Comment: 5 pages, 4 figures, accepted for publication in A&A Letters

Planets with masses larger than about 0.1 Earth-masses undergo rapid inward migration (type I migration) in a standard protoplanetary disk. Recent magnetohydrodynamical simulations revealed the presence of magnetically driven disk winds, which would alter the disk profile and the type I migration in the close-in region. We investigate orbital evolution of planetary embryos in disks that viscously evolve under the effects of disk winds. The aim is to discuss effects of altered disk profiles on type I migration. In addition, we aim to examine whether observed distributions of close-in super-Earths can be reproduced by simulations that include effects of disk winds. We perform N-body simulations of super-Earth formation from planetary embryos, in which a recent model for disk evolution is used. We explore a wide range of parameters and draw general trends. We also carry out N-body simulations of close-in super-Earth formation from embryos in such disks under various conditions. We find that the type I migration is significantly suppressed in many cases. Even in cases in which inward migration occurs, the migration timescale is lengthened to 1 Myr, which mitigates the type I migration problem. This is because the gas surface density is decreased and has a flatter profile in the close-in region due to disk winds. We find that when the type I migration is significantly suppressed, planets undergo late orbital instability during the gas depletion, leading to a non-resonant configuration. We also find that observed distributions of close-in super-Earths (e.g., period ratio, mass ratio) can be reproduced. In addition, we show that in some results of simulations, systems with a chain of resonant planets, like the TRAPPIST-1 system, form., Comment: 18 pages, 19 figures, accepted for publication in A&A

Using numerical simulations, we investigate the onset and suppression of parametric decay instability (PDI) in the solar wind, focusing on the suppression effect by the wind acceleration and expansion. We also discuss the role of PDI in generating the density fluctuation in the solar wind and the cross-helicity evolution up to $1 {\rm \ au}$. Numerical simulations from the coronal base to 1 au are conducted in a self-consistent manner; we take into account the feedback of wave energy and pressure in the background. Monochromatic waves with various injection frequencies $f_0$ are injected to discuss the suppression of PDI, while broadband waves are applied to compare the numerical results with observation. We find, via the simulations with monochromatic waves, that high-frequency ($f_0 \gtrsim10^{-3} {\rm \ Hz}$) Alfv\'en waves are subject to PDI. Meanwhile, the maximum growth rate of the PDI of low-frequency ($f_0 \lesssim10^{-4} {\rm \ Hz}$) Alfv\'en waves becomes negative due to acceleration and expansion effects. Medium-frequency ($f_0 \approx 10^{-3.5} {\rm \ Hz}$) Alfv\'en waves have a positive growth rate but do not show the signature of PDI up to 1 au because the growth rate is too small. The medium-frequency waves experience neither PDI nor reflection so they propagate through the solar wind most efficiently. This is called the frequency-filtering mechanism of the solar wind. The simulations with broadband waves indicate that the observed trend of the density fluctuation is well explained by the evolution of PDI while the observed cross-helicity evolution is in agreement with low-frequency wave propagation., Comment: Accepted for publication in The Astrophysical Journal

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

We propose a novel one-dimensional model that includes both shock and turbulence heating and qualify how these processes contribute to heating the corona and driving the solar wind. Compressible MHD simulations allow us to automatically consider shock formation and dissipation, while turbulent dissipation is modeled via a one-point closure based on Alfv\'en wave turbulence. Numerical simulations were conducted with different photospheric perpendicular correlation lengths $\lambda_0$, which is a critical parameter of Alfv\'en wave turbulence, and different root-mean-square photospheric transverse-wave amplitudes $\delta v_0$. For the various $\lambda_0$, we obtain a low-temperature chromosphere, high-temperature corona, and supersonic solar wind. Our analysis shows that turbulence heating is always dominant when $\lambda_0 \lesssim 1{\rm \ Mm}$. This result does not mean that we can ignore the compressibility because the analysis indicates that the compressible waves and their associated density fluctuations enhance the Alfv\'en wave reflection and therefore the turbulence heating. The density fluctuation and the cross helicity are strongly affected by $\lambda_0$, while the coronal temperature and mass loss rate depend weakly on $\lambda_0$., Comment: Accepted for publicaiton in The Astrophysical Journal

We studied roles of the magnetic field on the gas dynamics in the Galactic bulge by a three-dimensional global magnetohydrodynamical simulation data, particularly focusing on vertical flows that are ubiquitously excited by magnetic activity. In local regions where the magnetic filed is stronger, it is frequently seen that fast down-flows slide along inclined magnetic field lines that are associated with buoyantly rising magnetic loops. The vertical velocity of these down-flows reaches ~ 100 km s$^{-1}$ near the foot-point of the loops by the gravitational acceleration toward the Galactic plane. The two footpoints of rising magnetic loops are generally located at different radial locations and the field lines are deformed by the differential rotation. The angular momentum is transported along the field lines, and the radial force balance breaks down. As a result, a fast downflow is often observed only at the one footpoint located at the inner radial position. The fast downflow compresses the gas to form a dense region near the footpoint, which will be important in star formation afterward. Furthermore, the horizontal components of the velocity are also fast near the foot-point because the down-flow is accelerated along the magnetic sliding slope. As a result, the high-velocity flow creates various characteristic features in a simulated position-velocity diagram, depending on the viewing angle., Comment: Accepted to MNRAS. 11 pages, 13 figures

We investigated stellar winds from zero/low-metallicity low-mass stars by magnetohydrodynamical simulations for stellar winds driven by \Alfven waves from stars with mass $M_{\star}=(0.6-0.8)M_{\odot}$ and metallicity $Z=(0-1)Z_{\odot}$, where $M_{\odot}$ and $Z_{\odot}$ are the solar mass and metallicity, respectively. \Alfvenic waves, which are excited by the surface convection, travel upward from the photosphere and heat up the corona by their dissipation. For lower $Z$, denser gas can be heated up to the coronal temperature because of the inefficient radiation cooling. The coronal density of Pop.II/III stars with $Z\le 0.01Z_{\odot}$ is 1-2 orders of magnitude larger than that of the solar-metallicity star with the same mass, and as a result, the mass loss rate, $\dot{M}$, is $(4.5-20)$ times larger. This indicates that metal accretion on low-mass Pop.III stars is negligible. The soft X-ray flux of the Pop.II/III stars is also expected to be $\approx (1-30)$ times larger than that of the solar-metallicity counterpart owing to the larger coronal density, even though the radiation cooling efficiency is smaller. A larger fraction of the input \Alfvenic wave energy is transmitted to the corona in low $Z$ stars because they avoid severe reflection owing to the smaller density difference between the photosphere and the corona. Therefore, a larger fraction is converted to the thermal energy of the corona and the kinetic energy of the stellar wind. From this energetics argument, we finally derived a scaling of $\dot{M}$ as $\dot{M}\propto L R_{\star}^{11/9}M_{\star}^{-10/9}T_{\rm eff}^{11/2}\left[\max (Z/Z_{\odot},0.01)\right]^{-1/5}$, where $L$, $R_{\star}$, and $T_{\rm eff}$ are stellar luminosity, radius, and effective temperature, respectively., Comment: 16 pages, 8 figures embedded, comments are welcome

We investigate roles of magnetic activity in the Galactic bulge region in driving large-scale outflows of size $\sim 10$ kpc. Magnetic buoyancy and breakups of channel flows formed by magnetorotational instability excite Poynting flux by the magnetic tension force. A three-dimensional global numerical simulation shows that the average luminosity of such \Alfvenic Poynting flux is $10^{40} - 10^{41}$ erg s$^{-1}$. We examine the energy and momentum transfer from the Poynting flux to the gas by solving time-dependent hydrodynamical simulations with explicitly taking into account low-frequency \Alfvenic waves of period of 0.5 Myr in a one-dimensional vertical magnetic flux tube. The \Alfvenic waves propagate upward into the Galactic halo, and they are damped through the propagation along meandering magnetic field lines. If the turbulence is nearly trans-Alfv\'{e}nic, the wave damping is significant, which leads to the formation of an upward propagating shock wave. At the shock front, the temperature $\gtrsim 5\times 10^6$ K, the density $\approx 6\times 10^{-4}$ cm$^{-3}$, and the outflow velocity $\approx 400-500$ km s$^{-1}$ at a height $\approx 10$ kpc, which reasonably explain the basic physical properties of the thermal component of the Fermi bubbles., Comment: 14 pages, 10 figures included; comments are welcome

Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We aim to investigate effects of global gas flows on type I migration and also examine planet formation. A simplified approach was taken to address this issue, and N-body simulations with isolation-mass planets were also performed. In previous studies, the effect of gas flow induced by turbulence-driven accretion has been taken into account for its desaturation effect of the corotation torque. If more rapid gas flows (e.g., wind-driven accretion) are considered, the desaturation effect can be modified. In MRI-inactive disks, in which the wind-driven accretion dominates the disk evolution, the gas flow at the midplane plays an important role. If this flow is fast, the corotation torque is efficiently desaturated. Then, the fact that the surface density slope can be positive in the inner region due to the wind torque can generate an outward migration region extended to super-Earth mass planets. In this case, we observe that no planets fall onto the central star in N-body simulations with migration forces imposed to reproduce such migration pattern. We also see that super-Earth mass planets can undergo outward migration. Relatively rapid gas flows affects type I migration and thus the formation of close-in planets., Comment: 9 pages, 9 figures, accepted for publication in A&A

There is a subclass of the X-ray jets from young stellar objects which are heated very close to the footpoint of the jets, particularly DG Tau jets. Previous models attribute the strong heating to shocks in the jets. However, the mechanism that localizes the heating at the footpoint remains puzzling. We presented a different model of such X-ray jets, in which the disk atmosphere is magnetically heated. Our disk corona model is based on the so-called nanoflare model for the solar corona. We show that the magnetic heating near the disks can result in the formation of a hot corona with a temperature of > 10^6 K even if the average field strength in the disk is moderately weak, > 1 G. We determine the density and the temperature at the jet base by considering the energy balance between the heating and cooling. We derive the scaling relations of the mass loss rate and terminal velocity of jets. Our model is applied to the DG Tau jets. The observed temperature and estimated mass loss rate are consistent with the prediction of our model in the case of the disk magnetic field strength of ~20 G and the heating region of < 0.1 au. The derived scaling relation of the temperature of X-ray jets could be a useful tool to estimate the magnetic field strength. We also found that the jet X-ray can have a significant impact on the ionization degree near the disk surface and the dead-zone size., Comment: 31 pages (single column), 5 figures, Accepted to be published in ApJ

The vertical flow structure in the galactic center region remains poorly understood. We analyzed the MHD simulation data by Suzuki et al. (2015) for better understanding. As a result, we found the fast downflows with a speed of ~100 km/s near the foot-points of magnetic loops. These downflows are flowing along a magnetic field line and accelerated by the gravity. The direction of the fast flows is changed by the magnetic loop geometry, as it moves. As a result, not only vertical motions but also azimuthal and radial motions are excited. This feature could be relevant to the observed high velocity dispersion in the position-velocity diagram., Comment: 3 pages, 2 figures, modified proceeding based on a poster session at IAU 322 (The Multi-Messenger Astrophysics of the Galactic Centre), correct just one part of coauthor's affiliation

By performing a global magnetohydrodynamical (MHD) simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic centre (GC) region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram shape of the central molecular zone. In addition, Parker instability (magnetic buoyancy) creates vertical magnetic structure, which would correspond to observed molecular loops, and frequently excited vertical flows. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time dependent, which would contribute to the outflow from the bulge., Comment: 4 pages, based on a contributed talk at IAU 322 (The Multi-Messenger Astrophysics of the Galactic Centre)

Aims: We investigate the evolution of protoplanetary discs (PPDs hereafter) with magnetically driven disc winds and viscous heating. Methods: We consider an initially massive disc with ~0.1 Msun to track the evolution from the early stage of PPDs. We solve the time evolution of surface density and temperature by taking into account viscous heating and the loss of the mass and the angular momentum by the disc winds within the framework of a standard alpha model for accretion discs. Our model parameters, turbulent viscosity, disc wind mass loss, and disc wind torque, which are adopted from local magnetohydrodynamical simulations and constrained by the global energetics of the gravitational accretion, largely depends on the physical condition of PPDs, particularly on the evolution of the vertical magnetic flux in weakly ionized PPDs. Results: Although there are still uncertainties concerning the evolution of the vertical magnetic flux remaining, surface densities show a large variety, depending on the combination of these three parameters, some of which are very different from the surface density expected from the standard accretion. When a PPD is in a "wind-driven accretion" state with the preserved vertical magnetic field, the radial dependence of the surface density can be positive in the inner region <1-10 au. The mass accretion rates are consistent with observations, even in the very low level of magnetohydrodynamical turbulence. Such a positive radial slope of the surface density gives a great impact on planet formation because (i)it inhibits the inward drift or even results in the outward drift of pebble/boulder-sized solid bodies, and (ii) it also makes the inward type-I migration of proto-planets slower or even reversed. Conclusions: The variety of our calculated PPDs should yield a wide variety of exoplanet systems., Comment: 16 pages, 11 figures embedded, accepted by A&A (comments are welcome)

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., Comment: 11 pages, 9 figures, re-submitted to ApJ with minor revision

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., Comment: 12 pages, 14 figures, accepted for publication in ApJ

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

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., Comment: 8 pages, 4 figures, accepted for publication in A&A

By performing a global magneto-hydrodynamical simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic center region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. In addition, the magnetic pressure-gradient force also drives radial flows in a similar manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram-shape of the central molecular zone. This is a natural extension into the central few 100 pc of the magnetic activity, which is observed as molecular loops at radii from a few 100 pc to 1 kpc. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time-dependent, which we discuss from a viewpoint of the outflow from the bulge., Comment: 12 pages, 13 figures, MNRAS in press

We investigate the effect of the axion cooling on the nucleosynthesis in a massive star with $16M_{\odot}$ by standard stellar evolution calculation. We find that the axion cooling suppresses the nuclear reactions in carbon, oxygen and silicon burning phases because of the extraction of the energy. As a result, larger amounts of the already synthesized neon and magnesium remain without being consumed to produce further heavier elements. Even in the case with the axion-photon coupling constant $g_{a\gamma}= 10^{-11}$ GeV$^{-1}$, which is six times smaller than the current upper limit, the amount of neon and magnesium that remain just before the core-collapse supernova explosion is considerably larger than the standard value. This implies that we could give a more stringent constraint on $g_{a\gamma}$ from the nucleosynthesis of heavy elements in massive stars., Comment: 5 pages, 5 figures, submitted to Physical Review D