Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

We examine the sensitivity of neutrino emissions to stellar evolution models for a 15$M_\odot$ progenitor, paying particular attention to a phase prior to the collapse. We demonstrate that the number luminosities in both electron-type neutrinos ($\nu_e$) and their anti-partners ($\bar{\nu}_e$) differ by more than an order of magnitude by changing spatial resolutions and nuclear network sizes on stellar evolution models. We also develop a phenomenological model to capture the essential trend of the diversity, in which neutrino luminosities are expressed as a function of central density, temperature and electron fraction. In the analysis, we show that neutrino luminosity can be well characterized by these central quantities. This analysis also reveals that the most influential quantity to the time evolution of $\nu_e$ luminosity is matter density, while it is temperature for $\bar{\nu}_e$. These qualitative trends will be useful and applicable to constrain the physical state of progenitors at the final stages of stellar evolution from future neutrino observations, although more detailed systematic studies including various mass progenitors are required to assess the applicability.

We report, for the first time, the detection of the Mn-K$\alpha$ line in the Type IIb supernova (SN IIb) remnant, Cassiopeia A. Manganese ($^{55}$Mn after decay of $^{55}$Co), a neutron-rich element, together with chromium ($^{52}$Cr after decay of $^{52}$Fe), is mainly synthesized at the explosive incomplete Si burning regime. Therefore, the Mn/Cr mass ratio with its neutron excess reflects the neutronization at the relevant burning layer during the explosion. Chandra's archival X-ray data of Cassiopeia A indicate a low Mn/Cr mass ratio with values in the range 0.10--0.66, which, when compared to one-dimensional SN explosion models, requires that the electron fraction be 0.4990 $\lesssim Y_{\rm e} \lesssim$ 0.5 at the incomplete Si burning layer. An explosion model assuming a solar-metallicity progenitor with a typical explosion energy ($1 \times 10^{51}$ erg) fails to reproduce such a high electron fraction. In such models, the explosive Si-burning regime extends only to the Si/O layer established during the progenitor's hydrostatic evolution; the $Y_e$ in the Si/O layer is lower than the value required by our observational constraints. We can satisfy the observed Mn/Cr mass ratio if the explosive Si-burning regime were to extend into the O/Ne hydrostatic layer, which has a higher $Y_{\rm e}$. This would require an energetic ($> 2 \times 10^{51}$ erg) and/or asymmetric explosion of a sub-solar metallicity progenitor ($Z \lesssim 0.5Z_{\odot}$) for Cassiopeia A. The low initial metallicity can be used to rule out a single-star progenitor, leaving the possibility of a binary progenitor with a compact companion (white dwarf, neutron star or black hole). We discuss the detectability of X-rays from Bondi accretion onto such a compact companion around the explosion site. We also discuss other possible mass-loss scenarios for the progenitor system of Cassiopeia A.

Type Ibn supernovae (SNe Ibn) are intriguing stellar explosions whose spectra exhibit narrow helium lines with little hydrogen. They trace the presence of circumstellar material (CSM) formed via pre-SN eruptions of their stripped-envelope progenitors. Early work has generally assumed that SNe Ibn come from massive Wolf-Rayet (WR) stars via single star evolution. In this paper, we report ultraviolet (UV) and optical observations of two nearby Type Ibn SNe 2006jc and 2015G conducted with the Hubble Space Telescope (HST) at late times. A point source is detected at the position of SN 2006jc, and we confirm the conclusion of Maund et al. that it is the progenitor's binary companion. Its position on the Hertzsprung-Russell (HR) diagram corresponds to a star that has evolved off the main sequence (MS); further analysis implies a low initial mass for the companion star ($M_2$ $\le$ 11.9$^{+1.2}_{-0.8}$ $M_\odot$) and a secondary-to-primary initial mass ratio very close to unity ($q$ = $M_2/M_1$ $\sim$ 1); the SN progenitor's hydrogen envelope had been stripped through binary interaction. We do not detect the binary companion of SN 2015G. For both SNe, the surrounding stellar populations have relatively old ages and argue against any massive WR stars as their progenitors. These results suggest that SNe Ibn may have lower-mass origins in interacting binaries. As a result, they also provide evidence that the giant eruptions commonly seen in massive luminous blue variables (LBVs) can also occur in much lower-mass, stripped-envelope stars just before core collapse.

Motivated by the eruptive mass loss inferred from Luminous Blue Variable (LBV) stars, we present 1D hydrodynamical simulations of the response from sudden energy injection into the interior of a very massive ($100 \, M_\odot$) star. For a fiducial case with total energy addition set to a factor $f=0.5$ of the net stellar binding energy, and applied within the stellar envelope, we detail the dynamical response that leads to ejection of the outermost $7.2 \, M_\odot$. We find that the ejecta's variations in time $t$ and radius $r$ for the velocity $v$, density $\rho$, and temperature $T$ are quite well fit by similarity forms in the variable $r/t \approx v$. Specifically the scaled density follows a simple exponential decline $\rho t^{3} \sim \exp (-r/v_{\rm o} t)$. This `exponential similarity' leads to analytic scaling relations for total ejecta mass $\Delta M$ and kinetic energy $\Delta K$ that agree well with the hydrodynamical simulations, with the specific-energy-averaged speed related to the exponential scale speed $v_{\rm o}$ through ${\bar v} \equiv \sqrt{2 \Delta K/\Delta M} = \sqrt{12} \, v_{\rm o}$, and a value comparable to the star's surface escape speed, $v_{\rm esc}$. Models with energy added in the core develop a surface shock breakout that propels an initial, higher-speed ejecta ($>$5000km s$^{-1}$), but the bulk of the ejected material still follows the same exponential similarity scalings with ${\bar v} \approx v_{\rm esc}$. A broader parameter study examines how the ejected mass and energy depends on the energy-addition factor $f$, for three distinct model series that locate the added energy in either the core, envelope, or near-surface. We conclude by discussing the relevance of these results for understanding LBV outbursts and other eruptive phenomena, such as failed supernovae and pulsational pair instability events.

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method.

We carry out a comprehensive study of supernova ejecta-companion interaction in massive binary systems. We aim to physically understand the kinematics of the interaction and predict observational signatures. To do this we perform simulations over a vast parameter space of binary configurations, varying the masses of the progenitor and companion, structure of the companion, explosion energy, and orbital separation. Our results were not so consistent with classical models by Wheeler et al. 1975, sometimes deviating by an order of magnitude. We construct an alternative simple model which explains the simulated results reasonably well and can be used to estimate impact velocities for arbitrary explosion profiles and companion star structures. We then investigate the long term evolution after the supernova, where the companion can be inflated by the energy injected into the star. We find that the companion can become more than an order of magnitude overluminous straight after the supernova, but quickly fades away after $\sim10$ years and returns to its original luminosity in about a thermal timescale of the star. Finally, we also discuss the possible surface contamination of heavy elements from the slower ejecta.

This letter complements a formation scenario of the progenitor of the supernova iPTF13bvn proposed in Hirai (2017). Although the scenario was successful in reproducing various observational features of the explosion and pre-explosion photometry by assuming that the progenitor had a relatively large black hole companion, it lacked an explanation for the origin of the black hole itself. We now explore the possible evolutionary paths towards this binary with a relatively large black hole companion. We found that the black hole was probably produced by a very massive star that experienced common envelope evolution. According to our MESA stellar models, the primary mass should have been 70M to reproduce the required remnant mass and final separation. This indicates that iPTF13bvn was likely a rare case and normal type Ib supernovae originate from different paths.

We study formation of stellar mass binary black holes (BBHs) originating from Population III (PopIII) stars, performing stellar evolution simulations for PopIII binaries with MESA. We find that a significant fraction of PopIII binaries form massive BBHs through stable mass transfer between two stars in a binary, without experiencing common envelope phases. We investigate necessary conditions required for PopIII binaries to form BBHs coalescing within the Hubble time with a semi-analytical model calibrated by the stellar evolution simulations. The formation efficiency of coalescing PopIII BBHs is estimated for two different initial conditions for PopIII binaries with large and small separations, respectively. Consequently, in both models, $\sim 10\%$ of the total PopIII binaries form BBHs only through stable mass transfer and $\sim 10\%$ of these BBHs merge due to gravitational wave emission within the Hubble time. Furthermore, the chirp mass of merging BBHs has a flat distribution over $15\lesssim M_{\rm chirp}/M_\odot \lesssim 35$. This formation pathway of PopIII BBHs is presumably robust because stable mass transfer is less uncertain than common envelope evolution, which is the main formation channel for Population II BBHs. We also test the hypothesis that the BBH mergers detected by LIGO originate from PopIII stars using our result and the total number of PopIII stars formed in the early universe as inferred from the optical depth measured by Planck. We conclude that the PopIII BBH formation scenario can explain the mass-weighted merger rate of the LIGO's O1 events with the maximal PopIII formation efficiency inferred from the Planck measurement, even without BBHs formed by unstable mass transfer or common envelope phases.

The formation scenario of the progenitor of iPTF13bvn has been revisited. iPTF13bvn is unique in the sense that a corresponding pre-supernova image has been identified. This has enabled us to strongly constrain the nature of its progenitor. From the pre-supernova image, light curve and the latest observations, it is currently widely accepted that the progenitor of iPTF13bvn was in a binary system. The fast decline in the light curve suggests a progenitor mass of $\sim3.5M_\odot$, and the upper limit on the remaining companion is $\sim20M_\odot$. Recent works suggest that binary evolution models involving common envelope episodes can satisfy all the observational constraints. We have examined the common envelope scenario based on latest knowledges on common envelope evolution. We have found that the common envelope scenario for the progenitor of iPTF13bvn seems not to be suitable since it can not explain the pre-supernova radius. We also propose an alternative model with a black hole companion. Stellar evolution calculations with a large black hole companion were carried out and succeeded in satisfying all observational constraints. However, more studies should be carried out to explain the origin of the large black hole companion.

A new computationally efficient method has been introduced to treat self-gravity in mesh based hydrodynamical simulations. It is applied simply by slightly modifying the Poisson equation into an inhomogeneous wave equation. This roughly corresponds to the weak field limit of the Einstein equations in general relativity, and as long as the gravitation propagation speed is taken to be larger than the hydrodynamical characteristic speed, the results agree with solutions for the Poisson equation. The solutions almost perfectly agree if the domain is taken large enough, or appropriate boundary conditions are given. Our new method can not only significantly reduce the computational time compared with existent methods, but is also fully compatible with massive parallel computation, nested grids and adaptive mesh refinement techniques, all of which can accelerate the progress in computational astrophysics and cosmology.

We investigate the possible effects of the supernova ejecta hitting the companion star in iPTF 13bvn, focusing on the observable features when it becomes visible. iPTF 13bvn is a type Ib supernova that may become the first case that its progenitor is identified as a binary by near future observations. According to calculations by Bersten et al. (2014), the progenitor should have a mass $\approx3.5M_\odot$ to reproduce the supernova light curve, and such compact stars could only be produced via binary evolution. This is one of the reasons that we expect the progenitor to be a binary, but it should be confirmed by observing the remaining companion after the supernova. Their evolutionary calculations suggest that the companion star will be an overluminous OB star at the moment of supernova. With a combination of hydrodynamical and evolutionary simulations, we find that the secondary star will be heated by the supernova ejecta and expand to have larger luminosities and lower surface effective temperatures. The star will look rather like a red super giant, and this should be taken into account when searching for the companion star in the supernova ejecta in future observations.

The majority of massive stars are formed in binary systems. It is hence reasonable to expect that most core-collapse supernovae (CCSNe) take place in binaries and the existence of a companion star may leave some imprints in observed features. Having this in mind, we have conducted two-dimensional hydrodynamical simulations of the collisions of CCSNe ejecta with the companion star in an almost-equal-mass ($\sim 10M_\odot$) binary to find out possible consequences of such events. In particular we pay attention to the amount of mass removed and its dependence on the binary separation. In contrast to the previous surmise, we find that the companion mass is stripped not by momentum transfer but by shock heating. Up to $25\%$ of the original mass can be removed for the closest separations and the removed mass decreases as $M_{ub} \propto a^{-4.3}$ with the binary separation $a$. By performing some experimental computations with artificially-modified densities of incident ejecta, we show that if the velocity of ejecta is fixed, the density of incident ejecta is the single important parameter that actually determines the removed mass as $M_{ub} \propto \rho_{ej} ^{1.4} $. On the other hand, another set of simulations with modified velocities of incident ejecta demonstrate that the strength of the forward shock, which heats up the stellar material and causes the mass loss of the companion star, is actually the key parameter for the removed mass.