BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Masaru Shibata

Masaru Shibata contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcastro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
3topics
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

17 published item(s)

preprint2022arXiv

A low-mass binary neutron star: long-term ejecta evolution and kilonovae with weak blue emission

We study the long-term evolution of ejecta formed in a binary neutron star (BNS) merger that results in a long-lived remnant NS by performing a hydrodynamics simulation with the outflow data of a numerical relativity simulation as the initial condition. At the homologously expanding phase, the total ejecta mass reaches $\approx0.1\,M_\odot$ with an average velocity of $\approx0.1\,c$ and lanthanide fraction of $\approx 0.005$. We further perform the radiative transfer simulation employing the obtained ejecta profile. We find that, contrary to a naive expectation from the large ejecta mass and low lanthanide fraction, the optical emission is not as bright as that in GW170817/AT2017gfo, while the infrared emission can be brighter. This light curve property is attributed to preferential diffusion of photons toward the equatorial direction due to the prolate ejecta morphology, large opacity contribution of Zr, Y, and lanthanides, and low specific heating rate of the ejecta. Our results suggest that these light curve features could be used as an indicator for the presence of a long-lived remnant NS. We also found that the bright optical emission broadly consistent with GW170817/AT2017gfo is realized for the case that the high-velocity ejecta components in the polar region are suppressed. These results suggest that the remnant in GW170817/AT2017gfo is unlikely to be a long-lived NS, but might have collapsed to a black hole within ${\cal O}(0.1)$ s.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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

Implementation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its application to a binary neutron star merger

We implement advanced Riemann solvers HLLC and HLLD \cite{Mignone:2005ft,MUB:2009} together with an advanced constrained transport scheme \cite{Gardiner:2007nc} in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We validate our implementation by performing a series of one- and multi-dimensional test problems for relativistic hydrodynamics and magnetohydrodynamics in both Minkowski spacetime and a static black hole spacetime. We find that the numerical solutions with the advanced Riemann solvers are more accurate than those with the HLLE solver \cite{DelZanna:2002rv}, which was originally implemented in our code. As an application to numerical relativity, we simulate an asymmetric binary neutron star merger leading to a short-lived massive neutron star both with and without magnetic fields. We find that the lifetime of the rotating massive neutron star formed after the merger and also the amount of the tidally-driven dynamical ejecta are overestimated when we employ the diffusive HLLE solver. We also find that the magnetorotational instability is less resolved when we employ the HLLE solver because of the solver's large numerical diffusivity. This causes a spurious enhancement both of magnetic winding resulting from large scale poloidal magnetic fields, and also of the energy of the outflow induced by magnetic pressure.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analytic properties of electromagnetic field of binary compact stars and electromagnetic precursors to gravitational waves

We analytically study the properties of the electromagnetic field in vacuum around close binary compact stars containing at least one neutron star. We show that the orbital motion of the neutron star induces high multipole modes of the electromagnetic field just before the merger. These modes are superimposed to form a spiral arm configuration and its edge is found to be a likely site for magnetic reconnection. These modes also enhance the total Poynting flux from neutron star binaries by a factor of 2--4. We also indicate that the electric field induced by the orbital motion lead to a magnetosphere around binaries and estimate its plasma density, which has a different parameter dependence than Goldreich-Julian density. With these properties, we discuss possible electromagnetic counterparts to gravitational wave events, and identify radio precursors, such as fast radio bursts, as the most promising observational targets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraint on the ejecta mass for a black hole-neutron star merger event candidate S190814bv

We derive the upper limit to the ejecta mass of S190814bv, a black hole-neutron star merger candidate, through the radiative transfer simulations for kilonovae with the realistic ejecta density profile as well as the detailed opacity and heating rate models. The limits to the ejecta mass strongly depend on the viewing angle. For the face-on observations ($\le45^\circ$), the total ejecta mass should be smaller than $0.1\,M_\odot$ for the average distance of S190814bv ($D=267$ Mpc), while larger mass is allowed for the edge-on observations. We also derive the conservative upper limits of the dynamical ejecta mass to be $0.02\,M_\odot$, $0.03\,M_\odot$, and $0.05\,M_\odot$ for the viewing angle $\le 20^\circ$, $\le 50^\circ$, and for $\le 90^\circ$, respectively. We show that the {\it iz}-band observation deeper than $22$ mag within 2 d after the GW trigger is crucial to detect the kilonova with the total ejecta mass of $0.06\,M_\odot$ at the distance of $D=300$ Mpc. We also show that a strong constraint on the NS mass-radius relation can be obtained if the future observations put the upper limit of $0.03\,M_\odot$ to the dynamical ejecta mass for a BH-NS event with the chirp mass smaller than $\lesssim 3\,M_\odot$ and effective spin larger than $\gtrsim 0.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extreme mass ratio inspirals on the equatorial plane in the adiabatic order

We compute gravitational waves from inspiraling stellar-mass compact objects on the equatorial plane of a massive spinning black hole (BH). Our inspiral orbits are computed by taking into account the adiabatic change of orbital parameters due to gravitational radiation in the lowest order in mass ratio. We employ an interpolation method to compute the adiabatic change at arbitrary points inside the region of orbital parameter space computed in advance. Using the obtained inspiral orbits and associated gravitational waves, we compute power spectra of gravitational waves and the signal-to-noise ratio (SNR) for several values of the BH spin, the masses of the binary, and the initial orbital eccentricity during a hypothetical three-yrs LISA observation before final plunge. We find that (i) the SNR increases as the BH spin and the mass of the compact object increase for the BH mass $M \agt 10^6M_\odot$, (ii) the SNR has a maximum for $M \approx 10^6M_\odot$, and (iii) the SNR increases as the initial eccentricity increases for $M=10^6M_\odot$. We also show that incorporating the contribution from the higher multipole modes of gravitational waves is crucial for enhancing the detection rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetosphere of an orbiting neutron star

We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed within its surrounding magnetosphere. A detailed study of the resulting magnetospheric configurations in our setting is presented, incorporating as well the effects due to neutron star spin and the misalignment of the magnetic and orbital axes. We find many features resembling those of pulsar magnetospheres for the orbiting neutron star --even in the absence of spin--, being of particular interest the development of a spiral current sheet that extends beyond the light cylinder. Then, we use recent advances in pulsar theory to estimate electromagnetic emissions produced at the reconnection regions of such current sheets.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Properties of neutrino transfer in a deformed remnant of neutron star merger

We study properties of neutrino transfer in a remnant of neutron star merger, consisting of a massive neutron star and a surrounding torus. We perform numerical simulations of the neutrino transfer by solving the Boltzmann equation with momentum-space angles and energies of neutrinos for snapshots of the merger remnant having elongated shapes. The evaluation of the neutrino distributions in the multi-dimensions enable us to provide the detailed information of angle and energy spectra and neutrino reaction rates. We demonstrate features of asymmetric neutrino fluxes from the deformed remnant and investigate the neutrino emission region by determining the neutrinosphere for each energy. We examine the emission and absorption of neutrinos to identify important ingredients of heating rates through neutrino irradiation. We show that the contributions of $μ$- and $τ$-types neutrinos are important for the heating in the region above the massive neutron star. We also examine the angle moments and the Eddington tensor calculated directly by the neutrino distribution functions and compare them with those obtained by a moment closure approach, which is often used in the study of neutrino-radiation hydrodynamics. We show that the components of the Eddington tensor have non-monotonic behaviors and the approximation of the closure relation may become inaccurate for high energy neutrinos, whose fluxes are highly aspherical due to the extended merger remnant.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Diversity of kilonova light curves

We perform radiative transfer simulations for kilonova in various situations, including the cases of prompt collapse to a black hole from neutron-star mergers, high-velocity ejecta possibly accelerated by magnetars, and a black hole-neutron star merger. Our calculations are done employing ejecta profiles predicted by numerical-relativity simulations and a new line list for all the r-process elements. We found that (i) the optical emission for binary neutron stars promptly collapsing to a black hole would be fainter by $\gtrsim1$-$2\,{\rm mag}$ than that found in GW170817, while the infrared emission could be as bright as that in GW170817 if the post-merger ejecta is as massive as $\approx0.01\,M_\odot$; (ii) the kilonova would be brighter than that observed in GW170817 for the case that the ejecta is highly accelerated by the electromagnetic energy injection from the remnant, but it would decline rapidly and the magnitude would become fainter than in GW170817 within a few days; (iii) the optical emission from a black hole-neutron star merger ejecta could be as bright as that observed in GW170817 for the case that sufficiently large amount of matter is ejected ($\gtrsim0.02\,M_\odot$), while the infrared brightness would be brighter by $1$-$2\,{\rm mag}$ at the same time. We show that the difference in the ejecta properties would be imprinted in the differences in the peak brightness and time of peak. This indicates that we may be able to infer the type of the central engine for kilonovae by observation of the peak in the multiple band.

0 citations0 reviewsNo code linkedOpen access