Astrophysics of Galaxies (astro-ph.GA)

Abstract: We explore properties of stellar kinematics and ionized gas in a sample of 1106 local [U]LIRGs from the AKARI telescope. We combine data from $Wide-field\ Infrared\ Survey\ Explorer$ (WISE) and Sloan Digital Sky Survey (SDSS) Data Release 13 (DR13) to fit the spectral energy distribution (SED) of each source to constrain the contribution of AGN to the total IR luminosity and estimate physical parameters such as stellar mass and star-formation rate (SFR). We split our sample into AGNs and weak/non-AGNs. We find that our sample is considerably above the main sequence. The highest SFRs and stellar masses are associated with ULIRGs. We also fit the H$\beta$ and H$\alpha$ regions to characterize the outflows. We find that the incidence of ionized gas outflows in AGN [U]LIRGs ($\sim$ 72\%) is much higher than that in weak/non-AGN ones ($\sim$ 39\%). The AGN ULIRGs have extreme outflow velocities (up to $\sim$ 2300 km s$^{-1}$) and high mass outflow rates (up to $\sim$ 60 \solarm~yr$^{-1}$). Our results suggest that starbursts are insufficient to produce such powerful outflows. We explore the correlations of SFR and specific SFR (sSFR) with ionized gas outflows. We find that AGN hosts with the highest SFRs exhibit a negative correlation between outflow velocity and sSFR. Therefore, in AGNs containing large amounts of gas, the negative feedback scenario might be suggested.

Abstract: The hydrogen Balmer decrement (e.g., $\rm H\alpha/H\beta$) is widely adopted as an indicator of the internal reddening of active galactic nuclei (AGNs). This is challenged by some low-luminosity AGNs (LLAGNs) and changing-look AGNs (CLAGNs), which have steep Balmer decrement but without strong evidence for absorption. We compile a sample of normal AGNs and CLAGNs with a wider distribution of bolometric Eddington ratio ($\lambda_{\rm Edd}=L_{\rm bol}/L_{\rm Edd}$) and find a strong negative correlation between $\rm H\alpha/H\beta$ and $\lambda_{\rm Edd}$, which suggests that the Balmer decrement is also accretion-rate dependent. We further explore the Balmer decrement based on the photoionization model using the Cloudy code by considering spectral energy distribution (SED) from the accretion disk with different accretion rates (e.g., disk/corona and truncated disk at high and low Eddington ratios, respectively). Both the standard disk and truncated disk predict a negative correlation of $\rm H\alpha/H\beta-\lambda_{\rm Edd}$, where the relation is steeper in the case of the truncated disk. The negative correlations are also explored in two single CLAGNs. The measured negative correlation of $\rm H\alpha/H\beta$ -- $\lambda_{\rm Edd}$ is mainly caused by the lower responsivity $({\rm dlog}L_{\rm line}/{\rm dlog}L_{\rm cont})$ in $\rm H\alpha$ relative to that in $\rm H\beta$, due to the larger optical depth in the former. We propose that the steep Balmer decrements in low-Eddington-ratio AGNs (e.g., some Seyferts 1.5-1.9 and CLAGNs) are not simply caused by absorption but mainly caused by the relatively low flux of ionizing photons.

Abstract: The feedback from accretion of central supermassive black holes (SMBHs) is a hot topic in the co-evolution of the SMBHs and their host galaxies. By tracing the large scale outflow by the line profile and bulk velocity shift of $[ \rm O~{\scriptsize III}]~ \lambda 5007$, the evolutionary role of outflow is studied here on a large sample of 221 type 2 quasars (QSO2s) extracted from Reyes et al. By following our previous study on local Seyfert 2 galaxies, the current spectral analysis on the SDSS spectroscopic database enables us to arrive at following results: (1) by using the Lick indices, we confirm that QSO2s are on average associated with younger stellar populations than Seyfert galaxies; (2) QSO2s with a stronger outflow are tend to be associated with a younger stellar population, which implies a coevolution between the feedback from SMBH and the host in QSO2s; (3) although an occupation at the high $L_{\rm bol}/L_{\rm Edd}$ end, the QSO2s follow the $L_{\rm bol}/L_{\rm Edd}$-$D_{n}(4000)$ sequence established from local, less-luminous Seyfert galaxies, which suggests a decrease of accretion activity of SMBH and feedback as the circumnuclear stellar population continuously ages.

Abstract: In the context of a project aimed at characterizing the properties of star clusters in the Galactic bulge, here we present the determination of the internal kinematics and structure of the massive globular cluster NGC 6569. The kinematics has been studied by means of an unprecedented spectroscopic dataset acquired in the context of the ESO-VLT Multi-Instrument Kinematic Survey (MIKiS) of Galactic globular clusters, combining the observations from four different spectrographs. We measured the line-of-sight velocity of a sample of almost 1300 stars distributed between ~0.8" and 770" from the cluster center. From a sub-sample of high-quality measures, we determined the velocity dispersion profile of the system over its entire radial extension (from ~ 5" to ~ 200" from the center), finding the characteristic behavior usually observed in globular clusters, with a constant inner plateau and a declining trend at larger radii. The projected density profile of the cluster has been obtained from resolved star counts, by combining high-resolution photometric data in the center, and the Gaia EDR3 catalog radially extended out to ~20' for a proper sampling of the Galactic field background. The two profiles are properly reproduced by the same King model, from which we estimated updated values of the central velocity dispersion, main structural parameters (such as the King concentration, the core, half-mass, and tidal radii), total mass, and relaxation times. Our analysis also reveals a hint of ordered rotation in an intermediate region of the cluster (40"<r<90", corresponding to $ 2 r_c<r<4.5 r_c$), but additional data are required to properly assess this possibility.

Abstract: In this paper we aim to constrain for the first time the dust emission in the mid-to-far infrared domain, in the LMC, with the use of the Spitzer IRS and MIPS SED data, combined with Herschel data. We also consider UV extinction predictions derived from modeling. We selected 10 regions observed as part of the SAGE-Spec program, to probe dust properties in various environments (diffuse, molecular and ionized regions). All data were smoothed to the 40arcsec angular resolution. The SEDs were modeled with DustEM models, using the standard Mathis RF, as well as three additional RFs, with stellar clusters ages ranging from 4 Myr to 600 Myr. Standard dust models used to reproduce the Galactic diffuse medium are clearly not able to reproduce the dust emission in the MIR wavelength domain. This analysis evidences the need of adjusting parameters describing the dust size distribution and shows a clear distinct behavior according to the type of environments. In addition, whereas the small grain emission always seems to be negligible at long wavelengths in our Galaxy, the contribution of this small dust component could be more important than expected, in the submm-mm range, in the LMC averaged SED. Properties of the small dust component of the LMC are clearly different from those of our Galaxy. Its abundance, significantly enhanced, could be the result of large grains shattering due to strong shocks or turbulence. In addition, this grain component in the LMC systematically shows smaller grain size in the ionized regions compared to the diffuse medium. Predictions of extinction curves show significantly distinct behaviors depending on the dust models but also from one region to another. Comparison of model predictions with the LMC mean extinction curve shows that no model gives satisfactory agreement using the Mathis radiation field while using a harder radiation field tends to improve the agreement

Abstract: Low mass star formation regions are unlikely to fully populate their initial mass functions, leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 10^2 to 10^7 M_sun, using a binary stellar population synthesis code. For clusters with stellar masses less than 10^5 M_sun, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the initial mass function appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities which directly measure the impact of the most massive stars, such as N_{ion}, xi_{ion} and beta_{UV}, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density.

Abstract: The study of the elusive hot component ($T \gtrsim 10^7$ K) of the Milky Way circumgalactic medium (CGM) is a novel topic to understand Galactic formation and evolution. In this work, we use the stacking technique through 46 lines of sight with Chandra ACIS-S HETG totaling over 10Ms of exposure time and 9 lines of sight with ACIS-S LETG observations totaling over 1Ms of exposure time, to study in absorption the presence of highly ionized metals arising from the super-virial temperature phase of the CGM. Focusing in the spectral range $4 - 8$ $\r{A}$, we were able to confirm the presence of this hot phase with high significance. We detected transitions of Si XIV K$\alpha$ (with total significance of 6.0$\sigma$) and, for the first time, SXVI K (total significance 4.8$\sigma$) in the rest frame of our own Galaxy. For S XVI K$\alpha$ we found a column density of $1.50^{+0.44}_{-0.38} \times 10^{16} \mathrm{cm}^{-2}$. For Si XIV K$\alpha$ we measured a column density of $0.87\pm{0.16} \times 10^{16} \mathrm{cm}^{-2}$. The lines of sight used in this work are spread across the sky, probing widely separated regions of the CGM. Therefore, our results indicate that this newly discovered hot medium extends throughout the halo, and is not related only to the Galactic Bubbles. The hot gas location, distribution, and covering factor, however, remain unknown. This component might contribute significantly to the missing baryons and metals in the Milky Way.

Abstract: Combining new HI data from a synergetic survey of ASKAP WALLABY and FAST with the ALFALFA data, we study the effect of ram-pressure and tidal interactions in the NGC 4636 group. We develop two parameters to quantify and disentangle these two effects on gas stripping in HI-bearing galaxies: the strength of external forces at the optical-disk edge, and the outside-in extents of HI-disk stripping. We find that gas stripping is widespread in this group, affecting 80% of HI-detected non-merging galaxies, and that 34% are experiencing both types of stripping. Among the galaxies experiencing both effects, the strengths (and extents) of ram-pressure and tidal stripping are independent of each other. Both strengths are correlated with HI-disk shrinkage. The tidal strength is related to a rather uniform reddening of low-mass galaxies ($M_*<10^9\,\text{M}_\odot$) when tidal stripping is the dominating effect. In contrast, ram pressure is not clearly linked to the color-changing patterns of galaxies in the group. Combining these two stripping extents, we estimate the total stripping extent, and put forward an empirical model that can describe the decrease of HI richness as galaxies fall toward the group center. The stripping timescale we derived decreases with distance to the center, from $\mathord{\sim}1\,\text{Gyr}$ around $R_{200}$ to $\mathord{\lesssim}10\,\text{Myr}$ near the center. Gas-depletion happens $\mathord{\sim}3\,\text{Gyr}$ since crossing $2R_{200}$ for HI-rich galaxies, but much quicker for HI-poor ones. Our results quantify in a physically motivated way the details and processes of environmental-effects-driven galaxy evolution, and might assist in analyzing hydrodynamic simulations in an observational way.

Abstract: We present results from the first 22 GHz space very-long-baseline interferometric (VLBI) imaging observations of M87 by RadioAstron. As a part of the Nearby AGN Key Science Program, the source was observed in Feb 2014 at 22 GHz with 21 ground stations, reaching projected $(u,v)$-spacings up to $\sim11\,$G$\lambda$. The imaging experiment was complemented by snapshot RadioAstron data of M87 obtained during 2013--2016 from the AGN Survey Key Science Program. Their longest baselines extend up to $\sim25\,$G$\lambda$. For all these measurements, fringes are detected only up to $\sim$2.8 Earth Diameter or $\sim$3 G$\lambda$ baseline lengths, resulting in a new image with angular resolution of $\sim150\,\mu$as or $\sim20$ Schwarzschild radii spatial resolution. The new image not only shows edge-brightened jet and counterjet structures down to submilliarcsecond scales but also clearly resolves the VLBI core region. While the overall size of the core is comparable to those reported in the literature, the ground-space fringe detection and slightly super-resolved RadioAstron image suggest the presence of substructures in the nucleus, whose minimum brightness temperature exceeds $T_{\rm B, min}\sim10^{12}\,$K. It is challenging to explain the origin of this record-high $T_{\rm B, min}$ value for M87 by pure Doppler boosting effect with a simple conical jet geometry and known jet speed. Therefore, this can be evidence for more extreme Doppler boosting due to a blazar-like small jet viewing angle or highly efficient particle acceleration processes occurring already at the base of the outflow.

Abstract: We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 282 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.210 Jy km/s. Depending on the assumptions made, this value can be interpreted as either an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 7.30\times10^{10}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $\rho_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 2.02\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $\rho_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on these achieved sensitivities, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data.