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

Yuichiro Sekiguchi

Yuichiro Sekiguchi 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.COnucl-thastro-ph.GAhep-ph

Catalog footprint

What is connected

28works
7topics
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

28 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.

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
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

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

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
preprint2016arXiv

Alpha-constrained QSE Nucleosynthesis in High-entropy and Fast-expanding Material

We investigate the nucleosynthesis process in high-entropy ($s/k_{\rm B}\gtrsim100$) and very fast-expanding ($τ_{\rm exp}\sim10^{-3}\ {\rm s}$) materials. In such a material with the electron fraction near 0.5, an interesting nucleosynthesis process occurs. In this process, the abundance distribution of heavy-nuclei of $A>100$ achieve quasi-statistical equilibrium (QSE) at high temperature and the abundances are frozen at the end of the nucleosynthesis. We explain this abundance distribution using the "alpha-constrained QSE" abundances formulated in this paper. We demonstrate that this nucleosynthesis would occur in neutrino-driven winds from massive proto-neutron stars in hypernovae, where $A\sim140$ $p$-nuclei are synthesized.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical mass ejection from the merger of asymmetric binary neutron stars: Radiation-hydrodynamics study in general relativity

We perform neutrino radiation-hydrodynamics simulations for the merger of asymmetric binary neutron stars in numerical relativity. Neutron stars are modeled by soft and moderately stiff finite-temperature equations of state (EOS). We find that the properties of the dynamical ejecta such as the total mass, neutron richness profile, and specific entropy profile depend on the mass ratio of the binary systems for a given EOS in a unique manner. For the soft EOS (SFHo), the total ejecta mass depends weakly on the mass ratio, but the average of electron number per baryon ($Y_e$) and specific entropy ($s$) of the ejecta decreases significantly with the increase of the degree of mass asymmetry. For the stiff EOS (DD2), with the increase of the mass asymmetry degree, the total ejecta mass significantly increases while the average of $Y_e$ and $s$ moderately decreases. We find again that only for the soft EOS (SFHo), the total ejecta mass exceeds $0.01M_\odot$ irrespective of the mass ratio chosen in this paper. The ejecta have a variety of electron number per baryon with its average approximately between $Y_e \sim 0.2$ and $\sim 0.3$ irrespective of the EOS employed, which is well-suited for the production of the r-process heavy elements (second and third peaks), although its averaged value decreases with the increase of the degree of mass asymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of Overheated Regions and Truncated Disks around Black Holes; Three-dimensional General Relativistic Radiation-magnetohydrodynamics Simulations

Using three-dimensional general relativistic radiation magnetohydrodynamics simulations of accretion flows around stellar mass black holes, we report that the relatively cold disk ($\gtrsim 10^{7}$K) is truncated near the black hole. Hot and less-dense regions, of which the gas temperature is $ \gtrsim 10^9$K and more than ten times higher than the radiation temperature (overheated regions), appear within the truncation radius. The overheated regions also appear above as well as below the disk, and sandwich the cold disk, leading to the effective Compton upscattering. The truncation radius is $\sim 30 r_{\rm g}$ for $\dot{M} \sim L_{\rm Edd}/c^2$, where $r_{\rm g}, \dot M, L_\mathrm{Edd}, c$ are the gravitational radius, mass accretion rate, Eddington luminosity, and light speed. Our results are consistent with observations of very high state, whereby the truncated disk is thought to be embedded in the hot rarefied regions. The truncation radius shifts inward to $\sim 10 r_{\rm g}$ with increasing mass accretion rate $\dot{M} \sim 100 L_{\rm Edd}/c^2$, which is very close to an innermost stable circular orbit. This model corresponds to the slim disk state observed in ultra luminous X-ray sources. Although the overheated regions shrink if the Compton cooling effectively reduces the gas temperature, the sandwich-structure does not disappear at the range of $\dot{M} \lesssim 100L_{\rm Edd}/c^2$. Our simulations also reveal that the gas temperature in the overheated regions depends on black hole spin, which would be due to efficient energy transport from black hole to disks through the Poynting flux, resulting gas heating.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational waves from supermassive stars collapsing to a supermassive black hole

We derive the gravitational waveform from the collapse of a rapidly rotating supermassive star (SMS) core leading directly to a seed of a supermassive black hole (SMBH) in axisymmetric numerical-relativity simulations. We find that the peak strain amplitude of gravitational waves emitted during the black-hole formation is $\approx 5 \times 10^{-21}$ at the frequency $f \approx 5$\,mHz for an event at the cosmological redshift $z=3$, if the collapsing SMS core is in the hydrogen-burning phase. Such gravitational waves will be detectable by space laser interferometric detectors like eLISA with signal-to-noise ratio $\approx 10$, if the sensitivity is as high as LISA for $f=1$--10\,mHz. The detection of the gravitational-wave signal will provide a potential opportunity for testing the direct-collapse scenario for the formation of a seed of SMBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of rigidly rotating supermassive stars against gravitational collapse

We revisit secular stability against quasi-radial collapse for rigidly rotating supermassive stars (SMSs) in general relativity. We suppose that the SMSs are in a nuclear-burning phase and can be modeled by polytropic equations of state with the polytropic index $n_p$ slightly smaller than $3$. The stability is determined in terms of the turning-point method. We find a fitting formula of the stability condition for the plausible range of $n_p$ ($2.95 \alt n_p \alt 3$) for SMSs. This condition reconfirms that, while non-rotating SMSs with mass $\sim 10^5M_\odot$--$10^6M_\odot$ may undergo a general-relativistically induced quasi-radial collapse, rigidly rotating SMSs with a ratio of rotational to gravitational potential energy ($β$) of $\sim 10^{-2}$ are likely to be stable against collapse unless they are able to accrete $\sim 5$ times more mass during the (relatively brief) hydrogen-burning phase of their evolution. We discuss implications of our results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-resolution magnetohydrodynamics simulation of black hole-neutron star merger: Mass ejection and short gamma-ray burst

We report results of a high-resolution numerical-relativity simulation for the merger of black hole-magnetized neutron star binaries on Japanese supercomputer "K". We focus on a binary that is subject to tidal disruption and subsequent formation of a massive accretion torus. We find the launch of thermally driven torus wind, subsequent formation of a funnel wall above the torus and a magnetosphere with collimated poloidal magnetic field, and high Blandford-Znajek luminosity. We show for the first time this picture in a self-consistent simulation. The turbulence-like motion induced by the non-axisymmetric magnetorotational instability as well as the Kelvin-Helmholtz instability inside the accretion torus works as an agent to drive the mass accretion and converts the accretion energy to thermal energy, which results in the generation of a strong wind. By an in-depth resolution study, we reveal that high resolution is essential to draw such a picture. We also discuss the implication for the r-process nucleosynthesis, the radioactively-powered transient emission, and short gamma-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nucleosynthesis in neutrino-driven winds in hypernovae

We investigate the nucleosynthesis in the neutrino-driven winds blown off from a $3M_\odot$ massive proto-neutron star (mPNS) temporarily formed during the collapse of a $100M_\odot$ star. Such mPNSs would be formed in hypernovae. We construct steady and spherically symmetric wind solutions. We set large neutrino luminosities of $\sim 10^{53}\ {\rm erg\ s^{-1}}$ and average energies of electron neutrinos and antineutrinos in the ranges of $ε_{ν_e}=9-16\ {\rm MeV}$ and $ε_{\barν_e}=11-18\ {\rm MeV}$ based on a recent numerical relativity simulation. The wind solutions indicate much shorter temperature-decrease timescale than that of the winds from ordinary PNSs and, depending on $ε_ν$, the winds can be both neutron-rich and proton-rich. In the neutron-rich wind, the $r$-process occurs and the abundance distribution of a fiducial wind model of the mPNS gives an approximate agreement with the abundance pattern of metal-poor weak $r$ star HD~122563, although the third-peak elements are produced only when the $\barν_e$ energy is much larger than the $ν_e$ energy. In the proton-rich wind, the strong $νp$-process occurs and $A>100$ nuclides are synthesized. The synthesized nuclei can be neutron-rich in some cases because the large neutrino luminosity of the mPNS supplies a sufficient amount of neutrons.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dynamical mass ejection from binary neutron star mergers: Radiation-hydrodynamics study in general relativity

We perform radiation-hydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese "K" supercomputer, taking into account neutrino cooling and heating by an updated leakage-plus-transfer scheme for the first time. Neutron stars are modeled by three modern finite-temperature equations of state (EOS) developed by Hempel and his collaborators. We find that the properties of the dynamical ejecta of the merger such as total mass, average electron fraction, and thermal energy depend strongly on the EOS. Only for a soft EOS (the so-called SFHo), the ejecta mass exceeds $0.01M_{\odot}$. In this case, the distribution of the electron fraction of the ejecta becomes broad due to the shock heating during the merger. These properties are well-suited for the production of the solar-like $r$-process abundance. For the other stiff EOS (DD2 and TM1), for which a long-lived massive neutron star is formed after the merger, the ejecta mass is smaller than $0.01M_{\odot}$, although broad electron-fraction distributions are achieved by the positron capture and the neutrino heating.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High resolution numerical-relativity simulations for the merger of binary magnetized neutron stars

We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70\,m, which yields the highest-resolution results among those derived so far. By an in-depth resolution study, we clarify several amplification mechanisms of magnetic fields during the binary neutron star merger for the first time. First, the Kelvin-Helmholtz instability developed in the shear layer at the onset of the merger significantly amplifies the magnetic fields. A hypermassive neutron star (HMNS) formed after the merger is then subject to the nonaxisymmetric magnetorotational instability, which amplifies the magnetic field in the HMNS. These two amplification mechanisms cannot be found with insufficient-resolution runs. We also show that the HMNS eventually collapses to a black hole surrounded by an accretion torus which is strongly magnetized at birth.

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

Production of all the r-process nuclides in the dynamical ejecta of neutron star mergers

Recent studies suggest that binary neutron star (NS-NS) mergers robustly produce the heavy r-process nuclei above the atomic mass number A ~ 130 because of their ejecta consisting of almost pure neutrons (electron fraction of Y_e < 0.1). However, little production of the lighter r-process nuclei (A = 90-120) conflicts with the spectroscopic results of r-process-enhanced Galactic halo stars. We present, for the first time, the result of nucleosynthesis calculations based on the fully general-relativistic simulation of a NS-NS merger with approximate neutrino transport. It is found that the bulk of the dynamical ejecta are appreciably shock-heated and neutrino-processed, resulting in a wide range of Y_e (= 0.09-0.45). The mass-averaged abundance distribution of calculated nucleosynthesis yields is in reasonable agreement with the full-mass range (A = 90-240) of the solar r-process curve. This implies, if our model is representative of such events, that the dynamical ejecta of NS-NS mergers can be the origin of the Galactic r-process nuclei. Our result also shows that the radioactive heating after ~ 1 day from the merging, giving rise to r-process-powered transient emission, is dominated by the beta-decays of several species close to stability with precisely measured half-lives. This implies that the total radioactive heating rate for such an event can be well constrained within about a factor of two if the ejected material has a solar-like r-process pattern.

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

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
preprint2012arXiv

Current Status of Numerical-Relativity Simulations in Kyoto

We describe the current status of our numerical simulations for the collapse of a massive stellar core to a BH and the BNS mergers, performed in the framework of full general relativity incorporating finite-temperature EOS and neutrino cooling. For the stellar core collapse simulation, we present the latest numerical results. We employed a purely nucleonic EOS (Shen-EOS). As an initial condition, we adopted a 100 $M_{\odot}$ presupernova model calculated by Umeda and Nomoto. Changing the degree of rotation for the initial condition, we clarify the strong dependence of the outcome of the collapse on this. When the rotation is rapid enough, the shock wave formed at the core bounce is deformed to be a torus-like shape. Then, the infalling matter is accumulated in the central region due to the oblique shock at the torus surface, hitting the PNS and dissipating the kinetic energy there. As a result, outflows can be launched. The PNS eventually collapses to a BH and an accretion torus is formed around it. We also found that the evolution of the BH and torus depends strongly on the rotation initially given. In the BNS merger simulations, we in addition employ an EOS incorporating a degree of freedom for hyperons. The numerical simulations show that for the purely nucleonic EOS, a HMNS with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $M \lesssim 3.0M_{\odot}$. By contrast, the formed HMNS collapses to a BH in a shorter time scale with the hyperonic EOS for $M \gtrsim 2.7M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$10\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 2--$6 \times 10^{-22}$ at $f\approx 2$--2.5 kHz for a source distance of 100 Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Explicit-Implicit Scheme for Relativistic Radiation Hydrodynamics

We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion (FLD) approximation. For an hyperbolic term, of which the time scale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational waves, neutrino emissions, and effects of hyperons in binary neutron star mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling. It is found that for the nucleonic and hyperonic EOS, a hyper massive neutron star (HMNS) with a long lifetime $(t_{\rm life}\gtrsim 10 {\rm ms})$ is the outcome for the total mass $\approx 2.7 M_\odot$. For the total mass $\approx 3 M_\odot$, a long-lived (short-lived with $t_{\rm life}\approx 3 {\rm ms}$) HMNS is the outcome for the nucleonic (hyperonic) EOS. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$ -- $6 \times 10^{53} {\rm erg /s}$ and the effective amplitude of gravitational waves from the HMNS is 1 -- $4\times 10^{-22}$ at $f\approx 2$ -- $3.2 {\rm kHz}$ for a source of distance of 100 Mpc. During the HMNS phase, characteristic frequencies of gravitational waves shift to a higher frequency for the hyperonic EOS in contrast to the nucleonic EOS in which they remain constant approximately. Our finding suggests that the effects of hyperons are well imprinted in gravitational wave and its detection will give us a potential opportunity to explore the composition of the neutron star matter. We present the neutrino luminosity curve when a black hole is formed as well.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effects of hyperons in binary neutron star mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling for the first time. It is found that even for the hyperonic EOS, a hypermassive neutron star is first formed after the merger for the typical total mass $\approx$ 2.7M\bigodot, and subsequently collapses to a black hole (BH). It is shown that hyperons play a substantial role in the post-merger dynamics, torus formation around the BH, and emission of gravitational waves (GWs). In particular, the existence of hyperons is imprinted in GWs. Therefore, GW observations will provide a potential opportunity to explore the composition of the neutron star matter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Formation of black hole and accretion disk in a massive high-entropy stellar core collapse

We present the first numerical result of fully general relativistic axisymmetric simulations for the collapse of a rotating high-entropy stellar core to a black hole and an accretion disk. The simulations are performed taking into account the relevant microphysics. We adopt as initial condition a spherical core with constant electron fraction ($Y_e = 0.5$) and entropy per baryon $s$ = 8 $k_B$, and angular velocity is superimposed. In the early phase, the core collapses in a homologous manner. Then, it experiences a weak bounce due to the gas pressure of free nucleons. Because the bounce is weak, the core collapses eventually to a black hole. Subsequent evolution depends on initial angular velocity. When the rotation is not fast, a geometrically thin (but optically thick) accretion disk is formed, and shock waves are formed in the inner part of the disk. For the moderately rotating case, the thin accretion disk expands eventually to be a geometrically thick torus after sufficient accumulation of the thermal energy generated at the shocks. Furthermore, convection occurs inside the torus. Neutrino luminosities vary violently with time because of the convective motion. For the rapidly rotating case, by contrast, a geometrically thick torus is formed soon after the black hole formation, and convective activity is weak due to the presence of epicyclic mode.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational waves and neutrino emission from the merger of binary neutron stars

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $\alt 3.0M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$8\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 4--$6 \times 10^{-22}$ at $f=2.1$--2.5 kHz for a source distance of 100 Mpc. We also present the neutrino luminosity curve when a black hole is formed for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An implementation of the microphysics in full general relativity : General relativistic neutrino leakage scheme

Performing fully general relativistic simulations taking account of microphysical processes is one of long standing problems in numerical relativity. One of main difficulties in implementation of weak interactions in the general relativistic framework lies on the fact that the characteristic timescale of weak interaction processes (the WP timescale) in hot dense matters is much shorter than the dynamical timescale. Numerically this means that stiff source terms appears in the equations so that an implicit scheme is in general necessary to stably solve the relevant equations. Otherwise a very short timestep will be required to solve them explicitly. Furthermore, in the relativistic framework, the Lorentz factor is coupled with the rest mass density and the energy density. The specific enthalpy is also coupled with the momentum. Due to these couplings, it is very complicated to recover the primitive variables and the Lorentz factor from conserved quantities. At the current status, no implicit procedure have been proposed except for the case of the spherical symmetry. Therefore, an approximate, explicit procedure is developed in the fully general relativistic framework in this paper as an first implementation of the microphysics toward a more realistic sophisticated model. The procedure is based on the so-called neutrino leakage schemes which is based on the property that the characteristic timescale in which neutrinos leak out of the system (the leakage timescale) is much longer than the WP timescale. In this paper, I present a detailed neutrino leakage scheme and a simple and stable method for solving the equations explicitly in the fully general relativistic framework. I also perform a test simulation to check the validity of the present method, showing that it works fairly well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring binary-neutron-star-merger scenario of short-gamma-ray bursts by gravitational-wave observation

We elucidate the feature of gravitational waves (GWs) from binary neutron star merger collapsing to a black hole by general relativistic simulation. We show that GW spectrum imprints the coalescence dynamics, formation process of disk, equation of state for neutron stars, total masses, and mass ratio. A formation mechanism of the central engine of short $γ$-ray bursts, which are likely to be composed of a black hole and surrounding disk, therefore could be constrained by GW observation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar core collapse in full general relativity with microphysics - Formulation and Spherical collapse test -

One of the longstanding issues in numerical relativity is to enable a simulation taking account of microphysical processes (e.g., weak interactions and neutrino cooling). We develop an approximate and explicit scheme in the fully general relativistic framework as a first implementation of the microphysics toward a more realistic and sophisticated modeling. In this paper, we describe in detail a method for implementation of a realistic equation of state, the electron capture and the neutrino cooling in a multidimensional, fully general relativistic code. The procedure is based on the so-called neutrino leakage scheme. To check the validity of the code, we perform a two dimensional (2D) simulation of spherical stellar core collapse. Until the convective activities set in, our results approximately agree, or at least are consistent, with those in the previous so-called state-of-the-art simulations. In particular, the radial profiles of thermodynamical quantities and the time evolution of the neutrino luminosities agree quantitatively. The convection is driven by negative gradients of the entropy per baryon and the electron fraction as in the previous 2D Newtonian simulations. We clarify which gradient is more responsible for the convection. Gravitational waves from the convection are also calculated. We find that the characteristic frequencies of the gravitational-wave spectra are distributed for higher frequencies than those in Newtonian simulations due to the general relativistic effects.

0 citations0 reviewsNo code linkedOpen access