Researcher profile

Koutarou Kyutoku

Koutarou Kyutoku contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2022arXiv

Comprehensive study of mass ejection and nucleosynthesis in binary neutron star mergers leaving short-lived massive neutron stars

By performing general relativistic hydrodynamics simulations with an approximate neutrino-radiation transfer, the properties of ejecta in dynamical and post-merger phases are investigated for the cases in which the remnant massive neutron star collapses into a black hole in $\lesssim 20$ ms after the onset of the merger. The dynamical mass ejection is investigated in three-dimensional simulations. The post-merger mass ejection is investigated in two-dimensional axisymmetric simulations with viscosity using the three-dimensional post-merger systems as the initial conditions. We show that the typical neutron-richness of the dynamical ejecta is higher for the merger of more asymmetric binaries; hence, heavier $r$-process nuclei are dominantly synthesized. The post-merger ejecta are shown to have only a mild neutron-richness, which results in the production of lighter $r$-process nuclei, irrespective of binary mass ratios. Because of the larger disk mass, the post-merger ejecta mass is larger for more asymmetric binary mergers. Thus, the post-merger ejecta can compensate for the underproduced lighter $r$-process nuclei for asymmetric merger cases. As a result, by summing up both ejecta components, the solar residual $r$-process pattern is reproduced within the average deviation of a factor of three, irrespective of the binary mass ratio. Our result also indicates that the (about a factor of a few) light-to-heavy abundance scatter observed in $r$-process-enhanced stars can be attributed to variation in the binary mass ratio and total mass. Implications of our results associated with the mass distribution of compact neutron star binaries and the magnetar scenario of short gamma-ray bursts are discussed.

preprint2022arXiv

How many extra-Galactic stellar-mass binary black holes will be detected by space gravitational-wave interferometers?

On the basis of GWTC-3, we discuss the detection prospect of extra-Galactic binary black holes (BBHs) by space gravitational-wave interferometers. In particular, targeting BBHs with component masses around 5-100$M_\odot$, we directly incorporate the chirp mass distribution of the 62 BBHs detected at high significance. We find that, due to the reduction of both the comoving merger rate and a weighted average of chirp masses, the expected detection numbers are generally much smaller than the results obtained by the same authors immediately after the report of GW150914. For LISA, the total BBH detections are estimated to be $N_{\rm tot}\sim 2 (T/4{\rm yr})^{3/2}(ρ_{\rm thr}/10)^{-3}$, dominated by nearly monochromatic BBHs ($ρ_{\rm thr}$: the detection threshold, $T$: the observational period). TianQin will have a total detection number $N_{\rm tot}$ similar to LISA. Meanwhile, TianQin has potential to find $N_{\rm mer}\sim0.6 (T/4{\rm yr})^{7/4}(ρ_{\rm thr}/10)^{-3}$ BBHs that merge in the observational period. This number for merging BBHs is 4-5 times larger than that of LISA, because of the difference between the optimal bands. We also investigate prospects for joint operations of multiple detectors, finding that concurrent observations will be more advantageous than sequential ones.

preprint2021arXiv

Reducing orbital eccentricity in initial data of black hole--neutron star binaries in the puncture framework

We develop a method to compute low-eccentricity initial data of black hole--neutron star binaries in the puncture framework extending previous work on other types of compact binaries. In addition to adjusting the orbital angular velocity of the binary, the approaching velocity of a neutron star is incorporated by modifying the helical Killing vector used to derive equations of the hydrostationary equilibrium. The approaching velocity of the black hole is then induced by requiring the vanishing of the total linear momentum of the system, differently from the case of binary black holes in the puncture framework where the linear momentum of each black hole is specified explicitly. We successfully reduce the orbital eccentricity to <~0.001 by modifying the parameters iteratively using simulations of ~3 orbits both for nonprecessing and precessing configurations. We find that empirical formulas for binary black holes derived in the excision framework do not reduce the orbital eccentricity to ~0.001 for black hole--neutron star binaries in the puncture framework, although they work for binary neutron stars.

preprint2020arXiv

Mass ejection from disks surrounding a low-mass black hole: Viscous neutrino-radiation hydrodynamics simulation in full general relativity

New viscous neutrino-radiation hydrodynamics simulations are performed for accretion disks surrounding a spinning black hole with low mass $3M_\odot$ and dimensionless spin 0.8 or 0.6 in full general relativity, aiming at modeling the evolution of a merger remnant of massive binary neutron stars or low-mass black hole-neutron star binaries. We reconfirm the following results found by previous studies of other groups: 15-30% of the disk mass is ejected from the system with the average velocity of $\sim $5-10% of the speed of light for the plausible profile of the disk as merger remnants. In addition, we find that for the not extremely high viscous coefficient case, the neutron richness of the ejecta does not become very high, because weak interaction processes enhance the electron fraction during the viscous expansion of the disk before the onset of the mass ejection, resulting in the suppression of the lanthanide synthesis. For high-mass disks, the viscous expansion timescale is increased by a longer-term neutrino emission, and hence, the electron fraction of the ejecta becomes even higher. We also confirm that the mass distribution of the electron fraction depends strongly on the magnitude of the given viscous coefficient. This demonstrates that a first-principle magnetohydrodynamics simulation is necessary for black hole-disk systems with sufficient grid resolution and with sufficiently long timescale (longer than seconds) to clarify the nucleosynthesis and electromagnetic signals from them.

preprint2020arXiv

On the possibility of GW190425 being a black hole--neutron star binary merger

We argue that the kilonova/macronova associated with the gravitational-wave event GW190425 could have been bright enough to be detected if it was caused by the merger of a low-mass black hole and a neutron star. Although tidal disruption occurs for such a low-mass black hole as is generally expected, the masses of the dynamical ejecta are limited to <~10^{-3}M_sun, which is consistent with previous work in the literature. The remnant disk could be as massive as 0.05--0.1M_sun, and the disk outflow of ~0.01--0.03M_sun is likely to be driven by viscous or magnetohydrodynamic effects. The disk outflow may not be neutron-rich enough to synthesize an abundance of lanthanide elements, even in the absence of strong neutrino emitter, if the ejection is driven on the viscous time scale of >~0.3s. If this is the case, the opacity of the disk outflow is kept moderate, and a kilonova/macronova at the distance of GW190425 reaches a detectable brightness of 20--21mag at 1day after merger for most viewing angles. If some disk activity ejects the mass within ~0.1s, instead, lanthanide-rich outflows would be launched and the detection of emission becomes challenging. Future possible detections of kilonovae/macronovae from GW190425-like systems will disfavor the prompt collapse of binary neutron stars and a non-disruptive low-mass black hole--neutron star binary associated with a small neutron-star radius, whose mass ejection is negligible. The host-galaxy distance will constrain the viewing angle and deliver further information about the mass ejection.

preprint2020arXiv

Properties of the remnant disk and the dynamical ejecta produced in low-mass black hole-neutron star mergers

We systematically perform numerical-relativity simulations for low-mass black hole-neutron star mergers for the models with seven mass ratios $Q=M_{\rm BH}/M_{\rm NS}$ ranging from 1.5 to 4.4, and three neutron-star equations of state, focusing on properties of matter remaining outside the black hole and ejected dynamically during the merger. We pay particular attention to the dependence on the mass ratio of the binaries. It is found that the rest mass remaining outside the apparent horizon after the merger depends only weakly on the mass ratio for the models with low mass ratios. It is also clarified that the rest mass of the ejecta has a peak at $Q \sim 3$, and decreases steeply as the mass ratio decreases for the low mass-ratio case. We present a novel analysis method for the behavior of matter during the merger, focusing on the matter distribution in the phase space of specific energy and specific angular momentum. Then we model the matter distribution during and after the merger. Using the result of the analysis, we discuss the properties of the ejecta.

preprint2020arXiv

Sub-radian-accuracy gravitational waves from coalescing binary neutron stars II: Systematic study on the equation of state, binary mass, and mass ratio

We report results of numerical relativity simulations for {\it new} 26 non-spinning binary neutron star systems with 6 grid resolutions using an adaptive mesh refinement numerical re\ lativity code {\tt SACRA-MPI}. The finest grid spacing is $\approx 64$--$85$ m, depending on the systems. First, we derive long-term high-precision inspiral gravitational waveforms and show that the accumulated gravitational-wave phase error due to the finite grid resolution is less than $0.5$ rad during more than $200$ rad phase evolution irrespective of the systems. We also find that the gravitational-wave phase error for a binary system with a tabulated equation of state (EOS) is comparable to that for a piecewise polytropic EOS. Then we validate the SACRA inspiral gravitational waveform template, which will be used to extract tidal deformability from gravitational wave observation, and find that accuracy of \ our waveform modeling is $\lesssim 0.1$ rad in the gravitational-wave phase and $\lesssim 20 \%$ in the gravitational-wave amplitude up to the gravitational-wave frequency $1000$ Hz.\ Finally, we calibrate the proposed universal relations between a post-merger gravitational wave signal and tidal deformability/neutron star radius in the literature and show that th\ ey suffer from systematics and many relations proposed as universal are not very universal. Improved fitting formulae are also proposed.

preprint2020arXiv

Viscous evolution of a massive disk surrounding stellar-mass black holes in full general relativity

Long-term viscous neutrino-radiation hydrodynamics simulations in full general relativity are performed for a massive disk surrounding spinning stellar-mass black holes with mass $M_{\rm BH}=4$, $6$, and $10M_\odot$ and initial dimensionless spin $χ\approx 0.8$. The initial disk is chosen to have mass $M_{\rm disk}\approx 0.1$ or $3M_\odot$ as plausible models of the remnants for the merger of black hole-neutron star binaries or the stellar core collapse from a rapidly rotating progenitor, respectively. For $M_{\rm disk} \approx 0.1M_\odot$ with the outer disk edge initially located at $r_{\rm out} \sim 200$ km, we find that $15$%-$20$% of $M_{\rm disk}$ is ejected and the average electron fraction of the ejecta is $\langle Y_e \rangle = 0.30$-$0.35$ as found in the previous study. For $M_{\rm disk} \approx 3M_\odot$, we find that $\approx 10$%-$20$% of $M_{\rm disk}$ is ejected for $r_{\rm out}\approx 200$-$1000$ km. In addition, $\langle Y_e \rangle$ of the ejecta can be enhanced to be $\gtrsim 0.4$ because the electron fraction is increased significantly during the long-term viscous expansion of the disk with high neutrino luminosity until the mass ejection sets in. Our results suggest that not heavy $r$-process elements but light trans-iron elements would be synthesized in the matter ejected from a massive torus surrounding stellar-mass black holes. We also find that the outcomes of the viscous evolution for the high-mass disk case is composed of a rapidly spinning black hole surrounded by a torus with a narrow funnel, which appears to be suitable for generating gamma-ray bursts.

preprint2019arXiv

r-Process enrichment in the Galactic halo characterized by nucleosynthesis variation in the ejecta of coalescing neutron star binaries

A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is tightly correlated with [Eu/Fe], giving the lowest [Y/Eu] ratio that levels off at r-II stars. On the other hand, recent hydrodynamical simulations of coalescing double neutron stars (cNSNSs) have suggested that r-process sites may be separated into two classes providing different electron-fraction distributions: tidally-driven dynamical ejecta and (dynamical or postmerger) non-tidal ejecta. Here, we show that a widely spanning feature of [Eu/Fe] can be reproduced by models that consider the different masses of tidally-driven dynamical ejecta from both cNSNSs and coalescing black hole/neutron star binaries (cBHNSs). In addition, the observed [Y/Eu] trend is explained by the combined nucleosynthesis in two kinds of ejecta with varying mass asymmetry in double NS systems. Our scenario suggests that massive tidally-driven dynamical ejecta accompanied by massive non-tidal part from cNSNSs or cBHNSs could alone accommodate r-II abundances, including an actinide boost in some cases. The event rate for cNSNSs estimated from our study agrees with the latest result of ~1000 (90% confidence interval of 110-3840) Gpc$^{-3}$yr$^{-1}$ by gravitational-wave detection, and a few events per Gpc$^3$ per year of cBHNSs associated with r-process production are predicted to emerge.