Papers

We present the second data release (DR2) of the TNG50-SKIRT Atlas (TSA), a library of synthetic, spatially resolved galaxy observables. The atlas is constructed by post-processing a stellar-mass-complete ($10^{9.8}~{\text{M}}_\odot < M_\star < 10^{12}~{\text{M}}_\odot$) sample of 1154 $z=0$ galaxies from the TNG50 cosmological hydrodynamical simulation with the Monte Carlo radiative transfer code SKIRT. Compared to the first release, TSA DR2 extends the wavelength coverage from the ultraviolet to the submillimetre, including dust emission, and incorporates updated stellar population models together with an improved treatment of dust-enshrouded star-forming regions. The atlas provides spatially resolved spectral energy distributions, broadband images, and physical property maps for multiple viewing orientations, as well as a catalogue of integrated properties enabling direct comparison with unresolved observations. We validate the data products through extensive quality control, including an assessment of Monte Carlo noise, and demonstrate their internal consistency using diagnostic relations between luminosities and star formation rates. TSA DR2 provides a versatile resource for studies of dust attenuation and emission, star formation tracers, galaxy morphology, and multi-wavelength scaling relations across spatial scales. The atlas and associated data products are publicly released and are intended to support a wide range of observationally oriented studies of galaxy evolution.

Stellar shells are low-surface-brightness structures that typically appear as concentric arcs. They are observed in many giant elliptical and lenticular galaxies, as well as in some spiral and dwarf galaxies. They presumably result from minor and intermediate close-to-radial mergers of galaxies. An essential factor in determining the merger time are the distances of shells from the centre of the galaxy. Consequently, estimates of the merger time can be significantly affected if some shells remain undetected due to projection effects. In this study, we present the first systematic investigation of how measured shell radii depend on the orientation of a galaxy relative to the observer. Using the Illustris TNG50 simulation, we examine shell galaxies from nine selected lines of sight and measure the shell radii. We model shell evolution and calculate the merger time accordingly for each viewing angle and quantify the impact of orientation on the inferred merger age. Our results indicate that the line of sight can have a substantial effect on merger-time estimates derived from shell radii.

We present spectral measurements for 17 Little Red Dots (LRDs) and 14 blue broad-line active galactic nuclei (AGNs) at $2\lesssim z \lesssim 6$ using multi-epoch JWST NIRSpec MSA spectra from the NEXUS program, sampling rest-frame timescales of $\sim 1-3$ months. Overall, the LRD population shows significantly enhanced Balmer decrement compared with both blue JWST AGNs at similar redshifts and 56 low-redshift broad-line AGNs matched in H$\rm\alpha$ luminosity. The rest-optical continua of LRDs show little ensemble variability (rms $\lesssim 3\%$), and the total H$\rm\alpha$ emission also shows weaker ensemble variability compared with low-redshift AGNs matched in H$\rm\alpha$ luminosity and rest-frame timescales. Based on the flux uncertainties, we constrain the intrinsic H$\rm\alpha$ rms variability to be $\lesssim 4\%$ for the LRD population over these timescales. Combining our results with recent broad-line variability measurements of LRDs over yearly to decade timescales reveals a low-level white-noise pattern across all timescales, in stark contrast to the variability amplitude ($\sim 6\%$ over monthly timescales) and red-noise pattern observed in normal AGNs. These results add to the growing observational studies that suggest population-wise, LRDs have weak variability both in optical continuum and broad-line emission. Furthermore, the distinct white-noise broad-line variability pattern suggests different production mechanisms of broad-line emission in LRDs as opposed to normal AGNs, and/or different properties of the driving ionizing flux from the central engine.

This work investigates the long-term evolution of the vertical structure of superthin galaxies within the framework of Milgromian dynamics (MOND). By constructing an observationally constrained model of UGC 7321, a typical superthin galaxy, we test whether its disc can maintain an extremely flattened structure in a MOND gravitational field. We also construct models with different values of the MOND depth index D_M to study how the global MOND depth affects disc evolution. We perform three-dimensional hydrodynamical N-body simulations using the publicly available code Phantom of RAMSES. In the observationally constrained model of UGC 7321, the galaxy develops a strong bar and undergoes a buckling instability during the early stages of the simulation. The bar strength then decreases gradually, and the system eventually exhibits a weak bar structure. The vertical evolution reflects the combined effects of heating induced by non-axisymmetric structures and vertical confinement in the Milgromian potential. The stellar disc undergoes only limited vertical thickening, and the disc remains largely within the superthin regime, h_z/R_D<0.1, after 5.0 Gyr. The comparison of models with different D_M values suggests that models with lower D_M values, associated in our model suite with higher baryonic masses or more compact discs, exhibit stronger vertical heating and more significant disc thickening. By contrast, models with higher D_M values, corresponding to lower masses or more diffuse structures, tend to maintain a superthin structure. Overall, the simulation results indicate that superthin discs can remain vertically thin during long-term isolated evolution in MOND, and that the long-term maintenance of superthin structures is influenced, at least partly, by the degree to which a galaxy lies in the low-acceleration regime.

Protostellar outflow shocks reshape local dust properties and molecular chemistry. The L1157 outflow is an archetypal chemically rich shocked region, but the thermal dust associated with its successive shocks has remained unresolved because molecular-line contamination obscures the broadband continuum. We obtained new James Clerk Maxwell Telescope (825–906 $\mu$m) spectral-line observations and Submillimeter Array (1.1–1.4 mm) continuum observations toward L1157 B0-B1-B2, probing spatial scales from 0.4 pc to 1200 au. After removing molecular-line contamination on a pixel-by-pixel basis, we derived the dust temperature, column density, and dust opacity index from continuum data spanning 70 $\mu$m to 1.3 mm. The line-corrected continuum maps reveal the dust distribution across successive shocks. The dust opacity index ($\beta\approx1.8$–2.3) indicates that grains have not grown to millimeter sizes throughout the shocked regions. Combined with previous $\rm NH_3$ observations, we find that the dust emission resolves into compact clumps along the precessing jet, whereas gaseous $\rm NH_3$ peaks at the shock fronts, reaching abundances of $\sim10^{-5}$ relative to $\rm H_2$, even where the 0.85 and 1.3 mm dust emission is detected at only 3–5$\sigma$. Our newly developed physicochemical shock model shows that $\rm NH_3$ forms predominantly on grain surfaces and is released by shock-induced sputtering, with the highest abundances occurring where post-shock re-adsorption remains inefficient. These results establish spatially resolved dust continuum imaging as a direct observational probe of grain evolution and provide new observational constraints on dust-gas interactions in protostellar shocks.