Realistic Accretion Research Articles

Simulating X-Ray Reverberation in the Ultraviolet-emitting Regions of Active Galactic Nuclei Accretion Disks with Three-dimensional Multifrequency Radiation Magnetohydrodynamic Simulations

Active galactic nuclei (AGN) light curves observed with different wave bands show that the variability in longer wavelength bands lags the variability in shorter wavelength bands. Measuring these lags, or reverberation mapping, is used to measure the radial temperature profile and extent of AGN disks, typically with a reprocessing model that assumes X-rays are the main driver of the variability in other wavelength bands. To demonstrate how this reprocessing works with realistic accretion disk structures, we use 3D local shearing box multifrequency radiation magnetohydrodynamic simulations to model the UV-emitting region of an AGN disk, which is unstable to the magnetorotational instability and convection. At the same time, we inject hard X-rays (>1 keV) into the simulation box to study the effects of X-ray irradiation on the local properties of the turbulence and the resulting variability of the emitted UV light curve. We find that disk turbulence is sufficient to drive intrinsic variability in emitted UV light curves and that a damped random walk model is a good fit to this UV light curve for timescales >5 days. Meanwhile, X-ray irradiation has negligible impact on the power spectrum of the emitted UV light curve. Furthermore, the injected X-ray and emitted UV light curves are only correlated if there is X-ray variability on timescales >1 day, in which case we find a correlation coefficient r = 0.34. These results suggest that if the opacity for hard X-rays is scattering dominated as in the standard disk model, hard X-rays are not the main driver of reverberation signals.

Effects of accretion on the structure and rotation of forming stars

Context. Rotation period measurements of low-mass stars show that the spin distributions in young clusters do not exhibit the spin-up expected due to contraction in the phase when a large fraction of stars is still surrounded by accretion discs. Many physical models have been developed to explain this feature based on different types of star-disc interactions alone. In this phase, the stars accrete mass and angular momentum and may experience accretion-enhanced magnetised winds. The stellar structure and angular momentum content thus strongly depend on the properties of the accretion mechanism. At the same time, the accretion of mass and energy has a significant impact on the evolution of the stellar structure and the moment of inertia. Our understanding of the spin rates of young stars therefore requires a description of how accretion affects the stellar structure and angular momentum simultaneously. Aims. We aim to understand the role of accretion to explain the observed rotation-rate distributions of forming stars. Methods. We computed evolution models of accreting very young stars and determined in a self-consistent way the effect of accretion on stellar structure and the angular momentum exchanges between the stars and their disc. We then varied the deuterium content, the accretion history, the entropy content of the accreted material, and the magnetic field as well as the efficiency of the accretion-enhanced winds. Results. The models are driven alternatively both by the evolution of the momentum of inertia and by the star-disc interaction torques. Of all the parameters we tested, the magnetic field strength, the accretion history, and the deuterium content have the largest impact. The injection of heat plays a major role only early in the evolution. Conclusions. This work demonstrates the importance of the moment of inertia evolution under the influence of accretion for explaining the constant rotation-rate distributions that are observed during the star-disc interactions. When we account for rotation, the models computed with the recently calculated torque along with a consistent structural evolution of the accreting star are able to explain the almost constant spin evolution for the whole range of parameters we investigated, but it only reproduces a narrow range around the median of the observed spin rate distributions. Further development, including for example more realistic accretion histories based on dedicated disc simulations, are likely needed to reproduce the extremes of the spin rate distributions.

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

Relaxing CMB bounds on primordial black holes: The role of ionization fronts

The accretion of matter onto primordial black holes (PBHs) during the dark ages and the subsequent energy injection in the medium should have left imprints on the cosmic microwave background (CMB) anisotropies. Recent works have claimed stringent CMB limits on the PBH abundance, hardly compatible with a PBH interpretation of the gravitational-wave observations of binary BH mergers. By using a more realistic accretion model based on hydrodynamical simulations and conservative assumptions for the emission efficiency, we show that CMB limits on the PBH abundance are up to two orders of magnitude less stringent than previously estimated between $10$ and $10^4$ M$_\odot$. This reopens the possibility that PBHs might explain at the same time (at least a fraction of) the dark matter, some of the LIGO-Virgo-KAGRA binary BH mergers and the existence of super-massive BHs. More generally, we emphasize that PBH accretion can be a rather complex physical process with velocity dependences that are hard to assess, which introduces large uncertainties in accretion-based limits on the PBH abundance.

First minutes of a binary-driven hypernova

We simulate the first minutes of the evolution of a binary-driven hypernova (BdHN) event, with a special focus on the associated accretion processes of supernova (SN) ejecta onto the newborn neutron star ($\nu$NS) and the NS companion. We calculate the rotational evolution of the $\nu$NS and the NS under the torques exerted by the accreted matter and the magnetic field. We take into account general relativistic effects and use realistic hypercritical accretion rates obtained from three-dimensional smoothed-particle-hydrodynamics (SPH) numerical simulations of the BdHN for a variety of orbital periods. We show that the rotation power of the $\nu$NS has a unique double-peak structure while that of the NS has a single peak. These peaks are of comparable intensity and can occur very close in time or even simultaneously depending on the orbital period and the initial angular momentum of the stars. We outline the consequences of the above features in the early emission and their consequent observation in long gamma-ray bursts (GRBs).

The imprint of star formation on stellar pulsations

In the earliest phases of their evolution, stars gain mass through the acquisition of matter from their birth clouds. The widely accepted classical concept of early stellar evolution neglects the details of this accretion phase and assumes the formation of stars with large initial radii that contract gravitationally. In this picture, the common idea is that once the stars begin their fusion processes, they have forgotten their past. By analysing stellar oscillations in recently born stars, we show that the accretion history leaves a potentially detectable imprint on the stars’ interior structures. Currently available data from space would allow discriminating between these more realistic accretion scenarios and the classical early stellar evolution models. This opens a window to investigate the interior structures of young pulsating stars that will also be of relevance for related fields, such as stellar oscillations in general and exoplanet studies.

Grumpy: a simple framework for realistic forward modelling of dwarf galaxies

ABSTRACT We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass–halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in N-body simulations, such as colours, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colours and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude–halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultra-faint galaxies is sensitive to reionization redshift (zrei) and is predicted to be consistent with observations for zrei ≲ 9.

Jets in common envelopes: a low-mass main-sequence star in a red giant

ABSTRACTWe present small-scale 3D hydrodynamical simulations of the evolution of a 0.3 M⊙ main-sequence (MS) star that launches two perpendicular jets within the envelope of a 0.88 M⊙ red giant (RG). Based on previous large-scale simulations, we study the dynamics of the jets either when the secondary star is grazing, when it has plunged-in, or when it is well within the envelope of the RG (in each stage for ∼11 d). The dynamics of the jets through the common envelope (CE) depend on the conditions of the environment as well as on their powering. In the grazing stage and the commencement of the plunge self-regulated jets need higher efficiencies to break out of the envelope of the RG. Deep inside the CE, on the time-scales simulated, jets are choked independently of whether they are self-regulated or constantly powered. Jets able to break out of the envelope of the RG in large-scale simulations, are choked in our small-scale simulations. The accreted angular momentum on to the secondary star is not large enough to form a disc. The mass accretion on to the MS star is 1–10 per cent of the Bondi–Hoyle–Littleton rate (∼10−3–10−1 M⊙ yr−1). High-luminosity emission, from X-rays to ultraviolet and optical, is expected if the jets break out of the CE. Our simulations illustrate the need for inclusion of more realistic accretion and jet models in the dynamical evolution of the CEs.

Merging Criteria for Planetesimal Collisions

Abstract In the standard scenario of planet formation, terrestrial planets, ice giants, and cores of gas giants are formed by the accumulation of planetesimals. However, there are few N-body simulation studies of planetesimal accretion that correctly take into account the merging condition of planetesimals. In order to investigate a realistic accretion process of planetesimals, it is necessary to clarify the merging criteria of planetesimals at collision. We perform numerical collision experiments using smoothed particle hydrodynamics and obtain the merging criteria as a function of planetesimal mass and impact parameters for undifferentiated rocky and icy planetesimals and differentiated icy planetesimals. We vary the total mass of colliding planetesimals, their mass ratios, and the impact angle and obtain the critical impact velocity as the merging criteria distinguishing merging from hit-and-run collision. We find that the critical impact velocity normalized by the two-body surface escape velocity decreases with increasing impact angle. The critical impact velocity does not depend on the total mass, while it has a weak positive dependence on the mass ratio. These results barely depend on the composition and internal structure of the planetesimals.

Accretion of a Vlasov gas onto a black hole from a sphere of finite radius and the role of angular momentum

The accretion of a spherically symmetric, collisionless kinetic gas cloud on to a Schwarzschild black hole is analysed. Whereas previous studies have treated this problem by specifying boundary conditions at infinity, here the properties of the gas are given at a sphere of finite radius. The corresponding steady-state solutions are computed using four different models with an increasing level of sophistication, starting with the purely radial infall of Newtonian particles and culminating with a fully general relativistic calculation in which individual particles have angular momentum. The resulting mass accretion rates are analysed and compared with previous models, including the standard Bondi model for a hydrodynamic flow. We apply our models to the supermassive black holes Sgr A* and M87*, and we discuss how their low luminosity could be partially explained by a kinetic description involving angular momentum. Furthermore, we get results consistent with previous model-dependent bounds for the accretion rate imposed by rotation measures of the polarised light coming from Sgr A* and with estimations of the accretion rate of M87* from the Event Horizon Telescope collaboration. Our methods and results could serve as a first approximation for more realistic black hole accretion models in various astrophysical scenarios in which the accreted material is expected to be nearly collisionless.

Dual Component Passive Icephobic Coatings with Micron-Scale Phase-Separated 3D Structures.

A passive icephobic coating (τice < 20 kPa) is an enabling technology to many industries, including aerospace and energy and power generation, with recent efforts in materials research identifying strategies to achieve this low adhesion threshold. To better meet this need, we have combined low surface energy perfluoropolyether (PFPE) and hydrophilic poly(ethylene glycol) (PEG) species in a segmented polyurethane thermoplastic elastomer. Coating microstructure presents a segregated 3D morphology at the micron-scale (1-100 μm) with discrete PFPE and continuous PEG phases self-similar through the thickness. Spray application produces a solid, mechanically tough film free of additive fluids or sacrificial elements, demonstrating exceptional ice adhesion reduction up to 1000× lower versus aluminum (τice < 1 kPa), as measured under environmentally realistic accretion and centrifugal test shedding conditions. Finally, the modular nature of the synthetic system allows PEG and PFPE to be exchanged for poly(tetramethylene oxide) to investigate performance drivers.

Neutrino scattering off a black hole surrounded by a magnetized accretion disk

We study the neutrino scattering off a rotating black hole with a realistic accretion disk permeated by an intrinsic magnetic field. Neutrino trajectories in curved spacetime as well as the particle spin evolution in dense matter of an accretion disk and in the magnetic field are accounted for exactly. We obtain the fluxes of outgoing ultrarelativistic neutrinos taking into account the change of the neutrino polarization owing to spin oscillations. Using the conservative value of the neutrino magnetic moment and realistic radial distributions of the matter density and the magnetic field strength, we get that these fluxes are reduced by several percent compared to the case when no spin oscillations are accounted for. In some situations, there are spikes in the neutrino fluxes because of the neutrino interaction with the rotating plasma of an accretion disk. Taking into account the uncertainties in the astrophysical neutrino fluxes, the predicted effects turn out to be quite small to be observed with the current neutrino telescopes.

The Formation of Supermassive Black Holes at Early Universe

The high-redshift universe is now the frontier of modern astronomy, and is the key object for current and near-future telescopes. The high-redshift quasars put a tight constraint on existing BH growth model, as it is challenging to form a 1010 M⊙ BH at z∼6.In this work, we relax the widely-adopted Salpeter BH growth model, to consider a more realistic path. We consider the variation in the mass and angular momentum orientation of gas supply (i.e. through ∊L and n), as well as the change in accretion mode (hot versus cold). Moreover, the conventionally considered BH spin impact on the radiative efficiency ∊M is also taken into account. Our key results can be summarized as follows. Firstly, sufficient gas supply (larger in ∊L) is obviously a key factor to make the BH grow efficiently. Also, The BH spin (a*), through the radiative efficiency ∊M, has a dominant impact on the BH growth, i.e., those high-redshift quasars with MBH > 1010 M⊙ should be formed in a chaotic gas supply situation, where the angular momentum orientation of the gas is random. Finally, through analyzing the most realistic accretion model, we find that the only existing model of seed BH is through the direct collapse, driven by either dynamical processes or thermodynamics

Disc fragmentation and intermittent accretion on to supermassive stars

ABSTRACT Supermassive stars (SMSs) with ∼104–105 M⊙ are candidate objects for the origin of supermassive black holes observed at redshift z &amp;gt; 6. They are supposed to form in primordial-gas clouds that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and its quiescent periods are longer than the Kelvin–Helmholtz (KH) time-scales at the stellar surfaces. In this paper, we examine the role of the ionizing radiation feedback based on the accretion history in two possible SMS-forming clouds extracted from cosmological simulations, following their evolution with vertically integrated two-dimensional hydrodynamic simulations with detailed thermal and chemical models. The consistent treatment of the gas thermal evolution is crucial for obtaining the realistic accretion history, as we demonstrate by performing an additional run with a barotropic equation of state, in which the fluctuation of the accretion rate is artificially suppressed. We find that although the accretion becomes intermittent due to the formation of spiral arms and clumps in gravitationally unstable discs, the quiescent periods are always shorter than the KH time-scales, implying that SMSs can form without affected by the ionizing radiation.

Numerical investigation of plasma-driven superradiant instabilities

Photons propagating in a plasma acquire an effective mass μ, which is given by the plasma frequency and which scales with the square root of the plasma density. As noted previously in the literature, for electron number densities ne ∼ 10−3 cm−3 (such as those measured in the interstellar medium) the effective mass induced by the plasma is μ ∼ 10−12 eV. This would cause superradiant instabilities for spinning black holes of a few tens of solar masses. An obvious problem with this picture is that densities in the vicinity of black holes are much higher than in the interstellar medium because of accretion, and possibly also pair production. We have conducted numerical simulations of the superradiant instability in spinning black holes surrounded by a plasma with density increasing closer to the black hole, in order to mimic the effect of accretion. While we confirm that superradiant instabilities appear for plasma densities that are sufficiently low near the black hole, we find that astrophysically realistic accretion disks are unlikely to trigger such instabilities.

Spin oscillations of neutrinos scattered off a rotating black hole

Spin oscillations of neutrinos, gravitationally scattered off a black hole (BH), are studied. The cases of nonrotating and rotating BHs are analyzed. We derive the analytic expressions for the transition and survival probabilities of spin oscillations when neutrinos interact with these gravitational backgrounds. The obtained transition probabilities depend on the impact parameter, as well as the neutrino energy and the particle mass. We find that there is a possibility of spin oscillations of ultrarelativistic neutrinos scattering off a rotating BH. Then, considering the neutrino scattering off BH surrounded by background matter, we derive the effective Schrödinger equation for spin oscillations. The numerical solution of this equation is obtained in the case of a supermassive BH with a realistic accretion disk. Spin effects turn out to be negligible in the neutrino scattering in the Schwarzschild metric. In the Kerr metric, we find that the observed neutrino fluxes can be reduced almost 10% because of spin oscillations when ultrarelativistic neutrinos experience gravitational scattering. The neutrino interaction with an accretion disk results in the additional modification of the intensities of outgoing neutrino fluxes. We consider the applications of the obtained results for the neutrino astronomy.

Spin effects in neutrino gravitational scattering

We study spin oscillations of neutrinos gravitationally scattered off a nonrotating black hole (BH). We derive the transition and survival probabilities of spin oscillations in quadratures when neutrinos interacts with BH only. The dependence of the probabilities on the impact parameter is analyzed. Then, we obtain the effective Schr\"{o}dinger equation for neutrino spin oscillations in neutrino scattering off BH surrounded by background matter. This equation is solved numerically in the case of a supermassive BH with a realistic accretion disk. We find that the observed neutrino fluxes can be reduced almost 20% because of spin oscillations when neutrinos experience gravitational scattering. The neutrino interaction with an accretion disk results in the additional asymmetry in the intensities of outgoing fluxes depending on the neutrino trajectory.

X-ray reflection from cold white dwarfs in magnetic cataclysmic variables

We model X-ray reflection from white dwarfs (WD) in magnetic cataclysmic variables (mCVs) with a Monte Carlo simulation. A point source with a power-law spectrum or a realistic post-shock accretion column (PSAC) source irradiates a cool and spherical WD. The PSAC source emits thermal spectra of various temperatures stratified along the column according to the PSAC model. In the point source simulation, we confirm (1) a source harder and nearer to the WD enhances the reflection, (2) higher iron abundance enhances the equivalent widths (EWs) of fluorescent iron K${\alpha}_{1,2}$ lines and their Compton shoulder, and increases cut-off energy of a Compton hump, and (3) significant reflection appears from an area that is more than 90 deg apart from the position right under the point X-ray source because of the WD curvature. The PSAC simulation reveals that (1) a more massive WD basically enhances the intensities of the fluorescent iron K${\alpha}_{1,2}$ lines and the Compton hump, except for some specific accretion rate, because the more massive WD makes the hotter PSAC from which higher energy X-rays are preferentially emitted, (2) a larger specific accretion rate monotonically enhances the reflection because it makes the hotter and shorter PSAC, and (3) the intrinsic thermal component hardens by occultation of the cool base of the PSAC by the WD. We quantitatively evaluate influences of the parameters on the EWs and the Compton hump with both types of sources. We also calculate X-ray modulation profiles brought about by the WD spin. They depend on angles of the spin axis from the line of sight and from the PSAC, and whether the two PSACs can be seen. The reflection spectral model and the modulation model involve the fluorescent lines and the Compton hump and can directly be compared to the data which allow us to evaluate these geometrical parameters with unprecedented accuracy.

Simulations of the Fe Kα Energy Spectra from Gravitationally Microlensed Quasars

Abstract The analysis of the Chandra X-ray observations of the gravitationally lensed quasar RX J1131−1231 revealed the detection of multiple and energy-variable spectral peaks. The spectral variability is thought to result from the microlensing of the Fe Kα emission, selectively amplifying the emission from certain regions of the accretion disk with certain effective frequency shifts of the Fe Kα line emission. In this paper, we combine detailed simulations of the emission of Fe Kα photons from the accretion disk of a Kerr black hole with calculations of the effect of gravitational microlensing on the observed energy spectra. The simulations show that microlensing can indeed produce multiply peaked energy spectra. We explore the dependence of the spectral characteristics on black hole spin, accretion disk inclination, corona height, and microlensing amplification factor and show that the measurements can be used to constrain these parameters. We find that the range of observed spectral peak energies of QSO RX J1131−1231 can only be reproduced for black hole inclinations exceeding 70° and for lamppost corona heights of less than 30 gravitational radii above the black hole. We conclude by emphasizing the scientific potential of studies of the microlensed Fe Kα quasar emission and the need for more detailed modeling that explores how the results change for more realistic accretion disk and corona geometries and microlensing magnification patterns. A full analysis should furthermore model the signal-to-noise ratio of the observations and the resulting detection biases.

Retrograde binaries of massive black holes in circumbinary accretion discs

We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different disc's surface densities which alter the black hole's dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black hole's dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.