Search

The advent of large integral field spectroscopic surveys has found that elliptical galaxies (EGs) can be classified into two classes: the fast rotators (whose kinematics are dominated by rotation) and the slow rotators (which exhibit slow or no rotation pattern). It is often suggested that while the slow rotators typically have boxy isophotal shapes, have a high $\alpha$-to-iron abundance ratio, and are quite massive, the fast rotators often exhibit the opposite properties (that is, having disky isophotes, lower $\alpha$-to-iron ratio, and of typical masses). Whether the EGs consist of two distinct populations (i.e., a dichotomy exists), remains an unsolved issue. To examine the existence of the dichotomy, we used a sample of 1,895 EGs from the SDSS-IV MaNGA survey, and measured robustly the stellar kinematics, isophotal shapes, and [Mg/Fe] ratio. We confirmed the previous finding that the bulk of the EGs are disky (65%) and fast rotators (67%), but found no evidence supporting a dichotomy, based on a principal component analysis. The different classes (boxy/disky and slow/fast rotators) of EGs occupy slightly different loci in the principal component space. This may explain the observed trends that led to the premature support of a dichotomy based on small samples of galaxies., Comment: 27 pages, 19 figures, 5 tables. Submitted to ApJ. Comments are welcome!

Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only simulations and galaxies are strongly preferred, as producing mocks from full-fledged hydrodynamical simulations or semi-analytical models is too expensive. Here we present a Graph Neural Network-based model that is able to accurately predict key properties of galaxies such as stellar mass, $g-r$ color, star formation rate, gas mass, stellar metallicity, and gas metallicity, purely from dark matter properties extracted from halos along the full assembly history of the galaxies. Tests based on the TNG300 simulation of the IllustrisTNG project show that our model can recover the baryonic properties of galaxies to high accuracy, over a wide redshift range ($z = 0-5$), for all galaxies with stellar masses more massive than $10^9\,M_\odot$ and their progenitors, with strong improvements over the state-of-the-art methods. We further show that our method makes substantial strides toward providing an understanding of the implications of the IllustrisTNG galaxy formation model., Comment: 18 pages, 7 figures, 3 tables, 4 pages of appendices. Submitted to ApJ

Subhalo abundance matching (SHAM) has played an important role in improving our understanding of how galaxies populate their host dark matter halos. In essence, the SHAM framework is to find a dark matter halo property that best correlates with an attribute of galaxies, such as stellar mass. The peak value of the maximum circular velocity ($V_{\rm max}$) a halo/subhalo has ever attained throughout its lifetime, $V_{\rm peak}$, has been a popular choice for SHAM. A recent study by Tonnesen & Ostriker (2021) suggested that quantity $\phi$, which combines the present-day $V_{\rm max}$ and the peak value of halo dark matter mass, performs better in predicting stellar mass than $V_{\rm peak}$. Inspired by their approach, in this work, we find that further improvement can be achieved by a quantity $\psi_5$ that combines the 90th percentile of $V_{\rm max}$ a halo/subhalo has ever achieved with the 60th percentile of the dark matter halo time variation rate. Tests based on the simulation IllustrisTNG300 show that our new SHAM scheme, with just three free parameters, can improve the stellar mass prediction and mass-dependent clustering by 15% and 16% from $\phi$, respectively, over the redshift range $z=0-2$., Comment: 16 pages, 8 figures, 3 tables, Accepted by The Astrophysical Journal

The supernova remnant (SNR) B0532$-$67.5 in the Large Magellanic Cloud (LMC) was first diagnosed by its nonthermal radio emission and its SNR nature was confirmed by diffuse X-ray emission; however, no optical SNR shell is detected. The OB association LH75, or NGC 2011, is projected within the boundary of this SNR. We have analyzed the massive star population in and around SNR B0532$-$67.5: using optical photometric data to construct color-magnitude diagrams (CMDs), using stellar evolutionary tracks to estimate stellar masses, and using isochrones to assess the stellar ages. From these analyses, we find a 20-25 Myr population in LH75 and a younger population less than 10 Myr old to the southwest of LH75. The center of SNR B0532$-$67.5 is located closer to the core of LH75 than the massive stars to its southwest. We conclude that the SN progenitor was probably a member of LH75 with an initial mass $\sim$15 $M_\odot$. The SN exploded in an H I cavity excavated by the energy feedback of LH75. The low density of the ambient medium prohibits the formation of a visible nebular shell. Despite the low density in the ambient medium, physical properties of the hot gas within the SNR interior do not differ from SNRs with a visible shell by more than a factor of 2-3. The large-scale H I map shows that SNR B0532$-$67.5 is projected in a cavity that appears to be connected with the much larger cavity of the supergiant shell LMC-4., Comment: 14 pages, 2 tables, 9 figures, accepted for publication in The Astronomical Journal (AJ)

We have analyzed the ALMA archival data of the SO ($J_N=6_5-5_4$ and $J_N=7_6-6_5$), CO ($J=2-1$), and CCH ($N=3-2, J=7/2-5/2, F=4-3$) lines from the class 0 protobinary system, NGC1333 IRAS 4A. The images of SO ($J_N = 6_5-5_4$) and CO ($J=2-1$) successfully separate two northern outflow lobes connected to each protostar, IRAS 4A1 and IRAS 4A2. The outflow from IRAS 4A2 shows an S-shaped morphology, consisting of a flattened envelope around IRAS 4A2 with two outflow lobes connected to both edges of the envelope. The flattened envelope surrounding IRAS 4A2 has an opposite velocity gradient to that of the circumbinary envelope. The observed features are reproduced by the magnetohydrodynamic simulation of the collapsing core whose magnetic field direction is misaligned to the rotational axis. Our simulation shows that the intensity of the outflow lobes is enhanced on one side, resulting in the formation of S-shaped morphology. The S-shaped outflow can also be explained by the precessing outflow launched from an unresolved binary with a separation larger than 12 au (0.04arcsec). Additionally, we discovered a previously unknown extremely high velocity component at $\sim$45-90 km/s near IRAS 4A2 with CO. CCH ($J_{N,F}=7/2_{3,4}-5/2_{2,3}$) emission shows two pairs of blobs attaching to the bottom of shell like feature, and the morphology is significantly different from those of SO and CO lines. Toward IRAS 4A2, the S-shaped outflow shown in SO is overlapped with the edges of CCH shells, while CCH shells have the velocity gradients opposite to the flattened structure around IRAS 4A2., Comment: 55 pages, 20 figures