It is always interesting to investigate how well can a future experiment perform with respect to others (present or future ones). Cosmology is really an exciting field where a lot of puzzles are still unknown. In this article we consider a generalized dark energy (DE) scenario where anisotropic stress is present. We constrain this generalized cosmic scenario with an aim to investigate how gravitational waves standard sirens (GWSS) may constrain the anisotropic stress, which according to the standard cosmological probes, remains unconstrained. In order to do this, we generate the luminosity distance measurements from $\mathcal{O} (10^3)$ mock GW events which match the expected sensitivity of the Einstein Telescope. Our analyses report that, first of all, GWSS can give better constraints on various cosmological parameters compared to the usual cosmological probes, but the viscous sound speed appearing due to the dark energy anisotropic stress, is totally unconstrained even after the inclusion of GWSS.

Obtaining measurements of chromospheric and photometric activity of stars with near-solar fundamental parameters and rotation periods is important for a better understanding of solar-stellar connection. We select a sample of 2603 stars with near-solar fundamental parameters from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)-Kepler field and use LAMOST spectra to measure their chromospheric activity and Kepler light curves to measure their photospheric activity (i.e., the amplitude of the photometric variability). While the rotation periods of 1556 of these stars could not be measured due to the low amplitude of the photometric variability and highly irregular temporal profile of light curves, 254 stars were further identified as having near-solar rotation periods. We show that stars with near-solar rotation periods have chromospheric activities that are systematically higher than stars with undetected rotation periods. Furthermore, while the solar level of photospheric and chromospheric activity appears to be typical for stars with undetected rotation periods, the Sun appears to be less active than most stars with near-solar rotation periods (both in terms of photospheric and chromospheric activity).

We recently developed an automated method, auto-GMM to decompose simulated galaxies. It extracts kinematic structures in an accurate, efficient, and unsupervised way. We use auto-GMM to study the stellar kinematic structures of disk galaxies from the TNG100 run of IllustrisTNG. We identify four to five structures that are commonly present among the diverse galaxy population. Structures having strong to moderate rotation are defined as cold and warm disks, respectively. Spheroidal structures dominated by random motions are classified as bulges or stellar halos, depending on how tightly bound they are. Disky bulges are structures that have moderate rotation but compact morphology. Across all disky galaxies and accounting for the stellar mass within 3 half-mass radii, the kinematic spheroidal structures, obtained by summing up stars of bulges and halos, contribute ~45% of the total stellar mass, while the disky structures constitute 55%. This study also provides important insights about the relationship between kinematically and morphologically derived galactic structures. Comparing the morphology of kinematic structures with that of traditional bulge+disk decomposition, we conclude: (1) the morphologically decomposed bulges are composite structures comprised of a slowly rotating bulge, an inner halo, and a disky bulge; (2) kinematically disky bulges, akin to what are commonly called pseudo bulges in observations, are compact disk-like components that have rotation similar to warm disks; (3) halos contribute almost 30% of the surface density of the outer part of morphological disks when viewed face-on; and (4) both cold and warm disks are often truncated in central regions.

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

Using the $f\sigma_8(z)$ redshift space distortion (RSD) data, the $\sigma_8^0-\Omega_m^0$ tension is studied utilizing a parameterization of growth rate $f(z) = \Omega_m(z)^\gamma$. Here, $f(z)$ is derived from the expansion history $H(z)$ which is reconstructed from the observational Hubble data applying the Gaussian Process method. It is found that different priors of $H_0$ have a great influence on the evolution curve of $H(z)$ and the constraint of $\sigma_8^0-\Omega_m^0$. When using a larger $H_0$ prior, the low redshifts $H(z)$ deviate significantly from that of the $\Lambda$CDM model, which indicates that a dark energy model different from the cosmological constant can help to relax the $H_0$ tension problem. The tension between our best-fit values of $\sigma_8^0-\Omega_m^0$ and that of the \textit{Planck} 2018 $\Lambda$CDM (PLA) will disappear (less than $1\sigma$) when taking a prior for $H_0$ obtained from PLA. Moreover, the tension exceeds $2\sigma$ level when applying the prior $H_0 = 73.52 \pm 1.62$ resulted from the Hubble Space Telescope photometry.

The cosmographic approach is adopted to determine the spatial curvature (i.e., $\Omega_K$) combining the latest released cosmic chronometers data (CC), the Pantheon sample of type Ia supernovae observations, and the baryon acoustic oscillation measurements. We use the expanded transverse comoving distance $D_M(z)$ as a basic function for deriving $H(z)$ and the other cosmic distances. In this scenario, $\Omega_K$ can be constrained only by CC data. To overcome the convergence issues at high-redshift domains, two methods are applied: the Pad\'{e} approximants and the Taylor series in terms of the new redshift $y=z/(1+z)$. Adopting the Bayesian evidence, we find that there is positive evidence for the Pad\'{e} approximant up to order ($2,2$) and weak evidence for the Taylor series up to 3-rd order against $\Lambda\text{CDM}+\Omega_K$ model. The constraint results show that a closed universe is preferred by the present observations under all the approximants used in this study. And the tension level of the Hubble constant $H_0$ is less than $2\sigma$ significance between different approximants and the local distance ladder determination. For each assumed approximant, $H_0$ is anti-correlated with $\Omega_K$ and the sound horizon at the end of the radiation drag epoch, which indicates that the $H_0$ tension problem can be slightly relaxed by introducing $\Omega_K$ or any new physics which can reduce the sound horizon in the early universe.

Short (inner) bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scales, thus feeding central black holes. Recent numerical simulations have shown that the growth of central black holes in galaxies can destroy short bars, when the black hole reaches a mass of $\sim 0.1\%$ of the total stellar mass of the galaxy. We study $N$-body simulations of galaxies with single and double bars to track the long-term evolution of the central stellar mass distribution. We find that the destruction of the short bar contributes significantly to the growth of the bulge. The final bulge mass is roughly equal to the sum of the masses of the initial pseudo bulge and short bar. The initially boxy/peanut-shaped bulge of S\'ersic index $n\lesssim1$ is transformed into a more massive, compact structure that bears many similarities to a classical bulge, in terms of its morphology ($n \approx 2$), kinematics (dispersion-dominated, isotropic), and location on standard scaling relations (Kormendy relation, mass-size relation, and correlations between black hole mass and bulge stellar mass and velocity dispersion). Our proposed channel for forming classical bulges relies solely on the destruction of short bars without any reliance on mergers. We suggest that some of the less massive, less compact classical bulges were formed in this manner.

Galaxies host a wide array of internal stellar components, which need to be decomposed accurately in order to understand their formation and evolution. While significant progress has been made with recent integral-field spectroscopic surveys of nearby galaxies, much can be learned from analyzing the large sets of realistic galaxies now available through state-of-the-art hydrodynamical cosmological simulations. We present an unsupervised machine learning algorithm, named auto-GMM, based on Gaussian mixture models, to isolate intrinsic structures in simulated galaxies based on their kinematic phase space. For each galaxy, the number of Gaussian components allowed by the data is determined through a modified Bayesian information criterion. We test our method by applying it to prototype galaxies selected from the cosmological simulation IllustrisTNG. Our method can effectively decompose most galactic structures. The intrinsic structures of simulated galaxies can be inferred statistically by non-human supervised identification of galaxy structures. We successfully identify four kinds of intrinsic structures: cold disks, warm disks, bulges, and halos. Our method fails for barred galaxies because of the complex kinematics of particles moving on bar orbits.

The detection of gravitational waves (GW) by the LIGO and Virgo collaborations offers a whole new range of possible tests and opens up a new window which may shed light on the nature of dark energy and dark matter. In the present work we investigate how future gravitational waves data could help to constrain different dynamical dark energy models. In particular, we perform cosmological forecastings of a class of well known and most used dynamical dark energy models using the third-generation gravitational wave detector, the Einstein Telescope. We have considered 1000 simulated GW events in order to constrain the parameter space of the dynamical dark energy models. Our analyses show that the inclusion of the GW data from the Einstein Telescope, significantly improves the parameter space of the dynamical dark energy models compared to their constraints extracted from the standard cosmological probes, namely, the cosmic microwave observations, baryon acoustic oscillations distance measurements, Supernove type Ia, and the Hubble parameter measurements.

In this paper, we have constrained the neutrino mass and mass hierarchy in the $\Lambda$CDM cosmology with the neutrino mass hierarchy parameter $\Delta$, which represents different mass orderings, by using the {\it Planck} 2015 + BAO + SN + $H_{0}$ data set, together with the neutrino oscillation and neutrinoless double beta decay data. We find that the mass of the lightest neutrinos and the total neutrino mass are no more than $0.035$eV and $0.133$ eV at $95\%$ confidence level, respectively. Comparing the result of our joint analysis with that obtained using cosmological data alone, we find that, by adding the neutrino oscillation and neutrinoless double beta decay data, the tendency for normal hierarchy has increased a lot. By means of importance sampling, three other priors are taken into account, i.e., the flat logarithmic prior on the absolute value of the neutrino hierarchy parameter $\Delta$, the flat linear prior on the total neutrino mass $\Sigma m_{\nu}$, and the flat logarithmic prior on $\Sigma m_{\nu}$. We find that the preference for the normal hierarchy is in agreement whatever what kinds of priors we choose. Finally, we make a Bayesian model analysis about four priors and we find that flat-linear and the flat logarithmic priors on $\Sigma m_{\nu}$ are the most favored priors.

Compact ellipticals (cEs) are outliers from the scaling relations of early-type galaxies, particularly the mass-metallicity relation which is an important outcome of feedback. The formation of such low-mass, but metal-rich and compact, objects is a long-standing puzzle. Using a pair of high-resolution N-body+gas simulations, we investigate the evolution of a gas-rich low-mass galaxy on a highly radial orbit around a massive host galaxy. As the infalling low-mass galaxy passes through the host's corona at supersonic speeds, its diffuse gas outskirts are stripped by ram pressure, as expected. However, the compactness increases rapidly because of bursty star formation in the gas tidally driven to the centre. The metal-rich gas produced by supernovae and stellar winds is confined by the ram pressure from the surrounding environment, leading to subsequent generations of stars being more metal-rich. After the gas is depleted, tidal interactions enhance the metallicity further via the stripping of weakly bound, old, and metal-poor stars, while the size of the satellite is changed only modestly. The outcome is a metal-rich cE that is consistent with observations. These results argue that classical cEs are neither the stripped remnants of much more massive galaxies nor the merger remnants of normal dwarfs. We present observable predictions that can be used to test our model.

The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies.

Taking the neutrino oscillation data into consideration, a dimensionless parameter $\Delta = (m_3-m_1)/(m_3+m_1)$ is adopted to parameterize the three neutrino mass eigenstates and the normal (positive $\Delta$) or inverted (negative $\Delta$) mass hierarchies in three typical cosmological models. Using the currently available cosmic observational data, several Markov Chain Monte Carlo chains are obtained with uniform priors on the free parameters at first. Applying importance sampling the results are compared with three new priors, i.e., logarithmic prior on $|\Delta|$, linear and logarithmic priors on $\Sigma m_\nu$. It turns out that the three new priors increase the upper limits of neutrino mass, but do not change the tendency towards different model's preference for different hierarchies, i.e., the normal hierarchy tends to be favored by $\Lambda$CDM and $w$CDM, which, however, disappears in the $w_0 w_a$CDM model. In addition, the almost symmetrical contours in the $w-\Delta$, $w_0-\Delta$, $w_a-\Delta$ planes indicate that the normal and inverted hierarchy have strong degeneracy. Finally, we perform a Bayesian model comparison analysis, finding that flat linear prior on $\Delta$ and $w_0 w_a$CDM are the most preferred prior and model, respectively.

We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies.

The shape of LOSVDs (line-of-sight velocity distributions) carries important information about the internal dynamics of galaxies. The skewness of LOSVDs represents their asymmetric deviation from a Gaussian profile. Correlations between the skewness parameter ($h_3$) and the mean velocity (\vm) of a Gauss-Hermite series reflect the underlying stellar orbital configurations of different morphological components. Using two self-consistent $N$-body simulations of disk galaxies with different bar strengths, we investigate \hv\ correlations at different inclination angles. Similar to previous studies, we find anticorrelations in the disk area, and positive correlations in the bar area when viewed edge-on. However, at intermediate inclinations, the outer parts of bars exhibit anticorrelations, while the core areas dominated by the boxy/peanut-shaped (B/PS) bulges still maintain weak positive correlations. When viewed edge-on, particles in the foreground/background disk (the wing region) in the bar area constitute the main velocity peak, whereas the particles in the bar contribute to the high-velocity tail, generating the \hv\ correlation. If we remove the wing particles, the LOSVDs of the particles in the outer part of the bar only exhibit a low-velocity tail, resulting in a negative \hv\ correlation, whereas the core areas in the central region still show weakly positive correlations. We discuss implications for IFU observations on bars, and show that the variation of the \hv\ correlation in the disk galaxy may be used as a kinematic indicator of the bar and the B/PS bulge.

The growth of black holes (BHs) in disk galaxies lacking classical bulges, which implies an absence of significant mergers, appears to be driven by secular processes. Short bars of sub-kiloparsec radius have been hypothesized to be an important mechanism for driving gas inflows to small scale, feeding central BHs. In order to quantify the maximum BH mass allowed by this mechanism, we examine the robustness of short bars to the dynamical influence of BHs. Large-scale bars are expected to be robust, long-lived structures; extremely massive BHs, which are rare, are needed to completely destroy such bars. However, we find that short bars, which are generally embedded in large-scale outer bars, can be destroyed quickly when BHs of mass $M_{\rm bh}\sim0.05-0.2\%$ of the total stellar mass ($M_\star$) are present. In agreement with this prediction, all galaxies observed to host short bars have BHs with a mass fraction less than $0.2\%M_\star$. Thus the dissolution of short inner bars is possible, perhaps even frequent, in the universe. An important implication of this result is that inner bar-driven gas inflows may be terminated when BHs grow to $\sim0.1\%M_\star$. We predict that $0.2\%M_\star$ is the maximum mass of BHs allowed if they are fed predominately via inner bars. This value matches well the maximum ratio of BH-to-host galaxy stellar mass observed in galaxies with pseudo-bulges and most nearby AGN host galaxies. This hypothesis provides a novel explanation for the lower $M_{\rm bh}/M_\star$ in galaxies which have avoided significant mergers compared with galaxies with classical bulges.

High precision and long-lasting \emph{Kepler} data enabled us to estimate stellar properties with asteroseismology as an accurate tool. We performed asteroseismic analysis on six solar-like stars observed by the \emph{Kepler} mission: KIC 6064910, KIC 6766513, KIC 7107778, KIC 10079226, KIC 10147635 and KIC 12069127. The extraction of seismic information includes two parts. First, we obtained two global asteroseismic parameters, mean large separation $\Delta\nu$ and frequency of maximum power $\nu_{\rm{max}}$, with autocorrelation function and collapsed autocorrelation function. Second, we extracted individual oscillation modes $\nu_{nl}$ with low-$l$ degree using a least-squares fit. Stellar grid models were built with Yale Rotating Stellar Evolution Code (YREC) to analyze stellar properties. They covered the range of $M=0.8\sim1.8$ $M_{\odot}$ with a step of 0.02 $M_{\odot}$ and $\rm{[Fe/H]}=-0.3\sim0.4$ $\rm{dex}$ with a step of 0.1 $\rm{dex}$. We used a Bayesian approach to estimate stellar fundamental parameters of the six stars, under the constraints of asteroseismic parameters ($\Delta\nu$, $\nu_{\rm{max}}$) and non-asteroseismic parameters ($T_{\rm{eff}}$, $\rm{[Fe/H]}$). We discover that the six targets include five sub-giant stars with $1.2\sim1.5$ $M_{\odot}$ and one main-sequence star with $1.08M_{\odot}$, and with ages in the range of $3\sim5$ $\rm{Gyr}$.

In barred galaxies, the contours of stellar velocity dispersions ($\sigma$) are generally expected to be oval and aligned with the orientation of bars. However, many double-barred (S2B) galaxies exhibit distinct $\sigma$ peaks on the minor axis of inner bars, which we termed "$\sigma$-humps," while two local $\sigma$ minima are present close to the ends of inner bars, i.e., "$\sigma$-hollows." Analysis of numerical simulations shows that $\sigma_z$-humps or hollows should play an important role in generating the observed $\sigma$-humps+hollows in low-inclination galaxies. In order to systematically investigate the properties of $\sigma_z$ in barred galaxies, we apply the vertical Jeans equation to a group of well-designed three-dimensional bar+disk(+bulge) models. A vertically thin bar can lower $\sigma_z$ along the bar and enhance it perpendicular to the bar, thus generating $\sigma_z$-humps+hollows. Such a result suggests that $\sigma_z$-humps+hollows can be generated by the purely dynamical response of stars in the presence of a, sufficiently massive, vertically thin bar, even without an outer bar. Using self-consistent $N$-body simulations, we verify the existence of vertically thin bars in the nuclear-barred and S2B models which generate prominent $\sigma$-humps+hollows. Thus the ubiquitous presence of $\sigma$-humps+hollows in S2Bs implies that inner bars are vertically thin. The addition of a bulge makes the $\sigma_z$-humps more ambiguous and thus tends to somewhat hide the $\sigma_z$-humps+hollows. We show that $\sigma_z$ may be used as a kinematic diagnostic of stellar components that have different thickness, providing a direct perspective on the morphology and thickness of nearly face-on bars and bulges with integral field unit spectroscopy.

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.

Using high resolution $N$-body simulations, we recently reported that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady double-barred (S2B) structure. Here we study the kinematics of these S2B simulations, and compare them to integral-field observations from ATLAS$^{3D}$ and SAURON. We show that S2B galaxies exhibit several distinct kinematic features, namely: (1) significantly distorted isovelocity contours at the transition region between the two bars, (2) peaks in $\sigma_\mathrm{LOS}$ along the minor axis of inner bars, which we term "$\sigma$-humps", that are often accompanied by ring/spiral-like features of increased $\sigma_\mathrm{LOS}$, (3) $h_3-\bar{v}$ anti-correlations in the region of the inner bar for certain orientations, and (4) rings of positive $h_4$ when viewed at low inclinations. The most impressive of these features are the $\sigma$-humps, these evolve with the inner bar, oscillating in strength just as the inner bar does as it rotates relative to the outer bar. We show that, in cylindrical coordinates, the inner bar has similar streaming motions and velocity dispersion properties as normal large-scale bars, except for $\sigma_z$, which exhibits peaks on the minor axis, i.e., humps. These $\sigma_z$ humps are responsible for producing the $\sigma$-humps. For three well-resolved early-type S2Bs (NGC 2859, NGC 2950, and NGC 3941) and a potential S2B candidate (NGC 3384), the S2B model qualitatively matches the integral-field data well, including the "$\sigma$-hollows" previously identified. We also discuss the kinematic effect of a nuclear disk in S2Bs.

About one third of early-type barred galaxies host small-scale secondary bars. The formation and evolution of such double-barred galaxies remain far from being well understood. In order to understand the formation of such systems, we explore a large parameter space of isolated pure-disk simulations. We show that a dynamically cool inner disk embedded in a hotter outer disk can naturally generate a steady secondary bar while the outer disk forms a large-scale primary bar. The independent bar instabilities of inner and outer disks result in long-lived double-barred structures whose dynamical properties are comparable with observations. This formation scenario indicates that the secondary bar might form from the general bar instability, the same as the primary bar. Under some circumstances, the interaction of the bars and the disk leads to the two bars aligning or single, nuclear, bars only. Simulations that are cool enough of the center to experience clump instabilities may also generate steady double-barred galaxies. In this case, the secondary bars are "fast", i.e., the bar length is close to the co-rotation radius. This is the first time that double-barred galaxies containing a fast secondary bar are reported. Previous orbit-based studies had suggested that fast secondary bars are not dynamically possible.