This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper (Paper I), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass $10^{6}\rm{M}_{\odot}$. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars ($M_{\star} \lesssim 0.7~ {\rm M}_\odot$) is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor 2--2.5. The traditional estimate of the range of energies found in the debris is $\approx 1.4\times$ too large for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole $\sim O(10)r_{\rm t}$.

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses ($0.15 \leq M_*/M_\odot \leq 10$) and six black hole masses ($10^5 \leq M_{BH}/M_\odot \leq 5 \times 10^7$). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of $M_*$ for $M_*/M_\odot \lesssim 3$. Independent of $M_*$, relativistic effects increase the critical pericenter distance for full disruptions by up to a factor $\sim 3$ relative to the Newtonian prediction. The probability of a direct capture is also independent of $M_*$; at $M_{BH}/M_\odot \simeq 5 \times 10^6$ this probability is equal to that of a complete disruption. The width of the debris energy distribution $\Delta E$ can differ from the standard estimate by factors from 0.35 to 2, depending on $M_*$ and $M_{BH}$, implying a corresponding change in the characteristic mass-return timescale. The "frozen-in approximation" is inconsistent with $\Delta E$, and mass-loss continues over a long span of time. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly-universal relation between the star's angular momentum and the fraction of $M_*$ remaining. Within the "empty loss-cone" regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass ($10^{5}\leq M_{\rm BH}/\mathrm{M}_{\odot}\leq 5\times 10^{7}$) and stellar mass ($0.3 \leq M_{\star} /\mathrm{M}_{\odot}\leq 3$) as larger $M_{\rm BH}$ leads to stronger relativistic effects. For fixed $M_{\star}$, as $M_{\rm BH}$ increases, the ratio of the maximum pericenter distance yielding full disruptions ($\mathcal{R}_{\rm t}$) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for $M_{\rm BH} = 5\times10^{7}~{\rm M}_\odot$, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of two correction at $M_{\rm BH} = 5 \times 10^7~{\rm M}_\odot$. We find that for partial disruptions, the fractional remnant mass for a given ratio of the pericenter to $\mathcal{R}_{\rm t}$ is higher for larger $M_{\rm BH}$. These results have several implications. As $M_{\rm BH}$ increases above $\sim 10^7~{\rm M}_\odot$, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross section ratio for partial to complete disruptions depends only weakly on $M_{\rm BH}$. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For $M_{\rm BH} \gtrsim 10^7~{\rm M}_\odot$, the $M_{\rm BH}$-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Tidal disruption events (TDEs), events in which a star passes very close to a supermassive black hole, are generally imagined as leading either to the star's complete disruption or to its passage directly into the black hole. In the former case it is widely believed that in all cases the bound portion of the debris quickly "circularizes" due to relativistic apsidal precession, i.e., forms a compact accretion disk, and emits a flare of standardized lightcurve and spectrum. We show here that TDEs are more diverse and can be grouped into several distinct categories on the basis of stellar pericenter distance $r_p$; we calculate the relative frequency of these categories. In particular, because rapid circularization requires $r_p \lesssim 10r_g$ ($r_g \equiv GM_{\rm BH}/c^2$), it can happen in only a minority of total disruptions, $\lesssim 1/4$ when the black hole has mass $M_{\rm BH} = 10^6 M_\odot$. For larger pericenter distances, $10 < r_p/r_g < 27$ (for $M_{\rm BH}=10^6M_\odot$), main sequence stars are completely disrupted, but the bound debris orbits are highly eccentric and possess semimajor axes $\sim 100\times$ the scale of the expected compact disk. Partial disruptions with fractional mass-loss $\gtrsim 10\%$ should occur with a rate similar to that of total disruptions; for fractional mass-loss $\gtrsim 50\%$, the rate is $\approx 1/3$ as large. Partial disruptions---which must precede total disruptions when the stars' angular momenta evolve in the "empty loss-cone" regime---change the orbital energy by factors $\gtrsim O(1)$. Remnants of partial disruptions are in general far from thermal equilibrium. Depending on the orbital energy of the remnant and conditions within the stellar cluster surrounding the SMBH, it may return after hundreds or thousands of years and be fully disrupted, or it may rejoin the stellar cluster.

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the stellar mass dependence of the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than main sequence stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris, but of either sign with respect to the black hole potential, while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. The remnants most deeply bound to the black hole go through a second tidal disruption event upon their first return to pericenter; if they have not thermally relaxed, they will be completely disrupted.

In this paper, the third in this series, we continue our study of tidal disruption events of main-sequence stars by a non-spinning $10^{6}~\rm{M}_\odot$ supermassive black hole. Here we focus on the outcomes of partial disruptions. As the encounter becomes weaker, the debris mass is increasingly concentrated near the outer edges of the energy distribution. As a result, the mass fallback rate can deviate substantially from a $t^{-5/3}$ power-law, becoming more like a single peak with a tail declining as $t^{-p}$ with $p\simeq2-5$. Surviving remnants are spun-up in the prograde direction and are hotter than MS stars of the same mass. Their specific orbital energy is $\simeq10^{-3}\times$ that of the debris (but of either sign with respect to the black hole potential) while their specific angular momentum is close to that of the original star. Even for strong encounters, remnants have speeds at infinity relative to the black hole potential $\lesssim 300$ km s$^{-1}$, so they are unable to travel far out into the galactic bulge. Remnants bound to the black hole can possibly go through a second tidal disruption event.

A star is tidally disrupted by a supermassive black hole when their separation is shorter than the "tidal radius". This quantity is often estimated on an order-of-magnitude basis without reference to the star's internal structure. Using MESA models for main sequence stars and fully general relativistic dynamics, we find the physical tidal radius for complete disruption $\cal{R}_t$ for a $10^6M_\odot$ black hole (BH). We find that across a factor $\sim20$ in stellar mass $M_*$, i.e., $0.15M_{\odot}\leq M_*\leq3M_\odot$, $\cal{R}_t\sim27\times$(BH's gravitational radius). When comparing $\cal{R}_t$ with the commonly used order-of-magnitude estimate $r_t$, we find that $\cal{R}_t\sim1.05-1.45r_t$ for $0.15M_\odot\leq M_*\leq0.5M_\odot$, but between $0.5 M_\odot$ and $1 M_\odot$, $\cal{R}_t$ drops to $\sim 0.45r_t$, and it remains at this value up to $10 M_\odot$. The near-constancy of $\cal{R}_t$ implies a weaker dependence of the full disruption rate on $M_*$ than when predicted with $r_t$. The characteristic energy width of the debris $\Delta E$ ranges from $\sim1.2\Delta\cal{E}$ for low-mass stars to $\sim 0.35\Delta\cal{E}$ for higher-mass stars, where $\Delta\cal{E}=GM_{\rm BH}R_*/\cal{R}_t^{2}$. We present analytic fits for the $M_*$ dependence of $\cal{R}_t$ and $\Delta E$; these fits lead to analytic expressions for the time of peak mass fallback rate and the maximal mass fallback rate. Our results also bear on the fraction of events leading to fast or slow circularization, as well as on the character of the tidal event occurring when the remnant of a partial disruption returns to the black hole. Using a semi-analytic model, we show that $\cal{R}_t$ is primarily determined by the star's central density rather than its mean density. For high-mass stars, the full disruption rate is roughly 1/4 the partial disruption rate, while this ratio is close to unity for low-mass stars.

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. We consider eight different stellar masses, from $0.15~{\rm M}_\odot$ to $10~{\rm M}_\odot$. In the first paper (Ryu et al. 2019a), we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method and provide details about the outcomes of full disruptions. We find that, relative to the traditional order-of-magnitude estimate $r_{\rm t}$, the physical tidal radius of low-mass stars is larger by tens of percent, while for high-mass stars ($M_{\star} \gtrsim1~ {\rm M}_\odot$) it is smaller by a factor $2-2.5$. The traditional estimate of the range of energies found in the debris is approximately accurate for low-mass stars, but is a factor $\sim 2$ too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass-loss continues for a long time because the ${\it instantaneous}$ tidal radius, the distance out to which the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface, stays comparable to the distance to the black hole until the star has reached $O(10)~r_{\rm t}$. These findings indicate significant failings in the popular "frozen-in" approximation.

We have performed high resolution $3-$dimensional compressible hydrodynamics simulations to investigate the effects of shocks and turbulence on energy transport into hot Jupiter atmospheres, under a variety of shear gradients. We focus on a local atmospheric region to accurately follow the small-scale structures of turbulence and shocks. We find that the effects of turbulence above and below a shear layer are different in scale and magnitude: below the shear layer, the effects of turbulence on the vertical energy transfer are local, generally $\lesssim2~\times(\text{scale height})$. However, turbulence can have a spatially and thermally-large influence on almost the entire region above the shear layer. We also find that shock formation is local and transient. Once the atmosphere becomes steady, the time-averaged heat energy flux at $P\sim 1$ bar is insignificant, on the order of 0.001\% of the incoming stellar flux with a shear motion at $P\sim 1$ mbar, and 0.1\% with a deeper shear layer at $P\sim 100$ mbar. Accordingly, the diffusion coefficient is higher for the deeper shear layer. Therefore, our results suggest that turbulence near less dense regions ($P\sim 1$ mbar) does not cause a sufficiently deep and large penetration of thermal energy to account for radius inflation in hot Jupiters, regardless of how violent the turbulence is. However, as the shear layer gets deeper, heat energy transfer becomes more effective throughout the atmosphere (upwards and downwards) due to a larger kinetic energy budget. Therefore, it is more important how deep turbulence occurs in the atmosphere, than how unstable the atmosphere is for effective energy transfer.

We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-life for the distribution of disruption times. As in radioactive decay, the half-lives should provide a direct prediction for the relative fractions of each decay product. We test this prediction with simulated data and find good agreement with our hypothesis. We briefly discuss applications of this feature of the gravitational four-body problem to populations of black holes in globular clusters. This paper, the second in the series, builds on extending the remarkable similarity between gravitational chaos at the macroscopic scale and radioactive decay at the microscopic scale to larger-N systems.

We present the detection and follow-up observations of planetary candidates around low-mass stars observed by the K2 mission. Based on light-curve analysis, adaptive-optics imaging, and optical spectroscopy at low and high resolution (including radial velocity measurements), we validate 16 planets around 12 low-mass stars observed during K2 campaigns 5-10. Among the 16 planets, 12 are newly validated, with orbital periods ranging from 0.96-33 days. For one of the planets (K2-151b) we present ground-based transit photometry, allowing us to refine the ephemerides. Combining our K2 M-dwarf planets together with the validated or confirmed planets found previously, we investigate the dependence of planet radius $R_p$ on stellar insolation and metallicity [Fe/H]. We confirm that for periods $P\lesssim 2$ days, planets with a radius $R_p\gtrsim 2\,R_\oplus$ are less common than planets with a radius between 1-2$\,R_\oplus$. We also see a hint of the "radius valley" between 1.5 and 2$\,R_\oplus$ that has been seen for close-in planets around FGK stars. These features in the radius/period distribution could be attributed to photoevaporation of planetary envelopes by high-energy photons from the host star, as they have for FGK stars. For the M dwarfs, though, the features are not as well defined, and we cannot rule out other explanations such as atmospheric loss from internal planetary heat sources, or truncation of the protoplanetary disk. There also appears to be a relation between planet size and metallicity: those few planets larger than about 3 $R_\oplus$ are found around the most metal-rich M dwarfs.

Understanding of the role of X-rays for driving the thermal evolution of the intergalactic medium (IGM) at high redshifts is one of important questions in astrophysics. High-mass X-ray binaries (HMXBs) in early stellar populations are prime X-ray source; however, their formation efficiency is not well understood. Using $N$-body simulations, we estimate the HMXB formation rate via mutual gravitational interactions of nascent, small groups of the Population~III stars. We find that HMXBs form at a rate of one per $\gtrsim 10^{4}M_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} {\rm erg}~{\rm s}^{-1}$ in the $2-10$ keV band per star formation rate (SFR). This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed \textit{in situ} or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long gamma-ray bursts from Population~III stars and the direct collapse channel for massive black hole formation.

Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.

We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $\mu$ Col and binary $\iota$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $\mu$ Col.

In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $\mu$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.

We report extensive 3-yr multiwavelength observations of the WZ Sge-type dwarf nova SSS J122221.7-311525 during its unusual double superoutburst, the following decline and in quiescence. The second segment of the superoutburst had a long duration of 33 d and a very gentle decline with a rate of 0.02 mag/d, and it displayed an extended post-outburst decline lasting at least 500 d. Simultaneously with the start of the rapid fading from the superoutburst plateau, the system showed the appearance of a strong near-infrared excess resulting in very red colours, which reached extreme values (B-I~1.4) about 20 d later. The colours then became bluer again, but it took at least 250 d to acquire a stable level. Superhumps were clearly visible in the light curve from our very first time-resolved observations until at least 420 d after the rapid fading from the superoutburst. The spectroscopic and photometric data revealed an orbital period of 109.80 min and a fractional superhump period excess <0.8 per cent, indicating a very low mass ratio q<0.045. With such a small mass ratio the donor mass should be below the hydrogen-burning minimum mass limit. The observed infrared flux in quiescence is indeed much lower than is expected from a cataclysmic variable with a near-main-sequence donor star. This strongly suggests a brown-dwarf-like nature for the donor and that SSS J122221.7-311525 has already evolved away from the period minimum towards longer periods, with the donor now extremely dim.

We perform a large suite of $N=4$ numerical scattering experiments between two identical binaries consisting of identical point particles in a (continuous) background potential. For investigative purposes, albeit without loss of generality, we assume that the potential corresponds to a uniform (natal or star-forming) gas medium. We explore a range of constant gas densities, from $n=10~ {\rm cm}^{-3}$ to $10^{5}~ {\rm cm}^{-3}$. These densities are relevant for various astrophysical environments, including molecular clouds (i.e., star-forming regions) and denser, fragmented cores within these clouds. Our primary goal is to characterize the effects of the background potential on the subsequent stellar dynamics. We consider the outcome probabilities as well as the properties of any binaries formed during the binary-binary encounters, such as the distributions of binary binding energies and eccentricities. We also present the final velocity distributions of the ejected single stars. The background potential has two important effects on the stellar dynamics: 1) The potential acts to reset the zero-point of the total system energy, which in turn affects the types and properties of the products of the encounter; 2) For higher $n$ and weakly bound systems (i.e., large semimajor axes), the stellar dynamics are significantly affected when stars become trapped in the potential, oscillating around the system centre of mass (CM). This, in turn, increases the number of scattering events between stars (single, binary or triple) near the CM and makes it harder for single stars to escape to infinity. This ultimately leads to the preferential ionization of triples and wide binaries and the survival of compact binaries, with the single stars escaping at very high ejection velocities.

We present three-band simultaneous observations of a weak-line T-Tauri star CVSO~30 (PTFO~8-8695), which is one of the youngest objects having a candidate transiting planet. The data were obtained with the Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT) on the 188 cm telescope at Okayama Astrophysical Observatory in Japan. We observed the fading event in the $g^{\prime}_2$-, $r^{\prime}_2$-, and $z_{\rm s,2}$-bands simultaneously. As a result, we find a significant wavelength dependence of fading depths of about 3.1\%, 1.7\%, 1.0\% for the $g^{\prime}_2$-, $r^{\prime}_2$-, and $z_{\rm s,2}$-bands, respectively. A cloudless H/He dominant atmosphere of a hot Jupiter cannot explain this large wavelength dependence. Additionally, we rule out a scenario by the occultation of the gravity-darkened host star. Thus our result is in favor of the fading origin as circumstellar dust clump or occultation of an accretion hotspot.

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 study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an updated model of stellar evolution specific for Pop III stars. We find that the merger rate of these Pop III BH-BH systems is relatively small (< 0.1 Gpc^-3 yr^-1) at low redshifts (z<2), where it can be compared with the LIGO empirical estimate of 9-240 Gpc^-3 yr^-1 (Abbott et al. 2016). The predicted rates are even smaller for Pop III double neutron star and black hole neutron star mergers. Our rates are compatible with those of Hartwig et al. (2016), but significantly smaller than those found in previous work (Bond & Carr 1984; Belczynski et al. 2004; Kinugawa et al. 2014, 2016). We explain the reasons for this discrepancy by means of detailed model comparisons and point out that (i) identification of Pop III BH-BH mergers may not be possible by advanced LIGO, and (ii) the level of stochastic gravitational wave background from Pop III mergers may be lower than recently estimated (Kowalska et al. 2012; Inayoshi et al. 2016; Dvorkin et al. 2016). We further estimate gravitational wave detection rates for third-generation interferometric detectors. Our calculations are relevant for low to moderately rotating Pop III stars. We can now exclude significant (> 1 per cent) contribution of these stars to low-redshift BH-BH mergers. However, it remains to be tested whether (and at what level) rapidly spinning Pop III stars (homogeneous evolution) can contribute to BH-BH mergers in the local Universe.

We report the first ground-based transit observation of K2-3d, a 1.5 R_Earth planet supposedly within the habitable zone around a bright M-dwarf host star, using the Okayama 188 cm telescope and the multi(grz)-band imager MuSCAT. Although the depth of the transit (0.7 mmag) is smaller than the photometric precisions (1.2, 0.9, and 1.2 mmag per 60 s for the g, r, and z bands, respectively), we marginally but consistently identify the transit signal in all three bands, by taking advantage of the transit parameters from K2, and by introducing a novel technique that leverages multi-band information to reduce the systematics caused by second-order extinction. We also revisit previously analyzed Spitzer transit observations of K2-3d to investigate the possibility of systematic offsets in transit timing, and find that all the timing data can be explained well by a linear ephemeris. We revise the orbital period of K2-3d to be 44.55612 \pm 0.00021 days, which corrects the predicted transit times for 2019, i.e., the era of the James Webb Space Telescope, by \sim80 minutes. Our observation demonstrates that (1) even ground-based, 2 m class telescopes can play an important role in refining the transit ephemeris of small-sized, long-period planets, and that (2) a multi-band imager is useful to reduce the systematics of atmospheric origin, in particular for bluer bands and for observations conducted at low-altitude observatories.

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data.

We report spectroscopic and photometric follow-up of the peculiar nova V5852~Sgr (discovered as OGLE-2015-NOVA-01), which exhibits a combination of features from different nova classes. The photometry shows a flat-topped light curve with quasi-periodic oscillations, then a smooth decline followed by two fainter recoveries in brightness. Spectroscopy with the Southern African Large Telescope shows first a classical nova with an Fe II or Fe IIb spectral type. In the later spectrum, broad emissions from helium, nitrogen and oxygen are prominent and the iron has faded which could be an indication to the start of the nebular phase. The line widths suggest ejection velocities around $1000\,{\rm km\,s^{-1}}$. The nova is in the direction of the Galactic bulge and is heavily reddened by an uncertain amount. The $V$ magnitude 16 days after maximum enables a distance to be estimated and this suggests that the nova may be in the extreme trailing stream of the Sagittarius dwarf spheroidal galaxy. If so it is the first nova to be detected from that, or from any dwarf spheroidal galaxy. Given the uncertainty of the method and the unusual light curve we cannot rule out the possibility that it is in the bulge or even the Galactic disk behind the bulge.

We report on the discovery and characterization of the transiting planet K2-39b (EPIC 206247743b). With an orbital period of 4.6 days, it is the shortest-period planet orbiting a subgiant star known to date. Such planets are rare, with only a handful of known cases. The reason for this is poorly understood, but may reflect differences in planet occurrence around the relatively high-mass stars that have been surveyed, or may be the result of tidal destruction of such planets. K2-39 is an evolved star with a spectroscopically derived stellar radius and mass of $3.88^{+0.48}_{-0.42}~\mathrm{R_\odot}$ and $1.53^{+0.13}_{-0.12}~\mathrm{M_\odot}$, respectively, and a very close-in transiting planet, with $a/R_\star = 3.4$. Radial velocity (RV) follow-up using the HARPS, FIES and PFS instruments leads to a planetary mass of $50.3^{+9.7}_{-9.4}~\mathrm{M_\oplus}$. In combination with a radius measurement of $8.3 \pm 1.1~\mathrm{R_\oplus}$, this results in a mean planetary density of $0.50^{+0.29}_{-0.17}$ g~cm$^{-3}$. We furthermore discover a long-term RV trend, which may be caused by a long-period planet or stellar companion. Because K2-39b has a short orbital period, its existence makes it seem unlikely that tidal destruction is wholly responsible for the differences in planet populations around subgiant and main-sequence stars. Future monitoring of the transits of this system may enable the detection of period decay and constrain the tidal dissipation rates of subgiant stars.

We report the formation of intermediate-mass black holes (IMBHs) in suites of numerical $N$-body simulations of Population III remnant black holes (BHs) embedded in gas-rich protogalaxies at redshifts $z\gtrsim10$. We model the effects of gas drag on the BHs' orbits, and allow BHs to grow via gas accretion, including a mode of hyper-Eddington accretion in which photon trapping and rapid gas inflow suppress any negative radiative feedback. Most initial BH configurations lead to the formation of one (but never more than one) IMBH in the center of the protogalaxy, reaching a mass of $10^{3-5}\mathrm{M}_{\odot}$ through hyper-Eddington growth. Our results suggest a viable pathway to forming the earliest massive BHs in the centers of early galaxies. We also find that the nuclear IMBH typically captures a stellar-mass BH companion, making these systems observable in gravitational waves as extreme mass-ratio inspirals (EMRIs) with \textit{eLISA}.

We report on the detection and early characterization of a hot Jupiter in a 3-day orbit around K2-34 (EPIC~212110888), a metal-rich F-type star located in the K2 Cycle 5 field. Our follow-up campaign involves precise radial velocity (RV) measurements and high-contrast imaging using multiple facilities. The absence of a bright nearby source in our high-contrast data suggests that the transit-like signals are not due to light variations from such a companion star. Our intensive RV measurements show that K2-34b has a mass of $1.773\pm0.086M_J$, confirming its status as a planet. We also detect the Rossiter-McLaughlin effect for K2-34b and show that the system has a good spin-orbit alignment ($\lambda=-1_{-9}^{+10}$ degrees). High-contrast images obtained by the HiCIAO camera on the Subaru 8.2-m telescope reveal a faint companion candidate ($\Delta m_H=6.19\pm 0.11$ mag) at a separation of $0\farcs36$. Follow-up observations are needed to confirm that the companion candidate is physically associated with K2-34. K2-34b appears to be an example of a typical "hot Jupiter," albeit one which can be precisely characterized using a combination of K2 photometry and ground-based follow-up.

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $\gamma$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($\iota$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $\iota$ Dra b, HD 5608 b, and HD 14067 b.

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.

We validate a $R_p=2.32\pm 0.24R_\oplus$ planet on a close-in orbit ($P=2.260455\pm 0.000041$ days) around K2-28 (EPIC 206318379), a metal-rich M4-type dwarf in the Campaign 3 field of the K2 mission. Our follow-up observations included multi-band transit observations from the optical to the near infrared, low-resolution spectroscopy, and high-resolution adaptive-optics (AO) imaging. We perform a global fit to all the observed transits using a Gaussian process-based method and show that the transit depths in all passbands adopted for the ground-based transit follow-ups ($r'_2, z_\mathrm{s,2}, J, H, K_\mathrm{s}$) are within $\sim 2\sigma$ of the K2 value. Based on a model of the background stellar population and the absence of nearby sources in our AO imaging, we estimate the probability that a background eclipsing binary could cause a false positive to be $< 2\times 10^{-5}$. We also show that K2-28 cannot have a physically associated companion of stellar type later than M4, based on the measurement of almost identical transit depths in multiple passbands. There is a low probability for a M4 dwarf companion ($\approx 0.072_{-0.04}^{+0.02}$), but even if this were the case, the size of K2-28b falls within the planetary regime. K2-28b has the same radius (within $1\sigma$) and experiences a similar irradiation from its host star as the well-studied GJ~1214b. Given the relative brightness of K2-28 in the near infrared ($m_\mathrm{Kep}=14.85$ mag and $m_H=11.03$ mag) and relatively deep transit ($0.6-0.7\%$), a comparison between the atmospheric properties of these two planets with future observations would be especially interesting.

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.

The first astrophysical objects shaped the cosmic environment by reionizing and heating the intergalactic medium (IGM). In particular, X-rays are very efficient at heating the IGM before it became completely ionized, an effect that can be measured through the 21 cm line of neutral hydrogen. High-mass X-ray binaries (HMXBs), known to be prolific X-ray sources in star-forming galaxies at lower redshifts, are prime candidates for driving the thermal evolution of the IGM at redshifts $z > 20$. Despite their importance, the formation efficiency of HMXBs from the first stellar populations is not well understood---as such, their collective X-ray emission and the subsequent imprint on the 21 cm signature are usually evaluated using free parameters. Using $N$-body simulations, we estimate the rate of HMXB formation via mutual gravitational interactions of nascent, small groups of the first stars (Pop III stars). We run two sets of calculations: one in which stars form in small groups of five in nearly Keplerian initial orbits, and another in which two such groups collide (expected from mergers of host protogalaxies). We find that HMXBs form at a rate of one per $~10^{4}~{\rm M}_{\odot}$ in newly born stars, and that they emit with a power of $\sim 10^{41} ~{\rm erg}~{\rm s} ^{-1}$ in the $2-10$ keV band per solar mass per year of star formation. This value is a factor $\sim 10^{2}$ larger than what is observed in star forming galaxies at lower redshifts; the X-ray production from early HMXBs would have been even more copious, if they also formed $in$ $situ$ or via migration in protostellar disks. Combining our results with earlier studies suggests that early HMXBs were highly effective at heating the IGM and leaving a strong 21 cm signature. We discuss broader implications of our results, such as the rate of long GRBs from PopIII stars and the direct collapse black hole formation.

We report a development of a multi-color simultaneous camera for the 188cm telescope at Okayama Astrophysical Observatory in Japan. The instrument, named MuSCAT, has a capability of 3-color simultaneous imaging in optical wavelength where CCDs are sensitive. MuSCAT is equipped with three 1024x1024 pixel CCDs, which can be controlled independently. The three CCDs detect lights in $g'_2$ (400--550 nm), $r'_2$ (550--700 nm), and $z_{s,2}$ (820--920 nm) bands using Astrodon Photometrics Generation 2 Sloan filters. The field of view of MuSCAT is 6.1x6.1 arcmin$^2$ with the pixel scale of 0.358 arcsec per pixel. The principal purpose of MuSCAT is to perform high precision multi-color transit photometry. For the purpose, MuSCAT has a capability of self autoguiding which enables to fix positions of stellar images within ~1 pix. We demonstrate relative photometric precisions of 0.101%, 0.074%, and 0.076% in $g'_2$, $r'_2$, and $z_{s,2}$ bands, respectively, for GJ436 (magnitudes in $g'$=11.81, $r'$=10.08, and $z'$=8.66) with 30 s exposures. The achieved precisions meet our objective, and the instrument is ready for operation.

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.