Search

Non-axisymmetric structures, such as bars and spiral arms, are known to concentrate molecular gas and star formation in galaxy centers, actively building up the pseudo-bulges. However, a direct link between the neutral (i.e., molecular and atomic) gas distribution and the exerted torque forces over a broader radial range of galactic disks still remains to be explored. In the present work, we investigate this link by carefully evaluating the torque force field using the $3.6\, \mathrm{\mu m}$ images for 17 The H I Nearby Galaxy Survey (THINGS) galaxies, and measuring neutral gas distribution on resolved atomic and molecular line maps. We find that galaxies with stronger torque forces show a more concentrated neutral gas distribution over the disk-scale, defined as half the isophotal radius at $25.5\, \mathrm{mag\, arcsec^{-2}}$. The correlation holds regardless of whether the neutral gas fraction, or the effective stellar mass surface density is controlled for. In addition, $\mathrm{kpc}$-scale neutral gas over-densities tend to be located close to the local maxima of torque forces. Most of these correlations involving the torque forces are comparatively stronger than those using the traditional Fourier amplitudes to quantify the non-axisymmetric structures. These results are consistent with the scenario that non-axisymmetric structures exert torque forces, and trigger dissipative processes to transport gas inward, not only to build the pseudo-bulges, but also fuel the inner disk growth. In this regard, non-axisymmetric structures inducing stronger torque forces appear to be more efficient in these processes., Comment: Accepted for publication in ApJ, 16 pages (main text), 14 figures, 4 tables

Type-C quasi-periodic oscillations (QPOs) in black hole X-ray transients typically manifest in the low-hard and hard-intermediate states. This study presents a detailed spectral and temporal analysis of the black hole candidate Swift J1727.8-1613 using NICER observations from August and September 2023, with a focus on the first flare period. The time-averaged spectra, along with the rms and phase-lag spectra of the QPOs, were jointly fitted using the time-dependent Comptonization model \text{vkompthdk} to examine the geometry of the corona during this flare. The results provide a comprehensive view of the QPO, where we detected type-C QPOs with a centroid frequency increasing from 0.32 Hz to 2.63 Hz, while it elevated when it entered the flare state, and its energy spectral properties as they evolved during the first flare period. Correlations between spectral and temporal properties suggest that type-C QPOs are primarily modulated by Lense-Thirring precession. Based on simultaneous radio observations indicating discrete jet ejections, we propose, for the first time, a scenario where the temporarily extended corona contracts vertically from approximately 2714 km to less than 900 km, overlying the inner accretion disc, with a transient jet seemingly being launched. The corona then recovers to nearly 2000 km by the end of the first flare period of Swift J1727.8-1613, rather than undergoing horizontal changes. A phenomenological analysis of the corona scenario during this flare period was also conducted., Comment: 16 pages, 6 figures, 3 tables, Submitted to ApJ

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

We present a catalog of clouds identified from the $^{12}$CO (1--0) data of M83, which was observed using Atacama Large Millimeter/submillimeter Array (ALMA) with a spatial resolution of $\sim$46 pc and a mass sensitivity of $\sim$10$^4$ $M_{\odot}$ (3 $\sigma$). The almost full-disk coverage and high sensitivity of the data allowed us to sample 5724 molecular clouds with a median mass of $\sim1.9$ $\times$ $10^5$ $M_{\odot}$, which is comparable to the most frequently sampled mass of Giant Molecular Clouds by surveys in the Milky Way. About 60 percent of the total CO luminosity in M83's disk arises from clouds more massive than 10$^6$ $M_{\odot}$. Such massive clouds comprise 16 percent of the total clouds in number and tend to concentrate toward the arm, bar, and center, while smaller clouds are more prevalent in inter-arm regions. Most $>10^6$ $M_{\odot}$ clouds have peak brightness temperatures $T_{\mathrm{peak}}$ above 2 K with the current resolution. Comparing the observed cloud properties with the scaling relations determined by Solomon et al. 1987 (S87), $T_{\mathrm{peak}}$$>2$ K clouds follow the relations, but $T_{\mathrm{peak}}$$<2$ K clouds, which are dominant in number, deviate significantly. Without considering the effect of beam dilution, the deviations would suggest modestly high virial parameters and low surface mass densities for the entire cloud samples, which are similar to values found for the Milky Way clouds by Rice et al. (2016) and Miville-Desch{\^e}nes et al. (2017). However, once beam dilution is taken into account, the observed $\alpha_{\mathrm{vir}}$ and $\Sigma$ for a majority of the clouds (mostly $T_{\mathrm{peak}}$ $<2$ K) can be potentially explained with intrinsic $\Sigma$ of $\sim$100 $M_{\mathrm{\odot}}\ \mathrm{pc}^{-2}$ and $\alpha_{\mathrm{vir}}$ of $\sim$1, which are similar to the clouds of S87., Comment: Accepted for publication in ApJ

We search for the galaxies associated with the intervening MgII absorbers over a redshift range of $0.4 \le z \le 1.0$ using imaging data from DESI Legacy Imaging Surveys and measure the redshift based on direct detection of nebular emission in the background quasar spectra from SDSS survey. We find 270 MgII absorbers associated with strong [O II] nebular emissions, at $2.5\sigma$ level. Among them, for 213 MgII absorbers, we detect an absorber host galaxy at impact parameters of $4 \le \rho \le 16$ kpc, including a galaxy pair associated with three absorbers, with best-fit galaxy SED model based on multi-passband photometric data from DESI Legacy Imaging surveys, supplemented with the infrared VISTA and unWISE imaging surveys. The detection rate of MgII absorber host with strong [O II] nebular emission in the finite SDSS fiber of 2--3 arcsec diameter increases from 0.2% to $\sim$3% with increasing equivalent width from $0.3$\r{A} to $\sim$ 3.5\r{A}, which remains near-constant across the probed redshift range. The associated MgII host galaxies exhibit a wide range of stellar mass from 7.94 $\le \rm log(M_{\star}/M_{\odot}) \le $ 11.11, with an average star-formation rate (SFR) of $ 5.81 \rm M_{\odot}\ yr^{-1}$. The MgII absorber hosts selected based on [O II] nebular emission mostly exhibit active star-forming systems including 13.4% starburst systems, but 10.2% with suppressed SFR. The near-constant absorption strength at low-impact parameters suggests a high gas covering fraction. We find that the MgII equivalent width ($W_{2796}$) positively correlates with SFR and specific-SFR, likely indicate their wind origin. The average velocity offset between the host and absorber suggests that the MgII gas is bound within the dark matter halo., Comment: 9 Pages, 8 Figures, submitted to A&A

We investigate the kinematics and dynamics of the molecular and ionized gas in the host galaxies of three Palomar-Green quasars at low redshifts, benefiting from the archival millimeter-wave interferometric and optical integral field unit data. We study the kinematics of both cold molecular and hot ionized gas by analyzing the CO and H$\alpha$ data cubes, and construct the mass distributions of our sample through gas dynamics, utilizing a priori knowledge regarding the galaxy light distribution. We find no systematic offset between the stellar mass derived from our dynamical method and that from the broad-band photometry and mass-to-light ratio, suggesting the consistency of both methods. We then study the kinetic pressure and the weight of the interstellar medium using our dynamical mass model. By studying the relationship between kinetic pressure and gravitational pressure of the quasar host galaxies, we find an equivalence in the hydrostatic equilibrium states of ISM in the quasar host galaxies, similar to the result of gas equilibrium in normal star-forming galaxies, suggesting minimal quasar feedback. Regarding non-circular motion as indicative of quasar-driven outflows, we observe an exceptionally low coupling efficiency between molecular gas outflow and AGN bolometric luminosities. These results demonstrate the marginal influence of the central engine on the properties of cold molecular gas in quasar host galaxies., Comment: 27 pages, 12 figures, accepted for publication in ApJ

The properties of slim accretion disks, while crucial for our understanding of black hole growth, have yet to be studied extensively observationally. We analyze the multi-epoch broad-band spectral energy distribution of the changing-look active galactic nucleus 1ES 1927+654 to derive the properties of its complex, time-dependent accretion flow. The accretion rate decays as $\dot{M} \propto t^{-1.53}$, consistent with the tidal disruption of a $1.1\, M_\odot$ star. Three components contribute to the spectral energy distribution: a central overheated zone resembling a slim disk, an outer truncated thin disk, and a hot corona. Photon trapping in the slim disk triggered by the high initial $\dot{M}$ was characterized by a low radiation efficiency ($3\%$), which later more than doubled ($8\%$) after $\dot{M}$ dropped sufficiently low for the disk to transition to a geometrically thin state. The blackbody temperature profile $T \propto R^{-0.60}$ for the inner overheated zone matches the theoretical expectations of a slim disk, while the effective temperature profile of $T \propto R^{-0.69}$ for the outer zone is consistent with the predictions of a thin disk. Both profiles flatten toward the inner boundary of the disk as a result of Compton cooling in the corona. Our work presents compelling observational evidence for the existence of slim accretion disks and elucidates the key parameters governing their behavior, paving the way for further exploration in this area., Comment: Accepted for publication in The Astrophysical Journal. 27 pages

Time-domain studies of active galactic nuclei (AGNs) offer a powerful tool for understanding black hole accretion physics. Prior to the optical outburst on 23 December 2017, 1ES 1927+654 was classified as a "true" type~2 AGN, an unobscured source intrinsically devoid of broad-line emission in polarized spectra. Through our three-year monitoring campaign spanning X-ray to ultraviolet/optical wavelengths, we analyze the post-outburst evolution of the spectral energy distribution (SED) of 1ES 1927+654. Examination of the intrinsic SED and subsequent modeling using different models reveal that the post-outburst spectrum is best described by a combination of a disk, blackbody, and corona components. We detect systematic SED variability and identify four distinct stages in the evolution of these components. During the event the accretion rate is typically above the Eddington limit. The correlation between ultraviolet luminosity and optical to X-ray slope ($\alpha_\mathrm{OX}$) resembles that seen in previous studies of type 1 AGNs, yet exhibits two distinct branches with opposite slopes. The optical bolometric correction factor ($\kappa_{5100}$) is $\sim 10$ times higher than typical AGNs, again displaying two distinct branches. Correlations among the corona optical depth, disk surface density, and $\alpha_\mathrm{OX}$ provide compelling evidence of a disk-corona connection. The X-ray corona showcases systematic variation in the compactness-temperature plot. Between 200 and 650 days, the corona is "hotter-when-brighter", whereas after 650 days, it becomes "cooler-when-brighter". This bimodal behavior, in conjunction with the bifurcated branches of $\alpha_\mathrm{OX}$ and $\kappa_{5100}$, offers strong evidence of a transition from a slim disk to thin disk $\sim 650$ days after the outburst., Comment: Accepted for publication in The Astrophysical Journal. 36 pages. For companion paper on the broad-line region (Li et al. 2022, ApJ, 933, 70), see https://iopscience.iop.org/article/10.3847/1538-4357/ac714a

We present 0$.\!\!^{\prime\prime}{22}$-resolution CO(2$-$1) observations of the circumnuclear gas disk in the local compact galaxy NGC 384 with the Atacama Large Millimeter/submillimeter Array (ALMA). While the majority of the disk displays regular rotation with projected velocities rising to $370$ km s$^{-1}$, the inner $\sim$0\farcs{5} exhibits a kinematic twist. We develop warped disk gas-dynamical models to account for this twist, fit those models to the ALMA data cube, and find a stellar mass-to-light ratio in the $H$-band of \mlabstract\ and a supermassive black hole (BH) mass ($M_{\mathrm{BH}}$) of $M_{\mathrm{BH}}$ $= (7.26^{+0.43}_{-0.48}$ [$1\sigma$ statistical] $^{+0.55}_{-1.00}$ [systematic])$\times 10^8$ $M_\odot$. In contrast to most previous dynamical $M_{\mathrm{BH}}$ measurements in local compact galaxies, which typically found over-massive BHs compared to the local BH mass$-$bulge luminosity and BH mass$-$bulge mass relations, NGC 384 lies within the scatter of those scaling relations. NGC 384 and other local compact galaxies are likely relics of $z\sim2$ red nuggets, and over-massive BHs in these relics indicate BH growth may conclude before the host galaxy stars have finished assembly. Our NGC 384 results may challenge this evolutionary picture, suggesting there may be increased scatter in the scaling relations than previously thought. However, this scatter could be inflated by systematic differences between stellar- and gas-dynamical measurement methods, motivating direct comparisons between the methods for NGC 384 and the other compact galaxies in the sample., Comment: 19 pages, 8 figures, accepted for publication in ApJ

Most galaxies follow well-defined scaling relations of metallicity and stellar mass; however, some outliers at the low mass end of the observed galaxy population exhibit unusually high metallicity for their mass. Understanding how these objects get to be so metal-rich is vital for understanding the role of feedback in galaxy formation. Using the TNG50 simulation, we explore the origins of this phenomenon. We identify 227 metal-rich, Compact Stellar Systems (CSSs) that deviate significantly from this scaling relation. These CSSs are satellites located in the vicinity of massive host galaxies, with stellar masses ranging from $10^{8} M_{\odot}$ to $10^{10} M_{\odot}$ (including six systems that are close analogs of the M31-M32 system). Contrary to the previously assumed scenario that such objects are predominantly products of tidal stripping, our results suggest a more prevalent role for ram pressure in their formation. Indeed, 76\% (173) of these CSSs are formed through a burst of star formation occurring around the time of the first pericentric passage, typically at redshifts $z\lesssim1$, aided by strong ram pressure and tidal forces. The high ram pressure, resulting from the CSSs' rapid motion near the halo center, facilitates metal enrichment, producing high-metallicity CSSs by confining the metal-rich gas from bursty star formation, which leads to distinct stellar populations characterized by enhanced metallicity as well as high $\alpha$-abundance. Only the remaining 24\% (54) of metal-rich CSSs are generated through the tidal stripping of massive progenitors. Our results further indicate that M32 is more likely to have formed through intense star formation events rather than through gradual, tidal stripping, thereby providing crucial insights into the nature of low mass, compact galaxy formation., Comment: 28 pages, 13 figures. Submitted

The cold interstellar medium (ISM) as the raw material for star formation is critical to understanding galaxy evolution. It is generally understood that galaxies stop making stars when, in one way or another, they run out of gas. However, here we provide evidence that central spiral galaxies remain rich in atomic gas even if their star formation rate and molecular gas fraction have dropped significantly compared to "normal" star-forming galaxies of the same mass. Since HI is sensitive to external processes, here we investigate central spiral galaxies using a combined sample from SDSS, ALFALFA, and xGASS surveys. After proper incompleteness corrections, we find that the key HI scaling relations for central spirals show significant but regular systematic dependence on stellar mass. At any given stellar mass, the HI gas mass fraction is about constant with changing specific star formation rate (sSFR), which suggests that HI reservoir is ubiquitous in central spirals with any star formation status down to M* ~ 10^9 Msun. Together with the tight correlation between the molecular gas mass fraction and sSFR for galaxies across a wide range of different properties, it suggests that the decline of SFR of all central spirals in the local universe is due to the halt of H2 supply, though there is plenty of HI gas around. These hence provide critical observations of the dramatically different behavior of the cold multi-phase ISM, and a key to understand the star formation process and quenching mechanism., Comment: 18 pages, 7 figures; Accepted for publication in the ApJL; This is the fourth paper in the "From Haloes to Galaxies" series

Over the past three decades, multiple reverberation mapping (RM) campaigns conducted for the quasar PG 2130+099 have exhibited inconsistent findings with time delays ranging from $\sim$10 to $\sim$200 days. To achieve a comprehensive understanding of the geometry and dynamics of the broad-line region (BLR) in PG 2130+099, we continued an ongoing high-cadence RM monitoring campaign using the Calar Alto Observatory 2.2m optical telescope for an extra four years from 2019 to 2022. We measured the time lags of several broad emission lines (including He II, He I, H$\beta$, and Fe II) with respect to the 5100 {\AA} continuum, and their time lags continuously vary through the years. Especially, the H$\beta$ time lags exhibited approximately a factor of two increase in the last two years. Additionally, the velocity-resolved time delays of the broad H$\beta$ emission line reveal a back-and-forth change between signs of virial motion and inflow in the BLR. The combination of negligible ($\sim$10%) continuum change and substantial time-lag variation (over two times) results in significant scatter in the intrinsic $R_{\rm H\beta}-L_{\rm 5100}$ relationship for PG 2130+099. Taking into account the consistent changes in the continuum variability time scale and the size of the BLR, we tentatively propose that the changes in the measurement of the BLR size may be affected by 'geometric dilution'., Comment: 21 pages, 13 figures, 7 tables; accepted for publication in ApJ

In $\Lambda$CDM cosmology, galaxies form and evolve in their host dark matter (DM) halos. Halo mass is crucial for understanding the halo-galaxy connection. The abundance matching (AM) technique has been widely used to derive the halo masses of galaxy groups. However, quenching of the central galaxy can decouple the coevolution of its stellar mass and DM halo mass. Different halo assembly histories can also result in significantly different final stellar mass of the central galaxies. These processes can introduce substantial uncertainties in the halo masses derived from the AM method, particularly leading to a systematic bias between groups with star-forming centrals (blue groups) and passive centrals (red groups). To improve, we developed a new machine learning (ML) algorithm that accounts for these effects and is trained on simulations. Our results show that the ML method eliminates the systematic bias in the derived halo masses for blue and red groups and is, on average, $\sim1/3$ more accurate than the AM method. With careful calibration of observable quantities from simulations and observations from SDSS, we apply our ML model to the SDSS Yang et al. groups to derive their halo masses down to $10^{11.5}\mathrm{M_\odot}$ or even lower. The derived SDSS group halo mass function agrees well with the theoretical predictions, and the derived stellar-to-halo mass relations for both red and blue groups matches well with those obtained from direct weak lensing measurements. These new halo mass estimates enable more accurate investigation of the galaxy-halo connection and the role of the halos in galaxy evolution., Comment: 21 pages, 9 figures. Submitted to ApJ, comments welcome!

The spectral energy distribution (SED) of galaxies is essential for deriving fundamental properties like stellar mass and star formation history (SFH). However, conventional methods, including both parametric and non-parametric approaches, often fail to accurately recover the observed cosmic star formation rate (SFR) density due to oversimplified or unrealistic assumptions about SFH and their inability to account for the complex SFH variations across different galaxy populations. To address this issue, we introduce a novel approach that improves galaxy broad-band SED analysis by incorporating physical priors derived from hydrodynamical simulations. Tests using IllustrisTNG simulations demonstrate that our method can reliably determine galaxy physical properties from broad-band photometry, including stellar mass within 0.05 dex, current SFR within 0.3 dex, and fractional stellar formation time within 0.2 dex, with a negligible fraction of catastrophic failures. When applied to the SDSS main photometric galaxy sample with spectroscopic redshift, our estimates of stellar mass and SFR are consistent with the widely-used MPA-JHU and GSWLC catalogs. Notably, using the derived SFHs of individual SDSS galaxies, we estimate the cosmic SFR density and stellar mass density with remarkable consistency to direct observations up to $z \sim 6$. This marks the first time SFHs derived from SEDs can accurately match observations. Consequently, our method can reliably recover observed spectral indices such as $\rm D_{\rm n}(4000)$ and $\rm H\delta_{\rm A}$ by synthesizing the full spectra of galaxies using the estimated SFHs and metal enrichment histories, relying solely on broad-band photometry as input. Furthermore, this method is extremely computationally efficient compared to conventional approaches., Comment: 22 pages, 12 figures. Submitted to ApJ, comments welcome!

Galactic-wide outflows driven by active galactic nuclei (AGNs) is a routinely invoked feedback mechanism in galaxy evolution models. Hitherto, the interplay among the interstellar gas on galactic scales, the propagation of AGN outflows and the fundamental AGN parameters during evolution remains elusive. Powerful nuclear outflows are found to favorably exist at early AGN stages usually associated with high accretion rates and weak narrow emission lines. In a sample of quasars emitting Mg II narrow absorption lines (NALs) from the Sloan Digital Sky Survey, we discover an unprecedented phenomenon where galaxy-scale inflow-dominated transforming into outflow-dominated gas accompanied by an increasing strength of the narrow [O III] line, at a confidence level of 6.7{\sigma}. The fact that nuclear outflows diminish while galaxy-wide outflows intensifies as AGNs evolve implies that early-stage outflows interact with interstellar medium on galactic scales and trigger the gradual transformation into galaxy-wide outflows, providing observational links to the hypothetical multi-stage propagation of AGN outflows that globally regulates galaxy evolution., Comment: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy; 15 pages, 4 figures and 6 appendix figures

High-redshift radio(-loud) galaxies (H$z$RGs) are massive galaxies with powerful radio-loud active galactic nuclei (AGNs) and serve as beacons for protocluster identification. However, the interplay between H$z$RGs and the large-scale environment remains unclear. To understand the connection between H$z$RGs and the surrounding obscured star formation, we investigated the overdensity and spatial distribution of submillimeter-bright galaxies (SMGs) in the field of 4C\,23.56, a well-known H$z$RG at $z=2.48$. We used SCUBA-2 data ($\sigma\,{\sim}\,0.6$\,mJy) to estimate the $850\,{\rm \mu m}$ source number counts and examine the radial and azimuthal overdensities of the $850\,{\rm \mu m}$ sources in the vicinity of the H$z$RG. The angular distribution of SMGs is inhomogeneous around the H$z$RG 4C\,23.56, with fewer sources oriented along the radio jet. We also find a significant overdensity of bright SMGs (${\rm S}_{850\rm\,\mu m}\geq5\,$mJy). Faint and bright SMGs exhibit different spatial distributions. The former are concentrated in the core region, while the latter prefer the outskirts of the H$z$RG field. High-resolution observations show that the seven brightest SMGs in our sample are intrinsically bright, suggesting that the overdensity of bright SMGs is less likely due to the source multiplicity., Comment: 19 pages, 17 figures, 5 tables, accepted to A&A

Measuring the environments of massive galaxies at high redshift is crucial to understanding galaxy evolution and the conditions that gave rise to the distribution of matter we see in the Universe today. While high-$z$ radio galaxies (H$z$RGs) and quasars tend to reside in protocluster-like systems, the environments of their radio-quiet counterparts are relatively unexplored, particularly in the submillimetre, which traces dust-obscured star formation. In this study we search for 850 $\mu$m-selected submillimetre galaxies in the environments of massive ($M_{\star} > 10^{11} M_{\odot}$), radio-quiet ($L_{500 {\rm MHz}} \lesssim 10^{25}$ W Hz$^{-1}$) galaxies at $z \sim 1\text{--}3$ using S2COSMOS data. By constructing number counts in circular regions of radius 1--6 arcmin and comparing with blank-field measurements, we find no significant overdensities of SMGs around massive radio-quiet galaxies at any of these scales, despite being sensitive down to overdensities of $\delta \sim 0.4$. To probe deeper than the catalogue we also examine the distribution of peaks in the SCUBA-2 SNR map, which reveals only tentative signs of any difference in the SMG densities of the radio-quiet galaxy environments compared to the blank field, and only on smaller scales (1$^{\prime}$ radii, corresponding to $\sim0.5$ Mpc) and higher SNR thresholds. We conclude that massive, radio-quiet galaxies at cosmic noon are typically in environments with $\delta\lesssim0.4$, which are either consistent with the blank field or contain only weak overdensities spanning sub-Mpc scales. The contrast between our results and studies of H$z$RGs with similar stellar masses and redshifts implies an intrinsic link between the wide-field environment and radio AGN luminosity at high redshift., Comment: 13 pages, 5 figures. Published in MNRAS

We update the HI surface density measurements for a subset of 17 THINGS galaxies by dealing with the short-spacing problem of the original VLA HI images. It is the same sample that Bigiel et al. (2010) used to study the relation between HI surface densities and star formation rate surface densities in galaxy outer disks, which are beyond the optical radius r25. For ten galaxies, the update is based on combining original THINGS VLA HI images with HI images taken by the single-dish FAST in the FEASTS program. The median increment of HI surface densities in outer disks is 0.15 to 0.4 dex at a given new HI surface density. Several galaxies change significantly in the shape of radial profiles HI surface densities, and seven galaxies are now more than 1-$\sigma$ below the HI size-mass relation. We update the HI star formation laws in outer disks. The median relation between HI surface densities and star formation rate surface densities based on pixelwise measurements shifts downward by around 0.15 dex because the HI surface density values shift rightward, and the scatter increases significantly. The scatter of the relation, indicating the star forming efficiency, exhibits a much stronger positive correlation with the stellar mass surface density than before. Thus, detecting the previously missed, diffuse HI due to short-spacing problem of the VLA observation is important in revealing the true condition and variation of star formation possibly regulated by stellar feedbacks in localized environment of outer disks., Comment: 27 pages, the second paper of "FEASTS Combined with Interferometry". Data is available on the page http://kavli.pku.edu.cn/~jwang/FEASTS_data.html

The James Webb Space Telescope (JWST) observations have revolutionized extragalactic research, particularly with the discovery of little red dots (LRD), which we propose are dust-reddened broad-line active galactic nuclei (AGNs). Their unique v-shape spectral feature observed through JWST/NIRCam challenges us to discern the relative contributions of the galaxy and AGN. We study a spectral energy distribution (SED) model for LRDs from rest-frame UV to infrared bands. We hypothesize that the incident radiation from an AGN, characterized by a typical SED, is embedded in an extended dusty medium with an extinction law similar to those seen in dense regions such as Orion Nebula or certain AGN environments. The UV-optical spectrum is described by dust-attenuated AGN emission, featuring a red optical continuum at $\lambda>4000$ A and a flat UV spectral shape established through a gray extinction curve at $\lambda<3000$ A, due to the absence of small-size grains. There is no need for additional stellar emission or AGN scattered light. In the infrared, the SED is shaped by an extended dust and gas distribution ($\gamma<1$; $\rho\propto r^{-\gamma}$) with a characteristic gas density of $\simeq 10-10^3~{\rm cm}^{-3}$, which allows relatively cool dust temperatures to dominate the radiation, thereby shifting the energy peak from near- to mid-infrared bands. This model, unlike the typical AGN hot torus models, can produce an infrared SED flattening that is consistent with LRD observations through JWST MIRI. Such a density structure can arise from the coexistence of inflows and outflows during the early assembly of galactic nuclei. This might be the reason why LRDs emerge preferentially in the high-redshift universe younger than one billion years., Comment: 14 pages, 6 figures

We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among which Mrk 1239 and PDS 456 are for the first time spectroscopically monitored. We obtain multi-year monitoring data and perform multi-component spectral decomposition to extract the broad H$\beta$ profiles. We detect significant time lags between the H$\beta$ and continuum variations, generally obeying the previously established BLR size-luminosity relation. Velocity-resolved H$\beta$ time lags illustrate diverse, possibly evolving BLR kinematics. We further measure the H$\beta$ line widths from mean and rms spectra and the resulting virial products show good consistency among different seasons. Adopting a unity virial factor and the full width at half maximum of the broad H$\beta$ line from the mean spectrum as the measure of velocity, the obtained black hole mass averaged over seasons is $\log M_\bullet/M_\odot=8.02_{-0.14}^{+0.09}$, $6.92_{-0.12}^{+0.12}$, $8.01_{-0.25}^{+0.16}$, $7.44_{-0.14}^{+0.13}$, and $8.59_{-0.11}^{+0.07}$ for the five objects, respectively. The black hole mass estimations using other line width measures are also reported (up to the virial factors). For objects with previous RM campaigns, our mass estimates are in agreement with earlier results. In a companion paper, we will employ BLR dynamical modeling to directly infer the black hole mass and thereby determine the virial factors., Comment: 32 pages, 6 tables, 20 figures. To appear in ApJ

Dusty circumnuclear disks (CNDs) in luminous early-type galaxies (ETGs) show regular, dynamically cold molecular gas kinematics. For a growing number of ETGs, Atacama Large Millimeter/sub-millimeter Array (ALMA) CO imaging and detailed gas-dynamical modeling facilitate moderate-to-high precision black hole (BH) mass ($M_{BH}$) determinations. From the ALMA archive, we identified a subset of 26 ETGs with estimated $M_{BH}/M_{\odot} \gtrsim 10^8$ to a few $\times$10$^9$ and clean CO kinematics but that previously did not have sufficiently high angular resolution near-IR observations to mitigate dust obscuration when constructing stellar luminosity models. We present new optical and near-IR Hubble Space Telescope (HST) images of this sample to supplement the archival HST data, detailing the sample properties and data analysis techniques. After masking the most apparent dust features, we measure stellar surface brightness profiles and model the luminosities using the multi-Gaussian expansion (MGE) formalism. Some of these MGEs have already been used in CO dynamical modeling efforts to secure quality \mbh\ determinations, and the remaining ETG targets here are expected to significantly improve the high-mass end of the current BH census, facilitating new scrutiny of local BH mass-host galaxy scaling relationships. We also explore stellar isophotal behavior and general dust properties, finding these CNDs generally become optically thick in the near-IR ($A_H \gtrsim 1$ mag). These CNDs are typically well-aligned with the larger-scale stellar photometric axes with a few notable exceptions. Uncertain dust impact on the MGE often dominates the BH mass error budget, so extensions of this work will focus on constraining CND dust attenuation., Comment: 47 pages, 37 Figures, Accepted for publication in ApJ

The properties of the galaxies are tightly connected to their host halo mass and halo assembly history. Accurate measurement of the halo assembly history in observation is challenging but crucial to the understanding of galaxy formation and evolution. The stellar-to-halo mass ratio ($M_*/M_{\mathrm{h}}$) for the centrals has often been used to indicate the halo assembly time $t_{\mathrm{h,50}}$ of the group, where $t_{\mathrm{h,50}}$ is the lookback time at which a halo has assembled half of its present-day virial mass. Using mock data from the semi-analytic models, we find that $M_*/M_{\mathrm{h}}$ shows a significant scatter with $t_{\mathrm{h,50}}$, with a strong systematic difference between the group with a star-forming central (blue group) and passive central (red group). To improve the accuracy, we develop machine-learning models to estimate $t_{\mathrm{h,50}}$ for galaxy groups using only observable quantities in the mocks. Since star-formation quenching will decouple the co-growth of the dark matter and baryon, we train our models separately for blue and red groups. Our models have successfully recovered $t_{\mathrm{h,50}}$, within an accuracy of $\sim$ 1.09 Gyr. With careful calibrations of individual observable quantities in the mocks with SDSS observations, we apply the trained models to the SDSS Yang et al. groups and derive the $t_{\mathrm{h,50}}$ for each group for the first time. The derived SDSS $t_{\mathrm{h,50}}$ distributions are in good agreement with that in the mocks, in particular for blue groups. The derived halo assembly history, together with the halo mass, make an important step forward in studying the halo-galaxy connections in observation., Comment: 18 pages, 7 figures. Accepted by ApJ

Over three decades of reverberation mapping (RM) studies on local broad-line active galactic nuclei (AGNs) have measured reliable black-hole (BH) masses for $> 100$ AGNs. These RM measurements reveal a significant correlation between the Balmer broad-line region size and the AGN optical luminosity (the $R-L$ relation). Recent RM studies for AGN samples with more diverse BH accretion parameters (e.g., mass and Eddington ratio) reveal a substantial intrinsic dispersion around the average $R-L$ relation, suggesting variations in the overall spectral energy distribution shape as functions of accretion parameters. Here we perform a detailed photoionization investigation of expected broad-line properties as functions of accretion parameters, using the latest models for the AGN continuum implemented in {\tt qsosed}. We compare theoretical predictions with observations of a sample of 67 $z\lesssim0.5$ reverberation-mapped AGNs with both rest-frame optical and UV spectra in the moderate-accretion regime (Eddington ratio $\lambda_{\rm Edd}\equiv L/L_{\rm Edd}<0.5$). The UV/optical line strengths and their dependences on accretion parameters can be reasonably well reproduced by the locally-optimally-emitting cloud (LOC) photoionization models. We provide quantitative recipes that use optical/UV line flux ratios to infer the ionizing continuum, which is not directly observable. In addition, photoionization models with universal values of ionization parameter ($\log U_{\rm H}=-2$) and hydrogen density ($\log n({\rm H})=12$) can qualitatively reproduce the observed global $R-L$ relation for the current AGN sample. However, such models fail to reproduce the observed trend of decreasing BLR size with $L/L_{\rm Edd}$ at fixed optical luminosity, which may imply that the gas density increases with the accretion rate., Comment: 23 pages, 15 figures. Submitted to ApJ. Comments welcome!

The circumgalactic medium (CGM) is responsive to kinetic disruptions generated by nearby astrophysical events. In this work, we study the saturation and dissipation of turbulent hydrodynamics within the CGM through an extensive array of 252 numerical simulations with a large parameter space. These simulations are endowed with proper cooling mechanisms to consistently explore the parameter space spanned by the average gas density, metallicity, and turbulence driving strength. A dichotomy emerges in the dynamics dissipation behaviors. Disturbances that are hot and subsonic are characterized by weak compression and slow dissipation, resulting in density fluctuations typically $\lesssim 10^{-2}$. Conversely, warm supersonic turbulence, marked by significant compression shocks and subsequent rapid cooling, is associated with substantial clumping factors $\sim 10^0-10^1$. In the supersonic cases, the kinetic energy decay is divided into a rate-limiting phase of shock dissipation and a comparatively swift phase of thermal dissipation, predominantly occurring within the overdense regions. Upon turbulence driving turnoff, the strong density contrasts decay within a relatively brief timescale of $\sim 30 - 300~{\rm Myr}$, depending on the average gas density. Dense clouds are crushed on similar timescales of $ \sim 30 - 100 ~{\rm Myr} $, depending on turbulence driving strength but independent from average gas density. Results of this work also contribute a novel dataset of dissipation timescales that incorporates an understanding of kinematics and thermodynamics in addition to the traditional cooling rate tables, which may serve as a valuable asset for forthcoming simulations that aim to explore gas dynamics on galactic and cosmological scales., Comment: 24 pages, 20 figures

We investigate the optical variability of low-redshift ($0.15< z\leq0.4$) active galactic nuclei using the multi-epoch data from the Zwicky Transient Facility. We find that a damped random walk model well describes the ensemble structure function in the $g$ band. Consistent with previous studies, more luminous active galactic nuclei tend to have a steeper structure function at a timescale less than the break timescale and smaller variability amplitude. By comparing the structure functions in the optical with the mid-infrared obtained from the Wide-field Infrared Survey Explorer, we derive the size of the dusty torus using a toy model for the geometry of the torus. The size of the torus positively correlates with the luminosity of the active nucleus, following a relation that agrees well with previous studies based on reverberation mapping. This result demonstrates that the structure function method can be used as a powerful and highly efficient tool to examine the size of the torus., Comment: 7 pages, 5 figures. Accepted for publication in A&A

We present confusion-limited SCUBA-2 450-$\mu$m observations in the COSMOS-CANDELS region as part of the JCMT Large Program, SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). Our maps at 450 and 850 $\mu$m cover an area of 450 arcmin$^2$. We achieved instrumental noise levels of $\sigma_{\mathrm{450}}=$ 0.59 mJy beam$^{-1}$ and $\sigma_{\mathrm{850}}=$ 0.09 mJy beam$^{-1}$ in the deepest area of each map. The corresponding confusion noise levels are estimated to be 0.65 and 0.36 mJy beam$^{-1}$. Above the 4 (3.5) $\sigma$ threshold, we detected 360 (479) sources at 450 $\mu$m and 237 (314) sources at 850 $\mu$m. We derive the deepest blank-field number counts at 450 $\mu$m, covering the flux-density range of 2 to 43 mJy. These are in agreement with other SCUBA-2 blank-field and lensing-cluster observations, but are lower than various model counts. We compare the counts with those in other fields and find that the field-to-field variance observed at 450 $\mu$m at the $R=6^\prime$ scale is consistent with Poisson noise, so there is no evidence of strong 2-D clustering at this scale. Additionally, we derive the integrated surface brightness at 450 $\mu$m down to 2.1 mJy to be $57.3^{+1.0}_{-6.2}$~Jy deg$^{-2}$, contributing to (41$\pm$4)\% of the 450-$\mu$m extragalactic background light (EBL) measured by COBE and Planck. Our results suggest that the 450-$\mu$m EBL may be fully resolved at $0.08^{+0.09}_{-0.08}$~mJy, which extremely deep lensing-cluster observations and next-generation submillimeter instruments with large aperture sizes may be able to achieve., Comment: 29 pages, 14 figures, accepted for publication in ApJ

To investigate the short-term variability and determine the size of the optical continuum emitting region of intermediate-mass black holes (IMBHs), we carried out high-cadence, multiband photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Paloma spectroscopic observations, 23 years after the Sloan Digital Sky Survey, ruling out the supernovae origin of the broad Ha line. The photometric experiment was primarily conducted utilizing four-channel imagers MuSCAT 3 and 4 mounted on 2 m telescopes within the Las Cumbres Observatory Global Telescope Network. Despite the expectation of variability, we observed no significant variation (<1.4%) on timescales of 6 to 10 hr. This nondetection is likely due to substantial host galaxy light diluting the subtle active galactic nucleus (AGN) variability. Additionally, we cannot rule out that the target was in a relatively quiescent state without intranight variability during our monitoring, owing to the stochastic nature of AGN variations. To enhance the possibility of detecting subtle variability signals and lag in future IMBH reverberation campaigns, it may be beneficial to select targets with a higher AGN-to-host flux ratio, and conduct dual-band preliminary tests and tailored simulations., Comment: 11 pages, 6 figures

It is believed that dual active galactic nuclei (dual AGN) will form during galaxies merge. Studying dual-AGN emission can provide valuable insights into galaxy merging and evolution. To investigate parsec-scale radio emission properties, we observed eight radio components of four selected dual-AGN systems using the Very Long Baseline Array (VLBA) at 5 GHz in multiple-phase-center mode. Among them, two compact radio components, labeled J0051+0020B and J2300-0005A, were detected clearly on parsec scales for the first time. However, the radio emission of the other six components was resolved out in the high-resolution images. We provided the values or upper limits of the brightness temperature and radio emission power, and analyzed the emission origins in detail for each target. Based on their physical properties reported in this work and in the literature, we suggest the radio emission in J0051+0020B and J2300-0005A originates primarily from compact jets, while the other six sources show more complex emission mechanisms. In addition, our VLBA observations suggest the systematic X-ray deficit in our dual-AGN sample is likely attributed to the tidally induced effect and possible viewing angle effect., Comment: 15 pages, 4 figures

This study delves into the bulge-disc components and stellar mass distribution in the fast-rotating, highly massive spiral galaxy 2MASX~J23453268-0449256, distinguished by extraordinary radio jets extending to Mpc scales. Using high-resolution multi-wavelength Hubble Space Telescope (HST) observations and multi-parameter panchromatic spectral energy distribution (SED) fitting, we derive estimates for the star formation rate, total baryonic mass in stars, and warm dust properties. Our findings, validated at a spatial resolution of approximately 100 pc, reveal a pseudo-bulge rather than a classical bulge and a small nuclear bar and resonant ring, challenging conventional models of galaxy formation. Additionally, the lack of tidal debris and the highly symmetric spiral arms within a rotationally supported stellar disc indicate a tranquil coevolution of the galactic disc and its supermassive black hole (SMBH). Significantly, the galaxy exhibits suppressed star formation in its center, potentially influenced by feedback from the central accreting SMBH with powerful radio jets. Detailed multi-wavelength studies of potential star-forming gas disclose that, while hot X-ray gas cools down in the galaxy's halo, new stars do not form in the center, likely due to feedback effects. This study raises questions about the efficient fueling and sustained collimated jet ejection activity in J2345-0449, underscoring the imperative need for a comprehensive understanding of its central black hole engine properties, which are presently lacking. The exceptional rarity of galaxies like 2MASX~J23453268-0449256 presents intriguing challenges in unraveling the physical processes responsible for their unique characteristics., Comment: Submitted to MNRAS

Outflowing stars impinging upon ambient gas experience accelerations due to the gravitational feedback from the interaction morphology between the outflow and the ambient gas. Such ``negative dynamical friction'' (NDF), in contrast to the conventional ``dynamical friction'' (DF), is studied for its impact on the dynamics of open clusters (OCs) immersed in a uniform ambient gas. We modify the $N$-body integration code \texttt{REBOUND} with both NDF and DF implemented according to the outflow conditions of each star in a consistently constructed OC. The evolution of stars is also involved in determining the gas-star interactions throughout their stellar lives. Compared to DF-only and gas-free models with identical initial conditions, the NDF-affected cluster is puffier and evaporates faster, as indicated by various diagnostics, including lower velocity dispersions and larger half-mass and half-light radii. Neutron stars with fast winds are expelled from the cluster due to their intensive NDF effect, even without the ``kicks'' by asymmetric supernovae. Exploration of parameter space confirms that the NDF effect is generally enhanced with higher ambient gas densities, in qualitative agreement with the expression of acceleration. Outflow-ambient interactions should be considered for the proper interpretation of the stellar dynamics evolution in clusters., Comment: 20 pages, 18 figures

From the IllustrisTNG-50 simulation, a sample of 836 central disk galaxies with tiny stellar halos is chosen to study the inherent evolution of galaxies driven by nature. These galaxies are classified as compact, normal, or extended by referencing their locations on the mass-size ($M_\star-R_{\rm 1/2}$) diagram. This research demonstrates the distinctive evolutionary pathways of galaxies with different sizes in IllustrisTNG simulations, primarily driven by nature. It is confirmed that disk galaxies inherit the angular momentum of their parent dark matter halos. More compact galaxies form earlier within halos possessing lower specific angular momentum through heightened star formation during the early phase at redshifts above 2. During the later phase, the size of extended galaxies experiences more pronounced growth by accreting gas with high angular momentum. Additionally, we reveal that many key characteristics of galaxies are linked to their mass and size: (1) compact galaxies tend to exhibit higher metal content, proportional to the potential well $\frac{M_\star}{R_{\rm 1/2}}$, (2) compact galaxies host more massive bulges and black holes, and higher central concentration. Furthermore, our analysis indicates that galaxies of all types continue to actively engage in star formation, with no evident signs of quenching attributed to their varying sizes and angular momenta., Comment: 11 pages, 11 figures. Accepted for publication in A&A

We present a statistical study of the properties of diffuse HI in ten nearby galaxies, comparing the HI detected by the single-dish telescope FAST (FEASTS program) and the interferometer VLA (THINGS program), respectively. The THINGS' observation missed HI with a median of 23% due to the short-spacing problem of interferometry and limited sensitivity. We extract the diffuse HI by subtracting the dense HI, which is obtained from the THINGS data with a uniform flux-density threshold, from the total HI detected by FAST. Among the sample, the median diffuse-HI fraction is 34%, and more diffuse HI is found in galaxies exhibiting more prominent tidal-interaction signatures. The diffuse HI we detected seems to be distributed in disk-like layers within a typical thickness of $1\,\text{kpc}$, different from the more halo-like diffuse HI detected around NGC 4631 in a previous study. Most of the diffuse HI is cospatial with the dense HI and has a typical column density of $10^{17.7}$-$10^{20.1}\,\text{cm}^{-2}$. The diffuse and dense HI exhibits a similar rotational motion, but the former lags by a median of 25% in at least the inner disks, and its velocity dispersions are typically twice as high. Based on a simplified estimation of circum-galactic medium properties and assuming pressure equilibrium, the volume density of diffuse HI appears to be constant within each individual galaxy, implying its role as a cooling interface. Comparing with existing models, these results are consistent with a possible link between tidal interactions, the formation of diffuse HI, and gas accretion., Comment: 45 pages, 23 figures. In press at ApJ. Data will be released at the FEASTS site upon publication

We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$\beta$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of the 14 objects. These results include H$\beta$ time lags, supermassive black hole masses, and velocity-resolved time lags. The velocity-resolved lags allow us to investigate the kinematics of the broad-line region., Comment: 35 pages, 19 figures, accepted for publication in ApJ Supplement

The supermassive black holes ($M_{\rm BH} \sim 10^{6}$$-$$10^{10}~M_\odot$) that power luminous active galactic nuclei (AGNs), i.e., quasars, generally show a correlation between thermal disk emission in the ultraviolet (UV) and coronal emission in hard X-rays. In contrast, some "massive" black holes (mBHs; $M_{\rm BH} \sim 10^{5}$$-$$10^{6}~M_\odot$) in low-mass galaxies present curious X-ray properties with coronal radiative output up to 100$\times$ weaker than expected. To examine this issue, we present a pilot study incorporating Very Large Array radio observations of a sample of 18 high-accretion-rate (Eddington ratios $L_{\rm bol}/L_{\rm Edd} > 0.1$), mBH-powered AGNs ($M_{\rm BH} \sim 10^{6}~M_\odot$) with Chandra X-ray coverage. Empirical correlations previously revealed in samples of radio-quiet, high-Eddington AGNs indicate that the radio$-$X-ray luminosity ratio, $L_{\rm R}/L_{\rm X}$, is approximately constant. Through multiwavelength analysis, we instead find that the X-ray-weaker mBHs in our sample tend toward larger values of $L_{\rm R}/L_{\rm X}$ even though they remain radio-quiet per their optical$-$UV properties. This trend results in a tentative but highly intriguing correlation between $L_{\rm R}/L_{\rm X}$ and X-ray weakness, which we argue is consistent with a scenario in which X-rays may be preferentially obscured from our line of sight by a "slim" accretion disk. We compare this observation to weak emission-line quasars (AGNs with exceptionally weak broad-line emission and a significant X-ray-weak fraction) and conclude by suggesting that our results may offer a new observational signature for finding high-accretion-rate AGNs., Comment: Accepted for publication in ApJ. 26 pages, 8 figures, 6 tables

Utilizing a kinematic decomposition of simulated galaxies, we focus on galaxies with tiny kinematically inferred stellar halos, indicative of weak external influences. We investigate the intricate interplay between internal (natural) and external (nurture) processes in shaping the scaling relationships of specific angular momentum ($j_\star$), stellar mass ($M_\star$), and size of disk galaxies within the IllustrisTNG simulation. The correlation among mass, size, and angular momentum of galaxies is examined by comparing simulations with observations and the theoretical predictions of the exponential hypothesis. Galaxies with tiny stellar halos exhibit a large scatter in the $j_\star$-$M_\star$ relation, which suggests that it is inherently present in their initial conditions. The analysis reveals that the disks of these galaxies adhere to the exponential hypothesis, resulting in a tight fiducial $j_\star$-$M_\star$-scale length (size) relation that is qualitatively consistent with observations. The inherent scatter in $j_\star$ provides a robust explanation for the mass-size relation and its substantial variability. Notably, galaxies that are moderately influenced by external processes closely adhere to a scaling relation akin to that of galaxies with tiny stellar halos. This result underscores the dominant role of internal processes in shaping the overall $j_\star$-$M_\star$ and mass-size relation, with external effects playing a relatively minor role in disk galaxies. Furthermore, the correlation between galaxy size and the virial radius of the dark matter halo exists but fails to provide strong evidence of the connection between galaxies and their parent dark matter halos., Comment: 11 pages, 9 figures, accepted for publication in A&A

The formation of jets in black hole accretion systems is a long-standing problem. It has been proposed that a jet can be formed by extracting the rotation energy of the black hole ("BZ-jet") or the accretion flow ("disk-jet"). While both models can produce collimated relativistic outflows, neither has successfully explained the observed jet morphology. By employing general relativistic magnetohydrodynamic simulations, and considering nonthermal electrons accelerated by magnetic reconnection that is likely driven by magnetic eruption in the underlying accretion flow, we obtain images by radiative transfer calculations and compared them to millimeter observations of the jet in M87. We find that the BZ-jet originating from a magnetically arrested disk around a high-spin black hole can well reproduce the jet morphology, including its width and limb-brightening feature., Comment: 46 pages, 20 figures, 3 tables, published in Science Advances on 22 Mar 2024. arXiv admin note: substantial text overlap with arXiv:2206.05661

We present a series of studies on the connection between galaxy morphology and the structure of host dark-matter (DM) haloes using cosmological simulations. In this work, we introduce a new kinematic decomposition scheme that features physical identification of morphological components, enabling robust separation of thin and thick discs; and measure a wide range of halo properties, including their locations in the cosmic web, internal structures, and assembly histories. Our analysis of the TNG50 simulation reveals that the orbital-circularity threshold for disc differentiation varies across galaxies, with systematic trends in mass and redshift, so the widely used decomposition method with constant circularity cuts is oversimplified and underestimates thin disc at JWST redshifts. The energy threshold between the stellar halo and the inner galaxy is also a function of mass and redshift, minimizing at the sub-Galactic halo mass, where the circularity threshold peaks. Revisiting the issue of galaxy size predictor, we show that disc sizes in TNG50 exhibit correlations with three structural parameters besides virial mass and redshift: 1) a positive correlation with halo spin $\lambda$ across redshifts -- stronger than previously reported for zoom-in simulations but still weaker than the simple $r_{1/2}/R_{\rm vir} \propto \lambda$ scaling; 2) an anti-correlation with DM concentration $c$ that is well described by $r_{1/2}/R_{\rm vir} \propto c^{-0.7}$ even when $c$ is measured in the DM only run; 3) more actively accreting haloes having slightly larger discs, as well as more significant stellar haloes and lower thin-to-thick ratio. Disc mass fraction is higher in rounder haloes and in cosmic knots and filaments, implying that disc development needs both stable halo conditions and continuous material supply. Our methodology is public and adaptable to other simulations., Comment: 20 pages, 17 figures, submitted to MNRAS

The Galaxy inner parts are generally considered to be optically symmetric, as well as kinematically symmetric for most massive early-type galaxies. At the lower-mass end, many galaxies contain lots of small patches in their velocity maps, causing their kinematics to be nonsmooth in small scales and far from symmetry. These small patches can easily be mistaken for measurement uncertainties and have not been well discussed. We used the comparison of observations and numerical simulations to demonstrate the small patches existence beyond uncertainties. For the first time we have found that the fluctuation degrees have an approximate inverse loglinear relation with the galaxy stellar surface mass densities. This tight relation among galaxies that do not show obvious optical asymmetry that traces environmental perturbations indicates that stellar motion in galaxies has inherent asymmetry besides external environment influences. The degree of the kinetic asymmetry is closely related to and constrained by the intrinsic properties of the host galaxy., Comment: 13 pages,10 figures, accepted for publication in ApJL

We perform a comprehensive investigation of the PSF smoothing effect on the measurement of concentration-related parameters ($C$, Gini, $M_{20}$) of high redshift galaxies in the HST and JWST surveys. Our sample contains massive galaxies from the CANDELS/EGS survey (0 < z < 2), and the CEERS survey (1 < z < 3). The non-parametric concentration-related parameters ($C$, Gini, $M_{20}$) and the model-dependent parameters (n, Re) of these galaxies are derived from Statmorph and GALFIT, respectively. We try to evaluate the PSF smoothing effect by comparing the concentration-related parameters to the S\'{e}rsic index in both observations and mock images. We find that the concentration index is generally underestimated, especially for smaller galaxies with higher S\'{e}rsic index (eventually converging to the concentration index of the PSF). However, galaxies with lower S\'{e}rsic index ($n \leq 1$) or larger relative size are less affected by the PSF smoothing effect. The Gini coefficient and the absolute $M_{20}$ statistic also show similar behaviour as the concentration index. Caution should be taken for the possible correction of the concentration-related parameters, where both the relative size and the S\'{e}rsic index of the galaxy are important. Compared to the HST images, the PSF smoothing is much less severe for images in the CEERS survey due to the much higher spatial resolution. In fact, it is better to use the S\'{e}rsic index rather than the non-parametric morphology indicators to trace the light concentration for galaxies at high redshifts. From the single S\'{e}rsic modelling of the HST and JWST images, we also confirm that galaxies at higher redshift are more compact with smaller $R_e$. The lower mass galaxies are more disc-like ($n\sim1$) compared to the higher mass galaxies that are more spheroid dominated ($n\sim3$)., Comment: Accepted by A&A, 21 pages, 20 figures. Comments are welcome

We present molecular gas-dynamical mass measurements of the central black holes in the giant elliptical galaxies NGC 4786 and NGC 5193, based on CO(2$-$1) observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope near-infrared imaging. The central region in each galaxy contains a circumnuclear disk that exhibits orderly rotation with projected line-of-sight velocities of ${\sim} 270\, \mathrm{km}\,\mathrm{s^{-1}}$. We build gas-dynamical models for the rotating disk in each galaxy and fit them directly to the ALMA data cubes. At $0.31^{\prime \prime}$resolution, the ALMA observations do not fully resolve the black hole sphere of influence (SOI), and neither galaxy exhibits a central rise in rotation speed, indicating that emission from deep within the SOI is not detected. As a result, our models do not tightly constrain the central black hole mass in either galaxy, but they prefer the presence of a central massive object in both galaxies. We measure the black hole mass to be $(M_{\mathrm{BH}}/10^8\, M_{\odot}) = 5.0 \pm 0.2 \,[\mathrm{1\sigma \,statistical}] \,^{+1.4}_{-1.3} \,[\mathrm{systematic}]$ in NGC 4786 and $(M_{\mathrm{BH}}/10^8\, M_{\odot}) = 1.4 \pm 0.03 \, [\mathrm{1\sigma\,statistical}] ^{+1.5}_{-0.1} \,[\mathrm{systematic}]$ in NGC 5193. The largest component of each measurement's error budget is from the systematic uncertainty associated with the extinction correction in the host galaxy models. This underscores the importance of assessing the impact of dust attenuation on the inferred $M_{\mathrm{BH}}$., Comment: 24 pages, 8 figures. Accepted for publication in ApJ

The fundamental plane of black hole activity (BHFP) describes the correlation between radio luminosity ($L_R$), X-ray luminosity ($L_X$), and black hole mass. It reflects a disc-jet connection. However, dependence of BHFP on various physical properties of active galactic nuclei (AGNs) and host galaxies is unclear, especially for low-luminosity AGNs, which is important for understanding accretion physics in AGNs. Here we explore the dependence of BHFP on radio loudness, Eddington ratio ($\lambda_E$), redshift, and galaxy star formation properties at 0.1 < z < 4 for radio AGNs. Based on current deep and large surveys, our studies extend to low luminosities and high redshifts. From GOODS and COSMOS fields, we constructed a homogeneous radio AGN sample including 208 objects. We divided the radio AGN sample into 141 radio-quiet (RQ) AGNs and 67 radio-loud (RL) AGNs according to the radio loudness defined by the $L_R$-to-$L_X$ ratio. This ratio shows a bimodal distribution that is well described by two Gaussian models. Their cross point corresponds to a radio loudness threshold of $\log (L_R/L_X) = -2.73$. RQ and RL AGNs show different BHFPs, indicating a dependence of BHFP on the radio loudness. For both RQ and RL AGNs, BHFP shows a dependence on $\lambda_E$, but no dependence on redshift. BHFP shows a dependence on the galaxy star formation properties for RQ AGNs, while for RL AGNs this dependence disappears. BHFP sheds important light on the accretion physics and X-ray emission origins. X-ray emission of RQ AGNs originates from a combination of advection-dominated accretion flow (ADAF) and synchrotron radiation from jet at $0.01 < \lambda_E < 0.1$, and from jet at $0.1 < \lambda_E < 1$. X-ray emission in RL AGNs originates from a combination of ADAF and jet at $\lambda_E < 0.01$, from jet at $0.01 < \lambda_E < 0.1$, and from a combination of standard thin disc and jet at $\lambda_E > 0.1$., Comment: Accepted by A&A (in press), 13 pages, 6 figures, 2 tables; abstract slightly abridged

Resolved stellar morphology of z>1 galaxies was inaccessible before JWST. This limitation, due to the impact of dust on rest-frame UV light, had withheld major observational conclusions required to understand the importance of clumps in galaxy evolution. Essentially independent of this issue, we use the rest-frame near-IR for a stellar-mass dependent clump detection method and determine reliable estimations of selection effects. We exploit publicly available JWST/NIRCam and HST/ACS imaging data from CEERS, to create a stellar-mass based picture of clumps in a mass-complete sample of 418 galaxies within a wide wavelength coverage of 0.5-4.6${\mu}$m and a redshift window of 1

The inclination angle of substructures in active galaxies gives insights into physical components from scales of the vicinity of the central black hole to the entire host galaxy. We use the self-consistent reflection spectral model \textsc{RELXILL} to measure the inclination of the inner region of accretion disks with broadband ($0.3-78\,\rm keV$) X-ray observations, systematically studying the reliability of this methodology. To test the capability of the model to return statistically consistent results, we analyze multi-epoch, joint XMM-Newton and NuSTAR data of the narrow-line Seyfert~1 galaxy I\,Zwicky\,1 and the broad-line radio galaxy 3C\,382, which exhibit different degrees of spectral complexity and reflection features. As expected, we find that adding more data for analysis narrows the confidence interval and that multi-epoch, joint observations return optimal measurements; however, even single-epoch data can be well-fitted if the reflection component is sufficiently dominant. Mock spectra are used to test the capability of \textsc{RELXILL} to recover input parameters from typical single-epoch, joint observations. We find that inclination is well-recovered at 90\% confidence, with improved constraints at higher reflection fraction and higher inclination. Higher iron abundance and corona temperature tighten the constraints as well, but the effect is not as significant as a higher reflection fraction. The spin, however, have little effect in reflection-based inclination measurements. We conclude that broadband reflection spectroscopy can reliably measure inner accretion disk inclination., Comment: 27 pages, 11 figures, 11 tables; accepted by the Astrophysical Journal

In this work, we constrain the star-forming properties of all possible sites of incipient high-mass star formation in the Milky Way's Galactic Center. We identify dense structures using the CMZoom 1.3mm dust continuum catalog of objects with typical radii of $\sim$0.1pc, and measure their association with tracers of high-mass star formation. We incorporate compact emission at 8, 21, 24, 25, and 70um from MSX, Spitzer, Herschel, and SOFIA, catalogued young stellar objects, and water and methanol masers to characterize each source. We find an incipient star formation rate (SFR) for the CMZ of ~0.08 Msun yr^{-1} over the next few 10^5 yr. We calculate upper and lower limits on the CMZ's incipient SFR of ~0.45 Msun yr^{-1} and ~0.05 Msun yr^{-1} respectively, spanning between roughly equal to and several times greater than other estimates of CMZ's recent SFR. Despite substantial uncertainties, our results suggest the incipient SFR in the CMZ may be higher than previously estimated. We find that the prevalence of star formation tracers does not correlate with source volume density, but instead ~75% of high-mass star formation is found in regions above a column density ratio (N_{SMA}/N_{Herschel}) of ~1.5. Finally, we highlight the detection of ``atoll sources'', a reoccurring morphology of cold dust encircling evolved infrared sources, possibly representing HII regions in the process of destroying their envelopes., Comment: Accepted for publication in ApJ

We present a systematic study of the putative hot gas corona around late-type galaxies (LTGs) residing in the Virgo cluster, based on archival Chandra observations. Our sample consists of 21 nearly edge-on galaxies representing a star formation rate (SFR) range of ($0.2-3\rm~M_\odot~yr^{-1}$) a stellar mass ($M_*$) range of $(0.2-10) \times 10^{10}\rm~M_{\odot}$, the majority of which have not been explored with high-sensitivity X-ray observations so far. Significant extraplanar diffuse X-ray (0.5-2 keV) emission is detected in only three LTGs, which are also the three galaxies with the highest SFR. A stacking analysis is performed for the remaining galaxies without individual detection, dividing the whole sample into two subsets based on SFR, stellar mass, or specific SFR. Only the high-SFR bin yields a significant detection, which has a value of $L\rm_X \sim3\times10^{38}\rm~erg~s^{-1}$ per galaxy. The stacked extraplanar X-ray signals of the Virgo LTGs are consistent with the empirical $L\rm_X - SFR$ and $L\rm_X - M_*$ relations found among highly inclined disk galaxies in the field, but appear to be systematically lower than that of a comparison sample of simulated cluster star-formation galaxies identified from the Illustris-TNG100 simulation. The apparent paucity of hot gas coronae in the sampled Virgo LTGs might be understood as the net outcome of the long-lasting effect of ram pressure stripping exerted by the hot intra-cluster medium and in-disk star-forming activity acting on shorter timescales. A better understanding of the roles of environmental effects in regulating the hot gas content of cluster galaxies invites sensitive X-ray observations for a large galaxy sample., Comment: 18 pages, 7 figures. Accepted for publication in ApJ. Comments welcome

The nucleus of almost all massive galaxies contains a supermassive black hole (BH). The feedback from the accretion of these BHs is often considered to have crucial roles in establishing the quiescence of massive galaxies, although some recent studies show that even galaxies hosting the most active BHs do not exhibit a reduction in their molecular gas reservoirs or star formation rates. Therefore, the influence of BHs on galaxy star formation remains highly debated and lacks direct evidence. Here, based on a large sample of nearby galaxies with measurements of masses of both BHs and atomic hydrogen (HI), the main component of the interstellar medium, we show that the HI gas mass to stellar masses ratio ($\mu_{\rm HI} = M_{\rm HI}/M_{\star}$) is more strongly correlated with BH masses ($M_{\rm BH}$) than with any other galaxy parameters, including stellar mass, stellar mass surface density and bulge masses. Moreover, once the $\mu_{\rm HI}-M_{\rm BH}$ correlation is considered, $\mu_{\rm HI}$ loses dependence on other galactic parameters, demonstrating that $M_{\rm BH}$ serves as the primary driver of $\mu_{\rm HI}$. These findings provide important evidence for how the accumulated energy from BH accretion regulates the cool gas content in galaxies, by ejecting interstellar medium gas and/or suppressing gas cooling from the circumgalactic medium., Comment: Published in Nature. Updated to match the accepted version

Understanding the methodological robustness in identifying and quantifying high-redshift bars is essential for studying their evolution with the {\it James} {\it Webb} Space Telescope (JWST). We used nearby spiral galaxies to generate simulated images at various resolutions and signal-to-noise ratios, and obtained the simulated galaxy images observed in the Cosmic Evolution Early Release Science (CEERS) survey from Yu et al. Through a comparison of measurements before and after image degradation, we show that the bar measurements for massive galaxies remain robust against noise. While the measurement of the bar position angle remains unaffected by resolution, the measured bar ellipticity is significantly underestimated in low-resolution images. The size measurement is on average barely affected as long as the intrinsic bar size $a_{\rm bar,\,true}>2\times{\rm FWHM}$. To address these effects, correction functions are derived. We also find that the effectiveness of detecting bars remains at $\sim$\,100\% when the $a_{\rm bar,\,true}/{\rm FWHM}$ is above 2, below which the rate drops sharply, quantitatively validating the effectiveness of using $a_{\rm bar,\,true}>2\times {\rm FWHM}$ as a bar detection threshold. We analyze a set of simulated CEERS images, which take into account observational effects and plausible galaxy (and bar-size) evolution models, and show that a significant (and misleading) reduction in detected bar fraction with increasing redshift would apparently result even if the true bar fraction remained constant. Our results underscore the importance of disentangling the true bar fraction evolution from resolution effects and bar size growth., Comment: 17 pages, 15 figures; accepted for publication in Astronomy & Astrophysics (A&A)

We present 0$.\!\!^{\prime\prime}22$-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2$-$1) emission from the circumnuclear gas disk in the red nugget relic galaxy PGC 11179. The disk shows regular rotation, with projected velocities near the center of 400 km s$^{-1}$. We assume the CO emission originates from a dynamically cold, thin disk and fit gas-dynamical models directly to the ALMA data. In addition, we explore systematic uncertainties by testing the impacts of various model assumptions on our results. The supermassive black hole (BH) mass ($M_\mathrm{BH}$) is measured to be $M_\mathrm{BH} = (1.91\pm0.04$ [$1\sigma$ statistical] $^{+0.11}_{-0.51}$ [systematic])$\times 10^9$ $M_\odot$, and the $H$-band stellar mass-to-light ratio $M/L_H=1.620\pm0.004$ [$1\sigma$ statistical] $^{+0.211}_{-0.107}$ [systematic] $M_\odot/L_\odot$. This $M_\mathrm{BH}$ is consistent with the BH mass$-$stellar velocity dispersion relation but over-massive compared to the BH mass$-$bulge luminosity relation by a factor of 3.7. PGC 11179 is part of a sample of local compact early-type galaxies that are plausible relics of $z\sim2$ red nuggets, and its behavior relative to the scaling relations echoes that of three relic galaxy BHs previously measured with stellar dynamics. These over-massive BHs could suggest BHs gain most of their mass before their host galaxies do. However, our results could also be explained by greater intrinsic scatter at the high-mass end of the scaling relations, or by systematic differences in gas- and stellar-dynamical methods. Additional $M_\mathrm{BH}$ measurements in the sample, including independent cross-checks between molecular gas- and stellar-dynamical methods, will advance our understanding of the co-evolution of BHs and their host galaxies., Comment: 18 pages, 8 figures, 2 tables. Accepted in ApJ

The assembly of galaxies over cosmic time is tightly connected to the assembly of their host dark matter halos. We investigate the stellar mass growth history and the chemical enrichment history of central galaxies in SDSS-MaNGA. We find that the derived stellar metallicity of passive central galaxies is always higher than that of the star-forming ones. This stellar metallicity enhancement becomes progressively larger towards low-mass galaxies (at a given epoch) and earlier epochs (at a given stellar mass), which suggests strangulation as the primary mechanism for star formation quenching in central galaxies not only in the local universe, but also very likely at higher redshifts up to $z\sim3$. We show that at the same present-day stellar mass, passive central galaxies assembled half of their final stellar mass $\sim 2$ Gyr earlier than star-forming central galaxies, which agrees well with semi-analytic model. Exploring semi-analytic model, we find that this is because passive central galaxies reside in, on average, more massive halos with a higher halo mass increase rate across cosmic time. As a consequence, passive central galaxies are assembled faster and also quenched earlier than their star-forming counterparts. While at the same present-day halo mass, different halo assembly history also produces very different final stellar mass of the central galaxy within, and halos assembled earlier host more massive centrals with a higher quenched fraction, in particular around the "golden halo mass" at $10^{12}\mathrm{M_\odot}$. Our results call attention back to the dark matter halo as a key driver of galaxy evolution., Comment: 19 pages, 11 figures. Accepted by ApJ