We consider black hole formation in failed supernovae when a dense circumstellar medium (CSM) is present around the massive star progenitor. By utilizing radiation hydrodynamical simulations, we calculate the mass ejection of blue supergiants and Wolf-Rayet stars in the collapsing phase and the radiative shock occurring between the ejecta and the ambient CSM. We find that the resultant emission is redder and dimmer than normal supernovae (bolometric luminosity of $\sim 10^{40}-10^{41}\ {\rm erg\ s^{-1}}$, effective temperature of $\sim 5\times 10^3$ K, and timescale of 10-100 days) and shows a characteristic power-law decay, which may comprise a fraction of intermediate luminosity red transients (ILRTs) including AT 2017be. In addition to searching for the progenitor star in the archival data, we encourage X-ray follow-up observations of such ILRTs $\sim$ 1-10 yr after the collapse, targeting the fallback accretion disk.

We discuss the prospects of high precision pointing of our transmitter to habitable planets around Galactic main sequence stars. For an efficient signal delivery, the future sky positions of the host stars should be appropriately extrapolated with accuracy better than the beam opening angle $\Theta$ of the transmitter. Using the latest data release (DR2) of Gaia, we estimate the accuracy of the extrapolations individually for $4.7\times 10^7$ FGK stars, and find that the total number of targets could be $\sim 10^7$ for the accuracy goal better than 1". Considering the pairwise nature of communication, our study would be instructive also for SETI (Search for Extraterrestrial Intelligence), not only for sending signals outward.

WD J005311 is a newly identified white dwarf (WD) in a mid-infrared nebula. The spectroscopic observation indicates the existence of a neon-enriched carbon/oxygen wind with a terminal velocity of $v_{\infty,\rm obs}\sim 16,000\,\rm km\,s^{-1}$ and a mass loss rate of $\dot M_{\rm obs}\sim 3.5\times 10^{-6}\,M_\odot$ yr$^{-1}$. Here we consistently explain the properties of WD J005311 using a newly constructed wind solution, where the optically thick outflow is launched from the carbon burning shell on an oxygen-neon core and accelerated by the rotating magnetic field to become supersonic and unbound well below the photosphere. Our model implies that WD J005311 has a mass of $M_* \sim 1.1\mbox{-}1.3\,M_\odot$, a magnetic field of $B_* \sim (2\mbox{-}5)\times 10^7\,\rm G$, and a spin angular frequency of $\Omega \sim 0.2\mbox{-}0.5 \,\rm s^{-1}$. The large magnetic field and fast spin support the carbon-oxygen WD merger origin. WD J005311 will neither explode as a type Ia supernova nor collapse into a neutron star. If the wind continues to blow another few kyr, WD J005311 will spin down significantly and join to the known sequence of slowly-rotating magnetic WDs. Otherwise it may appear as a fast-spinning magnetic WD and could be a new high energy source.

We present results of a search for late-time radio emission and Fast Radio Bursts (FRBs) from a sample of type-I superluminous supernovae (SLSNe-I). We used the Karl G. Jansky Very Large Array to observe ten SLSN-I more than 5 years old at a frequency of 3 GHz. We searched fast-sampled visibilities for FRBs and used the same data to perform a deep imaging search for late-time radio emission expected in models of magnetar-powered supernovae. No FRBs were found. One SLSN-I, PTF10hgi, is detected in deep imaging, corresponding to a luminosity of $1.2\times10^{28}$ erg s$^{-1}$. This luminosity, considered with the recent 6 GHz detection of PTF10hgi in Eftekhari et al (2019), supports the interpretation that it is powered by a young, fast-spinning ($\sim$ ms spin period) magnetar with $\sim$ 15 Msun of partially ionized ejecta. Broadly, our observations are most consistent with SLSNe-I being powered by neutron stars with fast spin periods, although most require more free-free absorption than is inferred for PTF10hgi. We predict that radio observations at higher frequencies or in the near future will detect these systems and begin constraining properties of the young pulsars and their birth environments.

We present a bolometric light curve model of Type IIn supernovae powered by supernova ejecta colliding with a circumstellar medium. We estimate the conversion efficiency of the ejecta's kinetic energy to radiation at the reverse and forward shocks and find that a large density contrast makes a difference in the efficiency. The emission from the reverse shock can maintain high efficiency for a long time, and becomes important at the late phase of the light curve. We first construct a semi-analytical model that is applicable to the late phase of the light curve when the diffusion time of photons in the shocked region becomes negligible. We further develop radiation transfer simulations that incorporate these physical processes into the light curve. The numerical calculations predict light curves at early phases, which are testable by present and future short-cadence surveys. We compare our model with the bolometric light curve constructed from observations for a type IIn supernova 2005ip. Due to the reduced efficiency at the forward shock, we find from our model that the mass-loss rate of the progenitor star was $\approx 1\times 10^{-2}\ {\rm M_\odot \ yr^{-1}}$ for a wind velocity of $100\ {\rm km \ s^{-1}}$, an order of magnitude higher compared to previous work that used simple assumptions of the efficiency. This highlights the importance of taking these two components into account when extracting the physical parameters from observations.

We investigate effects of energetic pulsar wind nebulae (PWNe) on dust formation and evolution. Dust emission has been observed in many supernova remnants that also have neutron stars as compact remnants. We study the dependence of dust formation time and size on properties of the ejecta and central pulsar. We find that a pulsar with an initial spin period $P \sim 1\mbox{-}10\,\rm ms$ and a dipole magnetic field $B \sim 10^{12\mbox{-}15}\,\rm G$ can either accelerate or delay dust formation, with a timescale of several months to over ten years, and reduce the average size of dust by a factor of $\sim$ 10 or more compared to the non-pulsar case. We also find that infrared dust emission may be detectable in typical superluminous supernovae out to $\sim$ 100-1000 Mpc in 2-5 years after the explosion, although this depends sensitively on the spectral index of nonthermal emission from the nebula. We discuss implications to previous supernova observations. Some discrepancies between dust formation models and observations, such as the formation time in SN1987A or the dust size in the Crab Nebula, could be explained by the influence of a pulsar, and knowledge of the dust emission will be important for future ALMA observations of superluminous supernovae.

In the faint short gamma-ray burst sGRB 170817A followed by the gravitational waves (GWs) from a merger of two neutron stars (NSs) GW170817, the spectral peak energy is too high to explain only by canonical off-axis emission. We investigate off-axis appearance of an sGRB prompt emission scattered by a cocoon, which is produced through the jet-merger-ejecta interaction, with either sub-relativistic or mildly-relativistic velocities. We show that the observed properties of sGRB 170817A, in particular the high peak energy, can be consistently explained by the Thomson-scattered emission with a typical sGRB jet, together by its canonical off-axis emission, supporting that an NS-NS merger is the origin of sGRBs. The scattering occurs at $\lesssim 10^{10}$--$10^{12}\,{\rm cm}$ not far from the central engine, implying the photospheric or internal shock origin of the sGRB prompt emission. The boundary between the jet and cocoon is sharp, which could be probed by future observations of off-axis afterglows. The scattering model predicts a distribution of the spectral peak energy that is similar to the observed one but with a cutoff around $\sim$ MeV energy, and its correlations with the luminosity, duration, and time lag from GWs, providing a way to distinguish it from alternative models.

The flow of fallback matter being shocked and repelled back by an energy deposition from a central object is discussed by using newly found self-similar solutions. We show that there exists a maximum mass accretion rate if the adiabatic index of the flow is less than or equal to 4/3. Otherwise we can find a solution with an arbitrarily large accretion rate by appropriately shrinking the energy deposition region. Applying the self-similar solution to supernova fallback, we discuss how the fate of newborn pulsars or magnetars depends on the fallback accretion and their spin-down power. Combining the condition for the fallback accretion to bury the surface magnetic field into the crust, we argue that supernova fallback with a rate of dM_{fb}/dt ~ 10^{-(4-6)} M_sun /s could be the main origin of the diversity of Galactic young neutron stars, i.e., rotation-powered pulsars, magnetars, and central compact objects.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

We consider radio emission from a newborn black hole (BH), which is accompanied by a mini-disk with a mass of $\lesssim M_\odot$. Such a disk can be formed from an outer edge of the progenitor's envelope, especially for metal-poor massive stars and/or massive stars in close binaries. The disk accretion rate is typically super-Eddington and an ultrafast outflow with a velocity of $\sim 0.1\mbox{-}0.3\,c$ will be launched into the circumstellar medium. The outflow forms a collisionless shock, and electrons are accelerated and emit synchrotron emission in radio bands with a flux of $\sim 10^{26-30} \ \rm erg \ s^{-1} \ Hz^{-1}$ days to decades after the BH formation. The model predicts not only a fast UV/optical transient but also quasi-simultaneous inverse-Compton X-ray emission $\sim$ a few days after the BH formation, and the discovery of the radio counterpart with coordinated searches will enable us to identify this type of transients. The occurrence rate can be $0.1-10 \ \%$ of the core-collapse supernova rate, which makes them a promising target of dedicated radio observations such as the Jansky VLA Sky Survey.

We propose a long-term strategy for detecting thermal neutrinos from the remnant of binary-neutron-star mergers with a future M-ton water-Cherenkov detector such as Hyper-Kamiokande. Monitoring >~2500 mergers within <~200 Mpc, we may be able to detect a single neutrino with a human-time-scale operation of ~80 Mt years for the merger rate of 1 Mpc^{-3} Myr^{-1}, which is slightly lower than the median value derived by the LIGO-Virgo collaboration with GW170817. Although the number of neutrino events is minimal, contamination from other sources of neutrinos can be reduced efficiently to ~0.03 by analyzing only ~1 s after each merger identified with gravitational-wave detectors if Gadolinium is dissolved in the water. The contamination may be reduced further to ~0.01 if we allow the increase of waiting time by a factor of ~1.7. The detection of even a single neutrino can pin down the energy scale of thermal neutrino emission from binary-neutron-star mergers and could strongly support or disfavor formation of remnant massive neutron stars. Because the dispersion relation of gravitational waves are now securely constrained to that of massless particles with a corresponding limit on the graviton mass of <~10^{-22} eV/c^2 by binary-black-hole mergers, the time delay of a neutrino from gravitational waves can be used to put an upper limit of <~O(10) meV/c^2 on the absolute neutrino mass in the lightest eigenstate. Large neutrino detectors will enhance the detectability, and in particular, 5 Mt Deep-TITAND and 10 Mt MICA planned in the future will allow us to detect thermal neutrinos every ~16 and 8 years, respectively, increasing the significance.

We study the ejection of mass during stellar core-collapse when the stalled shock does not revive and a black hole forms. Neutrino emission during the protoneutron star phase causes a decrease in the gravitational mass of the core, resulting in an outward going sound pulse that steepens into a shock as it travels out through the star. We explore the properties of this mass ejection mechanism over a range of stellar progenitors using spherically-symmetric, time-dependent hydrodynamic simulations that treat neutrino mass loss parametrically and follow the shock propagation over the entire star. We find that all types of stellar progenitor can eject mass through this mechanism. The ejected mass is a decreasing function of the surface gravity of the star, ranging from several $M_\odot$ for red supergiants to $\sim 0.1M_\odot$ for blue supergiants and $\sim 10^{-3} M_\odot $ for Wolf-Rayet stars. We find that the final shock energy at the surface is a decreasing function of the core-compactness, and is $\lesssim 10^{47}-10^{48}$ erg in all cases. In progenitors with a sufficiently large envelope, high core-compactness, or a combination of both, the sound pulse fails to unbind mass. Successful mass ejection is accompanied by significant fallback accretion that can last from hours to years. We predict the properties of shock breakout and thermal plateau emission produced by the ejection of the outer envelope of blue supergiant and Wolf-Rayet progenitors in otherwise failed supernovae.

We examine the capability of pulsar timing arrays (PTAs) to detect very small-scale clumps of dark matter (DM), which are a natural outcome of the standard cold dark matter (CDM) paradigm. A clump streaming near the Earth or a pulsar induces an impulsive acceleration to encode residuals on pulsar timing data. We show that, assuming the standard abundance of DM clumps predicted by the CDM model, small-scale DM clumps with masses from $\sim 10^{-11} M_\odot$ to $\sim 10^{-8} \ M_\odot$ can be detectable by a PTA observation for a few decades with ${\cal O}(100)$ of pulsars with a timing noise of ${\cal O}(10)$ ns located at $\gtrsim 3$ kpc away from the Galactic center, as long as these mass scales are larger than the cutoff scale of the halo mass function that is determined by the particle nature of DM. Our result suggests that PTAs can provide a unique opportunity for testing one of the most fundamental predictions of the CDM paradigm. In addition, the detections and non-detections can constrain the cutoff mass scale inherent to the DM model.

The recent detection of gravitational waves and electromagnetic counterparts from the double neutron star merger event GW+EM170817, supports the standard paradigm of short gamma-ray bursts (SGRBs) and kilonovae/macronovae. It is important to reveal the nature of the compact remnant left after the merger, either a black hole or neutron star, and their physical link to the origin of the long-lasting emission observed in SGRBs. The diversity of the merger remnants may also lead to different kinds of transients that can be detected in future. Here we study the high-energy emission from the long-lasting central engine left after the coalescence, under certain assumptions. In particular, we consider the X-ray emission from a remnant disk and the non-thermal nebular emission from disk-driven outflows or pulsar winds. We demonstrate that late-time X-ray and high-frequency radio emission can provide useful constraints on properties of the hidden compact remnants and their connections to long-lasting SGRB emission, and we discuss the detectability of nearby merger events through late-time observations at $\sim30-100$ d after the coalescence. We also investigate the GeV-TeV gamma-ray emission that occurs in the presence of long-lasting central engines, and show the importance of external inverse-Compton radiation due to up-scattering of X-ray photons by relativistic electrons in the jet. We also search for high-energy gamma-rays from GW170817 in the Fermi-LAT data, and report upper limits on such long-lasting emission. Finally, we consider the implications of GW+EM170817 and discuss the constraints placed by X-ray and high-frequency radio observations.

We study optical counterparts of a new-born pulsar in a double neutron star system like PSR J0737-3039A/B. This system is believed to eject a small amount of mass of $\mathcal{O}(0.1M_{\odot})$ at the second core-collapse supernova. We argue that the initial spin of the new-born pulsar can be determined by the orbital period at the time when the second supernova occurs. The spin angular momentum of the progenitor is expected to be similar to that of the He-burning core, which is tidally synchronized with the orbital motion, and then the second remnant may be born as a millisecond pulsar. If the dipole magnetic field strength of the nascent pulsar is comparable to that inferred from the current spin-down rate of PSR J0737-3039B, the initial spin-down luminosity is comparable to the luminosity of super-luminous supernovae. We consider thermal emission arising from the supernova ejecta driven by the relativistic wind from such a new-born pulsar. The resulting optical light curves have a rising time $\sim 10$ days and peak luminosity $\sim 10^{44}$ erg/s. The optical emission may last for a month to several months, due to the reprocessing of X-rays and UV photons via photoelectric absorption. These features are broadly consistent with those of the rapidly-rising optical transients. The high spin-down luminosity and small ejecta mass are favorable for the progenitor of the repeating fast radio burst, FRB 121102. We discuss a possible connection between newborn double pulsars and fast radio bursts.

It has been widely argued that Type-I super-luminous supernovae (SLSNe-I) are driven by powerful central engines with a long-lasting energy injection after the core-collapse of massive progenitors. One of the popular hypotheses is that the hidden engines are fast-rotating pulsars with a magnetic field of $B\sim{10}^{13}-{10}^{15}$ G. Murase, Kashiyama & Meszaros (2016) showed that quasi-steady radio/submm emission from non-thermal electron-positron pairs in nascent pulsar wind nebulae can be used as a counterpart of such pulsar-driven supernovae (SNe). In this work, focusing on the nascent SLSN-I remnants, we examine constraints that can be placed by radio emission. We show that the Atacama Large Millimeter/submillimetre Array (ALMA) can detect the radio nebula from SNe at $D_{\rm L} \sim 1 \ \rm Gpc$ in a few years after the explosion, while the Jansky Very Large Array (VLA) can also detect the counterpart in a few decades. The proposed radio followup observation could solve the parameter degeneracy in the pulsar-driven SN model for optical/UV light curves, and could also give us clues to young neutron star scenarios for SLSNe-I and fast radio bursts.

We propose to use degree-scale angular clustering of fast radio bursts (FRBs) to identify their origin and the host galaxy population. We study the information content in autocorrelation of the angular positions and dispersion measures (DM) and in cross-correlation with galaxies. We show that the cross-correlation with Sloan Digital Sky Survey (SDSS) galaxies will place stringent constraints on the mean physical quantities associated with FRBs. If $\sim$10,000 FRBs are detected with $\lesssim \rm deg$ resolution in the SDSS field, the clustering analysis with the intrinsic DM scatter of $100\, {\rm pc}/{\rm cm}^3$ can constrain the global abundance of free electrons at $z\lt1$ and the large-scale bias of FRB host galaxies (the statistical relation between the distribution of host galaxies and cosmic matter density field) with fractional errors (with a $68\%$ confidence level) of $\sim10\%$ and $\sim20\%$, respectively. The mean near-source dispersion measure and the delay time distribution of FRB rates relative to the global star forming rate can be also determined by combining the clustering and the probability distribution function of DM. Our approach will be complementary to high-resolution ($\ll {\rm deg}$) event localization using e.g., VLA and VLBI for identifying the origin of FRBs and the source environment. We strongly encourage future observational programs such as CHIME, UTMOST, and HIRAX to survey FRBs in the SDSS field.

Recently a repeating fast radio burst (FRB) 121102 has been confirmed to be an extragalactic event and a persistent radio counterpart has been identified. While other possibilities are not ruled out, the emission properties are broadly consistent with Murase et al. (2016) that theoretically proposed quasi-steady radio emission as a counterpart of both FRBs and pulsar-driven supernovae. Here we constrain the model parameters of such a young neutron star scenario for FRB 121102. If the associated supernova has a conventional ejecta mass of $M_{\rm ej}\gtrsim{\rm a \ few}\ M_\odot$, a neutron star with an age of $t_{\rm age} \sim 10-100 \ \rm yrs$, an initial spin period of $P_{i} \lesssim$ a few ms, and a dipole magnetic field of $B_{\rm dip} \lesssim {\rm a \ few} \times 10^{13} \ \rm G$ can be compatible with the observations. However, in this case, the magnetically-powered scenario may be favored as an FRB energy source because of the efficiency problem in the rotation-powered scenario. On the other hand, if the associated supernova is an ultra-stripped one or the neutron star is born by the accretion-induced collapse with $M_{\rm ej} \sim 0.1 \ M_\odot$, a younger neutron star with $t_{\rm age} \sim 1-10$ yrs can be the persistent radio source and might produce FRBs with the spin-down power. These possibilities can be distinguished by the decline rate of the quasi-steady radio counterpart.

Tenuous wind bubbles, which are formed by the spin-down activity of central compact remnants, are relevant in some models of fast radio bursts (FRBs) and super-luminous supernovae. We study their high-energy signatures, focusing on the role of pair-enriched bubbles produced by young magnetars, rapidly-rotating neutron stars, and magnetized white dwarfs. (i) First, we study the nebular properties and the conditions allowing for escape of high-energy gamma-rays and radio waves, showing that their escape is possible for nebulae with ages of >10-100 yr. In the rapidly-rotating neutron star scenario, we find that radio emission from the quasi-steady nebula itself may be bright enough to be detected especially at sub-mm frequencies, which is relevant as a possible counterpart of pulsar-driven SNe and FRBs. (ii) Second, we consider the fate of bursting emission in the nebulae. We suggest that an impulsive burst may lead to a highly relativistic flow, which would interact with the nebula. If the shocked nebula is still relativistic, pre-existing non-thermal particles in the nebula can be significantly boosted by the forward shock, leading to short-duration (maybe millisecond or longer) high-energy gamma-ray flashes. Possible dissipation at the reverse shock may also lead to gamma-ray emission. (iii) After such flares, interactions with the baryonic ejecta may lead to afterglow emission with a duration of days to weeks. In the magnetar scenario, this burst-in-bubble model leads to the expectation that nearby (<10-100 Mpc) high-energy gamma-ray flashes may be detected by HAWC and CTA, and the subsequent afterglow emission may be seen by radio telescopes such as VLA. (iv) Finally, we discuss several implications specific to FRBs, including constraints on the emission regions and limits on soft gamma-ray counterparts.

The first direct detections of gravitational waves (GWs) from black hole (BH) mergers, GW150914, GW151226 and LVT151012, give a robust lower limit $\sim 70000$ on the number of merged, highly-spinning BHs in our Galaxy. The total spin energy is comparable to all the kinetic energy of supernovae that ever happened in our Galaxy. The BHs release the spin energy to relativistic jets by accreting matter and magnetic fields from the interstellar medium (ISM). By considering the distributions of the ISM density, BH mass and velocity, we calculate the luminosity function of the BH jets, and find that they can potentially accelerate TeV-PeV cosmic-ray particles in our Galaxy with total power $\sim 10^{37\pm 3}$ erg s$^{-1}$ as PeVatrons, positron factories and/or unidentified TeV gamma-ray sources. Additional $\sim 300$ BH jet nebulae could be detectable by CTA (Cherenkov Telescope Array). We also argue that the accretion from the ISM can evaporate and blow away cold material around the BH, which has profound implications for some scenarios to predict electromagnetic counterparts to BH mergers.

We analyze the early growth stage of direct-collapse black holes (DCBHs) with $\sim 10^{5} \ \rm M_\odot$, which are formed by collapse of supermassive stars in atomic-cooling halos at $z \gtrsim 10$. A nuclear accretion disk around a newborn DCBH is gravitationally unstable and fragments into clumps with a few $10 \ \rm M_\odot$ at $\sim 0.01-0.1 \ \rm pc$ from the center. Such clumps evolve into massive population III stars with a few $10-100 \ \rm M_\odot$ via successive gas accretion and a nuclear star cluster is formed. Radiative and mechanical feedback from an inner slim disk and the star cluster will significantly reduce the gas accretion rate onto the DCBH within $\sim 10^6 \ \rm yr$. Some of the nuclear stars can be scattered onto the loss cone orbits also within $\lesssim 10^6 \ \rm yr$ and tidally disrupted by the central DCBH. The jet luminosity powered by such tidal disruption events can be $L_{\rm j} \gtrsim 10^{50} \ \rm erg \ s^{-1}$. The prompt emission will be observed in X-ray bands with a peak duration of $\delta t_{\rm obs} \sim 10^{5-6} \ (1+z) \ \rm s$ followed by a tail $\propto t_{\rm obs}^{-5/3}$, which can be detectable by Swift BAT and eROSITA even from $z \sim 20$. Follow-up observations of the radio afterglows with, e.g., eVLA and the host halos with JWST could probe the earliest AGN feedback from DCBHs.

The recent discovery of the gravitational wave source GW150914 has revealed a coalescing binary black hole (BBH) with masses of $\sim 30~M_\odot$. Previous proposals for the origin of such a massive binary include Population III (PopIII) stars. PopIII stars are efficient producers of BBHs and of a gravitational wave background (GWB) in the $10-100$ Hz band, and also of ionizing radiation in the early Universe. We quantify the relation between the amplitude of the GWB ($\Omega_{\rm gw}$) and the electron scattering optical depth ($\tau_{\rm e}$), produced by PopIII stars, assuming that $f_{\rm esc}\approx 10\%$ of their ionizing radiation escapes into the intergalactic medium. We find that PopIII stars would produce a GWB that is detectable by the future O5 LIGO/Virgo if $\tau_{\rm e} \gtrsim 0.07$, consistent with the recent Planck measurement of $\tau_e=0.055 \pm 0.09$. Moreover, the spectral index of the background from PopIII BBHs becomes as small as ${\rm d}\ln \Omega_{\rm gw}/{\rm d}\ln f\lesssim 0.3$ at $f \gtrsim 30$ Hz, which is significantly flatter than the value $\sim 2/3$ generically produced by lower-redshift and less-massive BBHs. A detection of the unique flattening at such low frequencies by the O5 LIGO/Virgo will indicate the existence of a high-chirp mass, high-redshift BBH population, which is consistent with the PopIII origin. A precise characterization of the spectral shape near $30-50$ Hz by the Einstein Telescope could also constrain the PopIII initial mass function and star formation rate.

Recently, the direct detection of gravitational waves from black hole (BH) mergers was announced by the Advanced LIGO Collaboration. Multi-messenger counterparts of stellar-mass BH mergers are of interest, and it had been suggested that a small disk or celestial body may be involved in the binary of two BHs. To test such possibilities, we consider the fate of a wind powered by an active mini-disk in a relatively short, super-Eddington accretion episode onto a BH with ~10-100 solar masses. We show that its thermal emission could be seen as a fast optical transient with the duration from hours to days. We also find that the coasting outflow forms external shocks due to interaction with the interstellar medium, whose synchrotron emission might be expected in the radio band on a time scale of years. Finally, we also discuss a possible jet component and the associated high-energy neutrino emission as well as ultra-high-energy cosmic-ray acceleration.

Fast-spinning strongly magnetized newborn neutron stars, including nascent magnetars, are popularly implemented as the engine of luminous stellar explosions. Here, we consider the scenario that they power various stripped-envelope supernovae, not only super-luminous supernovae Ic but also broad-line supernovae Ibc and possibly some ordinary supernovae Ibc. This scenario is also motivated by the hypothesis that Galactic magnetars largely originate from fast-spinning neutron stars as remnants of stripped-envelope supernovae. By consistently modeling the energy injection from magnetized wind and Ni decay, we show that proto-neutron stars with >~ 10 ms rotation and B_dip >~ 5 x 10^14 G can be harbored in ordinary supernovae Ibc. On the other hand, millisecond proto-neuton stars can solely power broad-line supernovae Ibc if they are born with poloidal magnetic field of B_dip >~ 5 x 10^14 G, and superluminous supernovae Ic with B_dip >~ 10^13 G. Then, we study how multi-messenger emission can be used to discriminate such pulsar-driven supernova models from other competitive scenarios. First, high-energy x-ray and gamma-ray emission from embryonic pulsar wind nebulae is a promising smoking gun of the underlying newborn pulsar wind. Follow-up observations of stripped-envelope supernovae using NuSTAR ~ 50-100 days after the explosion is strongly encouraged for nearby objects. We also discuss possible effects of gravitational-waves on the spin-down of proto-neutron stars. If millisecond proto-neutron stars with B_dip <~ a few x 10^13 G emit gravitational waves through e.g., non-axisymmetric rotation deformed by the inner toroidal fields of B_t >~ 10^16 G, the gravitational wave signal can be detectable from ordinary supernova Ibc in the Virgo cluster by Advanced LIGO, Advanced Virgo, and KAGRA.

Extended emission is a mystery in short gamma-ray bursts (SGRBs). By making time resolved spectral analyses of brightest nine events observed by ${\it Swift}$ XRT, we obviously classify the early X-ray emission of SGRBs into two types. One is the extended emission with exponentially rapid decay, which shows significant spectral softening during hundreds seconds since the SGRB trigger and is also detected by ${\it Swift}$-BAT. The other is a dim afterglow only showing power-law decay over $10^4$ s. The correlations between the temporal decay and spectral indices of the extended emissions are inconsistent with the $\alpha$-$\beta$ correlation expected for the high-latitude curvature emission from a uniform jet. The observed too-rapid decay suggests the emission from a photosphere or a patchy surface, and manifests the stopping central engine via such as magnetic reconnection at the black hole.

We propose high-velocity collisions of protogalaxies as a new pathway to form supermassive stars (SMSs) with masses of ~ 10^5 Msun at high redshift (z > 10). When protogalaxies hosted by dark matter halos with a virial temperature of ~ 10^4 K collide with a relative velocity > 200 km/s, the gas is shock-heated to ~ 10^6 K and subsequently cools isobarically via free-free emission and He^+, He, and H line emission. Since the gas density (> 10^4 cm^{-3}) is high enough to destroy H_2 molecules by collisional dissociation, the shocked gas never cools below ~ 10^4 K. Once a gas cloud of ~ 10^5 Msun reaches this temperature, it becomes gravitationally unstable and forms a SMS which will rapidly collapse into a super massive black hole (SMBH) via general relativistic instability. We perform a simple analytic estimate of the number density of direct-collapse black holes (DCBHs) formed through this scenario (calibrated with cosmological N-body simulations) and find n_{DCBH} ~ 10^{-9} Mpc^{-3} (comoving) by z = 10. This could potentially explain the abundance of bright high-z quasars.

We examine the bright radio synchrotron counterparts of low-luminosity gamma-ray bursts (llGRBs) and relativistic supernovae (SNe) and find that they can be powered by spherical hypernova (HN) explosions. Our results imply that radio-bright HNe are driven by relativistic jets that are choked deep inside the progenitor stars or quasi-spherical magnetized winds from fast-rotating magnetars. We also consider the optical synchrotron counterparts of radio-bright HNe and show that they can be observed as precursors several days before the SN peak with an r-band absolute magnitude of M_r ~ -14 mag. While previous studies suggested that additional trans-relativistic components are required to power the bright radio emission, we find that they overestimated the energy budget of the trans-relativistic component by overlooking some factors related to the minimum energy of non-thermal electrons. If an additional trans-relativistic component exists, then a much brighter optical precursor with M_r ~ -20 mag can be expected. Thus, the scenarios of radio-bright HNe can be distinguished by using optical precursors, which can be detectable from < 100 Mpc by current SN surveys like the Kiso SN Survey, Palomar Transient Factory, and Panoramic Survey Telescope & Rapid Response System.

Newborn black holes in collapsing massive stars can be accompanied by a fallback disk. The accretion rate is typically super-Eddington and strong disk outflows are expected. Such outflows could be directly observed in some failed explosions of compact (blue supergiants or Wolf-Rayet stars) progenitors, and may be more common than long-duration gamma-ray bursts. Using an analytical model, we show that the fallback disk outflows produce blue UV-optical transients with a peak bolometric luminosity of ~10^(42-43) erg s^-1 (peak R-band absolute AB magnitudes of -16 to -18) and an emission duration of ~ a few to ~ 10 days. The spectra are likely dominated intermediate mass elements, but will lack much radioactive nuclei and iron-group elements. The above properties are broadly consistent with some of the rapid blue transients detected by Pan-STARRS and PTF. This scenario can be distinguished from alternative models using radio observations within a few years after the optical peak.

It has been suggested that some classes of luminous supernovae (SNe) and gamma-ray bursts (GRBs) are driven by newborn magnetars. Fast-rotating proto-neutron stars have also been of interest as potential sources of gravitational waves (GWs). We show that for a range of rotation periods and magnetic fields, hard X rays and GeV gamma rays provide us with a promising probe of pulsar-aided SNe. It is observationally known that young pulsar wind nebulae (PWNe) in the Milky Way are very efficient lepton accelerators. We argue that, if embryonic PWNe satisfy similar conditions at early stages of SNe (in ~1-10 months after the explosion), external inverse-Compton emission via upscatterings of SN photons is naturally expected in the GeV range as well as broadband synchrotron emission. To fully take into account the Klein-Nishina effect and two-photon annihilation process that are important at early times, we perform detailed calculations including electromagnetic cascades. Our results suggest that hard X-ray telescopes such as NuSTAR can observe such early PWN emission by followup observations in months-to-years. GeV gamma rays may also be detected by Fermi for nearby SNe, which serve as counterparts of these GW sources. Detecting the signals will give us an interesting probe of particle acceleration at early times of PWNe, as well as clues to driving mechanisms of luminous SNe and GRBs. Since the Bethe-Heitler cross section is lower than the Thomson cross section, gamma rays would allow us to study subphotospheric dissipation. We encourage searches for high-energy emission from nearby SNe, especially Type Ibc SNe including super-luminous objects.

We investigate the possible origin of extended emissions (EEs) of short gamma-ray bursts with an isotropic energy of ~ 10^(50-51) erg and a duration of a few 10 s to ~ 100 s, based on a compact binary (neutron star (NS)-NS or NS-black hole (BH)) merger scenario. We analyze the evolution of magnetized neutrino-dominated accretion disks of mass ~ 0.1 M_sun around BHs formed after the mergers, and estimate the power of relativistic outflows via the Blandford-Znajek (BZ) process. We show that a rotation energy of the BH up to > 10^52 erg can be extracted with an observed time scale of > 30 (1+z) s with a relatively small disk viscosity parameter of alpha < 0.01. Such a BZ power dissipates by clashing with non-relativistic pre-ejected matter of mass M ~ 10^-(2-4) M_sun, and forms a mildly relativistic fireball. We show that the dissipative photospheric emissions from such fireballs are likely in the soft X-ray band (1-10 keV) for M ~ 10^-2 M_sun possibly in NS-NS mergers, and in the BAT band (15-150 keV) for M ~ 10^-4 M_sun possibly in NS-BH mergers. In the former case, such soft EEs can provide a good chance of ~ 6 yr^-1 for simultaneous detections of the gravitational waves with a ~ 0.1 deg angular resolution by soft X-ray survey facilities like Wide-Field MAXI.

We investigate the shock acceleration of particles in massive galaxy mergers or collisions, and show that cosmic rays (CRs) can be accelerated up to the second knee energy ~0.1-1 EeV and possibly beyond, with a hard spectral index Gamma ~ 2. Such CRs lose their energy via hadronuclear interactions within a dynamical timescale of the merger shock, producing gamma rays and neutrinos as a by-product. If ~ 10 % of the shock dissipated energy goes into CR acceleration, some local merging galaxies will produce gamma-ray counterparts detectable by CTA. Also, based on the concordance cosmology, where a good fraction of the massive galaxies experience a major merger in a cosmological timescale, the neutrino counterparts can constitute ~ 20-60 % of the isotropic background detected by IceCube.

Recently, Thornton et al. reported the detection of four fast radio bursts (FRBs). The dispersion measures indicate that the sources of these FRBs are at cosmological distance. Given the large full sky event rate ~ 10^4 sky^-1 day^-1, the FRBs are a promising target for multi-messenger astronomy. Here we propose double degenerate, binary white-dwarf (WD) mergers as the source of FRBs, which are produced by coherent emission from the polar region of a rapidly rotating, magnetized massive WD formed after the merger. The basic characteristics of the FRBs, such as the energetics, emission duration and event rate, can be consistently explained in this scenario. As a result, we predict that some FRBs can accompany type Ia supernovae (SNe Ia) or X-ray debris disks. Simultaneous detection could test our scenario and probe the progenitors of SNe Ia, and moreover would provide a novel constraint on the cosmological parameters. We strongly encourage future SN and X-ray surveys that follow up FRBs.

We study a type of particle acceleration that operates via neutron-proton conversion in inelastic nuclear collisions. This mechanism can be expected for relativistic shocks at subphotospheres if relativistic outflows contain neutrons. Using a test-particle approximation, we numerically calculate the energy spectrum and the efficiency of accelerated particles, and show that a good energy fraction of the nucleons can be accelerated. This mechanism may especially be relevant if the shock is radiation-mediated, and it would enhance the detectability of GeV-TeV neutrinos.

Long GRBs (LGRBs) have typical duration of ~ 30 s and some of them are associated with hypernovae, like Type Ic SN 1998bw. Wolf-Rayet stars are the most plausible LGRB progenitors, since the free-fall time of the envelope is consistent with the duration, and the natural outcome of the progenitor is a Type Ic SN. While a new population of ultra-long GRBs (ULGRBs), GRB 111209A, GRB 101225A, and GRB 121027A, has a duration of ~ 10^4 s, two of them are accompanied by superluminous-supernova (SLSN) like bumps, which are <~ 10 times brighter than typical hypernovae. Wolf-Rayet progenitors cannot explain ULGRBs because of too long duration and too bright SN-like bump. A blue supergiant (BSG) progenitor model, however, can explain the duration of ULGRBs. Moreover, SLSN-like bump can be attributed to the so-called cocoon-fireball photospheric emissions (CFPEs). Since a large cocoon is inevitably produced during the relativistic jet piercing though the BSG envelope, this component can be a smoking-gun evidence of BSG model for ULGRBs. In this paper, we examine u, g, r, i, and J-band light curves of three ULGRBs and demonstrate that they can be fitted quite well by our BSG model with the appropriate choices of the jet opening angle and the number density of the ambient gas. In addition, we predict that for 121027A, SLSN-like bump could have been observed for ~ 20 - 80 days after the burst. We also propose that some SLSNe might be CFPEs of off-axis ULGRBs without visible prompt emission.

Relativistic outflows with neutrons inevitably lead to inelastic collisions, and resulting subphotospheric gamma rays may explain prompt emission of gamma-ray bursts. In this model, hadronuclear, quasithermal neutrinos in the 10-100 GeV range should be generated, and they may even have a high-energy tail by neutron-proton-converter or shock acceleration mechanisms. We demonstrate the importance of dedicated searches with DeepCore+IceCube, though such analyses have not been performed. Successful detections enable us to discriminate among prompt emission mechanisms, probe the jet composition, and see roles of relativistic neutrons as well as effects of cosmic-ray acceleration.

The IceCube collaboration discovery of 28 high-energy neutrinos over the energy range 30 TeV <~ E_nu <~ 1 PeV, a 4.3-sigma excess over expected backgrounds, represents the first high-confidence detection of cosmic neutrinos at these energies. In light of this discovery, we explore the possibility that some of the Sub-PeV cosmic neutrinos might originate in our Galaxy's TeV unidentified (TeV UnID) sources. While typically resolved at TeV energies, these sources lack prominent radio or X-ray counterparts, and so have been considered promising sites for hadron acceleration within our Galaxy. Modeling the TeV UnID sources as Galactic hypernova remnants, we predict Sub-PeV neutrino fluxes and spectra consistent with their contributing a minority of n_nu <~ 2 of the observed events. This is consistent with our analysis of the spatial distribution of the Sub-PeV neutrinos and TeV UnID sources, which finds that a best-fit of one, and maximum of 3.8 (at 90%-confidence), of the ~16 non-atmospheric Sub-PeV neutrinos may originate in TeV UnID sources, with the remaining 75% to 95% of events being drawn from an isotropic background. If our scenario is correct, we expect excess Sub-PeV neutrinos to accumulate along the Galactic plane, within |l| <~ 30 deg of the Galactic center and in the Cygnus region, as observations by IceCube and other high-energy neutrino facilities go forward. Our scenario also has implications for radio, X-ray, and TeV observations of the TeV UnID sources.

The recent gamma-ray burst GRB 130427A has an isotropic electromagnetic energy $E^{iso}\sim10^{54}$ erg, suggesting an ample supply of target photons for photo-hadronic interactions, which at its low redshift of $z\sim0.34$ would appear to make it a promising candidate for neutrino detection. However, the IceCube collaboration has reported a null result based on a search during the prompt emission phase. We show that this neutrino non-detection can provide valuable information about this GRB's key physical parameters such as the emission radius $R_{d}$, the bulk Lorentz factor $\Gamma$ and the energy fraction converted into cosmic rays $\epsilon_{p}$. The results are discussed both in a model-independent way and in the specific scenarios of an internal shock model (IS), a baryonic photospheric model (BPH) and magnetic photospheric model (MPH). We find that the constraints are most stringent for the magnetic photospheric model considered, but the constraints on the internal shock and the baryonic photosphere models are fairly modest.

Metal-poor massive stars may typically end up their lives as blue supergiants (BSGs). Gamma-ray bursts (GRBs) from such progenitors could have ultra-long duration of relativistic jets. For example Population III (Pop III) GRBs at z ~ 10-20 might be observable as X-ray rich events with a typical duration of T_90 ~ 10^4(1+z) sec. Recent GRB111209A at z = 0.677 has an ultra long duration of T_90 ~ 2.5*10^4 sec so that it have been suggested that the progenitor might be a metal-poor BSGs in the local universe. Here, we suggest luminous UV/optical/infrared emissions associated with such a new class of GRB from metal poor BSGs. Before the jet head breaks out the progenitor envelope, the energy injected by the jet is stored in a hot-plasma cocoon, which finally emerges and expands as a baryon-loaded fireball. We show that the photospheric emissions from the cocoon fireball could be intrinsically very bright (L_peak ~ 10^(42-44) erg/sec) in UV/optical bands (E_peak ~ 10 eV) with a typical duration of ~ 100 days in the rest frame. Such cocoon emissions from Pop III GRB might be detectable in infrared bands at ~ years after Pop III GRBs at up to z ~ 15 by up-coming facilities like JWST. We also suggest that GRB111209A might have been rebrightening in UV/optical bands up to an AB magnitude of < 26. The cocoon emissions from local metal-poor BSGs might have been already observed as luminous supernovae without GRB since they can be seen from the off-axis direction of the jet.

Relativistic shocks that accompany supernovae (SNe) produce X-ray burst emissions as they break out in the dense circumstellar medium around the progenitors. This phenomenon is sometimes associated with peculiar low-luminosity gamma-ray bursts (LL GRBs). Here, we investigate the high energy neutrino and gamma-ray counterparts of such a class of SNe. Just beyond the shock breakout radius, particle acceleration in the collisionless shock starts to operate in the presence of breakout photons. We show that protons may be accelerated to sufficiently high energies and produce high energy neutrinos and gamma rays via the photomeson interaction. These neutrinos and gamma rays may be detectable from 10 Mpc away by IceCube/KM3Net as multi-TeV transients almost simultaneously with the X-ray burst emission, and even from 100 Mpc away with follow-up observations by CTA using a wide-field sky monitor like Swift as a trigger. A statistical technique using a stacking approach could also be possible for the detection, with the aid of the SN optical/infrared counterparts. Such multi-messenger observations offer the possibility to probe the transition of relativistic shocks from radiation-mediated to collisionless ones, and would also constrain the mechanisms of particle acceleration and emission in LL GRBs.

Recent numerical simulations suggest that Population III (Pop III) stars were born with masses not larger than $\sim 100 M_{\odot}$ but typically $\sim 40M_{\odot}$. By self-consistently considering the jet generation and propagation in the envelope of these low mass Pop III stars, we find that a Pop III blue super giant star has the possibility to raise a gamma-ray burst (GRB) even though it keeps a massive hydrogen envelope. We evaluate observational characters of Pop III GRBs and predict that Pop III GRBs have the duration of $\sim 10^5$ sec in the observer frame and the peak luminosity of $\sim 5 \times 10^{50} {\rm erg} {\rm sec}^{-1}$. Assuming that the $E_p-L_p$ (or $E_p-E_{\gamma, \rm iso}$) correlation holds for Pop III GRBs, we find that the spectrum peak energy falls $\sim$ a few keV (or $\sim 100$ keV) in the observer frame. We discuss the detectability of Pop III GRBs by future satellite missions such as EXIST and Lobster. If the $E_p-E_{\gamma, \rm iso}$ correlation holds, we have the possibility to detect Pop III GRBs at $z \sim 9$ as long duration X-ray rich GRBs by EXIST. On the other hand, if the $E_p-L_p$ correlation holds, we have the possibility to detect Pop III GRBs up to $z \sim 19$ as long duration X-ray flashes by Lobster.

We investigate the capability of pulsar timing arrays (PTAs) as a probe of primordial black holes (PBHs), which might constitute the Galactic dark matter. A PBH passing nearby the Earth or a pulsar gives an impulse acceleration and induces residuals on otherwise orderly pulsar timing data. We show that the timing residuals induced at pulsars are optimal for searching heavier PBHs than those at the Earth, and the two probes are highly complemental. Future facilities like SKA could detect PBHs with masses around 10^(22-28) g even if only a small fraction (< 1%) of the Galactic dark matter consists of these PBHs.

Magnetic flares and induced oscillations of magnetars (super-magnetized neutron stars) are promising sources of gravitational waves (GWs). We suggest that the GW emission, if any, would last longer than the observed X-ray quasi-periodic oscillations (X-QPOs), calling for the longer-term GW analyses from a day to months than the current searches. Like the pulsar timing, the oscillation frequency would also evolve with time because of the decay or reconfiguration of magnetic field, which is crucial for the GW detection. With the observed GW frequency and its time-derivatives, we can probe the interior magnetic field strength and its evolution to open a new GW asteroseismology with the next generation interferometers like advanced LIGO, advanced VIRGO, LCGT and ET.

One of the most essential but uncertain processes for producing cosmic-rays (CRs) and their spectra is how accelerated particles escape into the interstellar space. We propose that the CR electron spectra at >~TeV energy can provide a direct probe of the CR escape complementary to the CR nuclei and gamma-rays. We calculate the electron spectra from a young pulsar embedded in the supernova remnant (SNR), like Vela, taking into account the energy-dependent CR escape. SNRs would accelerate and hence confine particles with energy up to 10^{15.5}eV. Only energetic particles can escape first, while the lower energy particles are confined and released later. Then the observed electron spectrum should have a low energy cutoff whose position marks the age of the pulsar/SNR. The low energy cutoff is observable in the $\gtrsim {\rm TeV}$ energy window, where other contaminating sources are expected to be few due to the fast cooling of electrons. The spectrum looks similar to a dark matter annihilation line if the low energy cutoff is close to the high energy intrinsic or cooling break. The future experiments such as CALET and CTA are capable of directly detecting the CR escape features toward revealing the origin of CRs.

We suggest that white dwarf (WD) pulsars can compete with neutron star (NS) pulsars for producing the excesses of cosmic ray electrons/positrons observed by the PAMELA, ATIC/PPB-BETS, Fermi and HESS experiments. A merger of two WDs leads to a rapidly spinning WD with a rotational energy comparable to the NS case. The birth rate is also similar, providing the right energy budget for the cosmic ray electrons/positrons. Applying the NS theory, we suggest that the WD pulsars can in principle produce electrons/positrons up to 10 TeV. In contrast to the NS model, the adiabatic and radiative energy losses of electrons/positrons are negligible since their injection continues after the expansion of the pulsar wind nebula, and hence it is enough that a fraction 1% of WDs are magnetized as observed. The long activity also increases the number of nearby sources, which reduces the Poisson fluctuation in the flux. The WD pulsars could dominate the quickly cooling electrons/positrons above TeV energy as a second spectral bump or even surpass the NS pulsars in the observing energy range 10 GeV - 1 TeV, providing a background for the dark matter signals and a nice target for the future AMS-02, CALET and CTA experiment.