Supermassive Primordial Black Holes for Nano-Hertz Gravitational Waves and High-redshift JWST Galaxies

Recently,observational hints for supermassive black holes have been accumulating,prompting the question:Can primordial black holes(PBHs)be supermassive,particularly with masses M■10^(9)M_(⊙)?A supercritical bubble,containing an inflating baby universe,that nucleated during inflation can evolve into a PBH in our observable universe.We find that when the inflaton slowly transitions past a neighboring vacuum,the nucleation rate of supercritical bubbles inevitably peaks,leading to a mass distribution of multiverse PBHs with a peak mass up to M■10^(11)M_(⊙).Thus,our mechanism naturally provides a primordial origin for supermassive black holes.

AGN Lifetimes in UV-selected Galaxies: A Clue to Supermassive Black Hole-galaxy Coevolution

The coevolution between supermassive black holes(SMBHs) and their host galaxies has been proposed for more than a decade,albeit with little direct evidence about black hole accretion activities regulating galaxy star formation at z> 1.In this paper,we study the lifetimes of X-ray active galactic nuclei(AGNs) in UV-selected red sequence(RS),blue cloud(BC) and green valley(GV) galaxies,finding that AGN accretion activities are most prominent in GV galaxies at z ~1.5-2,compared with RS and BC galaxies.We also compare AGN accretion timescales with typical color transition timescales of UV-selected galaxies.We find that the lifetime of GV galaxies at z~1.5-2 is very close to the typical timescale when the AGNs residing in them stay in the high-accretion-rate mode at these redshifts;for BC galaxies,the consistency between the color transition timescale and the black hole strong accretion lifetime is more likely to happen at lower redshifts(z <1).Our results support the scenario where AGN accretion activities govern UV color transitions of host galaxies,making galaxies and their central SMBHs coevolve with each other.

Do Radio-loud Active Galactic Nuclei really follow the same M_(BH)-σ_* Relation as Normal Galaxies?

In an examination of the relationship between the black hole mass MBH and stellar velocity dispersion σ＊ in radio-loud active galactic nuclei （AGNs）, we studied two effects which may cause uncertainties in the black hole mass estimates of radio-loud AGNs： the relativistic beaming effect on the observed optical continuum radiation and the orientation effect on the broad emission line width. After correcting these two effects, we re-examined the MBH- σ[OIII] relation for a sample of radio-loud and radio-quiet AGNs, and found the relation for radio-loud AGNs still deviated from that for nearby normal galaxies and radio-quiet AGNs. We also found there is no significant correlation between radio jet power and narrow [OIII] line width, indicating absence of strong interaction between radio jet and narrow line region. It may be that the deviation of the MBH-σ＊ relation of radio-loud AGNs is intrinsic, or that the [OIII] line width is not a good indicator of σ＊ for radio-loud AGNs.

The AGN Feedback in Compact Galaxies:On the Impact of a More Massive Central Black Hole

We conduct high-resolution hydrodynamical simulations using the MACER framework to investigate the interplay between the interstellar medium,active galactic nuclei(AGN)feedback and black hole(BH)feeding in a massive compact galaxy,with an emphasis on the impact of different central BH masses.We find that with a more massive central BH,high-speed outflows are more prominent,and the gas fraction in the compact galaxy is reduced.Due to the lower gas density and higher gas temperature,the compact galaxy with a more massive BH(MAS galaxy)remains predominantly single-phase with the cooling time t_(cool)■100t_(ff).In contrast,the compact galaxy with the reference BH mass(REF galaxy)maintains a higher gas fraction with a shorter cooling time,slightly more multiphase gas and less prominent outflows.We further demonstrate that the difference in gas thermal states and kinematics is caused by the stronger AGN feedback in the compact galaxy with a more massive BH,where the AGN wind power is twice as much as that with the reference BH.Since the AGN feedback efficiently suppresses the inflow rate and the BH feeding rate,the BH mass growth is significant in neither the compact galaxy with the reference BH nor that with the more massive BH,only by 24%and 11%of the initial BH mass,respectively,over the entire evolution time of 10 Gyr.We thus posit that without ex situ mass supply from mergers,the massive BHs in compact galaxies cannot grow significantly via gas accretion during the late phase,but might have already formed by the end of the rapid early phase of galaxy formation.

Observational signatures of close binaries of supermassive black holes in active galactic nuclei

Inspired by the General Relativity for many decades,experimental physicists and astronomers have a solid dream to detect gravitational waves(GWs)from mergers of black holes,which came true until the excellent performance of the Laser Interferometer Gravitational-Wave Observatory(LIGO)at hundreds Hz.Nano-Hz GWs are expected to be radiated by close-binaries of supermassive black holes(CB-SMBHs;defined as those with separations less than^0.1 pc)formed during galaxy mergers and detected through the Pulsar Timing Array(PTA)technique.As of the writing,there remains no nano-Hz GWs detection.Searching for CB-SMBHs is also observationally elusive though there exist a number of possible candidates.In this review,we focus on observational signatures of CB-SMBHs from theoretic expectations,simulations and observations.These signatures appear in energy distributions of multiwavelength continuum,long term variations of continuum,jet morphology,reverberation delay maps and spectroastrometry of broad emission lines,AGN type transitions between type-1 and type-2(changing-look),and gaseous dynamics of circumbinary disks,etc.Unlike hundred-Hz GWs from stellar mass black hole binaries,the waveform chirping of nano-Hz GWs is too slow to detect in a reasonable human timescale.We have to resort to electromagnetic observations to measure orbital parameters of CB-SMBHs to test nano-Hz GW properties.Reverberation mapping is a powerful tool for probing kinematics and geometry of ionized gas in the gravitational well of SMBHs(single or binary)and therefore provides a potential way to determine orbital parameters of CB-SMBHs.In particular,a combination of reverberation mapping with spectroastrometry(realized at the Very Large Telescope Interferometer)will further reinforce this capability.The Atacama Large Millimeter/submillimeter Array(ALMA)and the forthcoming Square Kilometre Array(SKA)are suggested to reveal dynamics of circumbinary disks through molecular emission lines.

Constraining the Earth-mass Primordial Black Hole Mergers Model of the Non-repeating FRBs Using the First CHIME/FRB Catalog

In this paper,we upgrade the constraints for the Earth-mass primordial black hole mergers model based on the first Canadian Hydrogen Intensity Mapping Experiment(CHIME)/fast radio burst(FRB)catalog.Assuming the null hypothesis that the observed non-repeating FRBs originate from Earth-mass primordial black hole mergers,we find that how the charges were distributed in the primordial black hole population is well described by a double powerlaw function with typical charge value of law function with typical charge value of q_(c)/10^(-5)=1.60_(-0.28)^(+0.28),where the power-law index α_(1)=2.33_(-0.18)^(+0.15) for q<q_(c) and α_(2)=4.56_(-0.26)^(+0.30)for q≥q_(c).Here,q represents the charge of the black hole in units of√GM,where M is the mass of the black hole.Furthermore,we infer the local event rate of the bursts is 8.8_(-2.1)^(+5.7)×10^(4)Gpc^(-3) yr^(-1),which indicates that an abundance of the primordial black hole population f■10^(-4) is needed to account for the observed FRBs by CHIME.The results of this paper lay the basis for further research on the electromagnetic radiation background generated by the merger of primordial black hole mergers.

The Origin of the Flat Rotation Curves in Spiral Galaxies: The Hidden Roles of Glitching SMDEOs and Emission of Gravitational Waves

Supermassive DEOs (SMDEOs) are cosmologically evolved objects made of irreducible incompressible supranuclear dense superfluids: The state we consider to govern the matter inside the cores of massive neutron stars. These cores are practically trapped in false vacua, rendering their detection by outside observers impossible. Based on massive parallel computations and theoretical investigations, we show that SMDEOs at the centres of spiral galaxies that are surrounded by massive rotating torii of normal matter may serve as powerful sources for gravitational waves carrying away roughly 1042 erg/s. Due to the extensive cooling by GWs, the SMDEO-Torus systems undergo glitching, through which both rotational and gravitational energies are abruptly ejected into the ambient media, during which the topologies of the embedding spacetimes change from curved into flatter ones, thereby triggering a burst gravitational energy of order 1059 erg. Also, the effects of glitches found to alter the force balance of objects in the Lagrangian-L1 region between the central SMDEO-Torus system and the bulge, enforcing the enclosed objects to develop violent motions, that may explain the origin of the rotational curve irregularities observed in the innermost part of spiral galaxies. Our study shows that the generated GWs at the centres of galaxies, which traverse billions of objects during their outward propagations throughout the entire galaxy, lose energy due to repeatedly squeezing and stretching the objects. Here, we find that these interactions may serve as damping processes that give rise to the formation of collective forces f∝m(r)/r, that point outward, endowing the objects with the observed flat rotation curves. Our approach predicts a correlation between the baryonic mass and the rotation velocities in galaxies, which is in line with the Tully-Fisher relation. The here-presented self-consistent approach explains nicely the observed rotation curves without invoking dark matter or modifying Newtonian gravitation in the low-field approximation.

The parameters of binary black hole system in PKS 1510-089

Observations of PKS 1510-089 indicate the existence of a deep flux minimum with a timescale of -35 min and an interval of about 336 ± 14 d. A binary black hole system is proposed to be at the nucleus of this object. The secondary black hole orbits around the primary black hole. The minimum is caused by the periodic eclipse of the primary black hole by the secondary black hole. Based on the observations of PKS 1510-089, we estimate the parameters of the binary black hole system. The masses for the primary and secondary black holes are 1.37 × 10^9M⊙（M⊙ is the solar mass） and 1.37 × 10^7M⊙, and the major axis for this pair being about 0.1 parsec（pc）.

The debeamed luminosity, sychrotron peak frequency and black hole mass of BL Lac objects

We estimate the intrinsic luminosities and synchrotron peak frequencies using the derived Doppler factor for a sample of 170 BL Lac objects, of which the synchrotron peak frequency is derived by fitting the SED constructed with the collected multi-band data from the literature. We find that the debeamed radio and optical core luminosities follow the same correlation found for FR I radio galaxies, which is in support of the unification of the BL Lac objects and the FR I galaxies based on orientation. For the debeamed luminosity at the synchrotron peak frequency, we find a significant positive correlation between the luminosity and intrinsic synchrotron peak frequency. This implies that the more powerful sources may have the majority of jet emission at higher frequency. At the synchrotron peak frequency, the intrinsic luminosity and black hole mass show strong positive correlation, while mild correlation is found in the case of jet power, indicating that the more powerful sources may have more massive black holes.

The relation between black hole masses and Lorentz factors of the jet components in blazars

We explore the relationship between black hole mass （MBH） and the motion of the jet components for a sample of blazars. The Very Long Baseline Array （VLBA） 2 cm Survey and its continuation： Monitoring of Jets in active galactic nuclei （AGNs） with VLBA Experiments （MOJAVE） have observed 278 radio-loud AGNs, of which 146 blazars have reliable measurements of their apparent velocities of the jet components. We calculate the minimal Lorentz factors for these sources from their measured apparent velocities, and their black hole masses are estimated with their broad-line widths. A significant intrinsic correlation is found between black hole masses and the minimal Lorentz factors of the jet components. The Eddington ratio is only weakly correlated with the minimal Lorentz factor, which may imply that the Blandford-Znajek （BZ） mechanism may dominate over the Blandford-Payne （BP） mechanism for the jet acceleration （at least） in blazars.

The Binary Black Hole Scenario for the BL Lacertae Object AO 0235+16

Recent analysis of the long term radio light curve of the extremely variableBL Lacertae object AO 0235+16 by Raiteri et al. have revealed the presence of recurrent outburstswith a period of ～ 5.7 +- 0.5 yr. Periodicity analysis of the optical light curve also showsevidence for a shorter period. Here we discuss whether such a behavior can be explained by a binaryblack hole model where the accretion disk of one of the supermassive black holes is precessing dueto the tidal effects of the companion. We estimate the mass of the accreting hole and analyzeconstraints on the secondary mass and the orbital parameters of the system. It is possible toprovide a viable interpretation of the available multiwavelength data.

Correlation between excitation index and Eddington ratio in radio galaxies

We use a sample of 111 radio galaxies with redshift z 〈 0.3 to investigate their nuclear properties. The black hole masses of the sources in this sample are estimated with the velocity dispersion/luminosity of the galaxies, or the width of the broad-lines. We find that the excitation index, the relative intensity of low and high excitation lines, is correlated with the Eddington ratio for this sample. The size of the narrow-line region （NLR） was found to vary with ionizing luminosity as RNLR ∝ Lion^0.25 （Liu et al. 2013）. Using this empirical relation, we find that the correlation between the excitation index and the Eddington ratio can be reproduced by photoionization models. We adopt two sets of spectral energy distributions （SEDs）, with or without a big blue bump in ultraviolet as the ionizing continuum, and infer that the modeled correlation between the excitation index and the Eddington ratio is insensitive to the applied SED. This means that the difference between high excitation galaxies and low excitation galaxies is not caused by the different accretion modes in these sources. Instead, it may be caused by the size of the NLR.

Radio Luminosity,Black Hole Mass and Eddington Ratio for Quasars from the Sloan Digital Sky Survey

We investigate the MBH-σ＊ relation for radio-loud quasars with redshifl z 〈 0.83 in Data Release 3 of the Sloan Digital Sky Survey （SDSS）. The sample consists of 3772 quasars with better models of the H/4 and [O Ⅲ] lines and available radio luminosity, including 306 radio-loud quasars, 3466 radio-quiet quasars with measured radio luminosity or upper-limit of radio luminosity （181 radio-quiet quasars with measured radio luminosity）. The virial supermassive black hole mass （MBH） is calculated from the broad Hβline, and the host stellar velocity dispersion （σ＊） is traced by the core [O Ⅲ] gaseous velocity dispersion. The radio luminosity and radio loudness are derived from the FIRST catalog. Our results are as follows： （1） For radio-quiet quasars, we confirm that there is no obvious deviation from the MBH-σ＊ relation defined for inactive galaxies when the uncertainties in ~IBH and the luminosity bias are concerned. （2） We find that the radio-loud quasars deviate more from the MBH-σ＊ relation than do the radio-quiet quasars. This deviation is only partly due to a possible cosmological evolution of the MBH-σ＊ relation and the luminosity bias. （3） The radio luminosity is proportional to MBH1.28＋0.23-0.16（LBol/LEdd） ^1.29＋0.31-0.24 for radio-quiet quasars and to -MBH3.10＋0.60-0.70（LBol/LEdd）^4.18＋1.40-1.10 - for radio-loud quasars. The weaker dependence of the radio luminosity on the mass and the Eddington ratio for radio-loud quasars shows that other physical effects would account for their radio luminosities, such as the spin of the black hole.

Does black hole spin play a key role in the FSRQ/BL Lac dichotomy?

Blazars are characterized by large intensity and spectral variations across the electromagnetic spectrum It is believed that jets emerging from them are almost aligned with the line-of-sight. The major- ity of identified extragalactic sources in γ-ray catalogs of EGRET and Fermi are blazars. Observationally, blazars can be divided into two classes： fiat spectrum radio quasars （FSRQs） and BL Lacs. BL Lacs usually exhibit lower γ-ray luminosity and harder power law spectra at γ-ray energies than FSRQs. We attempt to explain the high energy properties of FSRQs and BL Lacs from Fermi γ-ray space telescope observations. It was argued previously that the difference in accretion rates is mainly responsible for the large mismatch in observed luminosity in ＂7-ray. However, when intrinsic luminosities are derived by correcting for beaming effects, this difference in 7-ray luminosity between the two classes is significantly reduced. In order to ex- plain this difference in intrinsic luminosities, we propose that spin plays an important role in the luminosity distribution dichotomy of BL Lacs and FSRQs. As the outflow power of a blazar increases with increasing spin of a central black hole, we suggest that the spin plays a crucial role in making BL Lac sources low luminous and slow rotators compared to FSRQ sources.

On the Correlation between Radio Properties and Black Hole Mass of Quasars

The question whether the radio properties of quasars are related to the massof the central black hole or the accretion rate is important for our understanding of the formationof relativistic jets, but no consensus has been reached from statistical analyses. Using two largequasar samples, one radio-selected, one optical-selected, we re-examined these relations and findthat previous differences between radio- and optical- selected samples can be ascribed, at leastpartly, to the effect of the narrow line component. All previous claimed correlations are muchweaker, if exist at all.

A Multi-parameter Model for Radio Dichotomy of Active Galactic Nuclei and Jets

Based on the coexistence of the Blandford-Znajek and magnetic coupling processes in black hole （BH） accretion disc, a multi-parameter model for jet powers and radio loudnees of active galactic nuclei （AGNs） is studied. It turns out that radio-loudnees of AGNs couM be governed by five parameters： （i） the BH spin, （ii） a power-law index of the variation of the magnetic field on the disc; （iii） a parameter determining the position of the inner edge of the disc, （iv） the ratio of the pressure of the magnetic field on the horizon to the ram pressure of the innermost parts of an accretion flow, and （v） the ratio of the angular velocity of the open field lines to that of the horizon. The observed dichotomy between radio-loud and radio-quiet AGNs is well interpreted by the effects of the above parameters. Furthermore, we discuss the derivative of radio loudness of AGNs with respect to each parameter separately. In addition, the effect of the screw instability on radio loudness of AGNs is discussed.

Possibility of measuring spin precession of the nearest supermassive black hole by using S stars

The supermassive black hole （SMBH） with a mass of 4 million M⊙ inside the radio source Sgr A＊ in our Galactic center is the nearest SMBH. Once S stars with a shorter period are observed, relativistic precessions especially the Lense-Thirring effect can be measured by astronomical observations at the 10 ~tas level in the future. An interesting but so far unaddressed problem is that the SMBH not only has spin but also spin precession like similar objects. We study the effect of such spin precession on the orbital precessions of orbiting stars. Our results show that the spin precession can produce a periodic oscillation in the precession of the star＇s orbital plane, but has no obvious effect on the periapse shift. For stars with an orbital period of O（0.1） yr or less, such visible oscillations occur when the SMBH＇s spin-precession period ranges from about a few tens of years to hundreds of years. The period of oscillation is the same as the one of the spin precession. In principle, the precession of this oscillating orbital plane can be observed and then the spin and spin precession of the nearest SMBH can be determined.

A 2.5-dimensional viscous, resistive, advective magnetized accretion-outflow coupling in black hole systems: a higher order polynomial approximation

The correlated and coupled dynamics of accretion and outflow around black holes （BHs） are essentially governed by the fundamental laws of conservation as outflow extracts matter, momentum and energy from the accretion region. Here we analyze a robust form of 2.5-dimensional viscous, resistive, advective magnetized accretion-outflow coupling in BH systems. We solve the complete set of coupled MHD conservation equations self-consistently, through invoking a generalized polynomial expansion in two dimensions. We perform a critical analysis of the accretion-outflow region and provide a complete quasi-analytical family of solutions for advective flows. We obtain the physically plausible outflow solu- tions at high turbulent viscosity parameter a（〉0.3）, and at a reduced scale-height, as magnetic stresses compress or squeeze the flow region. We found that the value of the large-scale poloidal magnetic field Bp is enhanced with the increase of the geometrical thickness of the accretion flow. On the other hand, differential magnetic torque （-r2BBz） increases with the increase in M. Bp, -r2BBz as well as the plasma beta/3p get strongly augmented with the increase in the value of a, enhancing the transport of vertical flux outwards. Our solutions indicate that magnetocentrifugal acceleration plausibly plays a dominant role in effusing out plasma from the radial accretion flow in a moderately advective paradigm which is more centrifugally dominated. However in a strongly advective paradigm it is likely that the thermal pressure gradient would play a more contributory role in the vertical transport of plasma.

Do Very Short Gamma Ray Bursts Originate from Primordial Black Holes? Review

We present the state of current research of Very Short Gamma Ray Bursts (VSGRBs) from seven GRB detectors. We found that VSGRBs form distinct class of GRBs, which in our opinion, in most cases can originate from the evaporating Primordial Black Holes (PBHs). Arguments supporting our opinion: 1) GRBs with time duration (T90) < 100 ms form distinct class: VSGRBs. 2) We observe significant anisotropy in the galactic angular distribution of BATSE VSGRB events. 3) V/Vmax distribution for BATSE VSGRB events indicates the local distance production. 4) VSGBBs have more energetic γ-ray burst than other GRBs with longer duration (KONUS). 5) We observe small number of afterglows in SWIFT VSGRB sample (25%), in contrast with the noticeable afterglow frequency in SGRB sample (78%). 6) Time profile of rising part BATSE VSGRBs is in agreement with the evaporation PBH model.

Evidence of the Link between Broad Emission Line Regions and Accretion Disks in Active Galactic Nuclei

There is observational evidence that broad-line regions （BLRs） exist in most active galactic nuclei （AGNs）, but their origin is still unclear. One scenario is that the BLRs originate from winds accelerated from the hot coronae of the disks, and the winds are suppressed when the black hole is accreting at low rates. This model predicts a relation between rh （m ~ /~f//~fEdd） and the FWHM of broad emission lines. We estimate the central black hole masses for a sample of bright AGNs by using their broad Hβ line-widths and optical luminosities. The dimensionless accretion rates rh = M/MEdd are derived from the optical continuum luminosities by using two different models： using an empirical relation between the bolometric luminosity Lbol and the optical luminosity （rh = Lbol/LEdd, a fixed radiative efficiency is adopted）; and calculating the optical spectra of accretion disks as a function of rh. We find a slgmficant correlation between the denved of and the observed line width of Hβ, FWHMα m^-0.aT, which almost overlaps the disk-corona model calculations, if the viscosity a α≈ 0.1 - 0.2 is adopted. Our results provide strong evidence for the physical link between the BLRs and accretion disks in AGNs.