Showing 1–50 of 91 results for author: Suzuki, T K

Spiral Magnetic Field and Their Role on Accretion Dynamics in the Circumnuclear Disk of Sagittarius A*: Insight from λ = 850 μm Polarization Imaging

Authors: Kazuki Sato, Hiroko Shinnaga, Ray S. Furuya, Takeru K. Suzuki, Kensuke Kakiuchi, Jürgen Ott

Abstract: 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 λ = 850 μ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… ▽ More 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 λ = 850 μ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φ,Bz) = (0.4 \pm 0.1,-0.7 \pm 0.2,0.2 \pm 0.05) mG, respectively. A key finding is that the |φ| 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, λ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*. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Effect of Magnetic diffusion in the Chromosphere on the Solar Wind

Authors: Masato Matsuoka, Takeru K. Suzuki, Takato Tokuno, Kensuke Kakiuchi

Abstract: We investigate non-ideal magnetohydrodynamical (MHD) effects in the chromosphere on the solar wind by performing MHD simulations for Alfvén-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énic Poynting flux that reaches the corona is substantially reduced. As a resu… ▽ More We investigate non-ideal magnetohydrodynamical (MHD) effects in the chromosphere on the solar wind by performing MHD simulations for Alfvén-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énic 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én 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én waves in the chromosphere because the nonlinear generation of compressible perturbations is suppressed. △ Less

Submitted 8 May, 2024; v1 submitted 27 March, 2024; originally announced March 2024.

Comments: 17 pages, 11 figures, 1 table, submitted to ApJ

A Novel Method to Constrain Tidal Quality Factor from A Non-synchronized Exoplanetary System

Authors: Takato Tokuno, Akihiko Fukui, Takeru K. Suzuki

Abstract: We propose a novel method to constrain the tidal quality factor, $Q'$, from a non-synchronized star-planet system consisting of a slowly rotating low-mass star and a close-in Jovian planet, taking into account the tidal interaction and the magnetic braking. On the basis of dynamical system theory, the track of the co-evolution of angular momentum for such a system exhibits the existence of a forbi… ▽ More We propose a novel method to constrain the tidal quality factor, $Q'$, from a non-synchronized star-planet system consisting of a slowly rotating low-mass star and a close-in Jovian planet, taking into account the tidal interaction and the magnetic braking. On the basis of dynamical system theory, the track of the co-evolution of angular momentum for such a system exhibits the existence of a forbidden region in the $Ω_\mathrm{orb}$ -- $Ω_\mathrm{spin}$ plane , where $Ω_\mathrm{spin}$ and $Ω_\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 ($Ω_\mathrm{orb},Ω_\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. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 20 pages, 10 figures, 5 tables, submitted to ApJ

Dynamics Near the Inner Dead-Zone Edges in a Proprotoplanetary Disk

Authors: Kazunari Iwasaki, Kengo Tomida, Shinsuke Takasao, Satoshi Okuzumi, Takeru K. Suzuki

Abstract: 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ässer numbers of the Ohmic r… ▽ More 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ässer 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. △ Less

Submitted 8 June, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 43 pages, 37 figures, PASJ open access (https://doi.org/10.1093/pasj/psae036)

Evolution of Tidal Disruption Event Disks with Magnetically Driven Winds

Authors: T. Mageshwaran, Kimitake Hayasaki, Takeru K. Suzuki

Abstract: We present a time-dependent, one-dimensional magnetically driven disk-wind model based on magnetohydrodynamical (MHD) equations in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin and gas-pressure dominant and explicitly take into account magnetic braking as well as turbulent viscosity by an extended alpha viscosity prescription. We find a particular win… ▽ More We present a time-dependent, one-dimensional magnetically driven disk-wind model based on magnetohydrodynamical (MHD) equations in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin and gas-pressure dominant and explicitly take into account magnetic braking as well as turbulent viscosity by 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 demonstrates that the disk evolves with losing mass 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 if the wind is absent, which corresponds to the classical solution of Cannizzo et al. (1990). We find that the mass accretion rate is steeper than the $t^{-19/16}$ curve if the disk wind is present. This is because the wind removes a significant fraction of the mass and angular momentum. Mass accretion is further induced by magnetic braking, known as a wind-driven accretion mechanism, resulting in more rapid decay with time in both the mass accretion and loss rates. The ultraviolet (UV) luminosity is the highest among the optical, UV, and X-ray luminosities from early to late evolutionary phases, suggesting optical and X-ray emissions from the disk are observationally insignificant due to magnetically driven winds in TDEs. Our model predicts that late-time bolometric light curves steeper than $t^{-19/16}$ in UV-bright TDEs potentially serve as compelling indicators for magnetically driven winds. △ Less

Submitted 24 December, 2023; originally announced December 2023.

Comments: 20 pages, 13 figures, submitted to Physical Review D

MHD Simulation in Galactic Center Region with Radiative Cooling and Heating

Authors: Kensuke Kakiuchi, Takeru. K. Suzuki, Shu-ichiro Inutsuka, Tsuyoshi Inoue, Jiro Shimoda

Abstract: 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. Althou… ▽ More 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 μG}$ on average, and reaches several hundred ${\rm μ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. △ Less

Submitted 4 April, 2024; v1 submitted 27 June, 2023; originally announced June 2023.

Comments: Accepted to ApJ; 24 pages, 20 figures 3 tables

MHD in a cylindrical shearing box II: Intermittent Bursts and Substructures in MRI Turbulence

Authors: Takeru K. Suzuki

Abstract: 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 bo… ▽ More 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, $κ$, 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 $κ$ 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. △ Less

Submitted 3 November, 2023; v1 submitted 20 May, 2023; originally announced May 2023.

Comments: 26 pages, 21 figures, ApJ, in press

Journal ref: ApJ, 957, 99, 2003

Simulations of Protoplanetary Disk Dispersal: Stellar Mass Dependence of the Disk Lifetime

Authors: A. Komaki, S. Fukuhara, T. K. Suzuki, N. Yoshida

Abstract: Recent infrared and submillimeter observations suggest that the protoplanetary disk lifetime depends on the central stellar mass. The disk dispersal is thought to be driven by viscous accretion, magneto-hydrodynamics (MHD) winds, and photoevaporation by the central star. We perform a set of one-dimensional simulations of long-term disk evolution that include all the three processes. We vary the st… ▽ More Recent infrared and submillimeter observations suggest that the protoplanetary disk lifetime depends on the central stellar mass. The disk dispersal is thought to be driven by viscous accretion, magneto-hydrodynamics (MHD) winds, and photoevaporation by the central star. We perform a set of one-dimensional simulations of long-term disk evolution that include all the three processes. We vary the stellar mass in the range of 0.5-7M$_{\odot}$, and study the mass dependence of the disk evolution. We show that a significant fraction of the disk gas is lost by MHD winds in the early stage, but the later disk evolution is mainly governed by photoevaporation. The disk radius decreases as photoevaporation clears out the gas in the outer disk efficiently. The qualitative evolutionary trends of the disk mass are remarkably similar for the wide range of the central stellar mass we consider, and the time evolution of the disk mass can be well fitted by a simple function. The dispersal time is approximately ten million years for low mass stars with weak mass dependence, but gets as short as two million years around a 7M$_{\odot}$ star. In the latter case, a prominent inner hole is formed by the combined effect of accretion and MHD winds within about one million years. The strength of the MHD wind and viscous accretion controls the overall mass-loss rate, but does not alter the dependence of the dispersal timescale on the central stellar mass. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 14 pages, 10 figures, 1 table

Coronal Loops with Different Metallicities and Generalized RTV Scaling Laws

Authors: Haruka Washinoue, Takeru K. Suzuki

Abstract: 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 m… ▽ More 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. △ Less

Submitted 25 June, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: 14 pages, 12 figures, accepted to ApJ

Collisional Growth and Fragmentation of Dust Aggregates. II. Mass Distribution of Icy Fragments

Authors: Yukihiko Hasegawa, Takeru K. Suzuki, Hidekazu Tanaka, Hiroshi Kobayashi, Koji Wada

Abstract: 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 remn… ▽ More 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. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 30 pages, 23 figures, 4 tables, accepted for publication in ApJ

Transition of latitudinal differential rotation as a possible cause of weakened magnetic braking of solar-type stars

Authors: Takato Tokuno, Takeru K. Suzuki, Munehito Shoda

Abstract: 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 prope… ▽ More 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. △ Less

Submitted 28 February, 2023; v1 submitted 24 November, 2022; originally announced November 2022.

Comments: 20 pages, 16 figures, 3 table, accepted to MNRAS

Journal ref: MNRAS, 520, 1, p.418 (2023)

Role of Longitudinal Waves in Alfven-wave-driven Solar/Stellar Wind

Authors: Kimihiko Shimizu, Munehito Shoda, Takeru K. Suzuki

Abstract: 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 inc… ▽ More 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. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Appear in Proceedings of IAU 370 "Winds of Stars and Exoplanets" A. A. Vidotto, L. Fossati & J. Vink, eds

Three-dimensional Simulations of Magnetospheric Accretion in a T Tauri Star: Accretion and Wind Structures Just Around Star

Authors: Shinsuke Takasao, Kengo Tomida, Kazunari Iwasaki, Takeru K. Suzuki

Abstract: 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.… ▽ More 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. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 38 pages, 21 figures. Accepted for publication in ApJ

The Effect of the Chromospheric Temperature on Coronal Heating

Authors: Haruka Washinoue, Munehito Shoda, Takeru K. Suzuki

Abstract: 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én-wave heating of coronal loops, we study how… ▽ More 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é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. △ Less

Submitted 4 October, 2022; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: 18 pages, 19 figures, accepted for publication in ApJ

Role of Longitudinal Waves in Alfvén-wave-driven Solar Wind

Authors: Kimihiko Shimizu, Munehito Shoda, Takeru K. Suzuki

Abstract: We revisit the role of longitudinal waves in driving the solar wind. We study how the the $p$-mode-like vertical oscillation on the photosphere affects the properties of solar winds under the framework of Alfvén-wave-driven winds. We perform a series of one-dimensional magnetohydrodynamical numerical simulations from the photosphere to beyond several tens of solar radii. We find that the mass-loss… ▽ More We revisit the role of longitudinal waves in driving the solar wind. We study how the the $p$-mode-like vertical oscillation on the photosphere affects the properties of solar winds under the framework of Alfvén-wave-driven winds. We perform a series of one-dimensional magnetohydrodynamical numerical simulations from the photosphere to beyond several tens of solar radii. We find that the mass-loss rate drastically increases with the longitudinal wave amplitude at the photosphere up to $\sim 4$ times, in contrast to the classical understanding that the acoustic wave hardly affects the energetics of the solar wind. The addition of the longitudinal fluctuation induces the longitudinal-to-transverse wave mode conversion in the chromosphere, which results in the enhanced Alfvénic Poynting flux in the corona. Consequently, the coronal heating is promoted to give higher coronal density by the chromospheric evaporation, leading to the increased mass-loss rate. This study clearly shows the importance of the longnitudinal oscillation in the photosphere and the mode conversion in the chromosphere in determining the basic properties of the wind from solar-like stars. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 22 pages, 14 figures, Resubmitted to ApJ after referee comments

Rapid-then-slow migration reproduces mass distribution of TRAPPIST-1 system

Authors: Masahiro Ogihara, Eiichiro Kokubo, Ryuunosuke Nakano, Takeru K. Suzuki

Abstract: 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 fe… ▽ More 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. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 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

Journal ref: A&A 658, A184 (2022)

Coronal Properties of Low-mass Population III Stars and the Radiative Feedback in the Early Universe

Authors: Haruka Washinoue, Takeru K. Suzuki

Abstract: 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énic waves. Weaker $B_{\rm c}$ leads to cooler and less dense coronae… ▽ More 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é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. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 11 pages, 5 figures. Accepted for publication in MNRAS

Collisional Growth and Fragmentation of Dust Aggregates with Low Mass Ratios. I: Critical Collision Velocity for Water Ice

Authors: Yukihiko Hasegawa, Takeru K. Suzuki, Hidekazu Tanaka, Hiroshi Kobayashi, Koji Wada

Abstract: 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 co… ▽ More 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. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 33 pages, 17 figures, 3 tables, accepted for publication in ApJ

Photoevaporative Dispersal of Protoplanetary Disks around Evolving Intermediate-mass Stars

Authors: Masanobu Kunitomo, Shigeru Ida, Taku Takeuchi, Olja Panić, James M. Miley, Takeru K. Suzuki

Abstract: 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 (FU… ▽ More 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. △ Less

Submitted 13 March, 2021; originally announced March 2021.

Comments: 21 pages, 15 figures, published in ApJ. Evolutionary models of young 0.5-5 Msun stars are available at https://zenodo.org/record/4595807

Journal ref: ApJ (2021), 909, 109

Alfvén-wave driven magnetic rotator winds from low-mass stars I: rotation dependences of magnetic braking and mass-loss rate

Authors: Munehito Shoda, Takeru K. Suzuki, Sean P. Matt, Steven R. Cranmer, Aline A. Vidotto, Antoine Strugarek, Victor See, Victor Réville, Adam J. Finley, Allan Sacha Brun

Abstract: 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én-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force ex… ▽ More 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én-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 $τ_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( Ω_\ast / Ω_\odot \right)^{2.82}$, yielding a spin-down law $Ω_\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én 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 Ω_\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én radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfvén radius matches previous works which used thermally-driven winds. Our results suggest that the Alfvén-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence. △ Less

Submitted 22 May, 2020; v1 submitted 19 May, 2020; originally announced May 2020.

Comments: accepted for publication in The Astrophysical Journal

New growth mechanism of dust grains in protoplanetary disks with magnetically driven disk winds

Authors: Tetsuo Taki, Koh Kuwabara, Hiroshi Kobayashi, Takeru K. Suzuki

Abstract: 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 t… ▽ More 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} Σ_{\rm d}/Σ_{\rm g}\gtrsim η, \nonumber \end{equation} where $Σ_{\rm d}/Σ_{\rm g}$ is the dust-to-gas surface-density ratio and $η$ 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. △ Less

Submitted 29 December, 2020; v1 submitted 19 April, 2020; originally announced April 2020.

Comments: 24 pages, 10 figures, 1 table, accepted for publication in The Astrophysical Journal

Dispersal of protoplanetary disks by the combination of magnetically driven and photoevaporative winds

Authors: Masanobu Kunitomo, Takeru K. Suzuki, Shu-ichiro Inutsuka

Abstract: 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 P… ▽ More 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. △ Less

Submitted 14 November, 2021; v1 submitted 12 January, 2020; originally announced January 2020.

Comments: 11 pages, 6 figures; corrected Table 2 and Equations (17) and (29); see Erratum at https://doi.org/10.1093/mnras/stab2748

Three-dimensional simulation of the fast solar wind driven by compressible magnetohydrodynamic turbulence

Authors: Munehito Shoda, Takeru Ken Suzuki, Mahboubeh Asgari-Targhi, Takaaki Yokoyama

Abstract: Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of… ▽ More Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of power spectra and structure functions reveals that the turbulence is characterized by its imbalanced (in the sense of outward Alfvenic fluctuations) and anisotropic nature. The density fluctuation originates from the parametric decay instability of outwardly propagating Alfven waves and plays a significant role in the Alfven wave reflection that triggers turbulence. Our conclusion is that the fast solar wind is heated and accelerated by compressible MHD turbulence driven by parametric decay instability and resultant Alfven wave reflection. △ Less

Submitted 19 July, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: published from The Astrophysical Journal Letters

Alfvén wave-driven wind from RGB and AGB stars

Authors: Y. Yasuda, T. K. Suzuki, T. Kozasa

Abstract: We develop a magnetohydrodynamical model of Alfvén wave-driven wind in open magnetic flux tubes piercing the stellar surface of Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, and investigate the physical properties of the winds. The model simulations are carried out along the evolutionary tracks of stars with initial mass in the range of 1.5 to 3.0 $M_{\odot}$ and initial metallic… ▽ More We develop a magnetohydrodynamical model of Alfvén wave-driven wind in open magnetic flux tubes piercing the stellar surface of Red Giant Branch (RGB) and Asymptotic Giant Branch (AGB) stars, and investigate the physical properties of the winds. The model simulations are carried out along the evolutionary tracks of stars with initial mass in the range of 1.5 to 3.0 $M_{\odot}$ and initial metallicity $Z_{\rm ini}$=0.02. The surface magnetic field strength being set to be 1G, we find that the wind during the evolution of star can be classified into the following four types; the first is the wind with the velocity higher than 80 km s$^{-1}$ in the RGB and early AGB (E-AGB) phases; the second is the wind with outflow velocity less than 10 km s$^{-1}$ seen around the tip of RGB or in the E-AGB phase; the third is the unstable wind in the E-AGB and thermally pulsing AGB (TP-AGB) phases; the fourth is the stable massive and slow wind with the mass-loss rate higher than 10$^{-7} M_{\odot}$ yr$^{-1}$ and the outflow velocity lower than 20 km s$^{-1}$ in the TP-AGB phase. The mass-loss rates in the first and second types of wind are two or three orders of magnitude lower than the values evaluated by an empirical formula. The presence of massive and slow wind of the fourth type suggests the possibility that the massive outflow observed in TP-AGB stars could be attributed to the Alfvén wave-driven wind. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 17 pages, 15 figures, accepted for publication in ApJ

MSC Class: 85A04

Magnetohydrodynamics in a Cylindrical Shearing Box

Authors: Takeru K. Suzuki, Tetsuo Taki, Scott S. Suriano

Abstract: 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… ▽ More 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. △ Less

Submitted 1 July, 2019; v1 submitted 10 April, 2019; originally announced April 2019.

Comments: 27 pages, 10 figures, accepted by PASJ, movie is available at http://ea.c.u-tokyo.ac.jp/astro/Members/stakeru/research/cylshbx/

Coronae of Zero/Low-Metal Low-Mass Stars

Authors: Haruka Washinoue, Takeru K. Suzuki

Abstract: 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 o… ▽ More 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 △ Less

Submitted 26 September, 2019; v1 submitted 25 March, 2019; originally announced March 2019.

Comments: 11 pages, 12 figures

Giant protostellar flares: accretion-driven accumulation and reconnection-driven ejection of magnetic flux in protostars

Authors: Shinsuke Takasao, Kengo Tomida, Kazunari Iwasaki, Takeru K. Suzuki

Abstract: 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 di… ▽ More 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. △ Less

Submitted 21 May, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 9 pages, 5 figures. Accepted for publication in ApJ Letters

The formation of rings and gaps in wind-launching non-ideal MHD disks: three-dimensional simulations

Authors: Scott S. Suriano, Zhi-Yun Li, Ruben Krasnopolsky, Takeru K. Suzuki, Hsien Shang

Abstract: 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… ▽ More 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. △ Less

Submitted 16 December, 2018; v1 submitted 4 October, 2018; originally announced October 2018.

Comments: 19 pages, 13 figures, accepted by MNRAS

Hot Grain Dynamics by Electric Charging and Magnetic Trapping in Debris Disks

Authors: Hiroshi Kimura, Masanobu Kunitomo, Takeru K. Suzuki, Jan Robrade, Philippe Thebault, Ikuyuki Mitsuishi

Abstract: 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 s… ▽ More 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, ... △ Less

Submitted 9 August, 2018; originally announced August 2018.

Comments: 30 pages, 8 figures, 3 tables, to appear in Planetary and Space Science

Metal Pollution of Low-Mass Population III Stars through Accretion of Interstellar Objects like `Oumuamua

Authors: Ataru Tanikawa, Takeru K. Suzuki, Yasuo Doi

Abstract: 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… ▽ More 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^{-α}$, Pop.~III survivors can accrete ISO mass $\gtrsim 10^{-16}M_\odot$, or ISO iron mass $\gtrsim 10^{-17}M_\odot$, if $α< 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 $α> 4$, Pop.~III survivors have no chance at colliding with ISOs with $D \gtrsim 3$~km, and keep metal-free. If $3 < α< 4$, Pop.~III survivors would be most polluted by ISOs up to [Fe/H] $\sim -7$. If $α< 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 $α$ and the Galactocentric distance, we first show the importance of ISOs for the metal pollution of Pop.~III survivors. △ Less

Submitted 28 May, 2018; v1 submitted 22 April, 2018; originally announced April 2018.

Comments: 6 pages, 3 figures, accepted for publication in PASJ

Formation of the terrestrial planets in the solar system around 1 au via radial concentration of planetesimals

Authors: Masahiro Ogihara, Eiichiro Kokubo, Takeru K. Suzuki, Alessandro Morbidelli

Abstract: 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… ▽ More 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. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 5 pages, 4 figures, accepted for publication in A&A Letters

Formation of close-in super-Earths in evolving protoplanetary disks due to disk winds

Authors: Masahiro Ogihara, Eiichiro Kokubo, Takeru K. Suzuki, Alessandro Morbidelli

Abstract: 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… ▽ More 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. △ Less

Submitted 3 April, 2018; originally announced April 2018.

Comments: 18 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 615, A63 (2018)

Frequency-dependent Alfvén-wave propagation in the solar wind: Onset and suppression of parametric decay instability

Authors: Munehito Shoda, Takaaki Yokoyama, Takeru K. Suzuki

Abstract: 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… ▽ More 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én waves are subject to PDI. Meanwhile, the maximum growth rate of the PDI of low-frequency ($f_0 \lesssim10^{-4} {\rm \ Hz}$) Alfvén waves becomes negative due to acceleration and expansion effects. Medium-frequency ($f_0 \approx 10^{-3.5} {\rm \ Hz}$) Alfvén 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. △ Less

Submitted 1 May, 2018; v1 submitted 7 March, 2018; originally announced March 2018.

Comments: Accepted for publication in The Astrophysical Journal

A Three-dimensional Simulation of a Magnetized Accretion Disk: Fast Funnel Accretion onto a Weakly Magnetized Star

Authors: Shinsuke Takasao, Kengo Tomida, Kazunari Iwasaki, Takeru K. Suzuki

Abstract: 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 s… ▽ More 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. △ Less

Submitted 6 February, 2019; v1 submitted 22 January, 2018; originally announced January 2018.

Comments: 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

Journal ref: The Astrophysical Journal, 857, 4, 2018

A self-consistent model of the coronal heating and solar wind acceleration including compressible and incompressible heating processes

Authors: Munehito Shoda, Takaaki Yokoyama, Takeru K. Suzuki

Abstract: 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én wave turbulence. Numerical simulations wer… ▽ More 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én wave turbulence. Numerical simulations were conducted with different photospheric perpendicular correlation lengths $λ_0$, which is a critical parameter of Alfvén wave turbulence, and different root-mean-square photospheric transverse-wave amplitudes $δv_0$. For the various $λ_0$, we obtain a low-temperature chromosphere, high-temperature corona, and supersonic solar wind. Our analysis shows that turbulence heating is always dominant when $λ_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én wave reflection and therefore the turbulence heating. The density fluctuation and the cross helicity are strongly affected by $λ_0$, while the coronal temperature and mass loss rate depend weakly on $λ_0$. △ Less

Submitted 14 February, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: Accepted for publicaiton in The Astrophysical Journal

Magnetic Activity in the Galactic Centre Region -- Fast Downflows along Rising Magnetic Loops

Authors: Kensuke Kakiuchi, Takeru K. Suzuki, Yasuo Fukui, Kazufumi Torii, Rei Enokiya, Mami Machida, Ryoji Matsumoto

Abstract: 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 associate… ▽ More 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. △ Less

Submitted 6 March, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: Accepted to MNRAS. 11 pages, 13 figures

Stellar Winds and Coronae of Low-mass Pop. II/III Stars

Authors: Takeru K. Suzuki

Abstract: 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 fr… ▽ More 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. △ Less

Submitted 18 January, 2018; v1 submitted 12 October, 2017; originally announced October 2017.

Comments: 16 pages, 8 figures embedded, comments are welcome

Galactic Outflows by Alfvénic Poynting Flux: Application to Fermi Bubbles

Authors: Takeru K. Suzuki, Alex Lazarian

Abstract: 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… ▽ More 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é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. △ Less

Submitted 19 January, 2018; v1 submitted 11 October, 2017; originally announced October 2017.

Comments: 14 pages, 10 figures included; comments are welcome

Effects of global gas flows on type I migration

Authors: Masahiro Ogihara, Eiichiro Kokubo, Takeru K. Suzuki, Alessandro Morbidelli, Aurélien Crida

Abstract: 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 ai… ▽ More 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. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 9 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 608, A74 (2017)

A Theoretical Model of X-ray Jets from Young Stellar Objects

Authors: Shinsuke Takasao, Takeru K. Suzuki, Kazunari Shibata

Abstract: 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 he… ▽ More 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. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: 31 pages (single column), 5 figures, Accepted to be published in ApJ

Vertical flows and structures excited by magnetic activity in the Galactic center region

Authors: Kensuke Kakiuchi, Takeru K. Suzuki, Yasuo Fukui, Kazufumi Torii, Mami Machida, Ryoji Matsumoto

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2016; v1 submitted 13 October, 2016; originally announced October 2016.

Comments: 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

Investigating Magnetic Activity in the Galactic Centre by Global MHD Simulation

Authors: Takeru K. Suzuki, Yasuo Fukui, Kazufumi Torii, Mami Machida, Ryoji Matsumoto, Kensuke Kakiuchi

Abstract: 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 ampli… ▽ More 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. △ Less

Submitted 31 October, 2016; v1 submitted 7 September, 2016; originally announced September 2016.

Comments: 4 pages, based on a contributed talk at IAU 322 (The Multi-Messenger Astrophysics of the Galactic Centre)

Evolution of Protoplanetary Discs with Magnetically Driven Disc Winds

Authors: Takeru K. Suzuki, Masahiro Ogihara, Alessandro Morbidelli, Aurélien Crida, Tristan Guillot

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2016; v1 submitted 1 September, 2016; originally announced September 2016.

Comments: 16 pages, 11 figures embedded, accepted by A&A (comments are welcome)

Journal ref: A&A 596, A74 (2016)

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

Authors: Shigeo S. Kimura, Kenji Toma, Takeru K. Suzuki, Shu-ichiro Inutsuka

Abstract: 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… ▽ More 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. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 11 pages, 9 figures, re-submitted to ApJ with minor revision

Journal ref: Astrophys. J, 822, 88, 2016

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

Authors: Yuki A. Tanaka, Takeru K. Suzuki, Shu-ichiro Inutsuka

Abstract: 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 Alfveni… ▽ More 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. △ Less

Submitted 16 July, 2015; originally announced July 2015.

Comments: 12 pages, 14 figures, accepted for publication in ApJ

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

Authors: Tomoya Miyake, Takeru K. Suzuki, Shu-ichiro Inutsuka

Abstract: 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… ▽ More 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. △ Less

Submitted 12 February, 2016; v1 submitted 14 May, 2015; originally announced May 2015.

Comments: 8 pages, 6 figures, Submitted to ApJ

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

Authors: Masahiro Ogihara, Hiroshi Kobayashi, Shu-ichiro Inutsuka, Takeru K. Suzuki

Abstract: 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… ▽ More 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. △ Less

Submitted 5 May, 2015; originally announced May 2015.

Comments: 8 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A 579, A65 (2015)

Stochastic Noncircular Motion and Outflows Driven by Magnetic Activity in the Galactic Bulge Region

Authors: Takeru K. Suzuki, Yasuo Fukui, Kazufumi Torii, Mami Machida, Ryoji Matsumoto

Abstract: 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 o… ▽ More 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. △ Less

Submitted 22 September, 2015; v1 submitted 29 April, 2015; originally announced April 2015.

Comments: 12 pages, 13 figures, MNRAS in press

Effects of axions on Nucleosynthesis in massive stars

Authors: Shohei Aoyama, Takeru K. Suzuki

Abstract: 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 consu… ▽ More 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γ}= 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γ}$ from the nucleosynthesis of heavy elements in massive stars. △ Less

Submitted 8 February, 2015; originally announced February 2015.

Comments: 5 pages, 5 figures, submitted to Physical Review D

Connecting the Sun and the Solar Wind: The Self-consistent Transition of Heating Mechanisms

Authors: Takuma Matsumoto, Takeru K. Suzuki

Abstract: We have performed a 2.5 dimensional magnetohydrodynamic simulation that resolves the propagation and dissipation of Alfven waves in the solar atmosphere. Alfvenic fluctuations are introduced on the bottom boundary of the extremely large simulation box that ranges from the photosphere to far above the solar wind acceleration region. Our model is ab initio in the sense that no corona and no wind are… ▽ More We have performed a 2.5 dimensional magnetohydrodynamic simulation that resolves the propagation and dissipation of Alfven waves in the solar atmosphere. Alfvenic fluctuations are introduced on the bottom boundary of the extremely large simulation box that ranges from the photosphere to far above the solar wind acceleration region. Our model is ab initio in the sense that no corona and no wind are assumed initially.The numerical experiment reveals the quasi-steady solution that has the transition from the cool to the hot atmosphere and the emergence of the high speed wind. The global structure of the resulting hot wind solution fairly well agree with the coronal and the solar wind structure inferred from observations. The purpose of this study is to complement the previous paper by Matsumoto & Suzuki (2012) and describe the more detailed results and the analysis method. These results include the dynamics of the transition region and the more precisely measured heating rate in the atmosphere. Particularly, the spatial distribution of the heating rate helps us to interpret the actual heating mechanisms in the numerical simulation.Our estimation method of heating rate turned out to be a good measure for dissipation of Alfven waves and low beta fast waves. △ Less

Submitted 21 February, 2014; v1 submitted 3 February, 2014; originally announced February 2014.

Comments: 24 pages, 15 figures, MNRAS accepted