Formation of supermassive black hole seeds in nuclear star clusters via gas accretion and runaway collisions

Monthly Notices of the Royal Astronomical Society Oxford University Press 503:1 (2021) 1051-1069

Authors:

Arpan Das, Dominik RG Schleicher, Nathan WC Leigh, Tjarda CN Boekholt

Abstract:

More than 200 supermassive black holes (SMBHs) of masses ≳109M⊙≳109M⊙ have been discovered at z ≳ 6. One promising pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs) with masses ∼103−105M⊙∼103−105M⊙ into seed black holes which could grow upto few times 109M⊙109M⊙ SMBHs observed at z ∼ 7. In this paper, we explore how SMSs with masses ∼103−105M⊙∼103−105M⊙ could be formed via gas accretion and runaway stellar collisions in high-redshift, metal-poor nuclear star clusters (NSCs) using idealized N-body simulations. We explore physically motivated accretion scenarios, e.g. Bondi–Hoyle–Lyttleton accretion and Eddington accretion, as well as simplified scenarios such as constant accretions. While gas is present, the accretion time-scale remains considerably shorter than the time-scale for collisions with the most massive object (MMO). However, overall the time-scale for collisions between any two stars in the cluster can become comparable or shorter than the accretion time-scale, hence collisions still play a crucial role in determining the final mass of the SMSs. We find that the problem is highly sensitive to the initial conditions and our assumed recipe for the accretion, due to the highly chaotic nature of the problem. The key variables that determine the mass growth mechanism are the mass of the MMO and the gas reservoir that is available for the accretion. Depending on different conditions, SMSs of masses ∼103−105M⊙∼103−105M⊙ can form for all three accretion scenarios considered in this work.

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

Journal of Plasma Physics Cambridge University Press 87:1 (2021) 855870103

Authors:

Jm Garcia-Regana, M Barnes, I Calvo, Fi Parra, Ja Alcuson, R Davies, A Gonzalez-Jerez, A Mollen, E Sanchez, Jl Velasco, A Zocco

Abstract:

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. It is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

Multiple-scales approach to the averaging problem in cosmology

Journal of Cosmology and Astroparticle Physics IOP Publishing 2021:02 (2021) 049-049

On a new formulation for energy transfer between convection and fast tides with application to giant planets and solar type stars

Monthly Notices of the Royal Astronomical Society Royal Astronomical Society 503:4 (2021) 5789-5806

Abstract:

All the studies of the interaction between tides and a convective flow assume that the large scale tides can be described as a mean shear flow which is damped by small scale fluctuating convective eddies. The convective Reynolds stress is calculated using mixing length theory, accounting for a sharp suppression of dissipation when the turnover timescale is larger than the tidal period. This yields tidal dissipation rates several orders of magnitude too small to account for the circularization periods of late–type binaries or the tidal dissipation factor of giant planets. Here, we argue that the above description is inconsistent, because fluctuations and mean flow should be identified based on the timescale, not on the spatial scale, on which they vary. Therefore, the standard picture should be reversed, with the fluctuations being the tidal oscillations and the mean shear flow provided by the largest convective eddies. We assume that energy is locally transferred from the tides to the convective flow. Using this assumption, we obtain values for the tidal Q factor of Jupiter and Saturn and for the circularization periods of PMS binaries in good agreement with observations. The timescales obtained with the equilibrium tide approximation are however still 40 times too large to account for the circularization periods of late–type binaries. For these systems, shear in the tachocline or at the base of the convective zone may be the main cause of tidal dissipation.

Eccentric black hole mergers in active galactic nuclei

Astrophysical Journal Letters IOP Publishing 907:1 (2021) L20

Authors:

Hiromichi Tagawa, Bence Kocsis, Zoltan Haiman, Imre Bartos, Kazuyuki Omukai, Johan Samsing

Abstract:

The astrophysical origin of gravitational wave transients is a timely open question in the wake of discoveries by the Laser Interferometer Gravitational-Wave Observatory (LIGO)/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and the number of suspected hierarchical mergers observed by LIGO/Virgo. The binary eccentricity distribution can provide further information to distinguish between astrophysical models. Here we derive the eccentricity distribution of BH mergers in AGN disks. We find that eccentricity is mainly due to binary–single (BS) interactions, which lead to most BH mergers in AGN disks having a significant eccentricity at 0.01 Hz, detectable by the Laser Interferometer Space Antenna. If BS interactions occur in isotropic-3D directions, then 8%–30% of the mergers in AGN disks will have eccentricities at 10 Hz above e10 Hz ≳ 0.03, detectable by LIGO/Virgo/Kamioka Gravitational Wave Detector, while 5%–17% of mergers have e10 Hz ≥ 0.3. On the other hand, if BS interactions are confined to the AGN–disk plane due to torques from the disk, with 1–20 intermediate binary states during each interaction, or if BHs can migrate to ≲ 10−3 pc from the central supermassive BH, then 10%–70% of the mergers will be highly eccentric (e10 Hz ≥ 0.3), consistent with the possible high eccentricity in GW190521.