Schwarzschild Modeling of Barred S0 Galaxy NGC 4371

We apply the barred Schwarzschild method developed by Tahmasebzadeh et al. (2022) to a barred S0 galaxy, NGC 4371, observed by IFU instruments from the TIMER and ATLAS3D projects. We construct the gravitational potential by combining a fixed black hole mass, a spherical dark matter halo, and stellar mass distribution deprojected from $3.6$ $\mu$m S$^4$G image considering an axisymmetric disk and a triaxial bar. We create two sets of independent models, fitting the kinematic data derived from TIMER and ATLAS3D, separately. The models fit all the kinematic data remarkably well. We find a consistent bar pattern speed from the two sets of models, with $\Omega_{\rm p} = 23.6 \pm 2.8 \hspace{.08cm} \mathrm{km \hspace{.04cm} s^{-1} \hspace{.04cm} kpc^{-1} }$. The dimensionless bar rotation parameter is determined to be $ R_{\rm cor}/R_{\rm bar}=2.2 \pm 0.4$, indicating a slow bar in NGC 4371. Besides, we obtain a dark matter fraction $M_{\rm DM}/ M_{\rm total}$ of $\sim 0.51 \pm 0.06$ within the bar region. Our results support the scenario that bars may slow down in the presence of dynamical friction with a significant amount of dark matter in the disk regions. Based on our model, we further decompose the galaxy into multiple 3D orbital structures, including a BP/X bar, a classical bulge, a nuclear disk, and a main disk. The BP/X bar is not perfectly included in the input 3D density model, but BP/X-supporting orbits are picked through the fitting to the kinematic data. This is the first time a real barred galaxy has been modeled utilizing the Schwarzschild method including a 3D bar. Our model can be applied to a large number of nearby barred galaxies with IFU data, and it will significantly improve the previous models in which the bar is not explicitly included.