Resolved mass assembly and star formation in Milky Way Progenitors since $z=5$ from JWST/CANUCS: From clumps and mergers to well-ordered disks

We present a resolved study of 877 progenitors of Milky Way Analogs (MWAs) at $0.3<z<5$ selected with abundance matching in the ten fields of the Canadian NIRISS Unbiased Cluster Survey (CANUCS). Utilizing 18-21 bands of deep NIRCam, NIRISS, and HST photometry, we create resolved stellar mass maps and star formation rate maps via spectral energy distribution fitting with Dense Basis. We examine their resolved stellar mass and specific star formation rate (sSFR) profiles as a function of galactocentric radius, and find clear evidence for inside-out mass assembly. The total $\text{mstar}$ of the inner 2 kpc regions of the progenitors remain roughly constant (${10}^{9.3-9.4}\text{msun}$) at $2<z<5$, while the total $\text{mstar}$ of the regions beyond 2 kpc increases by 0.8 dex, from ${10}^{7.5}\text{msun}$ to ${10}^{8.3}\text{msun}$. Additionally, the sSFR of the outer regions increase with decreasing redshift, until $z\sim 2$. The median Sérsic index of the MWA progenitors stays nearly constant at $n\sim 1$ at $2<z<5$, while the half-mass radii of their stellar mass profiles double. We perform additional morphological measurements on the stellar mass maps via the Gini-M20 plane and asymmetry parameters. They show that the rate of double-peak mergers and disturbances to galaxy structure also increase with redshift, with $\sim 50\mathrm{\%}$ of galaxies at $4<z<5$ classified as disturbed, and $\sim 20\mathrm{\%}$ classified as ongoing mergers. Overall, the early evolution of MWAs is revealed as chaotic, with significant mergers and high SFRs. Mass growth is primarily inside-out and galaxies become more disk-like after $z=3$.