Using mock spectra based on Vazdekis/MILES library fitted within the wavelength region 3600-7350Å, we analyze the bias and scatter on the resulting physical parameters induced by the choice of fitting algorithms and observational uncertainties, but avoid effects of those model uncertainties. We consider two full-spectrum fitting codes: pPXF and STARLIGHT, in fitting for stellar population age, metallicity, mass-to-light ratio, and dust extinction. With pPXF we find that both the bias μ in the population parameters and the scatter σ in the recovered logarithmic values follows the expected trend $$\mu \propto \sigma \propto 1/(\rm S/N)$$. The bias increases for younger ages and systematically makes recovered ages older, M_*/L_r larger and metallicities lower than the true values. For reference, at S/N=30, and for the worst case (t = 10⁸yr), the bias is 0.06 dex in M/L_r, 0.03 dex in both age and [M/H]. There is no significant dependence on either E(B-V) or the shape of the error spectrum. Moreover, the results are consistent for both our 1-SSP and 2-SSP tests. With the STARLIGHT algorithm, we find trends similar to pPXF, when the input E(B-V)<0.2 mag. However, with larger input E(B-V), the biases of the output parameter do not converge to zero even at the highest S/N and are strongly affected by the shape of the error spectra. This effect is particularly dramatic for youngest age (t = 10⁸yr), for which all population parameters can be strongly different from the input values, with significantly underestimated dust extinction and [M/H], and larger ages and M_*/L_r. Results degrade when moving from our 1-SSP to the 2-SSP tests. The STARLIGHT convergence to the true values can be improved by increasing Markov Chains and annealing loops to the “slow mode”. For the same input spectrum, pPXF is about two order of magnitudes faster than STARLIGHT’s “default mode” and about three order of magnitude faster than STARLIGHT’s “slow mode”.

Cosmic hydrogen reionization and cosmic production of the first metals are major phase transitions of the Universe occurring during the first billion years after the Big Bang; however, these are still underexplored observationally. Using the JWST/NIRSpec prism spectroscopy, we report the discovery of a sub-L _* galaxy at z _spec = 8.1623 ± 0.0007, dubbed RX J2129–z8He II, via the detection of a series of strong rest-frame UV/optical nebular emission lines and the clear Lyman break. RX J2129–z8He II shows a pronounced UV continuum with an extremely steep (i.e., blue) spectral slope of β = − 2.53 − 0.07 + 0.06 , the steepest among all spectroscopically confirmed galaxies at z _spec ≳ 7, in support of its very hard ionizing spectrum that could lead to a significant leakage of its ionizing flux. Therefore, RX J2129–z8He II is representative of the key galaxy population driving the cosmic reionization. More importantly, we detect a strong He II λ1640 emission line in its spectrum, one of the highest redshifts at which such a line is robustly detected. Its high rest-frame equivalent width (EW = 21 ± 4 Å) and extreme flux ratios with respect to UV metal and Balmer lines raise the possibility that part of RX J2129–z8He II’s stellar population could be Pop III (Pop III)-like. Through careful photoionization modeling, we show that the physically calibrated phenomenological models of the ionizing spectra of Pop III stars with strong mass loss can successfully reproduce the emission line flux ratios observed in RX J2129–z8He II. Assuming the Eddington limit, the total mass of the Pop III stars within this system is estimated to be 7.8 ± 1.4 × 10⁵ M _⊙. To date, this galaxy presents the most compelling case in the early Universe where trace Pop III stars might coexist with metal-enriched populations.

The stellar mass-to-light ratio gradient in SDSS r-band ∇(M_*/L_r) of a galaxy depends on its mass assembly history, which is imprinted in its morphology and gradients of age, metallicity, and stellar initial mass function (IMF). Taking a MaNGA sample of 2051 galaxies with stellar masses ranging from 10⁹ to 10¹²M_⊙ released in SDSS DR15, we focus on face-on galaxies, without merger and bar signatures, and investigate the dependence of the 2D ∇(M_*/L_r) on other galaxy properties, including M_*/L_r-colour relationships by assuming a fixed Salpeter IMF as the mass normalization reference. The median gradient is ∇M_*/L_r ∼ −0.1 (i.e. the M_*/L_r is larger at the centre) for massive galaxies, becomes flat around M_* ∼ 10¹⁰M_⊙ and change sign to ∇M_*/L_r ∼ 0.1 at the lowest masses. The M_*/L_r inside a half-light radius increases with increasing galaxy stellar mass; in each mass bin, early-type galaxies have the highest value, while pure-disc late-type galaxies have the smallest. Correlation analyses suggest that the mass-weighted stellar age is the dominant parameter influencing the M_*/L_r profile, since a luminosity-weighted age is easily affected by star formation when the specific star formation rate (sSFR) inside the half-light radius is higher than 10⁻³ Gyr⁻¹. With increased sSFR gradient, one can obtain a steeper negative ∇(M_*/L_r). The scatter in the slopes of M_*/L-colour relations increases with increasing sSFR, for example, the slope for post-starburst galaxies can be flattened to 0.45 from the global value 0.87 in the M_*/L versus g − r diagram. Hence converting galaxy colours to M_*/L should be done carefully, especially for those galaxies with young luminosity-weighted stellar ages, which can have quite different star formation histories.

Three basic ingredients are required to generate a simple stellar population (SSP) library, i.e., an initial mass function (IMF), a stellar evolution model/isochrones, and an empirical/theoretical stellar spectral library. However, there are still some uncertainties to the determination and understanding of these ingredients. We perform the spectral fitting to test the relative parameter offsets between these uncertainties using two different stellar population models, two different empirical stellar libraries, two different isochrones, and the Salpeter and Chabrier IMFs. Based on these setups, we select five SSP libraries generated with the Galaxev/STELIB and Vazdekis/MILES models, and apply them to the pPXF full-spectrum fitting of both MaNGA and mock spectra. We find that: 1) Compared to the Galaxev/STELIB model, spectral fitting qualities with the Vazdekis/MILES model have significant improvements for those metal-rich (especially over-solar) spectra, which cause better reduced χ² distributions and more precisely fitted absorption lines. This might due to the lack of metal rich stars in the empirical STELIB library, or code improvement of the Vazdekis model. 2) When applying the Vazdekis/MILES model for spectral fitting, the IMF variation will lead to not only a systematic offset in M_*/L_r, but also offsets in age and metallicity, and these offsets increase with increasing stellar population ages. However, the IMF-variation caused metallicity offsets disappear in the case of Galaxev/STELIB based libraries. 3) The Padova2000 model provides a better match to the MaNGA galaxy spectra at [M/H]_L < −1.0, while the BaSTI model match the local galaxy spectra better at [M/H]_L > −1.0. Current tests suggest that spectral fitting with the Vazdekis/MILES+BaSTI combination would be a better choice for local galaxies.

Periodic quasars have been suggested to host supermassive binary black holes (BBHs) in their centers, and their optical/UV periodicities are interpreted as caused by either the Doppler-boosting (DB) effect of continuum emission from the disk around the secondary black hole (BH) or intrinsic accretion rate variation. However, no other definitive evidence has been found to confirm such a BBH interpretation(s). In this paper, we investigate the responses of broad emission lines (BELs) to the continuum variations for these quasars under two BBH scenarios and check whether they can be distinguished from each other and from that of a single BH system. We assume a simple circumbinary broad-line region (BLR) model, compatible with BLR size estimates, with a standard Γ distribution of BLR clouds. We find that BELs may change significantly and periodically under the BBH scenarios due to (1) the position variation of the secondary BH and (2) the DB effect, if significant, and/or intrinsic variation, which is significantly different from the case of a single BH system. For the two BBH scenarios, the responses of BELs to (apparent) continuum variations, caused by the DB effect or intrinsic rate variation, are also significantly different from each other, mainly because the DB effect has a preferred direction along the direction of motion of the secondary BH, while that due to intrinsic variation does not. Such differences in the responses of BELs from different scenarios may offer a robust way to distinguish different interpretations of periodic quasars and to identify BBHs, if any, in these systems.

Major mergers of galaxies are considered to be an efficient way to trigger Active Galactic Nuclei and are thought to be responsible for the phenomenon of quasars. This has however recently been challenged by observations of a large number of low luminosity Active Galactic Nuclei at low redshift ( z ≲ 1) without obvious major merger signatures. Minor mergers are frequently proposed to explain the existence of these Active Galactic Nuclei. In this paper, we perform nine high resolution hydrodynamical simulations of minor galaxy mergers, and investigate whether nuclear activities can be efficiently triggered by minor mergers, by setting various properties for the progenitor galaxies of those mergers. We find that minor galaxy mergers can activate the massive black hole in the primary galaxy with an Eddington ratio of f _Edd > 0.01 and > 0.05 (or a bolometric luminosity >10 ⁴³ and >10 ⁴⁴ erg s ⁻¹ ) with a duration of 2.71 and 0.49Gyr (or 2.69 and 0.19Gyr), respectively. The nuclear activity of the primary galaxy strongly depends on the nucleus separation, such that the nucleus is more active as the two nuclei approach each other. Dual Active Galactic Nuclei systems can still possibly be formed by minor mergers of galaxies, though the time duration for dual Active Galactic Nuclei is only ∼ 0.011 Gyr and ∼ 0.017 Gyr with Eddington ratio of f _Edd > 0.05 and bolometric luminosity >10 ⁴⁴ erg s ⁻¹ . This time period is typically shorter than that of dual Active Galactic Nuclei induced by major galaxy mergers.