The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constrained simulations in a volume of $(15\, h^{-1}{\rm Mpc})^3$ successfully recover the evolution of large-scale structure and the stellar and BH masses in the vicinity of a $\sim10^{12}\, M_{\odot}$ halo which we identified in BlueTides at $z\sim7$ hosting a $\sim10^9\, M_{\odot}$ SMBH. Among our constrained simulations, only the ones with a low-tidal field and high-density peak in the initial conditions induce the fastest BH growth required to explain the $z>6$ quasars. We run two sets of simulations with different BH seed masses of $5\times10^3$, $5\times10^4$, and $5\times10^5\, h^{-1}M_{\odot}$, (a) with the same ratio of halo to BH seed mass and (b) with the same halo threshold mass. At $z=6$, all the SMBHs converge in mass to $\sim10^9\, M_{\odot}$ except for the one with the smallest seed in (b) undergoing critical BH growth and reaching $10^8$ -- $10^9\, M_{\odot}$, albeit with most of the growth in (b) delayed compared to set (a). The finding of eight BH mergers in the small-seed scenario (four with masses $10^4$ -- $10^6\, M_{\odot}$ at $z>12$), six in the intermediate-seed scenario, and zero in the large-seed scenario suggests that the vast BHs in the small-seed scenario merge frequently during the early phases of the growth of SMBHs. The increased BH merger rate for the low-mass BH seed and halo threshold scenario provides an exciting prospect for discriminating BH formation mechanisms with the advent of multi-messenger astrophysics and next-generation gravitational wave facilities.

High-$z$ AGNs hosted in gas rich galaxies are expected to grow through significantly obscured accretion phases. This may limit or bias their observability. In this work, we use \textsc{BlueTides}, a large volume cosmological simulation of galaxy formation to examine quasar obscuration for the highest-redshift ($z \geq 7$) supermassive black holes residing in the center of galaxies. We find that for the bright quasars, most of the high column density gas ($>90\%$) resides in the innermost regions of the host galaxy, (typically within $< 10$ ckpc), while the gas in the outskirts is a minor contributor to the $N_\mathrm H$. The brightest quasars can have large angular variations in galactic obscuration, over 2 orders of magnitude, where the lines of sight with the lowest obscuration are those formed via strong gas outflows driven by AGN feedback. We find that for the overall AGN population, the mean $N_\mathrm H$ is generally larger for high luminosity and BH mass, while the $N_\mathrm H$ distribution is significantly broadened, developing a low $N_\mathrm H $ wing due to the angular variations driven by the AGN outflows/feedback. The obscured fraction P($N_{\rm H} > 10^{23} {\rm cm}^{-2}$) typically range from 0.6 to 1.0 for increasing $L_{X}$ (with $L_X > 10^{43} \rm{ergs/s}$), with no clear trend of redshift evolution. With respect to the galaxy host property, we find a linear relation between $N_{\rm H}$, $M_*$ and $M_{\rm H_2}$ with $\log N_{\rm H} = (0.24 \pm 0.03) \log M_{*} + (20.7 \pm 0.3)$ and $\log N_{\rm H} = (0.47 \pm 0.03) \log M_{\rm H_2} + (18.4 \pm 0.3)$. The dust optical depth in the UV band $\tau_{\mathrm UV}$ has tight positive correlation with $N_{\rm H}$. Our dust extincted UVLF is about 1.5 dex lower than the intrinsic UVLF, implying that more than 99\% of the $z \sim 7$ AGNs are heavily dust extincted and therefore would be missed by the UV band observation.

We examine the properties of the host galaxies of $z=7$ quasars using the large volume, cosmological hydrodynamical simulation BlueTides. We find that the most massive black holes and quasars are hosted by massive galaxies with stellar masses $\log(M_\ast/M_\odot)=10.8\pm0.2$, and $10.2\pm0.4$, which have a wide range of star formation rates, of $50\substack{+119 \\ -35}M_\odot/\rm{yr}$ and $19\substack{+28 \\ -12}M_\odot/\rm{yr}$, respectively. The hosts of the most massive black holes and quasars in BlueTides are generally bulge-dominated, with bulge-to-total mass ratio $B/T\simeq0.85\pm0.1$, however their morphologies are not biased relative to the overall $z=7$ galaxy sample. We find that the hosts of the most massive black holes and quasars are significantly more compact, with half-mass radii $R_{0.5}=0.29\substack{+0.15 \\ -0.10}$ kpc and $0.28\substack{+0.08 \\ -0.06}$ kpc respectively, relative to galaxies with similar masses, which have $R_{0.5}=0.5\substack{+0.3 \\ -0.2}$ kpc. We make mock James Webb Space Telescope (JWST) images of these quasars and their host galaxies. We find that distinguishing the host from the quasar emission will be possible but still challenging with JWST, due to the small sizes of quasar hosts. We find that currently observable quasar samples are biased tracers of the intrinsic black hole--stellar mass relations, following a relation that is 0.27 dex higher than that of the full galaxy sample. Finally, we find that black hole hosts are generally quite isolated. However, the most massive black holes are more likely to be found in denser environments than the typical $M_{\textrm{BH}}>10^{6.5}M_\odot$ black hole, indicating that minor mergers at least play some role in growing black holes in the early Universe.

We present the discovery and the photometric and spectroscopic study of H-rich Type II supernova (SN) KSP-SN-2016kf (SN2017it) observed in the KMTNet Supernova Program in the outskirts of a small irregular galaxy at $z\simeq0.043$ within a day from the explosion. Our high-cadence, multi-color ($BVI$) light curves of the SN show that it has a very long rise time ($t_\text{rise}\simeq 20$ days in $V$ band), a moderately luminous peak ($M_V\simeq -$17.6 mag), a notably luminous and flat plateau ($M_V\simeq -$17.4 mag and decay slope $s\simeq0.53$ mag per 100 days), and an exceptionally bright radioactive tail. Using the color-dependent bolometric correction to the light curves, we estimate the $^{56}$Ni mass powering the observed radioactive tail to be $0.10\pm0.01$ M$_\odot$, making it a H-rich Type II SN with one of the largest $^{56}$Ni masses observed to date. The results of our hydrodynamic simulations of the light curves constrain the mass and radius of the progenitor at the explosion to be $\sim$15 M$_\odot$ (evolved from a star with an initial mass of $\sim$ 18.8 M$_\odot$) and $\sim1040$ R$_\odot$, respectively, with the SN explosion energy of $\sim 1.3\times 10^{51}$ erg s$^{-1}$. The above-average mass of the KSP-SN-2016kf progenitor, together with its low metallicity $ Z/Z_\odot \simeq0.1-0.4$ obtained from spectroscopic analysis, is indicative of a link between the explosion of high-mass red supergiants and their low-metallicity environment. The early part of the observed light curves shows the presence of excess emission above what is predicted in model calculations, suggesting there is interaction between the ejecta and circumstellar material. We further discuss the implications of the high progenitor initial mass and low-metallicity environment of KSP-SN-2016kf on our understanding of the origin of Type II SNe.

During the last decades, rapid progress has been made in measurements of the rest-frame ultraviolet (UV) luminosity function (LF) for high-redshift galaxies ($z \geq 6$). The faint-end of the galaxy LF at these redshifts provides powerful constraints on different dark matter models that suppress small-scale structure formation. In this work we perform full hydrodynamical cosmological simulations of galaxy formation using an alternative DM model composed of extremely light bosonic particles ($m \sim 10^{-22}$ eV), also known as fuzzy dark matter (FDM), and examine the predictions for the galaxy stellar mass function and luminosity function at $z \geq 6$ for a range of FDM masses. We find that for FDM models with bosonic mass $m = 5\times10^{-22}$ eV, the number density of galaxies with stellar mass $\rm M_* \sim 10^7 M_{\odot}$ is suppressed by $\sim 40\%$ at z = 9, $\sim 20\%$ at z = 5, and the UV LFs within magnitude range of -16 < $M_{\rm UV}$ < -14 is suppressed by $\sim 60\%$ at $z = 9$, $\sim 20\%$ at $z = 5$ comparing to the CDM counterpart simulation. Comparing our predictions with current measurements of the faint-end LFs ($-18 \leqslant M_{\rm UV} \leqslant -14$), we find that FDM models with $m_{22} < 5\times10^{-22}$ are ruled out at $3\sigma$ confidence level. We expect that future LF measurements by James Webb Space Telescope (JWST), which will extend down to $M_{\rm UV} \sim -13$ for $z \lesssim 10$, with a survey volume that is comparable to the Hubble Ultra Deep Field (HUDF) would have the capability to constrain FDM models to $m\; \gtrsim 10^{-21}$ eV.

Many theoretical models predict that quasar driven outflows account for the observed quenching of star formation in massive galaxies. There is growing observational evidence for quasar-launched massive outflows at z>6, while the details of outflow-host galaxy interaction remain obscure. In this paper, we study the feedback around the highest redshift quasar in the BlueTides simulation, the largest volume cosmological hydrodynamic simulation so far carried out. We present predictions for gas outflows around the brightest $z = 7.54$ quasar which hosts the most massive black hole in the simulation volume, which has grown to black hole mass $6.7\times 10^{8}{\rm M}_\odot$ consistent with the current record holder for high-z quasars. We introduce a method to identify and trace the gas outflowing from the halo. We find that the total mass of the outflow gas is about $3.6\times 10^{9}{\rm M}_\odot$, constituting 6\% of the total gas in the halo. The outflow gas contains a cold, dense molecular component with mass about $2.6\times 10^{8}{\rm M}_\odot$, that originates from the inner region of the halo, within a few kpc of the central black hole. The velocities of the outflow gas reach thousands of km/s, within which the molecular component has mass averaged outward radial velocity of $1300$ km/s, consistent with observations. The averaged outflow rate is about $200-300 {\rm M}_\odot/yr$, with the outflowing gas mainly in a hotter ($T \sim 10^7$ K) and lower density state than the average of the host halo gas.

We present the discovery of a rapidly evolving transient by the Korean Microlensing Telescope Network Supernova Program (KSP). KSP is a novel high-cadence supernova survey that offers deep ($\sim21.5$ mag in $BVI$ bands) nearly continuous wide-field monitoring for the discovery of early and/or fast optical transients. KSP-OT-201509a, reported here, was discovered on 2015 September 27 during the KSP commissioning run in the direction of the nearby galaxy NGC~300, and stayed above detection limit for $\sim$ 22 days. We use our $BVI$ light-curves to constrain the ascent rate, $-3.7(7)$ mag day$^{-1}$ in $V$, decay time scale, $t^{V}_{2}=1.7(6)$ days, and peak absolute magnitude, $-9.65\leq M_{V}\leq -9.25$ mag. We also find evidence for a short-lived pre-maximum halt in all bands. The peak luminosity and lightcurve evolution make KSP-OT-201509a consistent with a bright, rapidly decaying nova outburst. We discuss constraints on the nature of the progenitor and its environment using archival HST/ACS images and conclude with a broad discussion on the nature of the system.

The spacetime geometry around astrophysical black holes is supposed to be well approximated by the Kerr metric, but deviations from the Kerr solution are predicted in a number of scenarios involving new physics. Broad iron K$\alpha$ lines are commonly observed in the X-ray spectrum of black holes and originate by X-ray fluorescence of the inner accretion disk. The profile of the iron line is sensitively affected by the spacetime geometry in the strong gravity region and can be used to test the Kerr black hole hypothesis. In this paper, we extend previous work in the literature. In particular: i) as test-metric, we employ the parametrization recently proposed by Konoplya, Rezzolla, and Zhidenko, which has a number of subtle advantages with respect to the existing approaches; ii) we perform simulations with specific X-ray missions, and we consider NuSTAR as a prototype of current observational facilities and eXTP as an example of the next generation of X-ray observatories. We find a significant difference between the constraining power of NuSTAR and eXTP. With NuSTAR, it is difficult or impossible to constrain deviations from the Kerr metric. With eXTP, in most cases we can obtain quite stringent constraints (modulo we have the correct astrophysical model).

Recently, a family of hairy black holes in 4-dimensional Einstein gravity minimally coupled to a complex, massive scalar field was discovered~\cite{hbh}. Besides the mass $M$ and spin angular momentum $J$, these objects are characterized by a Noether charge $Q$, measuring the amount of scalar hair, which is not associated to a Gauss law and cannot be measured at spatial infinity. Introducing a dimensionless scalar hair parameter $q$, ranging from 0 to 1, we recover (a subset of) Kerr black holes for $q=0$ and a family of rotating boson stars for $q=1$. In the present paper, we explore the possibility of measuring $q$ for astrophysical black holes with current and future X-ray missions. We study the iron K$\alpha$ line expected in the reflection spectrum of such hairy black holes and we simulate observations with Suzaku and eXTP. As a proof of concept, we point out, by analyzing a sample of hairy black holes, that current observations can already constrain the scalar hair parameter $q$, because black holes with $q$ close to 1 would have iron lines definitively different from those we observe in the available data. We conclude that a detailed scanning of the full space of solutions, together with data from the future X-ray missions, like eXTP, will be able to put relevant constraints on the astrophysical realization of Kerr black holes with scalar hair.

We investigate the multiple bifurcations in periodic orbit families in the potential field of a highly irregular-shaped celestial body. Topological cases of periodic orbits and four kinds of basic bifurcations in periodic orbit families are studied. Multiple bifurcations in periodic orbit families consist of four kinds of basic bifurcations. We found both binary period-doubling bifurcations and binary tangent bifurcations in periodic orbit families around asteroid 433 Eros. The periodic orbit family with binary period-doubling bifurcations is nearly circular, with almost zero inclination, and is reversed relative to the body of the asteroid 433 Eros. This implies that there are two stable regions separated by one unstable region for the motion around this asteroid. In addition, we found triple bifurcations which consist of two real saddle bifurcations and one period-doubling bifurcation. A periodic orbit family generated from an equilibrium point of asteroid 433 Eros has five bifurcations, which are one real saddle bifurcation, two tangent bifurcations, and two period-doubling bifurcations.