Based on systematic analyses of several popular non-relativistic energy density functionals with mean-field approximation, we estimate the value of the fourth-order symmetry energy $E_{\text{sym,4}}(\rho)$ at nuclear normal density $\rho_0$ and its density dependence, and explore the correlation between $E_{\text{sym,4}}(\rho_0)$ and other macroscopic quantities of nuclear matter properties. We use the empirical values of some nuclear macroscopic quantities to construct model parameter sets by Monte Carlo method for the conventional Skyrme-Hartree-Fock (SHF) model, the extended Skyrme-Hartree-Fock (eSHF) model, the Gogny-Hartree-Fock (GHF) model, and the momentum-dependent interaction (MDI) model. The value of $E_{\text{sym,4}}(\rho_0)$ is estimated to be $1.02\pm0.49$ MeV for the SHF model, $1.02\pm0.50$ MeV for the eSHF model, $0.70\pm0.60$ MeV for the GHF model, and $0.74\pm0.63$ MeV for the MDI model. Moreover, our results indicate that the density dependence of $E_{\text{sym,4}}(\rho)$ is model dependent, especially at higher densities. Furthermore, we find that the $E_{\text{sym},4}(\rho_0)$ has strong positive (negative) correlation with isoscalar (isovector) nucleon effective mass $m_{s,0}^*$ ($m_{v,0}^*$) at $\rho_0$. In particular, for the SHF and eSHF models, the $E_{\text{sym,4}}(\rho)$ is completely determined by the isoscalar and isovector nucleon effective masses $m_{s}^*(\rho)$ and $m_{v}^*(\rho)$, and the analytical expression is given. In the mean-field models, the magnitude of $ E_{\text{sym,4}}(\rho_0)$ is generally less than $2$ MeV, and its density dependence depends on models, especially at higher densities. $ E_{\text{sym,4}}(\rho_0)$ is strongly correlated with $m_{s,0}^*$ and $m_{v,0}^*$.

We investigated the dynamical structure of 53 elliptical galaxies, based on the {\it Chandra} archival X-ray data. In X-ray luminous galaxies, a temperature increases with radius and a gas density is systematically higher at the optical outskirts, indicating a presence of a significant amount of the group-scale hot gas. In contrast, X-ray dim galaxies show a flat or declining temperature profile against radius and the gas density is relatively lower at the optical outskirts. Thus it is found that X-ray bright and faint elliptical galaxies are clearly distinguished by the temperature and gas density profile. The mass profile is well scaled by a virial radius $r_{200}$ rather than an optical-half radius $r_e$, and is quite similar at $(0.001-0.03)r_{200}$ between X-ray luminous and dim galaxies, and smoothly connects to those of clusters of galaxies. At the inner region of $(0.001-0.01)r_{200}$ or $(0.1-1)r_e$, the mass profile well traces a stellar mass with a constant mass-to-light ratio of $M/L_{\rm B}=3-10(M_{\odot}/L_{\odot})$. $M/L_{\rm B}$ ratio of X-ray bright galaxies rises up steeply beyond $0.01r_{200}$, and thus requires a presence of massive dark matter halo. From the deprojection analysis combined with the {\it XMM-Newton} data, we found that X-ray dim galaxies, NGC 3923, NGC 720, and IC 1459, also have a high $M/L_{\rm B}$ ratio of 20--30 at 20 kpc, comparable to that of X-ray luminous galaxies. Therefore, dark matter is indicated to be common in elliptical galaxies, and their distribution almost follows the NFW profile, as well as galaxy clusters.