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

Ataru Tanikawa

Ataru Tanikawa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.GAastro-ph.COastro-ph.EPastro-ph.IM

Trust snapshot

Quick read

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

16 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 &#34;helium-ignited violent merger&#34; or &#34;dynamically-driven double-degenerate double-detonation&#34; (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&#39;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