Source author record

Takeru K. Suzuki

Takeru K. Suzuki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

26 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

Coronae of Zero/Low-Metal Low-Mass Stars

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

preprint2016arXiv

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

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

preprint2016arXiv

Evolution of Protoplanetary Discs with Magnetically Driven Disc Winds

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.

preprint2016arXiv

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

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

preprint2015arXiv

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

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

preprint2015arXiv

Effects of axions on Nucleosynthesis in massive stars

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.

preprint2015arXiv

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

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

preprint2015arXiv

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

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.

preprint2014arXiv

A Signature of Chromospheric Activity in Brown Dwarfs Revealed by 2.5-5.0 Micron AKARI Spectra

We propose that the 2.7 micron H_2O, 3.3 micron CH_4 and 4.6 micron CO absorption bands can be good tracers of chromospheric activity in brown dwarfs. In our previous study, we found that there are difficulties in explaining entire spectra between 1.0 and 5.0 microns with the Unified Cloudy Model (UCM), a brown dwarf atmosphere model. Based on simple radiative equilibrium, temperature in a model atmosphere usually decreases monotonically with height. However, if a brown dwarf has a chromosphere, as inferred by some observations, the temperature in the upper atmosphere is higher. We construct a simple model that takes into account heating due to chromospheric activity by setting a temperature floor in an upper atmosphere, and find that the model spectra of 3 brown dwarfs with moderate H-alpha emission, an indicator of chromospheric activity, are considerably improved to match the AKARI spectra. Because of the higher temperatures in the upper atmospheres, the amount of CH_4 molecules is reduced and the absorption band strengths become weaker. The strengths of the absorption bands of H_2O and CO also become weaker. On the other hand, other objects with weak H-alpha emission cannot be fitted by our treatment. We also briefly discuss magnetic heating processes which possibly operate in upper atmospheres, by extending our numerical simulations for the Sun and stars with surface convection to brown dwarf atmospheres.

preprint2014arXiv

Atmospheric escape by magnetically driven wind from gaseous planets

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

preprint2014arXiv

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

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.

preprint2014arXiv

MHD Simulations of Global Accretion Disks with Vertical Magnetic Fields

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

preprint2013arXiv

Driving Disk Winds and Heating up Hot Coronae by MRI Turbulence

We investigate the formation of hot coronae and vertical outflows in accretion disks by magneto-rotational turbulence. We perform local three-dimensional (3D) MHD simulations with the vertical stratification by explicitly solving an energy equation with various effective ratios of specific heats, gamma. Initially imposed weak vertical magnetic fields are effectively amplified by magnetorotational instability (MRI) and winding due to the differential rotation. In the isothermal case (gamma=1), the disk winds are mainly driven by the Poynting flux associated with the MHD turbulence and show quasi-periodic intermittency. On the other hand, in the non-isothermal cases with gamma >~ 1.1, the regions above 1-2 scale heights from the midplane are effectively heated up to form coronae with the temperature of ~ 50 times of the initial value, which are connected to the cooler midplane region through the pressure-balanced transition regions. As a result, the disk winds are mainly driven by the gas pressure with exhibiting more time-steady nature, although the nondimensional time-averaged mass loss rates are similar to that of the isothermal case. Sound-like waves are confined in the cool midplane region in these cases, and the amplitude of the density fluctuations is larger than that of the isothermal case.

preprint2013arXiv

Evolution of solar-type stellar wind

By extending our self-consistent MHD simulations for the solar wind, we study the evolution of stellar winds of solar-type stars from early main sequence stage to red giant phase. Young solar-type stars are active and the mass loss rates are larger by up to ~ 100 times than that of the present-day sun. We investigate how the stellar wind is affected when the magnetic field strength and fluctuation amplitude at the photosphere increase. While the mass loss rate sensitively depends on the input energy from the surface because of the global instability related to the reflection and nonlinear dissipation of Alfven waves, it saturates at ~ 100 times because most of the energy is used up for the radiative losses rather than the kinetic energy of the wind. After the end of the main sequence phase when the stellar radius expands by ~ 10 times, the steady hot corona with temperature $10^6$ K, suddenly disappears. Chromospheric materials, with hot bubbles embedded owing to thermal instability, directly stream out; the red giant wind is not a steady stream but structured outflow.

preprint2013arXiv

Saturation of Stellar Winds from Young Suns

We investigate mass losses via stellar winds from sun-like main sequence stars with a wide range of activity levels. We perform forward-type magnetohydrodynamical numerical experiments for Alfven wave-driven stellar winds with a wide range of the input Poynting flux from the photosphere. Increasing the magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, the mass loss rate rapidly increases at first owing to the suppression of the reflection of the Alfven waves. The surface materials are lifted up by the magnetic pressure associated with the Alfven waves, and the cool dense chromosphere is intermittently extended to 10 -- 20 % of the stellar radius. The dense atmospheres enhance the radiative losses and eventually most of the input Poynting energy from the stellar surface escapes by the radiation. As a result, there is no more sufficient energy remained for the kinetic energy of the wind; the stellar wind saturates in very active stars, as observed in Wood et al. The saturation level is positively correlated with B_{r,0}f_0, where B_{r,0} and f_0 are the magnetic field strength and the filling factor of open flux tubes at the photosphere. If B_{r,0}f_0 is relatively large >~ 5 G, the mass loss rate could be as high as 1000 times. If such a strong mass loss lasts for ~ 1 billion years, the stellar mass itself is affected, which could be a solution to the faint young sun paradox. We derive a Reimers-type scaling relation that estimates the mass loss rate from the energetics consideration of our simulations. Finally, we derive the evolution of the mass loss rates, \dot{M} t^{-1.23}, of our simulations, combining with an observed time evolution of X-ray flux from sun-like stars, which is shallower than \dot{M} t^{-2.33+/-0.55} in Wood et al.(2005).

preprint2011arXiv

Solar Wind and its Evolution

By using our previous results of magnetohydrodynamical simulations for the solar wind from open flux tubes, I discuss how the solar wind in the past is different from the current solar wind. The simulations are performed in fixed one-dimensional super-radially open magnetic flux tubes by inputing various types of fluctuations from the photosphere, which automatically determines solar wind properties in a forward manner. The three important parameters which determine physical properties of the solar wind are surface fluctuation, magnetic field strengths, and the configuration of magnetic flux tubes. Adjusting these parameters to the sun at earlier times in a qualitative sense, I infer that the quasi-steady-state component of the solar wind in the past was denser and slightly slower if the effect of the magneto-centrifugal force is not significant. I also discuss effects of magneto-centrifugal force and roles of coronal mass ejections.

preprint2010arXiv

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

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

preprint2010arXiv

Self-consistent Simulations of Alfven Wave Driven Winds from the Sun and Stars

We review our recent results of Alfven wave-driven winds. First, we present the result of a self-consistent 1D MHD simulations for solar winds from the photosphere to interplanetary region. Here, we emphasize the importance of the reflection of Alfven waves in the density stratified corona and solar winds. We also introduce the recent HINODE observation that might detect the reflection signature of transverse (Alfvenic) waves by Fujimura & Tsuneta (2009). Then, we show the results of Alfven wave-driven winds from red giant stars. We explain the change of the atmosphere properties from steady coronal winds to intermittent chromospheric winds and discuss how the wave reflection is affected by the decrease of the surface gravity with stellar evolution. We also discuss similarities and differences of accretion disk winds by MHD turbulence.

preprint2010arXiv

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

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

preprint2009arXiv

Thermal Response of A Solar-like Atmosphere to An Electron Beam from A Hot Jupiter: A Numerical Experiment

We investigate the thermal response of the atmosphere of a solar-type star to an electron beam injected from a hot Jupiter by performing a 1-dimensional magnetohydrodynamic numerical experiment with non-linear wave dissipation, radiative cooling, and thermal conduction. In our experiment, the stellar atmosphere is non-rotating and is modelled as a 1-D open flux tube expanding super-radially from the stellar photosphere to the planet. An electron beam is assumed to be generated from the reconnection site of the planet's magnetosphere. The effects of the electron beam are then implemented in our simulation as dissipation of the beam momentum and energy at the base of the corona where the Coulomb collisions become effective. When the sufficient energy is supplied by the electron beam, a warm region forms in the chromosphere. This warm region greatly enhances the radiative fluxes corresponding to the temperature of the chromosphere and transition region. The warm region can also intermittently contribute to the radiative flux associated with the coronal temperature due to the thermal instability. However, owing to the small area of the heating spot, the total luminosity of the beam-induced chromospheric radiation is several orders of magnitude smaller than the observed Ca II emissions from HD 179949.

preprint2007arXiv

Structured Red Giant Winds with Magnetized Hot Bubbles and the Corona/Cool Wind Dividing Line

By performing MHD simulations, we investigate the mass loss of intermediate- and low-mass stars from main sequence (MS) to red giant branch (RGB) phases. Alfven waves, which are excited by the surface convections travel outwardly and dissipate by nonlinear processes to accelerate and heat the stellar winds. We dynamically treat these processes in open magnetic field regions from the photospheres to 25 stellar radii. When the stars evolve to slightly blueward positions of the dividing line (Linsky & Haisch), the steady hot corona with temperature, ~ 1MK, suddenly disappears. Instead, many hot (~1MK) and warm (~10^5K) bubbles are formed in cool (T<~2x10^4K) chromospheric winds because of thermal instability; the red giant wind is not a steady stream but structured outflow. As a result, the mass loss rates, \dot{M}, largely vary in time by 3-4 orders or magnitude in the RGB stars. Supported by magnetic pressure, the density of hot bubbles can be kept low to reduce the radiative cooling and to maintain the high temperature long time. Even in the stars redward of the dividing line, hot bubbles intermittently exist, and they can be sources of UV/soft X-ray emissions from hybrid stars. Nearly static regions are formed above the photospheres of the RGB stars, and the stellar winds are effectively accelerated from several stellar radii. Then, the wind velocity is much smaller than the surface escape speed, because it is regulated by the slower escape speed at that location. We finally derive an equation that determines \dot{M} from the energetics of the simulated wave-driven winds in a forward manner. The relation explains \dot{M} from MS to RGB, and it can play a complementary role to the Reimers' formula, which is mainly for more luminous stars.

preprint2006arXiv

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

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

preprint2005arXiv

Alfven Wave-Driven Proto-Neutron Star Winds And R-Process Nucleosynthesis

We propose magnetic proto-neutron star (PNS) winds driven by \Alfven waves as well as the neutrino heating as an appropriate site for the r-process nucleosynthesis. \Alfven waves excited by surface motions of a PNS propagate outwardly, and they heat and accelerate the wind by dissipation. Compared with the wind purely driven by the neutrino heating, larger entropy per baryon and shorter dynamical time scale are achieved, which favors the r-process. We study reasonable cases that a wave amplitude is 10% of \Alfven speed at the surface to find that a PNS with surface field strength, $\gtrsim 5\times 10^{14}$G, gives suitable wind properties for the r-process, provided that dissipation length of the wave is at most $\sim$ 10 times of the PNS radius. We also compare properties of transcritical and subcritical winds in light of the r-process. We finally discuss possibilities of detections of $γ$-rays from radioactive nuclei and absorption lines due to Ba in supernova remnants which possess magnetars.