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

Ataru Tanikawa

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.SRastro-ph.IMphysics.comp-phastro-ph.COastro-ph.EPastro-ph

Catalog footprint

What is connected

30works
8topics
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

30 published item(s)

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Can Population III stars be major origins of both merging binary black holes and extremely metal poor stars?

Population (Pop) III stars, first stars, or metal-free stars are made of primordial gas. We have examined if they can be dominant origins of merging binary black holes (BHs) and extremely metal-poor stars. The abundance pattern of EMP stars is helpful to trace back the properties of Pop III stars. We have confirmed previous arguments that the observed BH merger rate needs Pop III star formation efficiency 10 times larger than theoretically predicted values, while the cosmic reionization history still permits such a high Pop III star formation efficiency. On the other hand, we have newly found that the elemental abundance pattern of EMP stars only allows the Pop III initial mass function with the minimum mass of $\sim 15 - 27$ $M_\odot$. In other words, the minimum mass must not deviate largely from the critical mass below and above which Pop III stars leave behind neutron stars and BHs, respectively. Pop III stars may be still a dominant origin of merging binary BHs but our study has reduced the allowed parameter space under a hypothesis that EMP stars are formed from primordial gas mixed with Pop III supernova ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compact Object Mergers in Hierarchical Triples from Low-Mass Young Star Clusters

A binary star orbited by an outer companion constitutes a hierarchical triple system. The outer body may excite the eccentricity of the inner binary through the von~Zeipel-Lidov-Kozai (ZLK) mechanism, triggering the gravitational wave (GW) coalescence of the inner binary when its members are compact objects. Here, we study a sample of hierarchical triples with an inner black hole (BH) -- BH binary, BH -- neutron star (NS) binary, and BH -- white dwarf (WD) binary, formed via dynamical interactions in low-mass young star clusters. Our sample of triples was obtained self-consistently from direct N-body simulations of star clusters which included up-to-date stellar evolution. We find that the inner binaries in our triples cannot merge via GW radiation alone, and the ZLK mechanism is essential to trigger their coalescence. Contrary to binaries assembled dynamically in young star clusters, binary BHs merging in triples have preferentially low mass ratios (q ~ 0.3) and higher primary masses (m_p > 40 MSun). We derive a local merger rate density of 0.60, 0.11 and 0.5 yr^-1 Gpc^-3 for BH-BH, BH-NS and BH-WD binaries, respectively. Additionally, we find that merging binaries have high eccentricities across the GW spectrum, including the LIGO-Virgo-KAGRA (LVK), LISA, and DECIGO frequencies. About 7% of BH-BH and 60% of BH-NS binaries will have detectable eccentricities in the LVK band. Our results indicate that the eccentricity and the mass spectrum of merging binaries are the strongest features for the identification of GW mergers from triples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Detectablity of Black Hole Binaries with Gaia: Dependence on Binary Evolution Models

Astrometric satellite Gaia is expected to observe non-interacting black hole (BH) binaries with luminous companions (LCs) (hereafter BH-LC binaries), a different population from BH X-ray binaries previously discovered. The detectability of BH-LC binaries with Gaia might be dependent on binary evolution models. We investigated the Gaia's detectability of BH-LC binaries formed through isolated binary evolution by means of binary population synthesis technique, and examined its dependence on single and binary star models: supernova models, common envelope (CE) ejection efficiency $α$, and BH natal kick models. We estimated that $1.1$ -- $46$ BH-LC binaries can be detected within five-year observation, and found that $α$ has the largest impacts on the detectable number. In each model, observable and intrinsic BH-LC binaries have similar distributions. Therefore, we found three important implications: (1) if the lower BH mass gap is not intrinsic (i.e. $3$ -- $5 M_\odot$ BHs exist), Gaia will observe $\leq 5 M_\odot$ BHs, (2) we may observe short orbital period binaries with light LCs if CE efficiency is significantly high, and (3) we may be able to identify the existence of natal kick from eccentricity distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational wave of intermediate-mass black holes in Population III star clusters

Previous theoretical studies suggest that the Population III (Pop3) stars tend to form in extremely metal poor gas clouds with approximately $10^5 M_\odot$ embedded in mini dark matter halos. Very massive stars can form via multiple collisions in Pop3 star clusters and eventually evolve to intermediate-mass black holes (IMBHs). In this work, we conduct star-by-star $N$-body simulations for modelling the long-term evolution of Pop3 star clusters. We find that if the mini dark matter halos can survive today, these star clusters can avoid tidal disruption by the galactic environment and can efficiently produce IMBH-BH mergers among a wide range of redshift from 0 to 20. The average gravitational wave event rate is estimated to be $0.1-0.8~\mathrm{yr}^{-1} \mathrm{Gpc}^{-3}$, and approximately $40-80$ percent of the mergers occur at high redshift ($z>6$). The characteristic strain shows that a part of low-redshift mergers can be detected by LISA, TianQin, and Taiji, whereas most mergers can be covered by DECIGO and advanced LIGO/VIRGO/Kagra. Mergers with pair-instability BHs have a rate of approximately $0.01-0.15$~yr$^{-1}$~Gpc$^{-3}$, which can explain the GW190521-like events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Merger rate density of binary black holes through isolated Population I, II, III and extremely metal-poor binary star evolution

We investigate the formation of merging binary black holes (BHs) through isolated binary evolution, performing binary population synthesis calculations covering an unprecedentedly wide metallicity range of Population (Pop) I, II, III, and extremely metal-poor (EMP) binary stars. We find that the predicted merger rate density and primary BH mass ($m_1$) distribution are consistent with the gravitational wave (GW) observations. Notably, Pop III and EMP ($< 10^{-2}$ $Z_\odot$) binary stars yield most of the pair instability (PI) mass gap events with $m_1 = 65$--$130$ $M_\odot$. Pop III binary stars contribute more to the PI mass gap events with increasing redshift, and all the PI mass gap events have the Pop III origin at redshifts $\gtrsim 8$. Our result can be assessed by future GW observations in the following two points. First, there are no binary BHs with $m_1=100$--$130$ $M_\odot$ in our result, and thus the $m_1$ distribution should suddenly drop in the range of $m_1=100$--$130$ $M_\odot$. Second, the PI mass gap event rate should increase toward higher redshift up to $\sim 11$, since those events mainly originate from the Pop III binary stars. We find that the following three assumptions are needed to reproduce the current GW observations: a top-heavy stellar initial mass function and the presence of close binary stars for Pop III and EMP binary stars, and inefficient convective overshoot in the main-sequence phase of stellar evolution. Without any of the above, the number of PI mass gap events becomes too low to reproduce current GW observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

MOCCA-SURVEY Database -- I. Tidal disruption events of white dwarfs in globular clusters and young massive clusters

We exploit more than 1000 star cluster Monte Carlo models from the MOCCA-SURVEY Database I, to infer the local rate density of white dwarf (WD) tidal disruption events (TDEs) in globular clusters (GCs) and young massive clusters (YMCs). We suggest that the WD TDE rate for GCs and YMCs in the local universe is $\sim 90$-$500~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, with 90 % of WD TDEs occurring in GCs. The total WD TDE rate density is $\sim 9$-$50$ times larger than estimated before. Our results show that thermonuclear explosions induced by WD TDEs can be observed at a rate of $\sim 100$-$550~{\rm yr}^{-1}$ by the next generation optical surveys, such as the Legacy Survey of Space & Time by the Vera C. Rubin Observatory. We also find that massive WDs are preferentially disrupted due to mass segregation, and that 20 % of exploding WDs have $\gtrsim 1.0 M_\odot$ despite of the small population of such WDs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes

We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory. Our formalism can model the eccentricity during the common envelope inspiral, and it gives results consistent with smoothed particles hydrodynamical simulations. We apply our formalism to common envelope events from binary population synthesis models and find that the final eccentricity distribution resembles the observed distribution of post-common-envelope binaries. Our model can be used for time-resolved common-envelope evolution in population synthesis calculations or as part of binary interactions in direct N-body simulations of star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The double detonation of a double degenerate system, from Type Ia supernova explosion to its supernova remnant

Type Ia supernovae (SNe) are believed to be caused by the thermonuclear explosion of a white dwarf (WD), but the nature of the progenitor system(s) is still unclear. Recent theoretical and observational developments have led to renewed interest in double degenerate models, in particular the "helium-ignited violent merger" or "dynamically-driven double-degenerate double-detonation" (D$^6$). In this paper we take the output of an existing D$^6$ SN model and carry it into the supernova remnant (SNR) phase up to 4000 years after the explosion, past the time when all the ejecta have been shocked. Assuming a uniform ambient medium, we reveal specific signatures of the explosion mechanism and spatial variations intrinsic to the ejecta. The first detonation produces an ejecta tail visible at early times, while the second detonation leaves a central density peak in the ejecta that is visible at late times. The SNR shell is off-centre at all times, because of an initial velocity shift due to binary motion. The companion WD produces a large conical shadow in the ejecta, visible in projection as a dark patch surrounded by a bright ring. This is a clear and long-lasting feature that is localized, and its impact on the observed morphology is dependent on the viewing angle of the SNR. These results offer a new way to diagnose the explosion mechanism and progenitor system using observations of a Type Ia SNR.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spin distribution of binary black holes formed in open clusters

We performed direct N-body simulations of open clusters with four different metallicities. To investigate the effective spins of merging binary black holes (BBHs) originated from these open clusters, we calculated the spin evolution of Wolf-Rayet (WR) stars with close companion stars (BBH progenitors), taking into account stellar wind mass loss and tidal spin-up of the WR stars. We found that BBH progenitors with smaller semi-major axes evolve to merging BBHs with greater effective spins because of strong tidal forces. In the local Universe, about 16% of merging BBHs get effective spins larger than 0.1 even if BHs and their progenitors do not get spin angular momenta other than tidal forces exerted by their companion stars. If we assume that WR stars have flat and isotropic distribution of dimensionless spins just after common envelope phases, the effective spin distribution of merging BBHs is similar to that inferred from gravitational wave observations with LIGO and Virgo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fitting formulae for evolution tracks of massive stars under extreme metal poor environments for population synthesis calculations and star cluster simulations

We have devised fitting formulae for evolution tracks of massive stars with $8 \lesssim M/M_\odot \lesssim 160$ under extreme metal poor (EMP) environments for $\log (Z/Z_\odot) = -2, -4, -5, -6$, and $-8$, where $M_\odot$ and $Z_\odot$ are the solar mass and metallicity, respectively. Our fitting formulae are based on reference stellar models which we have newly obtained by simulating the time evolutions of EMP stars. Our fitting formulae take into account stars ending with blue supergiant (BSG) stars, and stars skipping Hertzsprung gap (HG) phases and blue loops, which are characteristics of massive EMP stars. In our fitting formulae, stars may remain BSG stars when they finish their core Helium burning (CHeB) phases. Our fitting formulae are in good agreement with our stellar evolution models. We can use these fitting formulae on the SSE, BSE, NBODY4, and NBODY6 codes, which are widely used for population synthesis calculations and star cluster simulations. These fitting formulae should be useful to make theoretical templates of binary black holes formed under EMP environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation rate of LB-1-like systems through dynamical interactions

We estimate formation rates of LB-1-like systems through dynamical interactions in the framework of the theory of stellar evolution before the discovery of the LB-1 system. The LB-1 system contains $\sim 70M_\odot$ black hole (BH), so-called pair instability (PI)-gap BH, and B-type star with solar metallicity, and has nearly zero eccentricity. The most efficient formation mechanism is as follows. In an open cluster, a naked helium (He) star (with $\sim 20M_\odot$) collides with a heavy main-sequence (MS) star (with $\sim 50M_\odot$) which has a B-type companion. The collision results in a binary consisting of the collision product and B-type star with a high eccentricity. The binary can be circularized through the dynamical tide with radiative damping of the collision-product envelope. Finally, the collision product collapses to a PI-gap BH, avoiding pulsational pair instability and pair instability supernovae because its He core is as massive as the pre-colliding naked He star. We find that the number of LB-1-like systems in the Milky Way galaxy is $\sim 0.01 (ρ_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$, where $ρ_{\rm oc}$ is the initial mass densities of open clusters. If we take into account LB-1-like systems with O-type companion stars, the number increases to $\sim 0.03 (ρ_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$. This mechanism can form LB-1-like systems at least 10 times more efficiently than the other mechanisms: captures of B-type stars by PI-gap BHs, stellar collisions between other type stars, and stellar mergers in hierarchical triple systems. We conclude that no dynamical mechanism can explain the presence of the LB-1 system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gaia's Detectability of Black Hole-Main Sequence Star Binaries Formed in Open Clusters

Black hole-main sequence star (BH-MS) binaries are one of the targets of the future data releases of the astrometric satellite {\it Gaia}. They are supposed to be formed in two main sites: a galactic field and star clusters. However, previous work has never predicted the number of BH-MS binaries originating in the latter site. In this paper, we estimate the number of BH-MS binaries formed in open clusters and detectable with {\it Gaia} based on the results of {\it N}-body simulations. By considering interstellar extinction in the Milky Way (MW) and observational constraints, we predict $\sim 10$ BH-MS binaries are observable. We also find that chemical abundance patterns of companion MSs will help us to identify the origin of the binaries as star clusters. Such MSs are not polluted by outflows of the BH progenitors, such as stellar winds and supernova ejecta. Chemical anomalies might be a good test to confirm the origin of binaries with relatively less massive MSs ($\lesssim 5M_{\odot}$), orbital periods ($\sim 1.5\;$year) and higher eccentricities ($e \gtrsim 0.1$).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Merger rate density of binary black holes formed in open clusters

Several binary black holes (BBHs) have been observed using gravitational wave detectors. For the formation mechanism of BBHs, two main mechanisms, isolated binary evolution and dynamical formation in dense star clusters, have been suggested. Future observations are expected to provide more information about BBH distributions, and it will help us to distinguish the two formation mechanisms. For the star cluster channel, globular clusters have mainly been investigated. However, recent simulations have suggested that BBH formation in open clusters is not negligible. We estimate a local merger rate density of BBHs originated from open clusters using the results of our $N$-body simulations of open clusters with four different metallicities. We find that the merger rate per cluster is the highest for our 0.1 solar metallicity model. Assuming a cosmic star formation history and a metallicity evolution with dispersion, we estimate the local merger rate density of BBHs originated from open clusters to be $\sim 70~{\rm yr}^{-1} {\rm Gpc}^{-3}$. This value is comparable to the merger rate density expected from the first and second observation runs of LIGO and Virgo. In addition, we find that BBH mergers obtained from our simulations can reproduce the distribution of primary mass and mass ratio of merging BBHs estimated from the LIGO and Virgo observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rapid Transients Originating from Thermonuclear Explosions in Helium White Dwarf Tidal Disruption Events

We study the emission properties of thermonuclear explosions in a helium white dwarf (WD) tidal disruption event (TDE). We consider a TDE where a 0.2 $M_{\odot}$ helium WD is disrupted by a $10^{2.5}\,M_{\odot}$ intermediate-mass black hole (IMBH). The helium WD is not only tidally disrupted but is also detonated by the tidal compression and by succeeding shocks. We focus on the emission powered by radioactive nuclei in the unbound TDE ejecta. We perform hydrodynamic simulations coupled with nuclear reactions, post-process detailed nucleosynthesis calculations, and then radiative transfer simulations. We thus derive multi-band light curves and spectra. The helium WD TDE shows rapid ($Δt_{1\mathrm{mag}}\simeq5\text{--}10$ days) and relatively faint ($L_{\mathrm{peak}}\simeq10^{42}\,\mathrm{erg}\,\mathrm{s}^{-1}$) light curves, because the ejecta mass and $^{56}$Ni mass are low ($0.12\,M_{\odot}$ and $0.03\,M_{\odot}$, respectively). The spectra show strong calcium and Fe-peak features and very weak silicon features, reflecting the peculiar elemental abundance. The key feature is the Doppler shift of the spectral lines up to $\simeq\pm12,000\,\mathrm{km}\,\mathrm{s}^{-1}$, depending on the viewing angle, due to the bulk motion of the ejecta. Our model matches well with two rapid and faint transients reported in Pursiainen et al. (2018). The particular model presented here does not match with observed SNe Iax, calcium-rich transients, or .Ia explosion candidates, either in the spectra or light curves. However, we expect a large variety of the observational signatures once a wide range of the WD/black hole masses and orbital parameters are considered. This study contributes to the search for WD TDEs with current and upcoming surveys, and to the identification of IMBHs as disrupters in the TDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An optical search for transients lasting a few seconds

Using a prototype of the Tomo-e Gozen wide-field CMOS mosaic camera, we acquire wide-field optical images at a cadence of 2 Hz and search them for transient sources of duration 1.5 to 11.5 seconds. Over the course of eight nights, our survey encompasses the equivalent of roughly two days on one square degree, to a fluence equivalent to a limiting magnitude about $V = 15.6$ in a 1-second exposure. After examining by eye the candidates identified by a software pipeline, we find no sources which meet all our criteria. We compute upper limits to the rate of optical transients consistent with our survey, and compare those to the rates expected and observed for representative sources of ephemeral optical light.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Implementation and performance of FDPS: A Framework Developing Parallel Particle Simulation Codes

We present the basic idea, implementation, measured performance and performance model of FDPS (Framework for developing particle simulators). FDPS is an application-development framework which helps the researchers to develop particle-based simulation programs for large-scale distributed-memory parallel supercomputers. A particle-based simulation program for distributed-memory parallel computers needs to perform domain decomposition, redistribution of particles, and gathering of particle information for interaction calculation. Also, even if distributed-memory parallel computers are not used, in order to reduce the amount of computation, algorithms such as Barnes-Hut tree method should be used for long-range interactions. For short-range interactions, some methods to limit the calculation to neighbor particles are necessary. FDPS provides all of these necessary functions for efficient parallel execution of particle-based simulations as "templates", which are independent of the actual data structure of particles and the functional form of the interaction. By using FDPS, researchers can write their programs with the amount of work necessary to write a simple, sequential and unoptimized program of O(N^2) calculation cost, and yet the program, once compiled with FDPS, will run efficiently on large-scale parallel supercomputers. A simple gravitational N-body program can be written in around 120 lines. We report the actual performance of these programs and the performance model. The weak scaling performance is very good, and almost linear speedup was obtained for up to the full system of K computer. The minimum calculation time per timestep is in the range of 30 ms (N=10^7) to 300 ms (N=10^9). These are currently limited by the time for the calculation of the domain decomposition and communication necessary for the interaction calculation. We discuss how we can overcome these bottlenecks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The critical mass ratio of double white dwarf binaries for violent merger-induced Type Ia supernova explosions

Mergers of carbon-oxygen (CO) white dwarfs (WDs) are considered as one of the potential progenitors of type Ia supernovae (SNe Ia). Recent hydrodynamical simulations showed that the less massive (secondary) WD violently accretes onto the more massive (primary) one, carbon detonation occurs, the detonation wave propagates through the primary, and the primary finally explodes as a sub-Chandrasekhar mass SN Ia. Such an explosion mechanism is called the violent merger scenario. Based on the smoothed particle hydrodynamics (SPH) simulations of merging CO WDs, we derived more stringent critical mass ratio (qcr) leading to the violent merger scenario than the previous results. We conclude that this difference mainly comes from the differences in the initial condition, synchronously spinning of WDs or not. Using our new results, we estimated the brightness distribution of SNe Ia in the violent merger scenario and compared it with previous studies. We found that our new qcr does not significantly affect the brightness distribution. We present the direct outcome immediately following CO WD mergers for various primary masses and mass ratios. We also discussed the final fate of the central system of the bipolar planetary nebula Henize 2-428, which was recently suggested to be a double CO WD system whose total mass exceeds the Chandrasekhar-limiting mass, merging within the Hubble time. Even considering the uncertainties in the proposed binary parameters, we concluded that the final fate of this system is almost certainly a sub-Chandrasekhar mass SN Ia in the violent merger scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A systematic study of carbon-oxygen white dwarf mergers: mass combinations for Type Ia supernovae

Mergers of two carbon-oxygen (CO) white dwarfs (WDs) have been considered as progenitors of Type Ia supernovae (SNe Ia). Based on smoothed particle hydrodynamics (SPH) simulations, previous studies claimed that mergers of CO WDs lead to an SN Ia explosion either in the dynamical merger phase or stationary rotating merger remnant phase. However, the mass range of CO WDs that lead to an SN Ia has not been clearly identified yet. In the present work, we perform systematic SPH merger simulations for the WD masses ranging from $0.5~M_{\odot}$ to $1.1~M_{\odot}$ with higher resolutions than the previous systematic surveys and examine whether or not carbon burning occurs dynamically or quiescently in each phase. We further study the possibility of SN Ia explosion and estimate the mass range of CO WDs that lead to an SN Ia. We found that when the both WDs are massive, i.e., in the mass range of $0.9~M_{\odot} {\le} M_{1,2} {\le} 1.1~M_{\odot}$, they can explode as an SN Ia in the merger phase. On the other hand, when the more massive WD is in the range of $0.7~M_{\odot} {\le} M_{1} {\le} 0.9~M_{\odot}$ and the total mass exceeds $1.38~M_{\odot}$, they can finally explode in the stationary rotating merger remnant phase. We estimate the contribution of CO WD mergers to the entire SN Ia rate in our galaxy to be of ${\lt} 9\%$. So, it might be difficult to explain all galactic SNe Ia by CO WD mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hydrodynamical evolution of merging carbon-oxygen white dwarfs: their pre-supernova structure and observational counterparts

We perform smoothed particle hydrodynamics (SPH) simulations for merging binary carbon-oxygen (CO) white dwarfs (WDs) with masses of $1.1$ and $1.0$ $M_\odot$, until the merger remnant reaches a dynamically steady state. Using these results, we assess whether the binary could induce a thermonuclear explosion, and whether the explosion could be observed as a type Ia supernova (SN Ia). We investigate three explosion mechanisms: a helium-ignition following the dynamical merger (`helium-ignited violent merger model'), a carbon-ignition (`carbon-ignited violent merger model'), and an explosion following the formation of the Chandrasekhar mass WD (`Chandrasekhar mass model'). An explosion of the helium-ignited violent merger model is possible, while we predict that the resulting SN ejecta are highly asymmetric since its companion star is fully intact at the time of the explosion. The carbon-ignited violent merger model can also lead to an explosion. However, the envelope of the exploding WD spreads out to $\sim 0.1R_\odot$; it is much larger than that inferred for SN 2011fe ($< 0.1R_\odot $) while much smaller than that for SN 2014J ($\sim 1R_\odot$). For the particular combination of the WD masses studied in this work, the Chandrasekhar mass model is not successful to lead to an SN Ia explosion. Besides these assessments, we investigate the evolution of unbound materials ejected through the merging process (`merger ejecta'), assuming a case where the SN Ia explosion is not triggered by the helium- or carbon-ignition during the merger. The merger ejecta interact with the surrounding interstellar medium, and form a shell. The shell has a bolometric luminosity of more than $2 \times 10^{35}$ ergs$^{-1}$ lasting for $\sim 2 \times 10^4$ yr. If this is the case, Milky Way should harbor about $10$ such shells at any given time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Merger Criteria of Multiple Massive Black Holes and the Impact on the Host Galaxy

We perform N-body simulations on a multiple massive black hole (MBH) system in a host galaxy to derive the criteria for successive MBH merger. The calculations incorporate the dynamical friction by stars and general relativistic effects as pericentre shift and gravitational wave recoil. The orbits of MBHs are pursed down to ten Schwarzschild radii (~ 1 AU). As a result, it is shown that about a half of MBHs merge during 1 Gyr in a galaxy with mass $10^{11}M_{\odot}$ and stellar velocity dispersion 240 km/s, even if the recoil velocity is two times as high as the stellar velocity dispersion. The dynamical friction allows a binary MBH to interact frequently with other MBHs, and then the decay of the binary orbits leads to the merger through gravitational wave radiation, as shown by Tanikawa & Umemura (2011). We derive the MBH merger criteria for the masses, sizes, and luminosities of host galaxies. It is found that the successive MBH mergers are expected in bright galaxies, depending on redshifts. Furthermore, we find that the central stellar density is reduced by the sling-shot mechanism and that high-velocity stars with ~ 1000 km/s are generated intermittently in extremely radial orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical evolution of stellar-mass black holes in dense stellar clusters: estimate for merger rate of binary black holes originating from globular clusters

We have performed N-body simulations of globular clusters (GCs) in order to estimate a detection rate of mergers of Binary stellar-mass Black Holes (BBHs) by means of gravitational wave (GW) observatories. For our estimate, we have only considered mergers of BBHs which escape from GCs (BBH escapers). BBH escapers merge more quickly than BBHs inside GCs because of their small semi-major axes. N-body simulation can not deal with a GC with the number of stars N ~ 10^6 due to its high calculation cost. We have simulated dynamical evolution of small-N clusters (10^4 <~ N <~ 10^5), and have extrapolated our simulation results to large-N clusters. From our simulation results, we have found the following dependence of BBH properties on N. BBHs escape from a cluster at each two-body relaxation time at a rate proportional to N. Semi-major axes of BBH escapers are inversely proportional to N, if initial mass densities of clusters are fixed. Eccentricities, primary masses, and mass ratios of BBH escapers are independent of N. Using this dependence of BBH properties, we have artificially generated a population of BBH escapers from a GC with N ~ 10^6, and have estimated a detection rate of mergers of BBH escapers by next-generation GW observatories. We have assumed that all the GCs are formed 10 or 12Gyrs ago with their initial numbers of stars N_i=5 x 10^5 -- 2 x 10^6 and their initial stellar mass densities inside their half-mass radii ρ_h,i=6 x 10^3 -- 10^6M_sun pc^-3. Then, the detection rate of BBH escapers is 0.5 -- 20 yr^-1 for a BH retention fraction R_BH=0.5. A few BBH escapers are components of hierarchical triple systems, although we do not consider secular perturbation on such BBH escapers for our estimate. Our simulations have shown that BHs are still inside some of GCs at the present day. These BHs may marginally contribute to BBH detection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Few-Body Modes of Binary Formation in Core Collapse

At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further dissecting the subsequent steps leading to hard binary formation. For this purpose, we add some analysis tools in order to make the study less subjective. We find that the conventional treatment does remain valid for direct three-body scattering, but fails for resonant three-body scattering. Especially democratic resonance scattering, which forms an important part of the analytical theory of three-body binary formation, takes too much space and time to be approximated as being isolated, in the context of a cluster core around core collapse. We conclude that, while three-body encounters can be analytically approximated as isolated, subsequent strong perturbations typically occur whenever those encounters give rise to democratic resonances. We present analytical estimates postdicting our numerical results. If we only had been a bit more clever, we could have predicted this qualitative behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phantom-GRAPE: numerical software library to accelerate collisionless $N$-body simulation with SIMD instruction set on x86 architecture

(Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed. Using an Intel Core i7--2600 processor, we measure the performance of our library for both the forces. In the case of Newton's forces, we achieve 2 x 10^9 interactions per second with 1 processor core, which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With 4 processor cores, we obtain the performance of 8 x 10^9 interactions per second. In the case of the arbitrarily shaped forces, we can calculate 1 x 10^9 and 4 x 10^9 interactions per second with 1 and 4 processor cores, respectively. The performance with 1 processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend weakly on the number of particles. It is good contrast with the fact that the performance of force calculations accelerated by GPUs depends strongly on the number of particles. Substantially weak dependence of the performance on the number of particles is suitable to collisionless N-body simulations, since these simulations are usually performed with sophisticated N-body solvers such as Tree- and TreePM-methods combined with an individual timestep scheme. Collisionless N-body simulations accelerated with our library have significant advantage over those accelerated by GPUs, especially on massively parallel environments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

N-body simulation for self-gravitating collisional systems with a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of 20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved 90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N 105 on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Successive Merger of Multiple Massive Black Holes in a Primordial Galaxy

Using highly-accurate $N$-body simulations, we explore the evolution of multiple massive black holes (hereafter, MBHs) in a primordial galaxy that is composed of stars and MBHs. The evolution is pursed with a fourth-order Hermite scheme, where not only three-body interaction of MBHs but also dynamical friction by stars are incorporated. Initially, ten MBHs with equal mass of $10^7M_\odot$ are set in a host galaxy with $10^{11}M_\odot$. It is found that 4 - 6 MBHs merge successively within 1 Gyr, emitting gravitational wave radiation. The key process for the successive merger of MBHs is the dynamical friction by field stars, which enhances three-body interactions of MBHs when they enter the central regions of the galaxy. The heaviest MBH always composes a close binary at the galactic center, which shrinks owing to the angular momentum transfer by the third MBH and eventually merges. The angular momentum transfer by the third MBH is due to the sling-shot mechanism. We find that the secular Kozai mechanism does not work for a binary to merge if we include the relativistic pericenter shift. The simulations show that a multiple MBH system can produce a heavier MBH at the galactic center purely through $N$-body process. This merger path can be of great significance for the growth of MBHs in a primordial galaxy. The merger of multiple MBHs may be a potential source of gravitational waves for the Laser Interferometer Space Antenna (LISA) and pulsar timing.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Unexpected Formation Modes of the First Hard Binary in Core Collapse

The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that both assumptions are often incorrect: 1) quite a few three-body interactions produce a hard binary from scratch; 2) and in many cases there are more than three bodies directly and simultaneously involved in the production of the first binary. The main reason for the discrepancies is that the core of a star cluster, at the first deep collapse, contains typically only five or so stars. Therefore, the homogeneous background assumption, which still would be reasonable for, say, 25 stars, utterly breaks down. There have been some speculations in this direction, but we demonstrate this result here explicitly, for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of Hardness of Primordial Binaries on Evolution of Star Clusters

We investigate effects of hardness of primordial binaries on whole evolution of star clusters by means of N-body simulations. Using newly developed code, GORILLA, we simulated eleven N=16384 clusters with primordial binaries whose binding energies are equal in each cluster in range of 1-300kT_0, where 1.5kT_0 is average stellar kinetic energy at the initial time. We found that, in both soft (< 3kT_0) and hard (> 300kT_0) limits, clusters experience deep core collapse. In the intermediate hardness (10-100kT_0), the core collapses halt halfway due to an energy releases of the primordial binaries. The core radii at the halt can be explained by their energy

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mass-Loss Timescale of Star Clusters in an External Tidal Field. II. Effect of Mass Profile of Parent Galaxy

We investigate the long-term dynamical evolution of star clusters in a steady tidal field produced by its parent galaxy. In this paper, we focus on the influence of mass profile of the parent galaxy. The previous studies were done with the simplification where the parent galaxy was expressed by point mass. We express different mass profiles of the parent galaxy by the tidal fields in which the ratios of the epicyclic frequency to the angular velocity are different. We compare the mass-loss timescale of star clusters whose tidal radii are identical but in parent galaxies with different mass profile, by means of orbits calculations in fixed cluster potential and N-body simulations. In this situation, a cluster rotates around the parent galaxy more rapidly as the parent galaxy has shallower mass profile. We found that the mass-loss timescale increase 20% and 50% for the cases that the mass density profile of the parent galaxies are proportional to R^-2 and R^-1.5 where R is the distance from the galaxy center, compared to the point-mass case, in moderately strong tidal field. Counterintuitively, a cluster which rotates around the parent galaxy more rapidly has a longer lifetime. The increase of lifetime is due to the fact that the fraction occupied by regular-like orbit increases in shallower profile. Finally, we derive an evaluation formula for the mass-loss timescale of clusters. Our formula can explain a property of the population of the observed galactic globular clusters that their half-mass radii become smaller as their distances from the galactic center become smaller.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Mass Loss Timescale of Star Clusters in External Tidal Field

We investigate evolution of star clusters in external tidal field by means of $N$-body simulations. We followed seven sets of cluster models whose central concentration and strength of the tidal field are different. We found that the mass loss timescale due to escape of stars, $t_{mloss}$, and its dependence on the two-body relaxation timescale, $t_{rh,i}$, are determined by the strength of the tidal field. The logarithmic slope [$= dln(t_{mloss})/dln(t_{rh,i})$] approaches to near unity for the cluster models in weaker tidal field. The timescale and the dependence are almost independent of the central concentration for clusters in the tidal field of the same strength. In our results, the scaling found by Baumgardt (2001) can be seen only in the cluster models with moderately strong tidal field.

0 citations0 reviewsNo code linkedOpen access