Recent detection of exoplanets with Earth-like insolation attracts growing interest in how common Earth-like aqua planets are beyond the solar system. While terrestrial planets are often assumed to capture icy or water-rich planetesimals, a primordial atmosphere of nebular origin itself can produce water through oxidation of the atmospheric hydrogen with oxidising minerals from incoming planetesimals or the magma ocean. Thermodynamically, normal oxygen buffers produce water comparable in mole number to or more than hydrogen. Thus, the primordial atmosphere would likely be highly enriched with water vapour; however, the primordial atmosphereshave been always assumed to have the solar abundances. Here we integrate the 1D structure of such an enriched atmosphere of sub-Earths embedded in a protoplanetary disc around an M dwarf of 0.3$M_\odot$ and investigate the effects of water enrichment on the atmospheric properties with focus on water amount. We find that the well-mixed, highly-enriched atmosphere is more massive by a few orders of magnitude than the solar-abundance atmosphere, and that even a Mars-mass planet can obtain water comparable to the present Earth's oceans. Although close-in Mars-mass planets likely lose the captured water via disc dispersal and photo-evaporation, these results suggest that there are more sub-Earths with Earth-like water contents than previously predicted. How much water terrestrial planets really obtain and retain against subsequent loss, however, depends on efficiencies of water production, mixing in the atmosphere and magma ocean, and photo-evaporation, detailed investigation for which should be made in the future.

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.

Context.Recent studies suggest that many giant exoplanets are highly enriched with heavy elements compared to their host star andcontain several tens of Earth masses or more of heavy elements. Such enrichment is considered to have been brought by accretionof planetesimals in late formation stages. Previous dynamical simulations, however, show that planets are unable to collect so muchheavy elements throughin situplanetesimal accretion. Aims.We investigate whether a giant planet migrating inward can capture planetesimals efficiently to significantly increase its metal-licity. Methods.We performed orbital integrations of a migrating giant planet and planetesimals in a protoplanetary gas disc to infer theplanetesimal mass that is accreted by the planet. Results.We find that the two shepherding processes of mean motion resonances trapping and aerodynamic gas drag inhibit plan-etesimal capture of a migrating planet. However, the amplified libration allows the highly-excited planetesimals in the resonances toescape from the resonance trap and be accreted by the planet. Consequently, we show that a migrating giant planet captures planetes-imals with total mass of several tens of Earth masses, if the planet forms at a few tens of AU in a relatively massive disc. We alsofind that planetesimal capture occurs efficiently in a limited range of semi-major axis, and that the total captured planetesimal massincreases with increasing migration distances. Our results have important implications for understanding the relation between giantplanet metallicity and mass, as we suggest that it reflects the formation location of the planet, or more precisely, the location whererunaway gas accretion occurred. We also suggest the observed metal-rich close-in Jupiters migrated to their present locations fromafar, where they formed.

Recent observations of protoplanets embedded in circumstellar disks have shed light on the planet formation process. In particular, detection of hydrogen Balmer-line (H{\alpha}) emission gives direct constraints on late-stage accretion onto gas giants. Very recently Haffert et al. (2019) measured the spectral line-widths, in addition to intensities, of H{\alpha} emission from the two protoplanets orbiting PDS 70. Here, we study these protoplanets by applying radiation-hydrodynamic models of the shock-heated accretion flow onto protoplanets that Aoyama et al. (2018) has recently developed. As a result, we demonstrate that H{\alpha} line-widths combined with intensities lead to narrowing down the possible ranges of the protoplanetary accretion rate and/or mass significantly. While the current spectral resolution is not high enough to derive a definite conclusion regarding their accretion process, high-resolution spectral imaging of growing protoplanets is highly promising.

Terrestrial planets covered globally with thick oceans (termed ocean planets) in the habitable zone were previously inferred to have extremely hot climates in most cases. This is because ${\rm H_2O}$ high-pressure (HP) ice on the seafloor prevents chemical weathering and, thus, removal of atmospheric CO$_2$. Previous studies, however, ignored melting of the HP ice and horizontal variation in heat flux from oceanic crusts. Here we examine whether high heat fluxes near the mid-ocean ridge melts the HP ice and thereby removes atmospheric ${\rm CO_2}$. We develop integrated climate models of an Earth-size ocean planet with plate tectonics for different ocean masses, which include the effects of HP ice melting, seafloor weathering, and the carbonate-silicate geochemical carbon cycle. We find that the heat flux near the mid-ocean ridge is high enough to melt the ice, enabling seafloor weathering. In contrast to the previous theoretical prediction, we show that climates of terrestrial planets with massive oceans lapse into extremely cold ones (or snowball states) with CO$_2$-poor atmospheres. Such extremely cold climates are achieved mainly because the HP ice melting fixes seafloor temperature at the melting temperature, thereby keeping a high weathering flux regardless of surface temperature. We estimate that ocean planets with oceans several tens of the Earth's ocean mass no longer maintain temperate climates. These results suggest that terrestrial planets with extremely cold climates exist even in the habitable zone beyond the solar system, given the frequency of water-rich planets predicted by planet formation theories.

Studies of internal structure of gas giant planets suggest that their envelopes are enriched with heavier elements than hydrogen and helium relative to their central stars. Such enrichment likely occurred by solid accretion during late formation stages of gas giant planets in which gas accretion dominates protoplanetary growth. Some previous studies performed orbital integration of planetesimals around a growing protoplanet with the assumption of an uniform circumstellar disc to investigate how efficiently the protoplanet captures planetesimals. However, not only planetesimals but also disc gas are gravitationally perturbed by the protoplanet in its late formation stages, resulting in gap opening in the circumstellar disc. In this study, we investigate the effects of such gap formation on the capture of planetesimals by performing dynamical simulations of planetesimals around a growing proto-gas giant planet. Gap formation reduces the surface density of disc gas, makes a steep pressure gradient and limits the growth rate of the protoplanet. We find that the first effect enhances the capture of planetesimals, while the others reduce it. Consequently the amount of planetesimals captured during the gas accretion is estimated to be at most ~3 Mearth. We conclude that the in-situ capture of planetesimals needs the initial solid surface density more than five times higher than that of the minimum mass solar nebula for explaining the inferred large amount of heavy element in Jupiter. For highly dense warm Jupiters, we would need additional processes enhancing the capture and/or supply of planetesimals.

Recent transmission spectroscopy has revealed that clouds and hazes are common in the atmospheres of close-in exoplanets. In this study, using the photochemical, microphysical, and transmission spectrum models for close-in warm ($\lesssim$ 1000 K) exoplanet atmospheres that we newly developed in our preceding paper (Kawashima & Ikoma 2018), we investigate the vertical distributions of haze particles and gaseous species and the resultant transmission spectra over wide ranges of the model parameters including UV irradiation intensity, metallicity, carbon-to-oxygen ratio (C/O), eddy diffusion coefficient, and temperature. The sensitivity to metallicity is of particular interest. We find that a rise in metallicity leads basically to reducing the photodissociation rates of the hydrocarbons and therefore the haze monomer production rates. This is due to an enhanced photon-shielding effect by the major photon absorbers such as $\mathrm{H_2O}$, $\mathrm{CO}$, $\mathrm{CO_2}$, and $\mathrm{O_2}$, existing at higher altitudes than the hydrocarbons. We also find that at relatively short wavelengths ($\lesssim$ 2-3 $\mu$m), the absorption features in transmission spectra are most pronounced for moderate metallicities such as 100 times the solar metallicity, whereas the lower the metallicity the stronger the absorption features at relatively long wavelengths ($\gtrsim$ 2-3 $\mu$m), where the contribution of haze is small. These are because of the two competing effects of reduced haze production rate and atmospheric scale height for higher metallicities. For the other model parameters, we show that stronger absorption features appear in transmission spectra of the atmospheres with lower UV irradiation, lower C/O ratio, higher eddy diffusion coefficient, and higher temperature.

Some of the exoplanets so far observed show featureless or flat transmission spectra, possibly indicating the existence of clouds and/or haze in their atmospheres. Thanks to its large aperture size and broad wavelength coverage, JWST is expected to enable detailed investigation of exoplanet atmospheres, which could provide important constraints on the atmospheric composition obscured by clouds/haze. Here, we use four warm ($\lesssim 1000$ K) planets suitable for atmospheric characterization via transmission spectroscopy, GJ 1214b, GJ 436b, HD 97658b, and Kepler-51b, as examples to explore molecular absorption features detectable by JWST even in the existence of hydrocarbon haze in the atmospheres. We simulate photochemistry, the growth of hydrocarbon haze particles, and transmission spectra for the atmospheres of these four planets. Among the planetary parameters considered, super-Earths with hazy, relatively hydrogen-rich atmospheres are mostly expected to produce detectable molecular absorption features such as a quite prominent $\mathrm{CH_4}$ feature at 3.3 ${\rm \mu}$m even for the extreme case of the most efficient production of photochemical haze. For a planet that has extremely low gravity, such as Kepler-51b, haze particles grow significantly large in the upper atmosphere due to the small sedimentation velocity, resulting in the featureless or flat transmission spectrum in a wide wavelength range. This investigation shows that the transmission spectra with muted features measured by HST in most cases do not preclude strong features at the longer wavelengths accessible by JWST.

The World Space Observatory for Ultraviolet (WSO-UV) is an orbital optical telescope with a 1.7 m-diameter primary mirror currently under development. The WSO-UV is aimed to operate in the 115-310 nm UV spectral range. Its two major science instruments are UV spectrographs and a UV imaging field camera with filter wheels. The WSO-UV project is currently in the implementation phase, with a tentative launch date in 2023. Recently, two additional instruments devoted to exoplanets have been proposed for WSO-UV, which are the focus of this paper. UVSPEX, a UV-Spectrograph for Exoplanets, aims to determine atomic hydrogen and oxygen abundance in the exospheres of terrestrial exoplanets. The spectral range is 115-130 nm which enables simultaneous measurement of hydrogen and oxygen emission intensities during an exoplanet transit. Study of exosphere transit photometric curves can help differentiate among different types of rocky planets. The exospheric temperature of an Earth-like planet is much higher than that of a Venus-like planet, because of the low mixing ratio of the dominant coolant (CO2) in the upper atmosphere of the former, which causes a large difference in transit depth at the oxygen emission line. Thus, whether the terrestrial exoplanet is Earth-like, Venus-like, or other can be determined. SCEDI, a Stellar Coronagraph for Exoplanet Direct Imaging is aimed to directly detect the starlight reflected from exoplanets orbiting their parent stars or from the stellar vicinity including circumstellar discs, dust, and clumps. SCEDI will create an achromatic (optimized to 420-700 nm wavelength range), high-contrast stellocentric coronagraphic image of a circumstellar vicinity. The two instruments: UVSPEX and SCEDI, share common power and control modules. The present communication outlines the science goals of both proposed instruments and explains some of their engineering features.

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.

We explore the origin of the trend of heavy elements in observed massive exoplanets. Coupling of better measurements of mass ($M_p$) and radius of exoplanets with planet structure models enables estimating the total heavy element mass ($M_Z$) in these planets. The corresponding relation is characterized by a power-law profile, $M_Z \propto M_p^{3/5}$. We develop a simplified, but physically motivated analysis to investigate how the power-law profile can be produced under the current picture of planet formation. Making use of the existing semi-analytical formulae of accretion rates of pebbles and planetesimals, our analysis shows that the relation can be reproduced well if it traces the final stage of planet formation. In the stage, planets accrete solids from gapped planetesimal disks and gas accretion is limited by disk evolution. We also find that dust accretion accompanying with gas accretion does not contribute to $M_Z$ for planets with $M_p < 10^3 M_{\oplus}$. Our findings are broadly consistent with that of previous studies, yet we explicitly demonstrate how planetesimal dynamics is crucial for better understanding the relation. While our approach is simple, we can also reproduce the trend of a correlation between planet metallicity and $M_p$ that is obtained by detailed population synthesis calculations, when the same assumption is adopted. Our analysis suggests that pebble accretion would not play a direct role at the final stage of planet formation, whereas radial drift of pebbles might be important indirectly for metal enrichment of planets. Detailed numerical simulations and more observational data are required for confirming our analysis.

The water content of magma oceans is widely accepted as a key factor that determines whether a terrestrial planet is habitable. Water ocean mass is determined as a result not only of water delivery and loss, but also of water partitioning among several reservoirs. Here we review our current understanding of water partitioning among the atmosphere, magma ocean, and solid mantle of accreting planetary embryos and protoplanets just after giant collisions. Magma oceans are readily formed in planetary embryos and protoplanets in their accretion phase. Significant amounts of water are partitioned into magma oceans, provided the planetary building blocks are water-rich enough. Particularly important but still quite uncertain issues are how much water the planetary building blocks contain initially and how water goes out of the solidifying mantle and is finally degassed to the atmosphere. Constraints from both solar-system explorations and exoplanet observations and also from laboratory experiments are needed to resolve these issues.

Recently, properties of exoplanet atmospheres have been constrained via multi-wavelength transit observation, which measures an apparent decrease in stellar brightness during planetary transit in front of its host star (called transit depth). Sets of transit depths so far measured at different wavelengths (called transmission spectra) are somewhat diverse: Some show steep spectral slope features in the visible, some contain featureless spectra in the near-infrared, some show distinct features from radiative absorption by gaseous species. These facts infer the existence of haze in the atmospheres especially of warm, relatively low-density super-Earths and mini-Neptunes. Previous studies that addressed theoretical modeling of transmission spectra of hydrogen-dominated atmospheres with haze used some assumed distribution and size of haze particles. In this study, we model the atmospheric chemistry, derive the spatial and size distributions of haze particles by simulating the creation, growth and settling of hydrocarbon haze particles directly, and develop transmission spectrum models of UV-irradiated, solar-abundance atmospheres of close-in warm ($\sim$ 500 K) exoplanets. We find that the haze is distributed in the atmosphere much more broadly than previously assumed and consists of particles of various sizes. We also demonstrate that the observed diversity of transmission spectra can be explained by the difference in the production rate of haze monomers, which is related to the UV irradiation intensity from host stars.

The present infrared brightness of a planet originates partly from the accretion energy that the planet gained during its formation and hence provides important constraints to the planet formation process. A planet cools down from a hot initial state to the present state by losing energy through radiative emission from its atmosphere. Thus, the atmospheric properties affect the planetary cooling rate. Previous theories of giant planet cooling assume that the atmospheric composition is unchanged throughout evolution. Planet formation theories, however, suggest that the atmospheres especially of ice giants are rich in heavy elements in the early stages. Those heavy elements include condensable species such as water, ammonia, and methane, which are expected to have a great impact on atmospheric temperature and, thus, radiative emission through latent heat release. In this study we investigate the effect of such condensation on the planetary emission flux and quantify the impact on the cooling timescale. We then demonstrate that the latent heat of those species keeps the atmosphere hot and thus the emission flux high for billions of years, resulting in acceleration of the cooling of ice giants. This sheds light on the long-standing problem that Uranus is much less bright than theoretically predicted and is different in brightness from Neptune in spite of similarity in mass and radius. Also, we find that young ice giants with highly enriched atmospheres are much brighter in mid-infrared than those with unenriched atmospheres. This provides important implication for future direct imaging of extrasolar ice giants.

The Earth was born in violence. Many giant collisions of protoplanets are thought to have occurred during the terrestrial planet formation. Here we investigated the giant impact stage by using a hybrid code that consistently deals with the orbital evolution of protoplanets around the Sun and the details of processes during giant impacts between two protoplanets. A significant amount of materials (up to several tens of percent of the total mass of the protoplanets) is ejected by giant impacts. We call these ejected fragments the giant-impact fragments (GIFs). In some of the erosive hit-and-run and high-velocity collisions, metallic iron is also ejected, which comes from the colliding protoplanets' cores. From ten numerical simulations for the giant impact stage, we found that the mass fraction of metallic iron in GIFs ranges from ~ 1wt% to ~ 25wt%. We also discussed the effects of the GIFs on the dynamical and geochemical characteristics of formed terrestrial planets. We found that the GIFs have the potential to solve the following dynamical and geochemical conflicts: (1) The Earth, currently in a near circular orbit, is likely to have had a highly eccentric orbit during the giant impact stage. The GIFs are large enough in total mass to lower the eccentricity of the Earth to its current value via their dynamical friction. (2) The concentrations of highly siderophile elements (HSEs) in the Earth's mantle are greater than what was predicted experimentally. Re-accretion of the iron-bearing GIFs onto the Earth can contribute to the excess of HSEs. In addition, the estimated amount of iron-bearing GIFs provides significant reducing agent that could transform primitive CO2-H2O atmosphere and ocean into more reducing H2-bearing atmosphere. Thus, GIFs are important for the origin of Earth's life and its early evolution.

The Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT) is an optical three-band (g'_2-, r'_2-, and z_{s,2}-band) imager that was recently developed for the 188cm telescope at Okayama Astrophysical Observatory with the aim of validating and characterizing transiting planets. In a pilot observation with MuSCAT we observed a primary transit of HAT-P-14b, a high-surface gravity (g_p=38 ms^{-2}) hot Jupiter around a bright (V=10) F-type star. From a 2.9 hr observation, we achieved the five-minute binned photometric precisions of 0.028%, 0.022%, and 0.024% in the g'_2, r'_2, and z_{s,2} bands, respectively, which provided the highest-quality photometric data for this planet. Combining these results with those of previous observations, we search for variations of transit timing and duration over five years as well as variations of planet-star radius ratio (R_p/R_s) with wavelengths, but can find no considerable variation in any parameters. On the other hand, using the transit-subtracted light curves we simulate achievable measurement error of R_p/R_s with MuSCAT for various planetary sizes, assuming three types of host stars: HAT-P-14, the nearby K dwarf HAT-P-11, and the nearby M dwarf GJ1214. Comparing our results with the expected atmospheric scale heights, we find that MuSCAT is capable of probing the atmospheres of planets as small as a sub-Jupiter (R_p ~6 R_Earth) around HAT-P-14 in all bands, a Neptune (~4R_Earth) around HAT-P-11 in all bands, and a super-Earth (~2.5R_Earth) around GJ1214 in r'_2 and z_{s,2} bands. These results promise that MuSCAT will produce fruitful scientific outcomes in the K2 and TESS era.

Context. Within the core accretion scenario of planetary formation, most simulations performed so far always assume the accreting envelope to have a solar composition. From the study of meteorite showers on Earth and numerical simulations, we know that planetesimals must undergo thermal ablation and disruption when crossing a protoplanetary envelope. Once the protoplanet has acquired an atmosphere, the primordial envelope gets enriched in volatiles and silicates from the planetesimals. This change of envelope composition during the formation can have a significant effect in the final atmospheric composition and on the formation timescale of giant planets. Aims. To investigate the physical implications of considering the envelope enrichment of protoplanets due to the disruption of icy planetesimals during their way to the core. Particular focus is placed on the effect on the critical core mass for envelopes where condensation of water can occur. Methods. Internal structure models are numerically solved with the implementation of updated opacities for all ranges of metallicities and the software CEA to compute the equation of state. CEA computes the chemical equilibrium for an arbitrary mixture of gases and allows the condensation of some species, including water. This means that the latent heat of phase transitions is consistently incorporated in the total energy budget. Results. The critical core mass is found to decrease significantly when an enriched envelope composition is considered in the internal structure equations. A particular strong reduction of the critical core mass is obtained for planets whose envelope metallicity is larger than Z=0.45 when the outer boundary conditions are suitable for condensation of water to occur in the top layers of the atmosphere. We show that this effect is qualitatively preserved when the atmosphere is out of chemical equilibrium.

Motivated by recent detection of transiting high-density super-Earths, we explore the detectability of hot rocky super-Earths orbiting very close to their host stars. In the environment hot enough for their rocky surfaces to be molten, they would have the atmosphere composed of gas species from the magma oceans. In this study, we investigate the radiative properties of the atmosphere that is in the gas/melt equilibrium with the underlying magma ocean. Our equilibrium calculations yield Na, K, Fe, Si, SiO, O, and O$_2$ as the major atmospheric species. We compile the radiative-absorption line data of those species available in literature, and calculate their absorption opacities in the wavelength region of 0.1--100~$\mathrm{\mu m}$. Using them, we integrate the thermal structure of the atmosphere. Then, we find that thermal inversion occurs in the atmosphere because of the UV absorption by SiO. In addition, we calculate the ratio of the planetary to stellar emission fluxes during secondary eclipse, and find prominent emission features induced by SiO at 4~$\mathrm{\mu m}$ detectable by Spitzer, and those at 10 and 100~$\mathrm{\mu m}$ detectable by near-future space telescopes.

WASP-80b is a warm Jupiter transiting a bright late-K/early-M dwarf, providing a good opportunity to extend the atmospheric study of hot Jupiters toward the lower temperature regime. We report multi-band, multi-epoch transit observations of WASP-80b by using three ground-based telescopes covering from optical (g', Rc, and Ic bands) to near-infrared (NIR; J, H, and Ks bands) wavelengths. We observe 5 primary transits, each of which in 3 or 4 different bands simultaneously, obtaining 17 independent transit light curves. Combining them with results from previous works, we find that the observed transmission spectrum is largely consistent with both a solar abundance and thick cloud atmospheric models at 1.7$\sigma$ discrepancy level. On the other hand, we find a marginal spectral rise in optical region compared to the NIR region at 2.9$\sigma$ level, which possibly indicates the existence of haze in the atmosphere. We simulate theoretical transmission spectra for a solar abundance but hazy atmosphere, finding that a model with equilibrium temperature of 600 K can explain the observed data well, having a discrepancy level of 1.0$\sigma$. We also search for transit timing variations, but find no timing excess larger than 50 s from a linear ephemeris. In addition, we conduct 43 day long photometric monitoring of the host star in the optical bands, finding no significant variation in the stellar brightness. Combined with the fact that no spot-crossing event is observed in the five transits, our results confirm previous findings that the host star appears quiet for spot activities, despite the indications of strong chromospheric activities.

Recent progress in transit photometry opened a new window to the interior of super-Earths. From measured radii and masses, we can infer planetary internal compositions. It has been recently revealed that super-Earths are diverse in composition. Such a diversity is thought to arise from diversity in volatile content. The stability of the volatile components is to be examined, because hot super-Earths undergo photo-evaporative mass loss. While several studies investigated the impact of photo-evaporative mass loss on hydrogen-helium envelopes, there are few studies as to the impact on water-vapor envelopes. To obtain theoretical prediction to future observations, we also investigate the relationships among masses, radii, and semimajor axes of water-rich sub/super-Earths that have undergone photo-evaporative mass loss. We simulate the interior structure and evolution of sub/super-Earths that consist of a rocky core surrounded by a water envelope, including mass loss due to the stellar XUV-driven energy-limited hydrodynamic escape. We find that the photo-evaporative mass loss has a significant impact on the evolution of hot sub/super-Earths. We then derive the threshold planetary mass and radius below which the planet loses its water envelope completely as a function of the initial water content, and find that there are minimums of the threshold mass and radius. We constrain the domain in the parameter space of planetary mass, radius, and semimajor axis in which sub/super-Earths never retain water envelopes in 1-10 Gyr. This would provide an essential piece of information for understanding the origin of close-in low-mass planets. The current uncertainties in stellar XUV flux and its heating efficiency, however, prevent us from deriving robust conclusions. Nevertheless, it seems to be a robust conclusion that Kepler planet candidates contain a significant number of rocky sub/super-Earths.

Violent variation of transit depths and an ingress-egress asymmetry of the transit light curve discovered in KIC 12557548 have been interpreted as evidences of a catastrophic evaporation of atmosphere with dust (M_p gtrsim 1 M_oplus/Gyr) from a close-in small planet. To explore what drives the anomalous atmospheric escape, we perform time-series analysis of the transit depth variation of Kepler archival data for ~ 3.5 yr. We find a ~ 30% periodic variation of the transit depth with P1 = 22.83 pm 0.21 days, which is within the error of the rotation period of the host star estimated using the light curve modulation, Prot = 22.91 pm 0.24 days. We interpret the results as evidence that the atmospheric escape of KIC 12557548b correlates with stellar activity. We consider possible scenarios that account for both the mass loss rate and the correlation with stellar activity. X-ray and ultraviolet (XUV)-driven evaporation is possible if one accepts a relatively high XUV flux and a high efficiency for converting the input energy to the kinetic energy of the atmosphere. Star-planet magnetic interaction is another possible scenario though huge uncertainty remains for the mass loss rate.

We present 5 new transit light curves of GJ 1214b taken in BJHKs-bands. Two transits were observed in B-band using the Suprime-Cam and the FOCAS instruments onboard the Subaru 8.2m telescope, and one transit was done in JHKs-bands simultaneously with the SIRIUS camera on the IRSF 1.4m telescope. MCMC analyses show that the planet-to-star radius ratios are, Rp/Rs = 0.11651 \pm 0.00065 (B-band, Subaru/Suprime-Cam), Rp/Rs = 0.11601 \pm 0.00117 (B-band, Subaru/FOCAS), Rp/Rs = 0.11654 \pm 0.00080 (J-band, IRSF/SIRIUS), Rp/Rs = 0.11550 ^{+0.00142}_{-0.00153} (H-band, IRSF/SIRIUS), and Rp/Rs = 0.11547 \pm 0.00127 (Ks-band, IRSF/SIRIUS). The Subaru Suprime-Cam transit photometry shows a possible spot-crossing feature. Comparisons of the new transit depths and those from previous studies with the theoretical models by Howe & Burrows (2012) suggest that the high molecular weight atmosphere (e.g., 1% H$_2$O + 99% N$_2$) models are most likely, however, the low molecular weight (hydrogen dominated) atmospheres with extensive clouds are still not excluded. We also report a long-term monitoring of the stellar brightness variability of GJ 1214 observed with the MITSuME 50cm telescope in g'-, Rc-, and Ic-bands simultaneously. The monitoring was conducted for 32 nights spanning 78 nights in 2012, and we find a periodic brightness variation with a period of Ps = 44.3 \pm 1.2 days and semi-amplitudes of 2.1% \pm 0.4% in g'-band, 0.56% \pm 0.08% in Rc-band, and 0.32% \pm 0.04% in Ic-band.

We present optical (g', R_c, and I_c) to near-infrared (J) simultaneous photometric observations for a primary transit of GJ3470b, a Uranus-mass transiting planet around a nearby M dwarf, by using the 50-cm MITSuME telescope and the 188-cm telescope, both at Okayama Astrophysical Observatory. From these data, we derive the planetary mass, radius, and density as 14.1 \pm 1.3 M_earth, 4.32^{+0.21}_{-0.10} R_earth, and 0.94 \pm 0.12 g cm^{-3}, respectively, thus confirming the low density that was reported by Demory et al. based on the Spitzer/IRAC 4.5-micron photometry (0.72^{+0.13}_{-0.12} g cm^{-3}). Although the planetary radius is about 10% smaller than that reported by Demory et al., this difference does not alter their conclusion that the planet possesses a hydrogen-rich envelope whose mass is approximately 10% of the planetary total mass. On the other hand, we find that the planet-to-star radius ratio (R_p/R_s) in the J band (0.07577^{+0.00072}_{-0.00075}) is smaller than that in the I_c (0.0802 \pm 0.0013) and 4.5-micron (0.07806^{+0.00052}_{-0.00054}) bands by 5.9% \pm 2.0% and 3.0% \pm 1.2%, respectively. A plausible explanation for the differences is that the planetary atmospheric opacity varies with wavelength due to absorption and/or scattering by atmospheric molecules. Although the significance of the observed R_p/R_s variations is low, if confirmed, this fact would suggest that GJ3470b does not have a thick cloud layer in the atmosphere. This property would offer a wealth of opportunity for future transmission-spectroscopic observations of this planet to search for certain molecular features, such as H2O, CH4, and CO, without being prevented by clouds.

We report high precision transit photometry of GJ1214b in JHKs bands taken simultaneously with the SIRIUS camera on the IRSF 1.4m telescope at Sutherland, South Africa. Our MCMC analyses show that the observed planet-to-star radius ratios in JHKs bands are R_{\rm p}/R_{\rm s,J} = 0.11833 \pm 0.00077, R_{\rm p}/R_{\rm s,H} = 0.11522 \pm 0.00079, R_{\rm p}/R_{\rm s,Ks} = 0.11459 \pm 0.00099, respectively. The radius ratios are well consistent with the previous studies by Bean et al. (2011) within 1\sigma, while our result in Ks band is shallower than and inconsistent at 4\sigma\ level with the previous measurements in the same band by Croll et al. (2011). We have no good explanation for this discrepancy at this point. Our overall results support a flat transmission spectrum in the observed bands, which can be explained by a water-dominated atmosphere or an atmosphere with extensive high-altitude clouds or haze. To solve the discrepancy of the radius ratios and to discriminate a definitive atmosphere model for GJ1214b in the future, further transit observations around Ks band would be especially important.

Motivated by recent discoveries of low-density super-Earths with short orbital periods, we have investigated in-situ accretion of H-He atmospheres on rocky bodies embedded in dissipating warm disks, by simulating quasi-static evolution of atmospheres that connect to the ambient disk. We have found that the atmospheric evolution has two distinctly different outcomes, depending on the rocky body's mass: While the atmospheres on massive rocky bodies undergo runaway disk-gas accretion, those on light rocky bodies undergo significant erosion during disk dispersal. In the atmospheric erosion, the heat content of the rocky body that was previously neglected plays an important role. We have also realized that the atmospheric mass is rather sensitive to disk temperature in the mass range of interest in this study. Our theory is applied to recently-detected super-Earths orbiting Kepler-11 to examine the possibility that the planets are rock-dominated ones with relatively thick H-He atmospheres. The application suggests that the in-situ formation of the relatively thick H-He atmospheres inferred by structure modeling is possible only under restricted conditions; namely, relatively slow disk dissipation and/or cool environments. This study demonstrates that low-density super-Earths provide important clues to understanding of planetary accretion and disk evolution.

We have investigated how envelope pollution by icy planetesimals affects the critical core mass for gas giant formation and the gas accretion time-scales. In the core-accretion model, runaway gas accretion is triggered after a core reaches a critical core mass. All the previous studies on the core-accretion model assumed that the envelope has the solar composition uniformly. In fact, the envelope is likely polluted by evaporated materials of icy planetesimals because icy planetesimals going through the envelope experience mass loss via strong ablation and most of their masses are deposited in the deep envelope. In this paper, we have demonstrated that envelope pollution in general lowers the critical core masses and hastens gas accretion on to the protoplanet because of the increase in the molecular weight and reduction of adiabatic temperature gradient. Widely- and highly-polluted envelopes allow smaller cores to form massive envelopes before disc dissipation. Our results suggest that envelope pollution in the course of planetary accretion has the potential to trigger gas giant formation with small cores. We propose that it is necessary to take into account envelope pollution by icy planetesimals when we discuss gas giant formation based on the core-accretion model.

Recent radial-velocity surveys for GK clump giants have revealed that planets also exist around ~1.5-3 Msun stars. However, no planets have been found inside 0.6 AU around clump giants, in contrast to solar-type main-sequence stars, many of which harbor short-period planets such as hot Jupiters. In this study we examine the possibility that planets were engulfed by host stars evolving on the red-giant branch (RGB). We integrate the orbital evolution of planets in the RGB and helium burning (HeB) phases of host stars, including the effects of stellar tide and stellar mass loss. Then we derive the critical semimajor axis (or the survival limit) inside which planets are eventually engulfed by their host stars after tidal decay of their orbits. Especially, we investigate the impact of stellar mass and other stellar parameters on the survival limit in more detail than previous studies. In addition, we make detailed comparison with measured semimajor axes of planets detected so far, which no previous study did. We find that the critical semimajor axis is quite sensitive to stellar mass in the range between 1.7 and 2.1 Msun, which suggests a need for careful comparison between theoretical and observational limits of existence of planets. Our comparison demonstrates that all those planets are beyond the survival limit, which is consistent with the planet-engulfment hypothesis. However, on the high-mass side (> 2.1 Msun), the detected planets are orbiting significantly far from the survival limit, which suggests that engulfment by host stars may not be the main reason for the observed lack of short-period giant planets. To confirm our conclusion, the detection of more planets around clump giants, especially with masses > 2.5 Msun, is required.

Around low- and intermediate-mass (1.5-3 M_sun) red giants, no planets have been found inside 0.6 AU. Such a paucity is not seen in the case of 1 M_sun main sequence stars. In this study, we examine the possibility that short-period planets were engulfed by their host star evolving off the main sequence. To do so, we have simulated the orbital evolution of planets, including the effects of stellar tide and mass loss, to determine the critical semimajor axis, a_crit, beyond which planets survive the RGB expansion of their host star. We have found that a_crit changes drastically around 2 M_sun: In the lower-mass range, a_crit is more than 1 AU, while a_crit is as small as about 0.2 AU in the higher-mass range. Comparison with measured semimajor axes of known planets suggests that there is a lack of planets that only planet engulfment never accounts for in the higher-mass range. Whether the lack is real affects our understanding of planet formation. Therefore, increasing the number of planet samples around evolved intermediate-mass stars is quite meaningful to confirm robustness of the lack of planets.

The principal Hugoniot for liquid hydrogen was obtained up to 55 GPa under laser-driven shock loading. Pressure and density of compressed hydrogen were determined by impedance-matching to a quartz standard. The shock temperature was independently measured from the brightness of the shock front. Hugoniot data of hydrogen provide a good benchmark to modern theories of condensed matter. The initial number density of liquid hydrogen is lower than that for liquid deuterium, and this results in shock compressed hydrogen having a higher compression and higher temperature than deuterium at the same shock pressure.

We compute grids of radiative-convective model atmospheres for Jupiter, Saturn, Uranus, and Neptune over a range of intrinsic fluxes and surface gravities. The atmosphere grids serve as an upper boundary condition for models of the thermal evolution of the planets. Unlike previous work, we customize these grids for the specific properties of each planet, including the appropriate chemical abundances and incident fluxes as a function of solar system age. Using these grids, we compute new models of the thermal evolution of the major planets in an attempt to match their measured luminosities at their known ages. Compared to previous work, we find longer cooling times, predominantly due to higher atmospheric opacity at young ages. For all planets, we employ simple "standard" cooling models that feature adiabatic temperature gradients in the interior H/He and water layers, and an initially hot starting point for the calculation of subsequent cooling. For Jupiter we find a model cooling age 10% longer than previous work, a modest quantitative difference. This may indicate that the hydrogen equation of state used here overestimates the temperatures in the deep interior of the planet. For Saturn we find a model cooling age 20% longer than previous work. However, an additional energy source, such as that due to helium phase separation, is still clearly needed. For Neptune, unlike in work from the 1980s and 1990s, we match the measured Teff of the planet with a model that also matches the planet's current gravity field. This is predominantly due to advances in the equation of state of water. This may indicate that the planet possesses no barriers to efficient convection in its deep interior. However, for Uranus, our models exacerbate the well-known problem that Uranus is far cooler than calculations predict, which could imply strong barriers to interior convective cooling.

We investigate the critical core mass and the envelope growth timescale, assuming grain-free envelopes, to examine how small cores are allowed to form gas giants in the framework of the core accretion model. This is motivated by a theoretical dilemma concerning Jupiter formation: Modelings of Jupiter's interior suggest that it contains a small core of < 10 Earth mass, while many core accretion models of Jupiter formation require a large core of > 10 Earth mass to finish its formation by the time of disk dissipation. Reduction of opacity in the accreting envelope is known to hasten gas giant formation. Almost all the previous studies assumed grain-dominated opacity in the envelope. Instead, we examine cases of grain-free envelopes in this study. Our numerical simulations show that an isolated core of as small as 1.7 Earth mass is able to capture disk gas to form a gas giant on a timescale of million years, if the accreting envelope is grain-free; that value decreases to 0.75 Earth mass, if the envelope is metal-free, namely, composed purely of hydrogen and helium. It is also shown that alkali atoms, which are known to be one of the dominant opacity sources near 1500 K in the atmospheres of hot Jupiters, have little contribution to determine the critical core mass. Our results confirm that sedimentation and coagulation of grains in the accreting envelope is a key to resolve the dilemma about Jupiter formation.

The discovery of CoRoT-7b, a planet of radius 1.68 +/- 0.09 R_E, mass 4.8 +/- 0.8 M_E and orbital period of 0.854 days demonstrates that small planets can orbit extremely close to their star. We use knowledge of hot Jupiters, mass loss estimates and models for the interior structure and evolution of planets to understand its composition, structure and evolution. The inferred mass and radius of CoRoT-7b are consistent with a rocky planet that would be depleted in iron relative to Earth. However, a one sigma increase in mass (5.6 M_E) and decrease in size (1.59 R_E) are compatible with an Earth-like composition (33% iron, 67% silicates). Alternatively, it is possible that CoRoT-7b contains a significant amount of volatiles. An equally good fit to the data is found for a vapor envelope equal to 3% (and up to 10%) by mass above an Earth-like nucleus. Because of its intense irradiation and small size, the planet cannot possess an envelope of H and He of more than 1/10,000 of its total mass. We show that the mass loss is significant (~ 10^11 g/s) and independent of planetary composition. This is because the hydrodynamical escape rate is independent of the atmosphere's mean molecular mass, and owing to the intense irradiation, even a bare rocky planet would possess an equilibrium vapor atmosphere thick enough to capture stellar UV photons. This escape rate rules out the possibility of a H-He envelope as it would escape in only ~1 Ma. A water vapor atmosphere would escape in ~ 1 Ga, and thus it is a plausible scenario. The origin of CoRoT-7b cannot be inferred from present observations: It may have formed rocky; or be the remnant of a Uranus-like ice giant, or a gas giant with a small core that was stripped of its gaseous envelope.

We report the detection of 3 new extrasolar planets from the precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. The host stars, namely, 18 Del (G6 III), xi Aql (K0 III) and HD 81688 (K0 III-IV), are located at the clump region on the HR diagram with estimated masses of 2.1-2.3 M_solar. 18 Del b has a minimum mass of 10.3 M_Jup and resides in a nearly circular orbit with period of 993 days, which is the longest one ever discovered around evolved stars. xi Aql b and HD 81688 b have minimum masses of 2.8 and 2.7 M_Jup, and reside in nearly circular orbits with periods of 137 and 184 days, respectively, which are the shortest ones among planets around evolved stars. All of the substellar companions ever discovered around possible intermediate-mass (1.7-3.9 M_solar) clump giants have semimajor axes larger than 0.68 AU, suggesting the lack of short-period planets. Our numerical calculations suggest that Jupiter-mass planets within about 0.5 AU (even up to 1 AU depending on the metallicity and adopted models) around 2-3 M_solar stars could be engulfed by the central stars at the tip of RGB due to tidal torque from the central stars. Assuming that most of the clump giants are post-RGB stars, we can not distinguish whether the lack of short-period planets is primordial or due to engulfment by central stars. Deriving reliable mass and evolutionary status for evolved stars is highly required for further investigation of formation and evolution of planetary systems around intermediate-mass stars.

The origin of the Earth's ocean has been discussed on the basis of deuterium/hydrogen ratios (D/H) of several sources of water in the solar system. The average D/H of carbonaceous chondrites (CC's) is known to be close to the current D/H of the Earth's ocean, while those of comets and the solar nebula are larger by about a factor of two and smaller by about a factor of seven, respectively, than that of the Earth's ocean. Thus, the main source of the Earth's ocean has been thought to be CC's or adequate mixing of comets and the solar nebula. However, those conclusions are correct only if D/H of water on the Earth has remained unchanged for the past 4.5 Gyr. In this paper, we investigate evolution of D/H in the ocean in the case that the early Earth had a hydrogen-rich atmosphere, the existence of which is predicted by recent theories of planet formation no matter whether the nebula remains or not. Then we show that D/H in the ocean increases by a factor of 2-9, which is caused by the mass fractionation during atmospheric hydrogen loss, followed by deuterium exchange between hydrogen gas and water vapor during ocean formation. This result suggests that the apparent similarity in D/H of water between CC's and the current Earth's ocean does not necessarily support the CC's origin of water and that the apparent discrepancy in D/H is not a good reason for excluding the nebular origin of water.

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.

The high density of the close-in extrasolar planet HD149026b suggests the presence of a huge core in the planet, which challenges planet formation theory. We first derive constraints on the amount of heavy elements and hydrogen/helium present in the planet: We find that preferred values of the core mass are between 50 and 80 M_E. We then investigate the possibility of subcritical core accretion as envisioned for Uranus and Neptune and find that the subcritical accretion scenario is unlikely in the case of HD149026b for at least two reasons: (i) Subcritical planets are such that the ratio of their core mass to their total mass is above ~0.7, in contradiction with constraints for all but the most extreme interior models of HD149026b; (ii) High accretion rates and large isolation mass required for the formation of a subcritical core of 30 M_E are possible only at specific orbital distances in a disk with a surface density of dust equal to at least 10 times that of the minimum mass solar nebula. This value climbs to 30 when considering a 50 M_E core. These facts point toward two main routes for the formation of this planet: (i) Gas accretion that is limited by a slow viscous inflow of gas in an evaporating disk; (ii) A significant modification of the composition of the planet after as accretion has stopped. These two routes are not mutually exclusive. Illustrating the second route, we show that for a wide range of impact parameters, giant impacts lead to a loss of the gas component of the planet and thus may lead to planets that are highly enriched in heavy elements. In the giant impact scenario, we expect an outer giant planet to be present. Observational studies by imaging, astrometry and long term interferometry of this system are needed to better narrow down the ensemble of possibilities.

From an astrobiological point of view, special attention has been paid to the probability of habitable planets in extrasolar systems. The purpose of this study is to constrain a possible range of the mass of a terrestrial planet that can get water. We focus on the process of water production through oxidation of the atmospheric hydrogen--the nebular gas having been attracted gravitationally--by oxide available at the planetary surface. For the water production to work well on a planet, a sufficient amount of hydrogen and enough high temperature to melt the planetary surface are needed. We have simulated the structure of the atmosphere that connects with the protoplanetary nebula for wide ranges of heat flux, opacity, and density of the nebular gas. We have found both requirements are fulfilled for an Earth-mass planet for wide ranges of the parameters. We have also found the surface temperature of planets of <= 0.3 Earth masses is lower than the melting temperature of silicate (~ 1500K). On the other hand, a planet of more than several Earth masses becomes a gas giant planet through runaway accretion of the nebular gas.