Search

We investigate how the stellar rotational support changes as a function of spatially resolved stellar population age ($\rm D_n4000$) and relative central stellar surface density ($\Delta \Sigma_1$) for MaNGA isolated/central disk galaxies. We find that the galaxy rotational support $\lambda_{R_\mathrm{e}}$ varies smoothly as a function of $\Delta \Sigma_1$ and $\rm D_n4000$. $\rm D_n4000$ vs. $\Delta \Sigma_1$ follows a "J-shape", with $\lambda_{R_\mathrm{e}}$ contributing to the scatters. In this "J-shaped" pattern rotational support increases with central $\rm D_n4000$ when $\Delta \Sigma_1$ is low but decreases with $\Delta \Sigma_1$ when $\Delta \Sigma_1$ is high. Restricting attention to low-$\Delta \Sigma_1$ (i.e, large-radius) galaxies, we suggest that the trend of increasing rotational support with $\rm D_n4000$ for these objects is produced by a mix of two different processes, a primary trend characterized by growth in $\lambda_{R_\mathrm{e}}$ along with mass through gas accretion, on top of which disturbance episodes are overlaid, which reduce rotational support and trigger increased star formation. An additional finding is that star forming galaxies with low $\Delta \Sigma_1$ have relatively larger radii than galaxies with higher $\Delta \Sigma_1$ at fixed stellar mass. Assuming that these relative radii rankings are preserved while galaxies are star forming then implies clear evolutionary paths in central $\rm D_n4000$ vs. $\Delta \Sigma_1$. The paper closes with comments on the implications that these paths have for the evolution of pseudo-bulges vs. classical-bulges. The utility of using $\rm D_n4000$-$\Delta \Sigma_1$ to study $\lambda_{R_\mathrm{e}}$ reinforces the notion that galaxy kinematics correlate both with structure and with stellar-population state, and indicates the importance of a multi-dimensional description for understanding bulge and galaxy evolution., Comment: 24 pages, 22 figures (including Appendix), accepted for publication in MNRAS

Observations with the James Webb Space Telescope have revealed a high abundance of bright galaxies at redshift, $z\gtrsim 12$, which has been widely interpreted as conflicting with the $\Lambda$CDM model. In Cowley et al. (2018) predictions were made - prior to the JWST observations - for the expected abundance of these galaxies using the Durham semi-analytic galaxy formation model, GALFORM, which is known to produce a realistic population of galaxies at lower redshifts including the present day. Key to this model is the assumption of a "top-heavy" initial mass function of stars formed in bursts (required to explain the number counts and redshift distribution of sub-millimetre galaxies). Here, we compare the rest-frame ultraviolet luminosity functions derived from JWST observations with those predicted by the Cowley et al. model up to $z=14$ and make further predictions for $z=16$. We find that below $z\sim 10$, the Cowley et al. predictions agree very well with observations, while agreement at $z\gtrsim12$ requires extending the model to take into account the timescale for the growth of obscuring dust grains and its dependence on gas metallicity. We trace the evolution of these galaxies from $z=14$ to $z=0$ and find that their descendants typically reside in halos with a median mass of $10^{13.6}\,h^{-1}\,\mathrm{M_{\odot}}$. The stellar masses of the descendants range from $10^{7}\,h^{-1}\,\mathrm{M_{\odot}}$ to $10^{11.5}\,h^{-1}\,\mathrm{M_{\odot}}$. Although these galaxies were all central galaxies at $z=14$, nearly half of their descendants end up as satellites in massive halos., Comment: 14 pages, 8 figures, submitted to MNRAS on 4 June, 2024

Conventional galaxy mass estimation methods suffer from model assumptions and degeneracies. Machine learning, which reduces the reliance on such assumptions, can be used to determine how well present-day observations can yield predictions for the distributions of stellar and dark matter. In this work, we use a general sample of galaxies from the TNG100 simulation to investigate the ability of multi-branch convolutional neural network (CNN) based machine learning methods to predict the central (i.e., within $1-2$ effective radii) stellar and total masses, and the stellar mass-to-light ratio $M_*/L$. These models take galaxy images and spatially-resolved mean velocity and velocity dispersion maps as inputs. Such CNN-based models can in general break the degeneracy between baryonic and dark matter in the sense that the model can make reliable predictions on the individual contributions of each component. For example, with $r$-band images and two galaxy kinematic maps as inputs, our model predicting $M_*/L$ has a prediction uncertainty of 0.04 dex. Moreover, to investigate which (global) features significantly contribute to the correct predictions of the properties above, we utilize a gradient boosting machine. We find that galaxy luminosity dominates the prediction of all masses in the central regions, with stellar velocity dispersion coming next. We also investigate the main contributing features when predicting stellar and dark matter mass fractions ($f_*$, $f_{\rm DM}$) and the dark matter mass $M_{DM}$, and discuss the underlying astrophysics., Comment: 17 pages, 12 figures, published on MNRAS

We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses of $10^{9-12}\ {\rm M_{\odot}}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes for the stellar $\overline{\gamma}_*$, dark $\overline{\gamma}_{\rm DM}$, and total $\overline{\gamma}_{\rm T}$ mass from the MaNGA DynPop project. The $\overline{\gamma}_{\rm T}$ approaches a constant value of 2.2 for high $\sigma_{\rm e}$ galaxies, and flattens for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\lesssim2.3$, reaching 1.5 for lg$(\sigma_{\rm e}/{\rm km\ s^{-1}})\approx1.8$. The total and stellar slopes track each other tightly, with $\overline{\gamma}_{\rm T}\approx\overline{\gamma}_*-0.174$ over the full $\sigma_{\rm e}$ range. This confirms the dominance of stellar matter within $R_{\rm e}$. We also show that there is no perfect conspiracy between baryonic and dark matter, as $\overline{\gamma}_*$ and $\overline{\gamma}_{\rm DM}$ do not vary inversely within the $\sigma_{\rm e}$ range. We find that the central galaxies from TNG50 and TNG100 simulations do not reproduce the observed galaxy mass distribution, which we attribute to the overestimated dark matter fraction, possibly due to a constant IMF and excessive adiabatic contraction effects in the simulations. Finally, we present the stacked dark matter density profiles and show that they are slightly steeper than the pure dark matter simulation prediction of $\overline{\gamma}_{\rm DM}\approx1$, suggesting moderate adiabatic contraction in the central region of galaxies. Our work demonstrate the power of stellar dynamics modelling for probing the interaction between stellar and dark matter and testing galaxy formation theories., Comment: 17 pages, 14 figures; accepted by MNRAS

The effect of angular momentum on galaxy formation and evolution has been studied for several decades. Our recent two papers using IllustrisTNG-100 simulation have revealed the acquisition path of the angular momentum from large-scale environment (satellites within hundreds of kpc) through the circum-galactic medium (CGM) to the stellar discs, putting forward the co-rotation scenario across the three distance scales. In real observations, although the rotation signature for the CGM and environmental three-dimensional (3d) angular momentum are difficult to obtain, line-of-sight kinematics of group member galaxies and stellar disc kinematics of central galaxies are available utilizing existing group catalogue data and integral field unit (IFU) data. In this paper, we use (1) the group catalogue of SDSS DR7 and MaNGA IFU stellar kinematic maps and (2) the group catalogue of GAMA DR4 data and SAMI IFU stellar kinematic maps, to test if the prediction above can be seen in real data. We found the co-rotation pattern between stellar discs and satellites can be concluded with 99.7 percent confidence level ($\sim 3\sigma$) when combining the two datasets. And the random tests show that the signal can be scarcely drawn from random distribution., Comment: 8 pages, 8 figures

Using the final MaNGA sample of 10K galaxies, we investigate the dark matter fraction $f_{\rm DM}$ within one half-light radius $R_{\rm e}$ for about 6K galaxies with good kinematics spanning a wide range of morphologies and stellar velocity dispersion. We employ two techniques to estimate $f_{\rm DM}$: (i) Jeans Anisotropic Modelling (JAM), which performs dark matter decomposition based on stellar kinematics and (ii) comparing the total dynamical mass-to-light ratios $(M/L)_{\rm JAM}$ and $(M_{\ast}/L)_{\rm SPS}$ from Stellar Population Synthesis (SPS). We find that both methods consistently show a significant trend of increasing $f_{\rm DM}$ with decreasing $\sigma_{\rm e}$ and low $f_{\rm DM}$ at larger $\sigma_{\rm e}$. For 235 early-type galaxies with the best models, we explore the variation of stellar initial mass function (IMF) by comparing the stellar mass-to-light ratios from JAM and SPS. We confirm that the stellar mass excess factor $\alpha_{\rm IMF}$ increases with $\sigma_{\rm e}$, consistent with previous studies that reported a transition from Chabrier-like to Salpeter IMF among galaxies. We show that the $\alpha_{\rm IMF}$ trend cannot be driven by $M_{\ast}/L$ or IMF gradients as it persists when allowing for radial gradients in our model. We find no evidence for the total $M/L$ increasing toward the centre. We detect weak positive correlations between $\alpha_{\rm IMF}$ and age, but no correlations with metallicity. We stack galaxy spectra according to their $\alpha_{\rm IMF}$ to search for differences in IMF-sensitive spectral features (e.g. the $\rm Na_{\rm I}$ doublet). We only find marginal evidence for such differences, which casts doubt on the validity of one or both methods to measure the IMF., Comment: 20 pages, 12 figures, 2 tables; Accepted for publication in MNRAS on 22 April 2024

The second moment of the stellar velocity within the effective radius, denoted by $\sigma_{\rm e}^2$, is a crucial quantity in galaxy studies as it provides insight into galaxy properties and their mass distributions. However, large spectroscopic surveys typically do not measure $\sigma_{\rm e}$ directly, instead providing $\sigma_{\rm aper}$, the second moment of the stellar velocity within a fixed fiber aperture. In this paper, we derive an empirical aperture correction formula, given by $\sigma_{\rm aper}/\sigma_{\rm e}=(R_{\rm aper}/R_{\rm e})^{\alpha}$, using spatially resolved stellar kinematics extracted from approximately 10,000 Sloan Digital Sky Survey-Mapping Nearby Galaxies at Apache Point Observatory (SDSS-MaNGA) integral field unit observations. Our analysis reveals a strong dependence of $\alpha$ on the $r$-band absolute magnitude $M_{\rm r}$, $g-i$ color, and Sersic index $n_{\rm Ser}$, where $\alpha$ values are lower for brighter, redder galaxies with higher Sersic indices. Our results demonstrate that the aperture correction derived from previous literature on early-type galaxies cannot be applied to predict the aperture corrections for galaxies with intermediate Sersic indices. We provide a lookup table of $\alpha$ values for different galaxy types, with parameters in the ranges of $-18>M_{\rm r}>-24$, $0.4

We analyse the integral-field spectroscopy data for the $\approx10,000$ galaxies in final data release of the MaNGA survey. We identify 188 galaxies for which the emission lines cannot be described by single Gaussian components. These galaxies can be classified into (1) 38 galaxies with broad $H\alpha$ and [OIII] $\lambda$5007 lines, (2) 101 galaxies with broad $H\alpha$ lines but no broad [OIII] $\lambda$5007 lines, and (3) 49 galaxies with double-peaked narrow emission lines. Most of the broad line galaxies are classified as Active Galactic Nuclei (AGN) from their line ratios. The catalogue helps us further understand the AGN-galaxy coevolution through the stellar population of broad-line region host galaxies and the relation between broad lines' properties and the host galaxies' dynamical properties. The stellar population properties (including mass, age and metallicity) of broad-line host galaxies suggest there is no significant difference between narrow-line Seyfert-2 galaxies and Type-1 AGN with broad $H\alpha$ lines. We use the broad-$H\alpha$ line width and luminosity to estimate masses of black hole in these galaxies, and test the $M_{BH}-\sigma_{e}$ relation in Type-1 AGN host galaxies. Furthermore we find three dual AGN candidates supported by radio images from the VLA FIRST survey. This sample may be useful for further studies on AGN activities and feedback processes., Comment: 21 pages, 17 figures, LaTeX. Accepted by MNRAS

We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from several kpc to Mpc scales. Intriguingly, we find that the excess surface density derived from stellar kinematics in the inner region cannot be explained by simply adding an NFW dark matter halo extrapolated from lensing measurement at a larger scale to a stellar mass component derived from the NASA-Sloan Atlas (NSA) catalogue. We find that a good fit to both data sets requires a stellar mass normalization about 3 times higher than that derived from the NSA catalogue, which would require an unrealistically too-heavy initial mass function for stellar mass estimation. If we keep the stellar mass normalization to that of the NSA catalogue but allow a varying inner dark matter density profile, we obtain an asymptotic slope of $\gamma_{\rm gnfw}$= $1.82_{-0.25}^{+0.15}$ and $\gamma_{\rm gnfw}$= $1.48_{-0.41}^{+0.20}$ for the group bin and the cluster bin respectively, significantly steeper than the NFW case. We also compare the total mass inner density slopes with those from TNG300 and find that the values from the simulation are lower than the observation by about $2\sigma$ level., Comment: 12 pages, 9 figures, accepted by MNRAS

We present dynamical scaling relations, combined with the stellar population properties, for a subsample of about 6000 nearby galaxies with the most reliable dynamical models extracted from the full MaNGA sample of 10K galaxies. We show that the inclination-corrected mass plane (MP) for both early-type galaxies (ETGs) and late-type galaxies (LTGs), which links dynamical mass, projected half-light radius $R_{\rm e}$, and the second stellar velocity moment $\sigma_{\rm e}$ within $R_{\rm e}$, satisfies the virial theorem and is even tighter than the uncorrected one. We find a clear parabolic relation between $\lg(M/L)_{\rm e}$, the total mass-to-light ratio within a sphere of radius $R_{\rm e}$, and $\lg\sigma_{\rm e}$, with the $M/L$ increasing with $\sigma_{\rm e}$ and for older stellar populations. However, the relation for ETGs is linear and the one for the youngest galaxies is constant. We confirm and improve the relation between average logarithmic total density slopes $\overline{\gamma_{\rm T}}$ and $\sigma_{\rm e}$: $\overline{\gamma_{\rm T}}$ become steeper with increasing $\sigma_{\rm e}$ until $\lg(\sigma_{\rm e}/{\rm km\,s^{-1}})\approx 2.2$ and then remain constant around $\overline{\gamma_{\rm T}}\approx -2.2$. The $\overline{\gamma_{\rm T}}-\sigma_{\rm e}$ variation is larger for LTGs than ETGs. At fixed $\sigma_{\rm e}$ the total density profiles steepen with galaxy age and for ETGs. We find generally low dark matter fractions, median $f_{\rm DM}(

This is the second paper of the MaNGA DynPop series, which analyzes the global stellar population, radial gradients, and non-parametric star-formation history of $\sim 10$K galaxies from the MaNGA Survey final data release 17 (DR17) and relates them with dynamical properties of galaxies. We confirm the correlation between the stellar population properties and the stellar velocity dispersion $\sigma_{\rm e}$, but also find that younger galaxies are more metal-poor at fixed $\sigma_{\rm e}$. Stellar age, metallicity, and mass-to-light ratio $M_{\ast}/L$ all decrease with increasing galaxy rotation, while their radial gradients become more negative (i.e. lower value at the outskirts). The exception is the slow rotators, which also appear to have significantly negative metallicity gradients, confirming the mass-metallicity gradient correlation. Massive disk galaxies in the green valley, on the $(\sigma_{\rm e},\rm Age)$ plane, show the most negative age and metallicity gradients, consistent with their old central bulges surrounded by young star-forming disks and metal-poor gas accretion. Galaxies with high $\sigma_{\rm e}$, steep total mass-density slope, low dark matter fraction, high $M_{\ast}/L$, and high metallicity have the highest star-formation rate at earlier times, and are currently quenched. We also discover a population of low-mass star-forming galaxies with low rotation but physically distinct from the massive slow rotators. A catalogue of these stellar population properties is provided publicly., Comment: 26 pages, 20 figures, 1 table. Accepted for publication in MNRAS. The catalogue is available at https://manga-dynpop.github.io/

This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically-aligned JAM$_{\rm cyl}$ or spherically-aligned JAM$_{\rm sph}$. We also make four assumptions for the models' dark vs. luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically-constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results., Comment: 30 pages, 15 figures, 4 tables, accepted by MNRAS, the catalogue is available here https://manga-dynpop.github.io/

This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys., Comment: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

The connection between halo gas acquisition through the circumgalactic medium (CGM) and galaxy star formation has long been studied. In this series of two papers, we put this interplay within the context of the galaxy environment on large scales (several hundreds of kpc), which, to a certain degree, maps out various paths for galaxy interactions. We use the IllustrisTNG-100 simulation to demonstrate that the large-scale environment modulates the circumgalactic gas angular momentum, resulting in either enhanced (Paper I) or suppressed (Paper II) star formation inside a galaxy. In this paper (Paper I), we show that the large-scale environment around a star-forming galaxy is often responsible for triggering new episodes of star formation. Such an episodic star formation pattern is well synced with a pulsating motion of the circumgalactic gas, which, on the one hand receives angular momentum modulations from the large-scale environment, yielding in-spiralling gas to fuel the star-forming reservoir, while, on the other hand, is affected by the feedback activities from the galaxy centre. As a result, a present-day star-forming galaxy may have gone through several cycles of star-forming and quiescent phases during its evolutionary history, with the circumgalactic gas carrying out a synchronized cadence of "breathing in and out" motions out to $\sim 100$ kpc., Comment: 15 pages, 12 figures, accepted for publication in MNRAS. GIFs of two samples can be found in https://astro.tsinghua.edu.cn/~dxu/HTM/Science/MetTemEvolutionExample_WXL21.htm

The gas needed to sustain star formation in galaxies is supplied by the circumgalactic medium (CGM), which in turn is affected by accretion from large scales. In a series of two papers, we examine the interplay between a galaxy's ambient CGM and central star formation within the context of the large-scale environment. We use the IllustrisTNG-100 simulation to show that the influence exerted by the large-scale galaxy environment on the CGM gas angular momentum results in either enhanced (Paper I) or suppressed (Paper II, this paper) star formation inside a galaxy. We find that for present-day quenched galaxies, both the large-scale environments and the ambient CGM have always had higher angular momenta throughout their evolutionary history since at least $z=2$, in comparison to those around present-day star-forming disk galaxies, resulting in less efficient gas inflow into the central star-forming gas reservoirs. A sufficiently high CGM angular momentum, as inherited from the larger-scale environment, is thus an important factor in keeping a galaxy quenched, once it is quenched. The process above naturally renders two key observational signatures: (1) a coherent rotation pattern existing across multiple distances from the large-scale galaxy environment, to the circumgalactic gas, to the central stellar disk; and (2) an anti-correlation between galaxy star-formation rates and orbital angular momenta of interacting galaxy pairs or groups., Comment: 13 pages, 11 figures, accepted for publication in MNRAS

Galaxy morphologies, kinematics, and stellar populations are thought to be linked to each other. However, both simulations and observations have pointed out mismatches therein. In this work, we study the nature and origin of the present-day quenched, bulge-dominated, but dynamically cold galaxies within a stellar mass range of $10.3\,\leqslant\,\log\,M_{\ast}/\mathrm{M_{\odot}}\,\leqslant\,11.2$ in the IllustrisTNG-100 Simulation. We compare these galaxies with a population of normal star-forming dynamically cold disc galaxies and a population of normal quenched dynamically hot elliptical galaxies within the same mass range. The populations of the present-day quenched and bulge-dominated galaxies (both being dynamically cold and hot) used to have significantly higher star-formation rates and flatter morphologies at redshift of $z\sim 2$. They have experienced more frequent larger mass-ratio mergers below $z \sim 0.7$ in comparison to their star-forming disc counterparts, which is responsible for the formation of their bulge-dominated morphologies. The dynamically cold populations (both being star-forming and quenched) have experienced more frequent prograde and tangential mergers especially below $z \sim 1$, in contrast to the dynamically hot ellipticals, which have had more retrograde and radial mergers. Such different merging histories can well explain the differences on the cold and hot dynamical status among these galaxies. We point out that the real-world counterparts of these dynamically cold and hot bulge-dominated quenched populations are the fast- and slow-rotating early-type galaxies, respectively, as seen in observations and hence reveal the different evolution paths of these two distinct populations of early-type galaxies., Comment: 14 pages, 12 figures; accepted for publication in MNRAS

From Rybicki's analysis using the Fourier slice theorem, mathematically it is possible to reproduce uniquely an edge-on axisymmetric galaxy's 3D light distribution from its 2D surface brightness. Utilizing galaxies from a cosmological simulation, we examine the ability of Syer and Tremaine's made-to-measure method and Schwarzschild's method for stellar dynamical modeling to do so for edge-on oblate axisymmetric galaxies. Overall, we find that the methods do not accurately recover the 3D distributions, with the made-to-measure method producing more accurate estimates than Schwarzschild's method. Our results have implications broader than just luminosity density, and affect other luminosity-weighted distributions within galaxies, for example, age and metallicity., Comment: Accepted for publication by Research in Astronomy and Astrophysics. 26 pages, 15 figures, 4 tables

A key feature of a large population of low-mass, late-type disk galaxies are star-forming disks with exponential light distributions. They are typically also associated with thin and flat morphologies, blue colors, and dynamically cold stars moving along circular orbits within co-planar thin gas disks. However, the latter features do not necessarily always imply the former, in fact, a variety of different kinematic configurations do exist. In this work, we use the cosmological hydrodynamical IllustrisTNG Simulation to study the nature and origin of dynamically hot, sometimes even counter-rotating, star-forming disk galaxies in the lower stellar mass range (between $5\times 10^9\,\mathrm{M_{\odot}}$ and $2\times 10^{10}\,\mathrm{M_{\odot}}$). We find that being dynamically hot arises in most cases as an induced transient state, for example due to galaxy interactions and merger activities, rather than as an age-dependent evolutionary phase of star-forming disk galaxies. The dynamically hot but still actively star-forming disks show a common feature of hosting kinematically misaligned gas and stellar disks, and centrally concentrated on-going star formation. The former is often accompanied by disturbed gas morphologies, while the latter is reflected in low gas and stellar spins in comparison to their dynamically cold, normal disk counterparts. Interestingly, observed galaxies from MaNGA with kinematic misalignment between gas and stars show remarkably similar general properties as the IllustrisTNG galaxies, and therefore are plausible real-world counterparts. In turn, this allows us to make predictions for the stellar orbits and gas properties of these misaligned galaxies., Comment: 18 pages, 17 figures, 2 tables. Accepted for publication in MNRAS

It is believed that massive galaxies have quenched their star formation because of active galactic nucleus feedback. However, recent studies have shown that some massive galaxies are still forming stars. We analyze the morphology of star formation regions for galaxies of stellar mass larger than 10$^{11.3}$ M$_{\odot}$ at around redshift $z_r=0.5$ using $u-z$ color images. We find that about $20\%$ of the massive galaxies are star-forming (SF) galaxies, and most of them ($\sim 85\%$) have asymmetric structures induced by recent mergers. Moreover, for these asymmetric galaxies, we find that the asymmetry of the SF regions becomes larger for bluer galaxies. Using the Illustris simulation, we can qualitatively reproduce the observed relation between asymmetry parameter and color. Furthermore, using the merger trees in the simulation, we find a correlation between the color of the main branch galaxies at $z_r=0.5$ and the sum of the Star Formation Rates (SFRs) of the recently accreted galaxies, which implies that star formation of the accreted galaxies has contributed to the observed star formation of the massive (host) galaxies (ex situ star formation). Furthermore, we find two blue and symmetric galaxies, candidates for massive blue disks, in our observed sample, which indicates that about $\sim 10\%$ of massive SF galaxies are forming stars in the normal mode of disk star formation (in situ star formation). With the simulation, we find that the disk galaxies at $z_r \approx 0.5$ should have experienced few major mergers during the last 4.3 Gyrs., Comment: 17 pages, 11 figures, Accepted for Publication in ApJ

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $\sigma_{\rm e}$. Here we look for {\em additional} trends at fixed $\sigma_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($\gamma_{\rm rms} \equiv \sigma(R_{\rm e}/4)/\sigma_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + \sigma^2}$), which we show traces the total mass density gradient $\gamma_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $\gamma_{\rm rms}$ increases with $\sigma_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $\sigma_{\rm e}$, where galaxies with larger $\gamma_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $\sigma_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $\sigma_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations., Comment: 8 pages, 6 figures; Accepted for publication in MNRAS

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto \sigma^b I^c$ with $R$, $\sigma$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log \sigma\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=\delta$, where $\delta \sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5\sigma^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations., Comment: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed., Comment: 12 pages, 6 figures, 3 tables, submitted to MNRAS. Comments are welcome!

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1., Comment: 11 pages, 10 figures, submitted to MNRAS