BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kenta Hotokezaka

Kenta Hotokezaka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRastro-ph.GAastro-ph.COhep-exhep-phnucl-exnucl-thphysics.atom-ph

Catalog footprint

What is connected

35works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

Electromagnetic counterparts of binary neutron star mergers leading to a strongly magnetized long-lived remnant neutron star

We explore the electromagnetic counterparts that will associate with binary neutron star mergers for the case that remnant massive neutron stars survive for $\gtrsim 0.5\,$s after the merger. For this study, we employ the outflow profiles obtained by long-term general-relativistic neutrino-radiation magneto-hydrodynamics simulations with a mean field dynamo effect. We show that a synchrotron afterglow with high luminosity can be associated with the merger event if the magnetic fields of the remnant neutron stars are significantly amplified by the dynamo effect. We also perform a radiative transfer calculation for kilonovae and find that for the highly amplified magnetic field cases, the kilonovae can be bright in the early epoch, while it shows the optical emission rapid declining in a few days and the long-lasting ($\sim 10\,{\rm d}$) emission very bright in the near-infrared wavelength. All these features have not been found in GW170817, indicating that the merger remnant neutron star formed in GW170817 might have collapsed to a black hole within several hundreds ms or magnetic-field amplification might be a minor effect.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lanthanide Features in Near-infrared Spectra of Kilonovae

The observations of GW170817/AT2017gfo have provided us with evidence that binary neutron star mergers are sites of $r$-process nucleosynthesis. However, the observed signatures in the spectra of GW170817/AT2017gfo have not been fully decoded especially in the near-infrared (NIR) wavelengths. In this paper, we investigate the kilonova spectra over the entire wavelength range with the aim of elemental identification. We systematically calculate the strength of bound-bound transitions by constructing a hybrid line list that is accurate for important strong transitions and complete for weak transitions. We find that the elements on the left side of the periodic table, such as Ca, Sr, Y, Zr, Ba, La, and Ce, tend to produce prominent absorption lines in the spectra. This is because such elements have a small number of valence electrons and low-lying energy levels, resulting in strong transitions. By performing self-consistent radiative transfer simulations for the entire ejecta, we find that La III and Ce III appear in the NIR spectra, which can explain the absorption features at $λ\sim 12000$-14000 A in the spectra of GW170817/AT2017gfo. The mass fractions of La and Ce are estimated to be $>2\times 10^{-6}$ and $\sim$ (1-100)$\times 10^{-5}$, respectively. An actinide element Th can also be a source of absorption as the atomic structure is analogous to that of Ce. However, we show that Th III features are less prominent in the spectra because of the denser energy levels of actinides compared to those of lanthanides.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tungsten vs Selenium as a potential source of kilonova nebular emission observed by Spitzer

Infrared emission lines arising from transitions between fine structure levels of heavy elements are expected to produce kilonova nebular emission. For the kilonova in GW170817, strong emission at 4.5 ${\rm μm}$ at late times was detected by the Spitzer Space Telescope but no source was detected at 3.6 ${\rm μm}$. This peculiar spectrum indicates that there exist strong line emitters around 4.5 ${\rm μm}$ and the absence of strong lines around 3.6 ${\rm μm}$. To model the spectrum we prepare a line list based on the selection rules in LS coupling from the experimentally calibrated energy levels in the NIST database. This method enables to generate the synthetic spectra with accurate line wavelengths. We find that the spectrum is sensitive to the abundance pattern whether or not the first r-process peak elements are included. In both cases, the synthetic spectra can match the observed data, leading to two possible interpretations. If the first peak elements are abundant a Se III line dominates the flux. If otherwise, W III with Os III, Rh III, and Ce IV can be the main sources. Observing nebular spectra for the future kilonovae in a wider wavelength range can provide more conclusive elemental identification.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A transient radio source consistent with a merger-triggered core collapse supernova

A core-collapse supernova occurs when exothermic fusion ceases in the core of a massive star, typically due to exhaustion of nuclear fuel. Theory predicts that fusion could be interrupted earlier, by merging of the star with a compact binary companion. We report a luminous radio transient, VT J121001+495647, found in the Very Large Array Sky Survey. The radio emission is consistent with supernova ejecta colliding with a dense shell of material, potentially ejected by binary interaction in the centuries prior to explosion. We associate the supernova with an archival X-ray transient, which implies a relativistic jet was launched during the explosion. The combination of an early relativistic jet and late-time dense interaction is consistent with expectations for a merger-driven explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-thermal neutrinos created by shock acceleration in successful and failed core-collapse supernova

We present a comprehensive study of neutrino shock acceleration in core-collapse supernova (CCSN). The leading players are heavy leptonic neutrinos, $ν_μ$ and $ν_τ$; the former and latter potentially gain the energy up to $\sim 100$ MeV and $\sim 200$ MeV, respectively, through the shock acceleration. Demonstrating the neutrino shock acceleration by Monte Carlo neutrino transport, we make a statement that it commonly occurs in the early post bounce phase ($\lesssim 50$ ms after bounce) for all massive stellar collapse experiencing nuclear bounce and would reoccur in the late phase ($\gtrsim 100$ ms) for failed CCSNe. This opens up a new possibility to detect high energy neutrinos by terrestrial detectors from Galactic CCSNe; hence, we estimate the event counts for Hyper(Super)-Kamiokande, DUNE, and JUNO. We find that the event count with the energy of $\gtrsim 80$ MeV is a few orders of magnitude higher than that of the thermal neutrinos regardless of the detectors, and muon production may also happen in these detectors by $ν_μ$ with the energy of $\gtrsim 100$ MeV. The neutrino signals provide a precious information on deciphering the inner dynamics of CCSN and placing a constraint on the physics of neutrino oscillation; indeed, the detection of the high energy neutrinos through charged current reaction channels will be a smoking gun evidence of neutrino flavor conversion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Being Careful with the Field Formation Interpretation of GW190412

The LIGO/Virgo Scientific Collaboration recently announced the detection of a compact object binary merger, GW190412, as the first asymmetric binary black hole (BBH) merger with mass ratio $q\approx0.25$. Other than the mass ratio, this BBH has shown to have a positive effective spin of around $χ_{\rm eff}\approx0.28$. Assuming a field formation channel, associating this effective spin to either the primary or the secondary BH each has its implications: If the spin of the BBH comes form the primary BH, it has consequences for the efficiency of angular momentum transport in the formation of the BH. If, on the other hand, the spin is due to the secondary BH through tidal spin-up processes, one has to note that (i) such processes have very short delay-times, and (ii) subsequently, their local merger rate is determined by local star formation rate at assumed formation metallicity of the BBH. We show that the predicted merger rate density from this channel is $\lesssim 0.3~\rm Gpc^{-3} yr^{-1}$ and in tension with the rather high local merger rate of such systems which we estimate from this single event to be $\sim 1.7^{+2.5}_{-1.4}~\rm Gpc^{-3} yr^{-1}$ (90\% confidence interval, and assuming 50 days of observing time). Large natal kicks ($v\gtrsim 500\,{\rm km/s}$) would be required to get such BBHs with an in-plane spin component to account for the marginal detection of precession in GW190412. However, this would only exacerbate the tension as the estimated local merger rate would be further decreased. Similarly, the formation of such systems through the dynamical assembly is exceedingly rare, leaving this system a dilemma hard to account for with the currently accepted paradigms of BBH formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining properties of neutron star merger outflows with radio observations

The jet opening angle and inclination of GW170817 -- the first detected binary neutron star merger -- were vital to understand its energetics, relation to short gamma-ray bursts, and refinement of the standard siren-based determination of the Hubble constant, $H_0$. These basic quantities were determined through a combination of the radio lightcurve and Very Long Baseline Interferometry (VLBI) measurements of proper motion. In this paper we discuss and quantify the prospects for the use of radio VLBI observations and observations of scintillation-induced variability to measure the source size and proper motion of merger afterglows, and thereby infer properties of the merger including inclination angle, opening angle and energetics. We show that these techniques are complementary as they probe different parts of the circum-merger density/inclination angle parameter space and different periods of the temporal evolution of the afterglow. We also find that while VLBI observations will be limited to the very closest events it will be possible to detect scintillation for a large fraction of events beyond the range of current gravitational wave detectors. Scintillation will also be detectable with next generation telescopes such as the Square Kilometre Array, 2000 antenna Deep Synoptic Array and the next generation Very Large Array, for a large fraction of events detected with third generation gravitational wave detectors. Finally, we discuss prospects for the measurement of the $H_0$ with VLBI observations of neutron star mergers and compare this technique to other standard siren methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distinguishing the nature of comparable-mass neutron star binary systems with multimessenger observations: GW170817 case study

The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implications of the search for optical counterparts during the second part of the Advanced LIGO's and Advanced Virgo's third observing run: lessons learned for future follow-up observations

Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a "MassGap" component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent mixing of r-process elements in the Milky Way

We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three independent observations to constrain the model parameters: (i) the scatter of radioactively stable $r$-process element abundances, (ii) the largest $r$-process enrichment values observed in any solar neighborhood stars and (iii) the isotope abundance ratios of different radioactive $r$-process elements ($^{244}$Pu/$^{238}$U and $^{247}$Cm/$^{238}$U) at the early solar system as compared to their formation. Our results indicate that the Galactic $r$-process rate and the diffusion coefficient are respectively $r<4\times 10^{-5}\mbox{ yr}^{-1}, D>0.1 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ ($r<4\times 10^{-6}\mbox{ yr}^{-1}, D>0.5 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ for collapsars or similarly prolific $r$-process sources) with allowed values satisfying an approximate anti-correlation such that $D\approx r^{-2/3}$, implying that the time between two $r$-process events that enrich the same location in the Galaxy, is $τ_{\rm mix}\approx 100-200\mbox{ Myr}$. This suggests that a fraction of $\sim 0.8$ ($\sim 0.5$) of the observed $^{247}$Cm ($^{244}$Pu) abundance is dominated by one $r$-process event in the early solar system. Radioactively stable element abundances are dominated by contributions from $\sim 10$ different events in the early solar system. For metal poor stars (with [Fe/H]$\lesssim -2$), their $r$-process abundances are dominated by either a single or several events, depending on the star formation history.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Velocity correction for Hubble constant measurements from standard sirens

Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. A first implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 68.3^{+ 4.6}_{-4.5}$ km/s/Mpc. While this revision is minor, it demonstrates that our method makes it possible for obtaining an unbiased and accurate measurements of $H_0$ at the precision required for the standard siren cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An ASKAP search for a radio counterpart to the first high-significance neutron star-black hole merger LIGO/Virgo S190814bv

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $ΔT$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

GROWTH on S190814bv: Deep Synoptic Limits on the Optical/Near-Infrared Counterpart to a Neutron Star-Black Hole Merger

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $κ< 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ χ< 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Implications of the search for optical counterparts during the first six months of the Advanced LIGO's and Advanced Virgo's third observing run: possible limits on the ejecta mass and binary properties

GW170817 showed that neutron star mergers not only emit gravitational waves but also can release electromagnetic signatures in multiple wavelengths. Within the first half of the third observing run of the Advanced LIGO and Virgo detectors, there have been a number of gravitational wave candidates of compact binary systems for which at least one component is potentially a neutron star. In this article, we look at the candidates S190425z, S190426c, S190510g, S190901ap, and S190910h, predicted to have potentially a non-zero remnant mass, in more detail. All these triggers have been followed up with extensive campaigns by the astronomical community doing electromagnetic searches for their optical counterparts; however, according to the released classification, there is a high probability that some of these events might not be of extraterrestrial origin. Assuming that the triggers are caused by a compact binary coalescence and that the individual source locations have been covered during the EM follow-up campaigns, we employ three different kilonova models and apply them to derive possible constraints on the matter ejection consistent with the publicly available gravitational-wave trigger information and the lack of a kilonova detection. These upper bounds on the ejecta mass can be related to limits on the maximum mass of the binary neutron star candidate S190425z and to constraints on the mass-ratio, spin, and NS compactness for the potential black hole-neutron star candidate S190426c. Our results show that deeper electromagnetic observations for future gravitational wave events near the horizon limit of the advanced detectors are essential.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Radioactive heating rate of r-process elements and macronova light curve

We study the heating rate of r-process nuclei and thermalization of decay products in neutron star merger ejecta and macronova (kilonova) light curves. Thermalization of charged decay products, i.e., electrons, $α$-particles, and fission fragments is calculated according to their injection energy. The $γ$-ray thermalization processes are also properly calculated by taking the $γ$-ray spectrum of each decay into account. We show that the $β$-decay heating rate at later times approaches a power-law decline as $\propto t^{-2.8}$, which agrees with the result of Waxman et al. (2019). We present a new analytic model to calculate macronova light curves, in which the density structure of the ejecta is accounted for. We demonstrate that the observed bolometric light curve and temperature evolution of the macronova associated with GW170817 are reproduced well by the $β$-decay heating rate with the solar r-process abundance pattern. We interpret the break in the observed bolometric light curve around a week as a result of the diffusion wave crossing a significant part of the ejecta rather than a thermalization break. We also show that the time-weighted integral of the bolometric light curve (Katz integral) is useful to provide an estimate of the total r-process mass from the observed data, which is independent of the highly uncertain radiative transfer. For the macronova in GW170817, the ejecta mass is robustly estimated as $\approx 0.05M_{\odot}$ for $A_{\rm min}\leq 72$ and $85\leq A_{\rm min}\leq 130$ with the solar r-process abundance pattern. The code for computation of the heating rate and light curve for given initial nuclear abundances is publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Are Ultra-Long Gamma-Ray Bursts Caused by Blue Supergiant Collapsars, Newborn Magnetars, or White Dwarf Tidal Disruption Events?

Ultra-long gamma-ray bursts (ulGRBs) are a new population of GRBs with extreme durations of $\sim 10^{4}$ s. Leading candidates for their origin are blue supergiant collapsars, magnetars, and white dwarf tidal disruption events (WD-TDEs) caused by massive black holes (BHs). Recent observations of supernova-like (SN-like) bumps associated with ulGRBs challenged both the WD-TDE and the blue supergiant models because of the detection of SNe and the absence of hydrogen lines, respectively. We propose that WD-TDEs can accommodate the observed SN-like bumps if the fallback WD matter releases energy into the unbound WD ejecta. The observed ejecta energy, luminosity, and velocity are explained by the gravitational energy, Eddington luminosity, and escape velocity of the formed accretion disk, respectively. We also show that the observed X-rays can ionize the ejecta, eliminating lines. The SN-like light curves (SN 2011kl) for the ulGRB 111209A are consistent with all three models, although a magnetar model is unnatural because the spin-down time required to power the SN-like bump is a hundred times longer than the GRB. Our results imply that TDEs are a possible energy source for SN-like events in general and for ulGRBs in particular.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effects of neutron-star dynamic tides on gravitational waveforms within the effective-one-body approach

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measurability of the tidal deformability by gravitational waves from coalescing binary neutron stars

Combining new gravitational waveforms derived by long-term (14--16 orbits) numerical-relativity simulations with waveforms by an effective-one-body (EOB) formalism for coalescing binary neutron stars, we construct hybrid waveforms and estimate the measurability for the dimensionless tidal deformability of the neutron stars, $Λ$, by advanced gravitational-wave detectors. We focus on the equal-mass case with the total mass $2.7M_\odot$. We find that for an event at a hypothetical effective distance of $D_{\rm eff}=200$ Mpc, the distinguishable difference in the dimensionless tidal deformability will be $\approx 100$, 400, and 800 at 1-$σ$, 2-$σ$, and 3-$σ$ levels, respectively, for advanced LIGO. If the true equation of state is stiff and the typical neutron-star radius is $R \gtrsim 13 $ km, our analysis suggests that the radius will be constrained within $\approx 1$ km at 2-$σ$ level for an event at $D_{\rm eff}=200$ Mpc. On the other hand, if the true equation of state is soft and the typical neutron-star radius is $R\lesssim 12$ km , it will be difficult to narrow down the equation of state among many soft ones, although it is still possible to discriminate the true one from stiff equations of state with $R\gtrsim 13$ km. We also find that gravitational waves from binary neutron stars will be distinguished from those from spinless binary black holes at more than 2-$σ$ level for an event at $D_{\rm eff}=200$ Mpc. The validity of the EOB formalism, Taylor-T4, and Taylor-F2 approximants as the inspiral waveform model is also examined.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Natal Kicks and Time Delays in Merging Neutron Star Binaries - Implications for r-process nucleosynthesis in Ultra Faint Dwarfs and in the Milky Way

Merging neutron star binaries are prime candidate sources for heavy r-process nucleosynthesis. The amount of heavy r-process material is consistent with the mass ejection and rates of mergers, and abundances of relic radioactive materials suggest that heavy r-process material is produced in rare events. Observations of possible macronovae provide further support for this model. Still, some concerns remain. One is the observation of heavy r-process elements in Ultra Faint Dwarf (UFD) galaxies. The escape velocities from UFDs are so small that the natal kicks, taking place at neutron stars birth, might eject such binaries from UFDs. Furthermore the old stellar populations of UFDs requires that r-process nucleosynthesis must have taken place very early on, while it may take several Gyr for compact binaries to merge. This last problem arises also within the Milky Way where heavy r-process materials has been observed in some low metallicity stars. We show here that since a significant fraction of neutron star binaries form with a very small proper motion those won't be ejected even from a UFD. Furthermore, a significant fraction of these binaries merge within a few hundred Myr and in some cases the merger time is even shorter than a hundred Myr. This population of "rapid mergers" explains the appearance of heavy r-process material in both UFDs and in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

r-process Production Sites as inferred from Eu Abundances in Dwarf Galaxies

Recent observations of $r$-process material in ultra-faint dwarf galaxies (UFDs) shed light on the sources of these elements. Strong upper limits on the Eu mass in some UFDs combined with detections of much larger masses in a UFD, Reticulum II, and other dwarf galaxies imply that Eu production is dominated by rare events, and that the minimal Eu mass observed in any UFD is approximately the amount of Eu mass produced per event. This is consistent with other independent observations in the Galaxy. We estimate, using a model independent likelihood analysis, the rate and Eu (Fe) mass produced per $r$-process (Fe production) event in dwarf galaxies including classical dwarfs and UFDs. The mass and rate of the Fe production events are consistent with the normal core-collapse supernova~(ccSN) scenario. The Eu mass per event is $3\times 10^{-5}M_{\odot}<\tilde{m}_{\rm Eu}<2\times 10^{-4}M_{\odot}$, corresponding to a total $r$-process mass per event of $6\times 10^{-3}M_{\odot}<\tilde{m}_{r-process}<4\times 10^{-2}M_{\odot}$. The rate of $r$-process events is $2.5\times 10^{-4}<R_{rp/SN}<1.4\times 10^{-3}$ as compared with the ccSNe rate. These values are consistent with the total Eu mass observed in our own Galaxy, suggesting that the same mechanism is behind the production of $r$-process events in both dwarf galaxies and the Milky Way, and that it may be the dominant mechanism for production of $r$-process elements in the Universe. The results are consistent with neutron star mergers estimates, but cannot rule out other rare core collapse scenarios, provided that they produce a significant amounts of $r$-process material per event.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radio Counterparts of Compact Binary Mergers detectable in Gravitational Waves: A Simulation for an Optimized Survey

Mergers of binary neutron stars and black hole-neutron star binaries produce gravitational-wave (GW) emission and outflows with significant kinetic energies. These outflows result in radio emissions through synchrotron radiation. We explore the detectability of these synchrotron generated radio signals by follow-up observations of GW merger events lacking a detection of electromagnetic counterparts in other wavelengths. We model radio light curves arising from (i) sub-relativistic merger ejecta and (ii) ultra-relativistic jets. The former produces radio remnants on timescales of a few years and the latter produces $γ$-ray bursts in the direction of the jet and orphan-radio afterglows extending over wider angles on timescales of weeks. Based on the derived light curves, we suggest an optimized survey at $1.4$ GHz with five epochs separated by a logarithmic time interval. We estimate the detectability of the radio counterparts of simulated GW-merger events to be detected by advanced LIGO and Virgo by current and future radio facilities. The detectable distances for these GW merger events could be as high as 1 Gpc. $20$--$60\%$ of the long-lasting radio remnants will be detectable in the case of the moderate kinetic energy of $3\cdot 10^{50}$ erg and a circum-merger density of $0.1 {\rm cm^{-3}}$ or larger, while $5$--$20\%$ of the orphan radio afterglows with kinetic energy of $10^{48}$ erg will be detectable. The detection likelihood increases if one focuses on the well-localizable GW events. We discuss the background noise due to radio fluxes of host galaxies and false positives arising from extragalactic radio transients and variable Active Galactic Nuclei and we show that the quiet radio transient sky is of great advantage when searching for the radio counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing the Magnetar Model via Late Time Radio Observations of Two Macronova Candidates

Compact binary mergers may have already been observed as they are the leading model for short gamma-ray bursts (sGRBs). Radioactive decay within the ejecta from these mergers is expected to produce an infra-red flare, dubbed macronova (or kilonova), on a time scale of a week. Recently two such macronova candidates were identified in followup observations of sGRBs, strengthening the possibility that those indeed arise from mergers. The same ejecta will also produce a long term (months to years) radio emission due to its interaction with the surrounding ISM. In search for this emission, we observed the two macronova candidates, GRB 130603B and GRB 060614 with the Jansky very large array (VLA) and the Australia Telescope Compact Array (ATCA). Our observations resulted in null-detections, putting strong upper limits on the kinetic energy and mass of the ejecta. A possible outcome of a merger is a highly magnetized neutron star (a magnetar), which has been suggested as the central engine for GRBs. Such a magnetar will deposit a significant fraction of its energy into the ejecta leading to a brighter radio flare. Our results, therefore, rule out magnetars in these two events.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Macronova in GRB 050709 and the GRB/macronova connection

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A possible Macronova in the late afterglow of the `long-short' burst GRB 060614

Long-duration ($>2$ s) $γ$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczyński macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploring tidal effects of coalescing binary neutron stars in numerical relativity II: Longterm simulations

We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mass ejection from neutron star mergers: different components and expected radio signals

In addition to producing a strong gravitational signal, a short gamma-ray burst (GRB), and a compact remnant, neutron star mergers eject significant masses at significant kinetic energies. This mass ejection takes place via dynamical mass ejection and a GRB jet but other processes have also been suggested: a shock-breakout material, a cocoon resulting from the interaction of the jet with other ejecta, and viscous and neutrino driven winds from the central remnant or the accretion disk. The different components of the ejected masses include up to a few percent of a solar mass, some of which is ejected at relativistic velocities. The interaction of these ejecta with the surrounding interstellar medium will produce a long lasting radio flare, in a similar way to GRB afterglows or to radio supernovae. The relative strength of the different signals depends strongly on the viewing angle. An observer along the jet axis or close to it will detect a strong signal at a few dozen days from the radio afterglow (or the orphan radio afterglow) produced by the highly relativistic GRB jet. For a generic observer at larger viewing angles, the dynamical ejecta, whose contribution peaks a year or so after the event, will generally dominate. Depending on the observed frequency and the external density, other components may also give rise to a significant contribution. We also compare these estimates with the radio signature of the short GRB 130603B. The radio flare from the dynamical ejecta might be detectable with the EVLA and the LOFAR for the higher range of external densities $n\gtrsim 0.5{\rm cm^{-3}}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radioactive decay products in neutron star merger ejecta: heating efficiency and $γ$-ray emission

The radioactive decay of the freshly synthesized $r$-process nuclei ejected in compact binary mergers power optical/infrared macronovae (kilonovae) that follow these events. The light curves depend critically on the energy partition among the different products of the radioactive decay and this plays an important role in estimates of the amount of ejected $r$-process elements from a given observed signal. We study the energy partition and $γ$-ray emission of the radioactive decay. We show that $20$-$50\%$ of the total radioactive energy is released in $γ$-rays on timescales from hours to a month. The number of emitted $γ$-rays per unit energy interval has roughly a flat spectrum between a few dozen keV and $1$ MeV so that most of this energy is carried by $\sim 1$ MeV $γ$-rays. However at the peak of macronova emission the optical depth of the $γ$-rays is $\sim 0.02$ and most of the $γ$-rays escape. The loss of these $γ$-rays reduces the heat deposition into the ejecta and hence reduces the expected macronova signals if those are lanthanides dominated. This implies that the ejected mass is larger by a factor of $2$-$3$ than what was previously estimated. Spontaneous fission heats up the ejecta and the heating rate can increase if a sufficient amount of transuranic nuclei are synthesized. Direct measurements of these escaping $γ$-rays may provide the ultimate proof for the macronova mechanisms and an identification of the $r$-process nucleosynthesis sites. However, the chances to detect these signals are slim with current X-ray and $γ$-ray missions. New detectors, more sensitive by at least a factor of ten, are needed for a realistic detection rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Jet collimation in the ejecta of double neutron star merger: new canonical picture of short gamma-ray bursts

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10°. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Possible Indirect Confirmation of the Existence of Pop III Massive Stars by Gravitational Wave

We perform population synthesis simulations for Population III (Pop III) coalescing binary neutron stars (NS-NSs), neutron star - black hole binaries (NS-BHs), and binary black holes (BH-BHs) which merge within the age of the universe. We found that the typical mass of Pop III BH-BHs is $\sim 30\ \rm{M}_{\odot}$ so that the inspiral chirp signal of gravitational waves can be detected up to z=0.28 by KAGRA, Adv. LIGO, Adv. Virgo and GEO network. Our simulations suggest that the detection rate of the coalescing Pop III BH-BHs is $140 (68) \cdot ({\rm SFR}_{\rm p}/10^{-2.5} \rm{M}_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}) \cdot {\rm Err}_{\rm sys} ~{\rm events} ~{\rm yr}^{-1}$ for the flat (Salpeter) initial mass function (IMF), respectively, where $\rm SFR_p$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate and the possible systematic errors due to the assumptions in Pop III population synthesis, respectively. $\rm Err_{sys}=1$ correspond to conventional parameters for Pop I stars. From the observation of the chirp signal of the coalescing Pop III BH-BHs, we can determine both the mass and the redshift of the binary for the cosmological parameters determined by Planck satellite. Our simulations suggest that the cumulative redshift distribution of the coalescing Pop III BH-BHs depends almost only on the cosmological parameters. We might be able to confirm the existence of Pop III massive stars of mass $\sim 30~\rm M_{\odot}$ by the detections of gravitational waves if the merger rate of the Pop III massive BH-BHs dominates that of Pop I BH-BHs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploring tidal effects of coalescing binary neutron stars in numerical relativity

We study gravitational waves emitted in the late inspiral stage of binary neutron stars by analyzing the waveform obtained in numerical-relativity simulations. For deriving the physical gravitational waveforms from the numerical results, the resolution extrapolation plays an essential role for our simulations. The extrapolated gravitational-wave phases are compared with those calculated in the post-Newtonian (PN) and effective-one-body (EOB) formalisms including corrections of tidal effects. We show that the extrapolated gravitational-wave phases in numerical relativity agree well with those by the PN and EOB calculations for most of the inspiral stage except for a tidally-dominated, final inspiral stage, in which the PN and EOB results underestimate the tidal effects. Nevertheless, the accumulated phase difference between our extrapolated results and the results by the PN/EOB calculations is at most 1--3 radian in the last 15 cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Progenitor models of the electromagnetic transient associated with the short GRB 130603B

An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiative Transfer Simulations for Neutron Star Merger Ejecta

The merger of binary neutron stars (NSs) is among the most promising gravitational wave (GW) sources. Next-generation GW detectors are expected to detect signals from the NS merger within 200 Mpc. Detection of electromagnetic wave (EM) counterpart is crucial to understand the nature of GW sources. Among possible EM emission from the NS merger, emission powered by radioactive r-process nuclei is one of the best targets for follow-up observations. However, prediction so far does not take into account detailed r-process element abundances in the ejecta. We perform radiative transfer simulations for the NS merger ejecta including all the r-process elements from Ga to U for the first time. We show that the opacity in the NS merger ejecta is about kappa = 10 cm^2 g^{-1}, which is higher than that of Fe-rich Type Ia supernova ejecta by a factor of ~ 100. As a result, the emission is fainter and longer than previously expected. The spectra are almost featureless due to the high expansion velocity and bound-bound transitions of many different r-process elements. We demonstrate that the emission is brighter for a higher mass ratio of two NSs and a softer equation of states adopted in the merger simulations. Because of the red color of the emission, follow-up observations in red optical and near-infrared (NIR) wavelengths will be the most efficient. At 200 Mpc, expected brightness of the emission is i = 22 - 25 AB mag, z = 21 - 23 AB mag, and 21 - 24 AB mag in NIR JHK bands. Thus, observations with wide-field 4m- and 8m-class optical telescopes and wide-field NIR space telescopes are necessary. We also argue that the emission powered by radioactive energy can be detected in the afterglow of nearby short gamma-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radioactively Powered Emission from Black Hole-Neutron Star Mergers

Detection of electromagnetic counterparts of gravitational wave (GW) sources is important to unveil the nature of compact binary coalescences. We perform three-dimensional, time-dependent, multi-frequency radiative transfer simulations for radioactively powered emission from the ejecta of black hole (BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of the BH, and equations of state of dense matter, BH-NS mergers can eject more material than NS-NS mergers. In such cases, radioactively powered emission from the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We show that, in spite of the expected larger distances to BH-NS merger events, observed brightness of BH-NS mergers can be comparable to or even higher than that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger ejecta are confined to a small solid angle, the emission from BH-NS merger ejecta tends to be bluer than that from NS-NS merger ejecta for a given total luminosity. Thanks to this property, we might be able to distinguish BH-NS merger events from NS-NS merger events by multi-band observations of the radioactively powered emission. In addition to the GW observations, such electromagnetic observations can potentially provide independent information on the nature of compact binary coalescences.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Remnant massive neutron stars of binary neutron star mergers: Evolution process and gravitational waveform

Massive (hypermassive and supramassive) neutron stars are likely to be often formed after the merger of binary neutron stars. We explore the evolution process of the remnant massive neutron stars and gravitational waves emitted by them, based on numerical-relativity simulations for binary neutron star mergers employing a variety of equations of state and choosing a plausible range of the neutron-star mass of binaries. We show that the lifetime of remnant hypermassive neutron stars depends strongly on the total binary mass and also on the equations of state. Gravitational waves emitted by the remnant massive neutron stars universally have a quasiperiodic nature of an approximately constant frequency although the frequency varies with time. We also show that the frequency and time-variation feature of gravitational waves depend strongly on the equations of state. We derive a fitting formula for the quasiperiodic gravitational waveforms, which may be used for the data analysis of a gravitational-wave signal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Binary Neutron Star Mergers: Dependence on the Nuclear Equation of State

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one σlevel. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.

0 citations0 reviewsNo code linkedOpen access