Stellar Mass-Size Relation and Morphology of Massive Quiescent Galaxies at $3 < z < 4$ with JWST

Abstract: We present the morphological properties of 17 spectroscopically confirmed massive quiescent galaxies ($10.2 < \log(M_{\ast}/M_{\odot}) < 11.2$) at $3.0 < z < 4.3$, observed with JWST/NIRSpec and NIRCam. Using S\'{e}rsic profile fits to the F277W and F444W imaging, we derive the size-mass relation and find typical sizes of $\sim$0.6-0.8 kpc at $M_{\ast} = 5 \times 10^{10}~M_{\odot}$, consistent with a factor of $\sim$7 growth from $z \sim 4$ to the present day, including a factor of $\sim$2 by $z \sim 2$. Galaxies with earlier formation times ($z_{50} \gtrsim 5$) tend to be larger (up to 0.3 dex above the relation), while late-forming galaxies are more compact. Using a random forest regressor, we identify local environmental density, quantified by $\log(1+\delta^{\prime}_{3})$ from the distances to the three nearest neighbors, as the strongest predictor of bulge-to-total ratio ($B/T$), which ranges from 0.25 to 1. Applying the same method to the IllustrisTNG simulation, the ex-situ stellar mass fraction ($f_{\ast,\mathrm{ex\text{-}situ}}$) -- a proxy for mergers -- is the dominant predictor of $B/T$. Galaxies with high $B/T$ in dense environments show bursty star formation and short quenching timescales ($\lesssim0.4$ Gyr), consistent with bulge growth via merger-induced starbursts; in simulations, such galaxies exhibit elevated ex-situ fractions ($\sim20\%$-$30\%$). In contrast, some high-$B/T$ galaxies in intermediate-density environments have low ex-situ fractions, suggesting that additional mechanisms -- such as violent disk instabilities -- may contribute. These findings point to multiple bulge growth pathways at high redshift, unified by rapid gas accretion, central starbursts, and AGN feedback, as predicted by cosmological simulations.