A substantial number of super-Earths have been discovered, and atmospheres of transiting super-Earths have also been observed by transmission spectroscopy. Several lines of observational evidence indicate that most super-Earths do not possess massive H$_2$/He atmospheres. However, accretion and retention of less massive atmospheres on super-Earths challenge planet formation theory. We consider the following three mechanisms: (i) envelope heating by pebble accretion, (ii) mass loss during giant impacts, and (iii) atmospheric loss by stellar X-ray and EUV photoevaporation. We investigate whether these mechanisms influence the amount of the atmospheres that form around super-Earths. We develop a code combining an N-body simulation of pebble-driven planetary formation and an atmospheric evolution simulation. We demonstrate that the observed orbital properties of super-Earths are well reproduced by the results of our simulations. However, (i) heating by pebble accretion ceases prior to disk dispersal, (ii) the frequency of giant impact events is too low to sculpt massive atmospheres, and (iii) many super-Earths having H$_2$/He atmospheres of $\gtrsim 10$ wt% survive against stellar irradiations for 1 Gyr. Therefore, it is likely that other mechanisms such as suppression of gas accretion are required to explain less massive atmospheres ($\lesssim 10$ wt%) of super-Earths.

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.

Recently, transmission spectroscopy in the atmospheres of the TRAPPIST-1 planets revealed flat and featureless absorption spectra, which rule out cloud-free hydrogen-dominated atmospheres. Earth-sized planets orbiting TRAPPIST-1 likely have either a clear or a cloudy/hazy hydrogen-poor atmosphere. In this paper, we investigate whether a proposed formation scenario is consistent with expected atmospheric compositions of the TRAPPIST-1 planets. We examine the amount of a hydrogen-rich gas that TRAPPIST-1-like planets accreted from the ambient disk until disk dispersal. Since TRAPPIST-1 planets are trapped into a resonant chain, we simulate disk gas accretion onto a migrating TRAPPIST-1-like planet. We find that the amount of an accreted hydrogen-rich gas is as small as 10$^{-2}$ wt% and 0.1 wt% for TRAPPIST-1b and 1c, 10$^{-2}$ wt% for 1d, 1 wt% for 1e, a few wt% for 1f and 1g and 1 wt% for 1h, respectively. We also calculate a long-term thermal evolution of TRAPPIST-1-like planets after disk dissipation and estimate the mass loss of their hydrogen-rich atmospheres driven by a stellar X-ray and UV irradiation. We find that all the accreted hydrogen-rich atmospheres can be lost via hydrodynamic escape. Therefore, we conclude that TRAPPIST-1 planets should have no primordial hydrogen-rich gases but secondary atmospheres such as a Venus-like one and water vapor, if they currently retain atmospheres.

We present $\kappa$ Andromeda b's photometry and astrometry taken with Subaru/SCExAO+HiCIAO and Keck/NIRC2, combined with recently published SCExAO/CHARIS low-resolution spectroscopy and published thermal infrared photometry to further constrain the companion's atmospheric properties and orbit. $\kappa$ And b's Y/Y-K colors are redder than field dwarfs, consistent with its youth and lower gravity. Empirical comparisons of its Y-band photometry and CHARIS spectrum to a large spectral library of isolated field dwarfs reaffirm the conclusion from Currie et al. (2018) that it likely has a low gravity but admit a wider range of most plausible spectral types (L0-L2). Our gravitational classification also suggests that the best-fit objects for $\kappa$ And b may have lower gravity than those previously reported. Atmospheric models lacking dust/clouds fail to reproduce its entire 1--4.7 $\mu m$ spectral energy distribution, cloudy atmosphere models with temperatures of $\sim$ 1700--2000 $K$ better match $\kappa$ And b data. Most well-fitting model comparisons favor 1700--1900 $K$, a surface gravity of log(g) $\sim$ 4--4.5, and a radius of 1.3--1.6\,$R_{\rm Jup}$; the best-fit model (DRIFT-Phoenix) yields the coolest and lowest-gravity values: $T_{\rm eff}$=1700 K and $\log g$=4.0. An update to $\kappa$ And b's orbit with ExoSOFT using new astrometry spanning seven years reaffirms its high eccentricity ($0.77\pm0.08$). We consider a scenario where unseen companions are responsible for scattering $\kappa$ And b to a wide separation and high eccentricity. If three planets, including $\kappa$ And b, were born with coplanar orbits and one of them was ejected by gravitational scattering, a potential inner companion with mass $\gtrsim10M_{\rm Jup}$ could be located at $\lesssim$ 25 au.

Planets occur most frequently around cool dwarfs, but only a handful of specific examples are known to orbit the latest-type M stars. Using TESS photometry, we report the discovery of two planets transiting the low-mass star called LP 791-18 (identified by TESS as TOI 736). This star has spectral type M6V, effective temperature 2960 K, and radius 0.17 R_Sun, making it the third-coolest star known to host planets. The two planets straddle the radius gap seen for smaller exoplanets; they include a 1.1 R_Earth planet on a 0.95 day orbit and a 2.3 R_Earth planet on a 5 day orbit. Because the host star is small the loss of light during these planets' transits is fairly large (0.4% and 1.7%). This has allowed us to detect both planets' transits from ground-based photometry, refining their radii and orbital ephemerides. In the future, radial velocity observations and transmission spectroscopy can both probe these planets' bulk interior and atmospheric compositions, and additional photometric monitoring would be sensitive to even smaller transiting planets.

In clustered environments, stellar encounters can liberate planets from their host stars via close encounters. Although the detection probability of planets suggests that the planet population in open clusters resembles that in the field, only a few dozen planet-hosting stars have been discovered in open clusters. We explore the survival rates of planets against stellar encounters in open clusters similar to the Pleiades, Hyades, and Praesepe and embedded clusters. We performed a series of N-body simulations of high-density and low-density open clusters, open clusters that grow via mergers of subclusters, and embedded clusters. We semi-analytically calculated the survival rate of planets in star clusters up to 1Gyr using relative velocities, masses, and impact parameters of intruding stars. Less than 1.5% of close-in planets within 1 AU and at most 7% of planets with 1-10 AU are ejected by stellar encounters in clustered environments after the dynamical evolution of star clusters. If a planet population from 0.01-100 AU in an open cluster initially follows the probability distribution function of exoplanets with semi-major axis ($a_p$) between 0.03-3 AU in the field discovered by RV surveys, the PDF of surviving planets beyond ~10 AU in open clusters can be slightly modified to $\propto a_p^{-0.76}$. The production rate of free-floating planets (FFPs) per star is 0.0096-0.18, where we have assumed that all the stars initially have one giant planet with a mass of 1--13 MJ in a circular orbit. The expected frequency of FFPs is compatible with the upper limit on that of FFPs indicated by recent microlensing surveys. Our survival rates of planets in open clusters suggest that planets within 10 AU around FGKM-type stars are rich in relatively-young (<~10-100 Myr for open clusters and ~1-10 Myr for embedded clusters), less massive open clusters, which are promising targets for planet searches.

We have newly developed a Parallelized Particle-Particle Particle-tree code for Planet formation, PENTACLE, which is a parallelized hybrid $N$-body integrator executed on a CPU-based (super)computer. PENTACLE uses a 4th-order Hermite algorithm to calculate gravitational interactions between particles within a cutoff radius and a Barnes-Hut tree method for gravity from particles beyond. It also implements an open-source library designed for full automatic parallelization of particle simulations, FDPS (Framework for Developing Particle Simulator) to parallelize a Barnes-Hut tree algorithm for a memory-distributed supercomputer. These allow us to handle $1-10$ million particles in a high-resolution $N$-body simulation on CPU clusters for collisional dynamics, including physical collisions in a planetesimal disc. In this paper, we show the performance and the accuracy of PENTACLE in terms of $\tilde{R}_{\rm cut}$ and a time-step $\Delta t$. It turns out that the accuracy of a hybrid $N$-body simulation is controlled through $\Delta t / \tilde{R}_{\rm cut}$ and $\Delta t / \tilde{R}_{\rm cut} \sim 0.1$ is necessary to simulate accurately accretion process of a planet for $\geq 10^6$ years. For all those who interested in large-scale particle simulations, PENTACLE customized for planet formation will be freely available from https://github.com/PENTACLE-Team/PENTACLE under the MIT lisence.

The core-accretion model predicts that planetary cores as massive as super-Earths undergo runaway gas accretion to become gas giants. However, the exoplanet census revealed the prevalence of super-Earths close to their host stars, which should have avoided runaway gas accretion. In fact, mass-radius relationships of transiting planets suggest that some close-in super-Earths possess H$_2$/He atmospheres of ~ 0.1-10% by mass. Previous studies indicated that properties of a disk gas such as metallicity and the inflow/outflow cycle of a disk gas around a super-Earth can regulate accumulation of a H$_2$/He atmosphere onto itself. In this paper, we propose a new mechanism that radial mass accretion in a disk can limit the gas accretion onto super-Earth cores. Recent magneto-hydrodynamic simulations found that magnetically driven disk winds can drive a rapid gas flow near the disk surface. Such a rapid gas flow may slip out of a planetary core and regulate gas supply to an accreting gas onto the core. We performed N-body simulations for formation of super-Earths with accretion of atmospheres in a viscous accretion disk including effects of wind-driven accretion. We found that even super-Earth cores can avoid triggering runaway gas accretion if the inflow of a disk gas toward the cores is limited by viscous accretion. Our model predicts that super-Earths having H$_2$/He atmosphere of ~ 0.1-10 wt % form within $\lesssim$ 1 au of the central star, whereas gas giants are born in the outer region. This mechanism can explain the radial dependence of observed giant planets beyond the solar system.

Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves, but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on their context. Characterization of temperate Earth-size planets in the coming years will focus on those around nearby late-type stars. JWST and later 30 meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the WFIRST mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables.

We report the photometry of six transits of the hot Jupiter HAT-P-29b obtained from 2013 October to 2015 January. We analyze the new light curves, in combination with the published photometric, and Doppler velocimetric, and spectroscopic measurements, finding an updated orbital ephemeris for the HAT-P-29 system, $T_{\rm C}[0]= 2456170.5494(15)\,[\rm{BJD_{TDB}}]$ and $P=5.723390(13) \,{\rm days}$. It is $17.63\,{\rm s}$ ($4.0\,\sigma$) longer than the previously published value, amounting to errors exceeding $2.5\,\rm{hrs}$ at the time of writing (on UTC 2018 June 1). The measured transit mid-times for HAT-P-29b show no compelling evidence of timing anomalies from a linear model, which rules out the presence of a perturbers with masses greater than $0.6$, $0.7$, $0.5$, and $0.4\,{\rm M_\oplus}$ near the $1:2$, $2:3$, $3:2$, and $2:1$ resonances with HAT-P-29b, respectively.

We report a measurement of the B-mode polarization power spectrum in the cosmic microwave background (CMB) using the POLARBEAR experiment in Chile. The faint B-mode polarization signature carries information about the Universe's entire history of gravitational structure formation, and the cosmic inflation that may have occurred in the very early Universe. Our measurement covers the angular multipole range 500 < l < 2100 and is based on observations of an effective sky area of 25 square degrees with 3.5 arcmin resolution at 150 GHz. On these angular scales, gravitational lensing of the CMB by intervening structure in the Universe is expected to be the dominant source of B-mode polarization. Including both systematic and statistical uncertainties, the hypothesis of no B-mode polarization power from gravitational lensing is rejected at 97.1% confidence. The band powers are consistent with the standard cosmological model. Fitting a single lensing amplitude parameter A_BB to the measured band powers, A_BB = 1.12 +/- 0.61 (stat) +0.04/-0.12 (sys) +/- 0.07 (multi), where A_BB = 1 is the fiducial WMAP-9 LCDM value. In this expression, "stat" refers to the statistical uncertainty, "sys" to the systematic uncertainty associated with possible biases from the instrument and astrophysical foregrounds, and "multi" to the calibration uncertainties that have a multiplicative effect on the measured amplitude A_BB.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

We search for signatures of a distant planet around the two-million-year-old classical T-Tauri star CI Tau hosting a hot Jupiter candidate (M_{p} sin{i} ~ 8.1 M_{Jupiter}) in an eccentric orbit (e ~0.3). To probe the existence of an outer perturber, we reanalyzed 1.3 mm dust continuum observations of the protoplanetary disk around CI Tau obtained by the Atacama Large Millimeter/submillimeter Array. We found a gap structure at ~0".8 in CI Tau's disk. Our visibility fitting assuming an axisymmetric surface brightness profile suggested that the gap is located at a deprojected radius of 104.5+/-1.6 au and has a width of 36.9+/-2.9 au. The brightness temperature around the gap was calculated to be ~2.3 K lower than that of the ambient disk. Gap-opening mechanisms such as secular gravitational instability and dust trapping can explain the gap morphology in the CI Tau disk. The scenario that an unseen planet created the observed gap structure cannot be ruled out, although the coexistence of an eccentric hot Jupiter and a distant planet around the young CI Tau would be challenging for gravitational scattering scenarios. The mass of the planet was estimated to be between ~0.25 M_{Jupiter} and ~0.8 M_{Jupiter} from the gap width and depth (0.41 +0.04/-0.06) in the modeled surface brightness image, which is lower than the current detection limits of high-contrast direct imaging. The young classical T-Tauri CI Tau may be a unique system to explore the existence of a potential distant planet as well as the origin of an eccentric hot Jupiter.

For the first time in human history, we will soon be able to apply the scientific method to the question "Are We Alone?" The rapid advance of exoplanet discovery, planetary systems science, and telescope technology will soon allow scientists to search for life beyond our Solar System through direct observation of extrasolar planets. This endeavor will occur alongside searches for habitable environments and signs of life within our Solar System. While the searches are thematically related and will inform each other, they will require separate observational techniques. The search for life on exoplanets holds potential through the great diversity of worlds to be explored beyond our Solar System. However, there are also unique challenges related to the relatively limited data this search will obtain on any individual world. This white paper reviews the scientific community's ability to use data from future telescopes to search for life on exoplanets. This material summarizes products from the Exoplanet Biosignatures Workshop Without Walls (EBWWW). The EBWWW was constituted by a series of online and in person activities, with participation from the international exoplanet and astrobiology communities, to assess state of the science and future research needs for the remote detection of life on planets outside our Solar System.

We present ten $R$-band photometric observations of eight different transits of the hot Jupiter HAT-P-33b, which has been targeted by our Transiting Exoplanet Monitoring Project (TEMP). The data were obtained by two telescopes at the Xinglong Station of National Astronomical Observatories of China (NAOC) from 2013 December through 2016 January, and exhibit photometric scatter of $1.6-3.0\,\rm{mmag}$. After jointly analyzing the previously published photometric data, radial-velocity (RV) measurements, and our new light curves, we revisit the system parameters and orbital ephemeris for the HAT-P-33b system. Our results are consistent with the published values except for the planet-to-star radius ratio ($R_{P}/R_{*}$), the ingress/egress duration ($\tau$) and the total duration ($T_{14}$), which together indicate a slightly shallower and shorter transit shape. Our results are based on more complete light curves, whereas the previously published work had only one complete transit light curve. No significant anomalies in Transit Timing Variations (TTVs) are found, and we place upper mass limits on potential perturbers, largely supplanting the loose constraints provided by the extant RV data. The TTV limits are stronger near mean-motion resonances, especially for the low-order commensurabilities. We can exclude the existence of a perturber with mass larger than 0.6, 0.3, 0.5, 0.5, and $0.3\,{\rm M_\oplus}$ near the 1:3, 1:2, 2:3, 3:2, and 2:1 resonances, respectively.

We report on the confirmation that the candidate transits observed for the star EPIC 211525389 are due to a short-period Neptune-sized planet. The host star, located in K2 campaign field 5, is a metal-rich ([Fe/H] = 0.26$\pm$0.05) G-dwarf (T_eff = 5430$\pm$70 K and log g = 4.48$\pm$0.09), based on observations with the High Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. High-spatial resolution AO imaging with HiCIAO on the Subaru telescope excludes faint companions near the host star, and the false positive probability of this target is found to be <$10^{-6}$ using the open source vespa code. A joint analysis of transit light curves from K2 and additional ground-based multi-color transit photometry with MuSCAT on the Okayama 1.88m telescope gives the orbital period of P = 8.266902$\pm$0.000070 days and consistent transit depths of $R_p/R_\star \sim 0.035$ or $(R_p/R_\star)^2 \sim 0.0012$. The transit depth corresponds to a planetary radius of $R_p = 3.59_{-0.39}^{+0.44} R_{\oplus}$, indicating that EPIC 211525389 b is a short-period Neptune-sized planet. Radial velocities of the host star, obtained with the Subaru HDS, lead to a 3\sigma\ upper limit of 90 $M_{\oplus} (0.00027 M_{\odot})$ on the mass of EPIC 211525389 b, confirming its planetary nature. We expect this planet, newly named K2-105 b, to be the subject of future studies to characterize its mass, atmosphere, spin-orbit (mis)alignment, as well as investigate the possibility of additional planets in the system.

We have developed a novel two-layer anti-reflection (AR) coating method for large-diameter infrared (IR) filters made of alumina, for the use at cryogenic temperatures in millimeter wave measurements. Thermally- sprayed mullite and polyimide foam (Skybond Foam) are used as the AR material. An advantage of the Skybond Foam is that the index of refraction is chosen between 1.1 and 1.7 by changing the filling factor. Combination with mullite is suitable for wide-band millimeter wave measurements with sufficient IR cutoff capability. We present the material properties, fabrication of a large-diameter IR filter made of alumina with this AR coating method, and characterizations at cryogenic temperatures. This technology can be applied to a low-temperature receiver system with a large-diameter focal plane for next-generation cosmic microwave background (CMB) polarization measurements, such as POLARBEAR-2 (PB-2).

POLARBEAR-2 (PB-2) is a cosmic microwave background (CMB) polarization experiment that will be located in the Atacama highland in Chile at an altitude of 5200 m. Its science goals are to measure the CMB polarization signals originating from both primordial gravitational waves and weak lensing. PB-2 is designed to measure the tensor to scalar ratio, r, with precision {\sigma}(r) < 0.01, and the sum of neutrino masses, {\Sigma}m{\nu}, with {\sigma}({\Sigma}m{\nu}) < 90 meV. To achieve these goals, PB-2 will employ 7588 transition-edge sensor bolometers at 95 GHz and 150 GHz, which will be operated at the base temperature of 250 mK. Science observations will begin in 2017.

We constrain anisotropic cosmic birefringence using four-point correlations of even-parity $E$-mode and odd-parity $B$-mode polarization in the cosmic microwave background measurements made by the POLARization of the Background Radiation (POLARBEAR) experiment in its first season of observations. We find that the anisotropic cosmic birefringence signal from any parity-violating processes is consistent with zero. The Faraday rotation from anisotropic cosmic birefringence can be compared with the equivalent quantity generated by primordial magnetic fields if they existed. The POLARBEAR nondetection translates into a 95% confidence level (C.L.) upper limit of 93 nanogauss (nG) on the amplitude of an equivalent primordial magnetic field inclusive of systematic uncertainties. This four-point correlation constraint on Faraday rotation is about 15 times tighter than the upper limit of 1380 nG inferred from constraining the contribution of Faraday rotation to two-point correlations of $B$-modes measured by Planck in 2015. Metric perturbations sourced by primordial magnetic fields would also contribute to the $B$-mode power spectrum. Using the POLARBEAR measurements of the $B$-mode power spectrum (two-point correlation), we set a 95% C.L. upper limit of 3.9 nG on primordial magnetic fields assuming a flat prior on the field amplitude. This limit is comparable to what was found in the Planck 2015 two-point correlation analysis with both temperature and polarization. We perform a set of systematic error tests and find no evidence for contamination. This work marks the first time that anisotropic cosmic birefringence or primordial magnetic fields have been constrained from the ground at subdegree scales.

The readout of transition-edge sensor (TES) bolometers with a large multiplexing factor is key for the next generation Cosmic Microwave Background (CMB) experiment, Polarbear-2, having 7,588 TES bolometers. To enable the large arrays, we have been developing a readout system with a multiplexing factor of 40 in the frequency domain. Extending that architecture to 40 bolometers requires an increase in the bandwidth of the SQUID electronics above 4 MHz. This paper focuses on cryogenic readout and shows how it affects cross talk and the responsivity of the TES bolometers. A series resistance, such as equivalent series resistance (ESR) of capacitors for LC filters, leads to non-linear response of the bolometers. A wiring inductance modulates a voltage across the bolometers and causes cross talk. They should be controlled well to reduce systematic errors in CMB observations. We have been developing a cryogenic readout with a low series impedance and have tuned bolometers in the middle of their transition at a high frequency (> 3 MHz).

We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out by SQUIDs with 40 channel frequency domain multiplexing. Refractive optical elements are made with high purity alumina to achieve high optical throughput. The receiver is designed to achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each frequency band. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz frequency bands for foreground control. The Simons Array will be able to constrain tensor-to-scalar ratio and sum of neutrino masses to $\sigma(r) = 6\times 10^{-3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

Atmosphere is one of the most important noise sources for ground-based cosmic microwave background (CMB) experiments. By increasing optical loading on the detectors, it amplifies their effective noise, while its fluctuations introduce spatial and temporal correlations between detected signals. We present a physically motivated 3d-model of the atmosphere total intensity emission in the millimeter and sub-millimeter wavelengths. We derive a new analytical estimate for the correlation between detectors time-ordered data as a function of the instrument and survey design, as well as several atmospheric parameters such as wind, relative humidity, temperature and turbulence characteristics. Using an original numerical computation, we examine the effect of each physical parameter on the correlations in the time series of a given experiment. We then use a parametric-likelihood approach to validate the modeling and estimate atmosphere parameters from the POLARBEAR-I project first season data set. We derive a new 1.0% upper limit on the linear polarization fraction of atmospheric emission. We also compare our results to previous studies and weather station measurements. The proposed model can be used for realistic simulations of future ground-based CMB observations.

K2-19 (EPIC201505350) is an interesting planetary system in which two transiting planets with radii ~ 7 $R_{Earth}$ (inner planet b) and ~ 4 $R_{Earth}$ (outer planet c) have orbits that are nearly in a 3:2 mean-motion resonance. Here, we present results of ground-based follow-up observations for the K2-19 planetary system. We have performed high-dispersion spectroscopy and high-contrast adaptive-optics imaging of the host star with the HDS and HiCIAO on the Subaru 8.2m telescope. We find that the host star is relatively old (>8 Gyr) late G-type star ($T_{eff}$ ~ 5350 K, $M_s$ ~ 0.9 $M_{Sun}$, and $R_{s}$ ~ 0.9 $R_{Sun}$). We do not find any contaminating faint objects near the host star which could be responsible for (or dilute) the transit signals. We have also conducted transit follow-up photometry for the inner planet with KeplerCam on the FLWO 1.2m telescope, TRAPPISTCAM on the TRAPPIST 0.6m telescope, and MuSCAT on the OAO 1.88m telescope. We confirm the presence of transit-timing variations, as previously reported by Armstrong and coworkers. We model the observed transit-timing variations of the inner planet using the synodic chopping formulae given by Deck & Agol (2015). We find two statistically indistinguishable solutions for which the period ratios ($P_{c}/P_{b}$) are located slightly above and below the exact 3:2 commensurability. Despite the degeneracy, we derive the orbital period of the inner planet $P_b$ ~ 7.921 days and the mass of the outer planet $M_c$ ~ 20 $M_{Earth}$. Additional transit photometry (especially for the outer planet) as well as precise radial-velocity measurements would be helpful to break the degeneracy and to determine the mass of the inner planet.

Mini-Neptunes and volatile-poor super-Earths coexist on adjacent orbits in proximity to host stars such as Kepler-36 and Kepler-11. Several post-formation processes have been proposed for explaining the origin of the compositional diversity: the mass loss via stellar XUV irradiation, degassing of accreted material, and in-situ accumulation of the disk gas. Close-in planets are also likely to experience giant impacts during the advanced stage of planet formation. This study examines the possibility of transforming volatile-rich super-Earths / mini-Neptunes into volatile-depleted super-Earths through giant impacts. We present the results of three-dimensional giant impact simulations in the accretionary and disruptive regimes. Target planets are modeled with a three-layered structure composed of an iron core, silicate mantle and hydrogen/helium envelope. In the disruptive case, the giant impact can remove most of the H/He atmosphere immediately and homogenize the refractory material in the planetary interior. In the accretionary case, the planet can retain more than half of the gaseous envelope, while a compositional gradient suppresses efficient heat transfer as its interior undergoes double-diffusive convection. After the giant impact, a hot and inflated planet cools and contracts slowly. The extended atmosphere enhances the mass loss via both a Parker wind induced by thermal pressure and hydrodynamic escape driven by the stellar XUV irradiation. As a result, the entire gaseous envelope is expected to be lost due to the combination of those processes in both cases. We propose that Kepler-36b may have been significantly devolatilized by giant impacts, while a substantial fraction of Kepler-36c's atmosphere may remain intact. Furthermore, the stochastic nature of giant impacts may account for the large dispersion in the mass--radius relationship of close-in super-Earths and mini-Neptunes.

We report detections of new exoplanets from a radial velocity (RV) survey of metal-rich FGK stars by using three telescopes. By optimizing our RV analysis method to long time-baseline observations, we have succeeded in detecting five new Jovian-planets around three metal-rich stars HD 1605, HD 1666, and HD 67087 with the masses of $1.3 M_{\odot}$, $1.5 M_{\odot}$, and $1.4 M_{\odot}$, respectively. A K1 subgiant star HD 1605 hosts two planetary companions with the minimum masses of $ M_p \sin i = 0.96 M_{\mathrm{JUP}}$ and $3.5 M_{\mathrm{JUP}}$ in circular orbits with the planets' periods $P = 577.9$ days and $2111$ days, respectively. HD 1605 shows a significant linear trend in RVs. Such a system consisting of Jovian planets in circular orbits has rarely been found and thus HD 1605 should be an important example of a multi-planetary system that is likely unperturbed by planet-planet interactions. HD 1666 is a F7 main sequence star which hosts an eccentric and massive planet of $ M_p \sin i = 6.4 M_{\mathrm{JUP}}$ in the orbit with $a_{\rm p} = 0.94$ AU and an eccentricity $e=0.63$. Such an eccentric and massive planet can be explained as a result of planet-planet interactions among Jovian planets. While we have found the large residuals of $\mathrm{rms} = 35.6\ \mathrm{m\ s^{-1}}$, the periodogram analysis does not support any additional periodicities. Finally, HD 67087 hosts two planets of $ M_p \sin i = 3.1 M_{\mathrm{JUP}}$ and $4.9 M_{\mathrm{JUP}}$ in orbits with $P=352.2$ days and $2374$ days, and $e=0.17$ and $0.76$, respectively. Although the current RVs do not lead to accurate determinations of its orbit and mass, HD 67087 c can be one of the most eccentric planets ever discovered in multiple systems.

POLARBEAR-2 is a next-generation receiver for precision measurements of the polarization of the cosmic microwave background (Cosmic Microwave Background (CMB)). Scheduled to deploy in early 2015, it will observe alongside the existing POLARBEAR-1 receiver, on a new telescope in the Simons Array on Cerro Toco in the Atacama desert of Chile. For increased sensitivity, it will feature a larger area focal plane, with a total of 7,588 polarization sensitive antenna-coupled Transition Edge Sensor (TES) bolometers, with a design sensitivity of 4.1 uKrt(s). The focal plane will be cooled to 250 milliKelvin, and the bolometers will be read-out with 40x frequency domain multiplexing, with 36 optical bolometers on a single SQUID amplifier, along with 2 dark bolometers and 2 calibration resistors. To increase the multiplexing factor from 8x for POLARBEAR-1 to 40x for POLARBEAR-2 requires additional bandwidth for SQUID readout and well-defined frequency channel spacing. Extending to these higher frequencies requires new components and design for the LC filters which define channel spacing. The LC filters are cold resonant circuits with an inductor and capacitor in series with each bolometer, and stray inductance in the wiring and equivalent series resistance from the capacitors can affect bolometer operation. We present results from characterizing these new readout components. Integration of the readout system is being done first on a small scale, to ensure that the readout system does not affect bolometer sensitivity or stability, and to validate the overall system before expansion into the full receiver. We present the status of readout integration, and the initial results and status of components for the full array.

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.

Gravitational lensing due to the large-scale distribution of matter in the cosmos distorts the primordial Cosmic Microwave Background (CMB) and thereby induces new, small-scale $B$-mode polarization. This signal carries detailed information about the distribution of all the gravitating matter between the observer and CMB last scattering surface. We report the first direct evidence for polarization lensing based on purely CMB information, from using the four-point correlations of even- and odd-parity $E$- and $B$-mode polarization mapped over $\sim30$ square degrees of the sky measured by the POLARBEAR experiment. These data were analyzed using a blind analysis framework and checked for spurious systematic contamination using null tests and simulations. Evidence for the signal of polarization lensing and lensing $B$-modes is found at 4.2$\sigma$ (stat.+sys.) significance. The amplitude of matter fluctuations is measured with a precision of $27\%$, and is found to be consistent with the Lambda Cold Dark Matter ($\Lambda$CDM) cosmological model. This measurement demonstrates a new technique, capable of mapping all gravitating matter in the Universe, sensitive to the sum of neutrino masses, and essential for cleaning the lensing $B$-mode signal in searches for primordial gravitational waves.

We reconstruct the gravitational lensing convergence signal from Cosmic Microwave Background (CMB) polarization data taken by the POLARBEAR experiment and cross-correlate it with Cosmic Infrared Background (CIB) maps from the Herschel satellite. From the cross-spectra, we obtain evidence for gravitational lensing of the CMB polarization at a statistical significance of 4.0$\sigma$ and evidence for the presence of a lensing $B$-mode signal at a significance of 2.3$\sigma$. We demonstrate that our results are not biased by instrumental and astrophysical systematic errors by performing null-tests, checks with simulated and real data, and analytical calculations. This measurement of polarization lensing, made via the robust cross-correlation channel, not only reinforces POLARBEAR auto-correlation measurements, but also represents one of the early steps towards establishing CMB polarization lensing as a powerful new probe of cosmology and astrophysics.

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.

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties.

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 report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

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.

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 report the detection of a large mass planet orbiting around the K0 metal-rich subgiant HD38801 ($V=8.26$) by precise radial velocity (RV) measurements from the Subaru Telescope and the Keck Telescope. The star has a mass of $1.36M_{\odot}$ and metallicity of [Fe/H]= +0.26. The RV variations are consistent with a circular orbit with a period of 696.0 days and a velocity semiamplitude of 200.0\mps, which yield a minimum-mass for the companion of $10.7\mjup$ and semimajor axis of 1.71 AU. Such super-massive objects with very low-eccentricities and hundreds of days period are uncommon among the ensemble of known exoplanets.

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.

We report the detection of a substellar companion orbiting the G5 dwarf HD 16760 from the N2K sample. Precise Doppler measurements of the star from Subaru and Keck revealed a Keplerian velocity variation with a period of 466.47+-0.35 d, a semiamplitude of 407.71+-0.84 m/s, and an eccentricity of 0.084+-0.003. Adopting a stellar mass of 0.78+-0.05 M_Sun, we obtain a minimum mass for the companion of 13.13+-0.56 M_JUP, which is close to the planet/brown-dwarf transition, and the semimajor axis of 1.084+-0.023 AU. The nearly circular orbit despite the large mass and intermediate orbital period makes this companion unique among known substellar companions.

We study the three-dimensional (3D) compact U(1) lattice gauge theory coupled with $N$-flavor Higgs fields by means of the Monte Carlo simulations. This model is relevant to multi-component superconductors, antiferromagnetic spin systems in easy plane, inflational cosmology, etc. It is known that there is no phase transition in the N=1 model. For N=2, we found that the system has a second-order phase transition line $\tilde{c}_1(c_2)$ in the $c_2$(gauge coupling)$-c_1$(Higgs coupling) plane, which separates the confinement phase and the Higgs phase. Numerical results suggest that the phase transition belongs to the universality class of the 3D XY model as the previous works by Babaev et al. and Smiseth et al. suggested. For N=3, we found that there exists a critical line similar to that in the N=2 model, but the critical line is separated into two parts; one for $c_2 < c_{2{\rm tc}}=2.4\pm 0.1$ with first-order transitions, and the other for $ c_{2{\rm tc}} < c_2$ with second-order transitions, indicating the existence of a tricritical point. We verified that similar phase diagram appears for the N=4 and N=5 systems. We also studied the case of anistropic Higgs coupling in the N=3 model and found that there appear two second-order phase transitions or a single second-order transition and a crossover depending on the values of the anisotropic Higgs couplings. This result indicates that an "enhancement" of phase transition occurs when multiple phase transitions coincide at a certain point in the parameter space.