We report a new evaluation of the accretion properties of PDS~70b obtained with VLT/MUSE. The main difference from previous studies in Haffert et al. (2019) and Aoyama & Ikoma (2019) is in the mass accretion rate. Simultaneous multiple line observations, such as H$\alpha$ and H$\beta$, can better constrain the physical properties of an accreting planet. While we clearly detected H$\alpha$ emissions from PDS~70b, no H$\beta$ emissions were detected. We estimate the line flux of H$\beta$ with a 3-$\sigma$ upper limit to be 2.3~$\times$~10$^{-16}$~erg~s$^{-1}$~cm$^{-2}$. The flux ratio $F_{\rm H\beta}$/$F_{\rm H\alpha}$ for PDS~70b is $<$~0.28. Numerical investigations by Aoyama et al. (2018) suggest that $F_{\rm H\beta}$/$F_{\rm H\alpha}$ should be close to unity if the extinction is negligible. We attribute the reduction of the flux ratio to the extinction, and estimate the extinction of H$\alpha$ ($A_{\rm H\alpha}$) for PDS~70b to be $>$~2.0~mag using the interstellar extinction value. %The expected $A_{\rm H\alpha}$ value in the gap of the protoplanetary disk at the PDS~70b location is 2.4~mag, which is consistent with the estimated extinction. By combining with the H$\alpha$ linewidth and the dereddening line luminosity of H$\alpha$, %we derive the PDS~70b dynamical mass and mass accretion rate to be \hashimotor{12~$\pm$~3~$M_{\rm Jup}$} and $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$, respectively. we derive the PDS~70b mass accretion rate to be $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$. The PDS~70b mass accretion rate is an order of magnitude larger than that of PDS~70. We found that the filling factor $f_{\rm f}$ (the fractional area of the planetary surface emitting H$\alpha$) is $\gtrsim$0.01, which is similar to the typical stellar value. The small value of $f_{\rm f}$ indicates that the H$\alpha$ emitting areas are localized at the surface of PDS~70b.

We herein develop a new simple model for giant planet formation, which predicts the final mass of a giant planet born in a given disk, by adding the disk mass loss due to photoevaporation and a new type II migration formula to our previous model. The proposed model provides some interesting results. First, it gives universal evolution tracks in the diagram of planetary mass and orbital radius, which clarifies how giant planets migrate at growth in the runaway gas accretion stage. Giant planets with a few Jupiter masses or less suffer only a slight radial migration in the runaway gas accretion stage. Second, the final mass of giant planets is approximately given as a function of only three parameters: the initial disk mass at the starting time of runaway gas accretion onto the planet, the mass loss rate due to photoevaporation, and the starting time. On the other hand, the final planet mass is almost independent of the disk radius, viscosity, and initial orbital radius. The obtained final planet mass is similar to or less than 10% of the initial disk mass. Third, the proposed model successfully explains properties in the mass distribution of giant exoplanets with the mass distribution of observed protoplanetary disks for a reasonable range of the mass loss rate due to photoevaporation.

When a planet forms a deep gap in a protoplanetary disk, dust grains cannot pass through the gap. As a consequence, the density of the dust grains can increase up to the same level of the density of the gas at the outer edge. The feedback on the gas from the drifting dust grains is not negligible, in such a dusty region. We carried out two-dimensional two-fluid (gas and dust) hydrodynamic simulations. We found that when the radial flow of the dust grains across the gap is halted, a broad ring of the dust grains can be formed because of the dust feedback and the diffusion of the dust grains. The minimum mass of the planet to form the broad dust ring is consistent with the pebble-isolation mass, in the parameter range of our simulations. The broad ring of the dust grains is good environment for the formation of the protoplanetary solid core. If the ring is formed in the disk around the sun-like star at $\sim 2\ \mbox{AU}$, a massive solid core ($\sim 50M_{\oplus}$) can be formed within the ring, which may be connected to the formation of Hot Jupiters holding a massive solid core such as HD 149026b. In the disk of the dwarf star, a number of Earth-sized planets can be formed within the dust ring around $\sim 0.5\ \mbox{AU}$, which potentially explain the planet system made of multiple Earth-sized planets around the dwarf star such as TRAPPIST-1.

The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-jupiter-mass planets are piled up at > 1 au, while jupiter/sub-jupiter-mass planets are broadly distributed from ~0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al. (2018), by comparing the migration, growth timescales of gas giants, and disk lifetime and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. (2018) proposed that type II migration speed is no other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent MHD simulations. Population synthesis simulations show that super-jupiter-mass planets remain at > 1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at > 1 au would be theoretically solved, if the new formula is confirmed and wind-driven disk accretion dominates.

Progress in understanding of giant planet formation has been hampered by a lack of observational constraints to growing protoplanets. Recently, detection of an H{\alpha}-emission excess via direct imaging was reported for the protoplanet LkCa15b orbiting the pre-main-sequence star LkCa15. However, the physical mechanism for the H{\alpha} emission is poorly understood. According to recent high-resolution three-dimensional hydrodynamic simulations of the flow accreting onto protoplanets, the disk gas flows down almost vertically onto and collides with the surface of a circum-planetary disk at a super-sonic velocity and thus passes through a strong shockwave. The shock-heated gas is hot enough to generate H{\alpha} emission. Here we develop a one-dimensional radiative hydrodynamic model of the flow after the shock by detailed calculations of chemical reactions and electron transitions in hydrogen atoms, and quantify the hydrogen line emission in the Lyman-, Balmer-, and Paschen-series from the accreting gas giant system. We then demonstrate that the H{\alpha} intensity is strong enough to be detected with current observational technique. Comparing our theoretical H{\alpha} intensity with the observed one from LkCa15b, we constrain the protoplanet mass and the disk gas density. Observation of hydrogen line emission from protoplanets is highly encouraged to obtain direct constraints of accreting gas giants, which will be key in understanding the formation of gas giants.

The ubiquity of super-Earths poses a problem for planet formation theory to explain how they avoided becoming gas giants. Rapid recycling of the envelope gas of planets embedded in a protoplanetary disc has been proposed to delay the cooling and following accretion of disc gas. We compare isothermal and non-isothermal 3D hydrodynamical simulations of the gas flow past a planet to investigate the influence on the feasibility of the recycling mechanism. Radiative cooling is implemented by using the beta cooling model. We find that, in either case, gas enters the Bondi sphere at high latitudes and leaves through the midplane regions, or vice versa when disc gas rotates sub-Keplerian. However, in contrast to the isothermal case where the recycling flow reaches the deeper part of the envelope, the inflow is inhibited from reaching the deep envelope in the non-isothermal case. Once the atmosphere starts cooling, buoyant force prevents the high-entropy disc gas from intruding the low-entropy atmosphere. We suggest that the buoyancy barrier isolates the lower envelope from the recycling and allows further cooling, which may lead runaway gas accretion onto the core.

A giant planet embedded in a protoplanetary disk creates a gap. This process is important for both theory and observations. Using results of a survey for a wide parameter range with two-dimensional hydrodynamic simulations, we constructed an empirical formula for the gap structure (i.e., the radial surface density distribution), which can reproduce the gap width and depth obtained by two-dimensional simulations. This formula enables us to judge whether an observed gap is likely to be caused by an embedded planet or not. The propagation of waves launched by the planet is closely connected to the gap structure. It makes the gap wider and shallower as compared with the case where an instantaneous wave damping is assumed. The hydrodynamic simulations shows that the waves do not decay immediately at the launching point of waves, even when the planet is as massive as Jupiter. Based on the results of hydrodynamic simulations, we also obtained an empirical model of wave propagation and damping for the cases of deep gaps. The one-dimensional gap model with our wave propagation model is able to well reproduce the gap structures in hydrodynamic simulations. In the case of a Jupiter-mass planet, we also found that the waves with smaller wavenumber (e.g., $m=2$) are excited and transport the angular momentum to the location far away from the planet. The wave with $m=2$ is closely related with a secondary wave launched by the site opposite from the planet.

We study the accelerating present universe in terms of the time evolution of the equation of state $w(z)$ (redshift $z$) due to thawing and freezing scalar potentials in the quintessence model. The values of $dw/da$ and $d^2w/da^2$ at scale factor of $a = 1$ are associated with two parameters of each potential. For five type scalar potentials, the scalar fields $Q$ and $w$ as function of time $t$ and/or $z$ are numerically calculated under the fixed boundary condition of $w(z=0)=-1+\Delta$. The observational constraint $w_{obs}$ (Planck 2015) is imposed to test whether the numerical $w(z)$ is in $w_{obs}$. Some solutions show the thawing features in the freezing potentials. Mutually exclusive allowed regions in $dw/da$ vs. $d^2w/da^2$ diagram are obtained in order to identify the likely scalar potential and even the potential parameters by future observational tests.

We present a near-infrared direct imaging search for accretion signatures of possible protoplanets around the young stellar object (YSO) TW Hya, a multi-ring disk exhibiting evidence of planet formation. The Pa$\beta$ line (1.282 $\mu$m) is an indication of accretion onto a protoplanet, and its intensity is much higher than that of blackbody radiation from the protoplanet. We focused on the Pa$\beta$ line and performed Keck/OSIRIS spectroscopic observations. Although spectral differential imaging (SDI) reduction detected no accretion signatures, the results of the present study allowed us to set 5$\sigma$ detection limits for Pa$\beta$ emission of $5.8\times10^{-18}$ and $1.5\times10^{-18}$ erg/s/cm$^2$ at 0\farcs4 and 1\farcs6, respectively. We considered the mass of potential planets using theoretical simulations of circumplanetary disks and hydrogen emission. The resulting masses were $1.45\pm 0.04$ M$_{\rm J}$ and $2.29 ^{+0.03}_{-0.04}$ M$_{\rm J}$ at 25 and 95 AU, respectively, which agree with the detection limits obtained from previous broadband imaging. The detection limits should allow the identification of protoplanets as small as $\sim$1 M$_{\rm J}$, which may assist in direct imaging searches around faint YSOs for which extreme adaptive optics instruments are unavailable.

A giant planet creates a gap in a protoplanetary disk, which might explain the observed gaps in protoplanetary disks. The width and depth of the gaps depend on the planet mass and disk properties. We have performed two--dimensional hydrodynamic simulations for various planet masses, disk aspect ratios and viscosities, to obtain an empirical formula for the gap width. The gap width is proportional to the square root of the planet mass, -3/4 power of the disk aspect ratio and -1/4 power of the viscosity. This empirical formula enables us to estimate the mass of a planet embedded in the disk from the width of an observed gap. We have applied the empirical formula for the gap width to the disk around HL~Tau, assuming that each gap observed by ALMA observations is produced by planets, and discussed the planet masses within the gaps. The estimate of planet masses from the gap widths is less affected by the observational resolution and dust filtration than that from the gap depth.

Firstly, we study the final masses of giant planets growing in protoplanetary disks through capture of disk gas, by employing an empirical formula for the gas capture rate and a shallow disk gap model, which are both based on hydrodynamical simulations. The shallow disk gaps cannot terminate growth of giant planets. For planets less massive than 10 Jupiter masses, their growth rates are mainly controlled by the gas supply through the global disk accretion, rather than their gaps. The insufficient gas supply compared with the rapid gas capture causes a depletion of the gas surface density even at the outside of the gap, which can create an inner hole in the protoplanetary disk. Our model can also predict the depleted gas surface density in the inner hole for a given planet mass. Secondly, our findings are applied to the formation of our solar system. For the formation of Jupiter, a very low-mass gas disk with a few or several Jupiter masses is required at the beginning of its gas capture because of the non-stopping capture. Such a low-mass gas disk with sufficient solid material can be formed through viscous evolution from an initially $\sim$10AU-sized compact disk with the solar composition. By the viscous evolution with a moderate viscosity of $\alpha \sim 10^{-3}$, most of disk gas accretes onto the sun and a widely spread low-mass gas disk remains when the solid core of Jupiter starts gas capture at $t \sim 10^7$ yrs. The depletion of the disk gas is suitable for explaining the high metallicity in giant planets of our solar system. A very low-mass gas disk also provides a plausible path where type I and II planetary migrations are both suppressed significantly. In particular, we also show that the type II migration of Jupiter-size planets becomes inefficient because of the additional gas depletion due to the rapid gas capture by themselves.

A giant planet embedded in a protoplanetary disk forms a gap. An analytic relationship among the gap depth, planet mass $M_{p}$, disk aspect ratio $h_p$, and viscosity $\alpha$ has been found recently, and the gap depth can be written in terms of a single parameter $K= (M_{p}/M_{\ast})^2 h_p^{-5} \alpha^{-1}$. We discuss how observed gap features can be used to constrain the disk and/or planet parameters based on the analytic formula for the gap depth. The constraint on the disk aspect ratio is critical in determining the planet mass so the combination of the observations of the temperature and the image can provide a constraint on the planet mass. We apply the formula for the gap depth to observations of HL~Tau and HD~169142. In the case of HL~Tau, we propose that a planet with $\gtrsim 0.3$ is responsible for the observed gap at $30$~AU from the central star based on the estimate that the gap depth is $\lesssim 1/3$. In the case of HD~169142, the planet mass that causes the gap structure recently found by VLA is $\gtrsim 0.4 M_J$. We also argue that the spiral structure, if observed, can be used to estimate the lower limit of the disk aspect ratio and the planet mass.

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

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

We investigate accretion of solid materials onto circumplanetary disks from heliocentric orbits rotating in protoplanetary disks, which is a key process for the formation of regular satellite systems. In the late stage of gas-capturing phase of giant planet formation, the accreting gas from protoplanetary disks forms circumplanetary disks. Since the accretion flow toward the circumplanetary disks affects the particle motion through gas drag force, we use hydrodynamic simulation data for the gas drag term to calculate the motion of solid materials. We consider wide range of size for the solid particles ($10^{-2}$-$10^6$m), and find that the accretion efficiency of the solid particles peaks around 10m-sized particles because energy dissipation of drag with circum-planetary disk gas in this size regime is most effective. The efficiency for particles larger than 10m size becomes lower because gas drag becomes less effective. For particles smaller than 10m, the efficiency is lower because the particles are strongly coupled with the back-ground gas flow, which prevent particles from accretion. We also find that the distance from the planet where the particles are captured by the circumplanetary disks is in a narrow range and well described as a function of the particle size.

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

We investigate gas accretion flow onto a circumplanetary disk from a protoplanetary disk in detail by using high-resolution three-dimensional nested-grid hydrodynamic simulations, in order to provide a basis of formation processes of satellites around giant planets. Based on detailed analyses of gas accretion flow, we find that most of gas accretion onto circumplanetary disks occurs nearly vertically toward the disk surface from high altitude, which generates a shock surface at several scale heights of the circumplanetary disk. The gas that has passed through the shock surface moves inward because its specific angular momentum is smaller than that of the local Keplerian rotation, while gas near the midplane in the protoplanetary disk cannot accrete to the circumplanetary disk. Gas near the midplane within the planet's Hill sphere spirals outward and escapes from the Hill sphere through the two Lagrangian points L$_1$ and L$_2$. We also analyze fluxes of accreting mass and angular momentum in detail and find that the distributions of the fluxes onto the disk surface are well described by power-law functions and that a large fraction of gas accretion occurs at the outer region of the disk, i.e., at about 0.1 times the Hill radius. The nature of power-law functions indicates that, other than the outer edge, there is no specific radius where gas accretion is concentrated. These source functions of mass and angular momentum in the circumplanetary disk would provide us with useful constraints on the structure and evolution of the circumplanetary disk, which is important for satellite formation.

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals' eccentricity, inclination, planet's radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals' kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with one to ten Earth-masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals.

We construct an analytic model for the rate of gas accretion onto a planet embedded in a protoplanetary disk as a function of planetary mass, disk viscosity, disk scale height, and unperturbed surface density in order to study the long-term accretion and final masses of gas giant planets. We first derive an analytical formula for surface density profile near the planetary orbit from considerations of the balance of force and the dynamical stability. Using it in the empirical formula linking surface density with gas accretion rate that is derived based on hydrodynamic simulations of Tanigawa and Watanabe (2002, ApJ 586, 506), we then simulate the mass evolution of gas giant planets in viscously-evolving disks. We finally determine the final mass as a function of semi-major axis of the planet. We find that the disk can be divided into three regions characterized by different processes by which the final mass is determined. In the inner region, the planet grows quickly and forms a deep gap to suppress the growth by itself before disk dissipation. The final mass of the planet in this region is found to increase with the semi-major axis in a similar way to the mass given by the viscous condition for gap opening, but the former is larger by a factor of approximately 10 than the latter. In the intermediate region, viscous diffusion of the disk gas limits the gas accretion before the planet form a deep gap. The final mass can be up to the disk mass, when disk viscous evolution occurs faster than disk evaporation. In the outer region, planets capture only tiny amounts of gas within the lifetime of the disk to form Neptune-like planets. We derive analytic formulae for the final masses in the different regions and the locations of the boundaries, which are helpful to gain a systematic understanding of the masses of gas giant planets.

We study the ^{1}S_{0} proton pairing gap in neutron star matter putting emphasis on influence of the Dirac effective mass and the proton fraction on the gap within the relativistic Hartree-Bogoliubov model. The gap equation is solved using the Bonn-B potential as a particle-particle channel interaction. It is found that the maximal pairing gap \Delta_{max} is 1-2 MeV, which has a strong correlation with the Dirac effective mass. Hence we suggest that it serves as a guide to narrow down parameter sets of the relativistic effective field theory. Furthermore, the more slowly protons increase with density in the core region of neutron stars, the wider the superfluid range and the slightly lower the peak of the gap become.

The Kelvin--Helmholtz (KH) and tearing instabilities are likely to be important for the process of fast magnetic reconnection that is believed to explain the observed explosive energy release in solar flares. Theoretical studies of the instabilities, however, typically invoke simplified initial magnetic and velocity fields that are not solutions of the governing magnetohydrodynamic (MHD) equations. In the present study, the stability of a reconnecting current sheet is examined using a class of exact global MHD solutions for steady state incompressible magnetic reconnection, discovered by Craig & Henton. Numerical simulation indicates that the outflow solutions where the current sheet is formed by strong shearing flows are subject to the KH instability. The inflow solutions where the current sheet is formed by a fast and weakly sheared inflow are shown to be tearing unstable. Although the observed instability of the solutions can be interpreted qualitatively by applying standard linear results for the KH and tearing instabilities, the magnetic field and plasma flow, specified by the Craig--Henton solution, lead to the stabilization of the current sheet in some cases. The sensitivity of the instability growth rate to the global geometry of magnetic reconnection may help in solving the trigger problem in solar flare research.

We calculate a \Lambda\Lambda pairing gap in binary mixed matter of nucleons and \Lambda hyperons within the relativistic Hartree-Bogoliubov model. Lambda hyperons to be paired up are immersed in background nucleons in a normal state. The gap is calculated with a one-boson-exchange interaction obtained from a relativistic Lagrangian. It is found that at background density \rho_{N}=2.5\rho_{0} the \Lambda\Lambda pairing gap is very small, and that denser background makes it rapidly suppressed. This result suggests a mechanism, specific to mixed matter dealt with relativistic models, of its dependence on the nucleon density. An effect of weaker \Lambda\Lambda attraction on the gap is also examined in connection with revised information of the \Lambda\Lambda interaction.