2025-07-17 | | Total: 29
Context. Dust is a fundamental component of the interstellar medium (ISM) and plays a critical role in galaxy evolution. Dust grains influence the ISM by cooling the gas, altering its chemistry, and absorbing stellar radiation, re-emitting it at longer wavelengths in the far-infrared (FIR) and sub-millimeter regimes. The cold dust component, which dominates the dust mass, is primarily heated by stellar radiation, including both young, massive stars and the diffuse emission from older stars. Understanding dust heating is essential to trace the connection between stellar populations and their environments. Aims. We aim to identify the dominant heating mechanisms of the cold dust in typical nearby spiral galaxies and explore the contributions of young and evolved stars to dust heating. Methods. Using 18 large, face-on spiral galaxies from the DustPedia project, we apply two complementary approaches: (1) correlation analysis between dust temperature (T_dust), SFR surface density (Sigma_SFR), and stellar mass surface density (Sigma_Mstar); and (2) study of the relationship between T_dust and dust mass surface density (Sigma_dust). Results. T_dust peaks at ~24 K in galaxy centers and drops to ~15 K at large radii. Galaxies with and without AGNs show similar T_dust profiles. For ~72% of the sample, both methods agree on the dominant heating source. Overall, we find that both young and evolved stars contribute to dust heating, with their relative roles varying between galaxies.
We present observations of the nearby extremely metal-poor galaxy I Zw 18 using the Keck Cosmic Web Imager (KCWI) and the JWST Mid-InfraRed Instrument (MIRI) Integral Field Spectrographs (IFS). From optical and mid-IR oxygen emission lines, we measured direct-method abundances for three ionic states of oxygen, including O$^{3+}$/H$^+$. In contrast to previous studies of I Zw 18, the high spatial resolution afforded by KCWI and MIRI/MRS revealed chemical inhomogeneities on 60 pc scales in the form of metal-poor pockets and metal-enriched gas. These are located outside I Zw 18's star-forming complexes having possibly been dispersed beyond these regions via stellar feedback effects. We found that metallicities derived using a single low-ionization density tracer, and T$_{\rm e}$([O II]) derived from a temperature relationship commonly used in high-$z$ galaxy studies, exhibited the largest scatter and underestimated the metallicity compared to those derived using multi-ion densities and estimated T$_{\rm e}$([N II]). Finally, we compared O$^{3+}$/H$^+$ abundances from a theoretical ionization correction factor (ICF) against observed values and found that the oxygen ICF underestimates the O$^{3+}$/H$^+$ abundance by a factor of 2, indicating that either additional ionizing sources are needed or standard stellar population models are unable to produce the requisite ionizing flux.
We propose an empirical model to describe and constrain the baryon cycle process during galaxy evolution. This model utilizes the evolution of star formation rate, derived from the stellar mass-halo mass relations (SHMRs) across different redshifts, and the cold gas content, derived from the NeutralUniverseMachine model, to constrain gas accretion and recycle of gas outflow in the model galaxy. Additionally, through detailed modeling of each cycle process, particularly the recycling process, and utilizing the gas-phase mass-metallicity relation (MZR) at $z=0$ as a constraint, our model establishes a relation between the recycle fraction and halo mass. It is found that the fraction of gas recycled from the outflow is a function of halo mass, with a value of $25\%$ in galaxies with halo mass $\sim10^{10.4}M_{\rm \odot}$, increasing to $75\%$ in halos with mass $\sim 10^{12}M_{\rm \odot}$. We also find that the mass loading factor from the FIRE-2 simulation matches well with the constraints from both observational data and our model. Furthermore, using the gas content and metallicity of the circumgalactic medium (CGM) obtained from hydrodynamical simulations as constraints, our model predicts that on average $70\%$ of universal baryon accretion is accreted to the halo and $80\%$ of the non-recycled gas in the outflow has escaped from the galaxy, entering the intergalactic medium (IGM). However, we note that future observational data are needed to finally constrain the mass and metal exchange between the CGM and the IGM.
The phenomenon of magnetic braking is one of the significant physical effects of the magnetic field in rotating molecular clouds. The physical characteristics of the core can affect on the core rotation rate. and one of the important parameter is the core density structure. According to observation, by regarding the power-law density distribution, $r^{-p}$, for molecular cloud cores, using smoothed particle hydrodynamics simulation, the results show that the increasing of density steepness (i.e., larger $p$) leads to the intensity of the toroidal components of the magnetic field and as a result larger $B_{\phi}$-components lead to more transfer of angular momentum to the outward. Thus, results show that the magnetic braking being stronger with increasing density slope in non-uniform molecular core. For example, the rotation of the system can approximately decrease by fifty percent from $p=0.2$ to $p=1.8$ for a non-uniform system.
We present a statistical analysis of the 175 SPARC galactic rotation curves to test the hypothesis of whether the Keplerian velocity tapering at large radii ($v(r)\propto 1/\sqrt{r}$) of the Navarro-Frenk-White (NFW) halo model agrees with observational data. The null hypothesis is Rubin's flat-rotation curve, $v(r)=\text{constant}$ -such as can be obtained from a spherical, isothermal-like density profile, or alternatively with a very prolate halo-. To decide whether we adopt the null (Rubin behaviour) or alternative (NFW behaviour) hypothesis, we evaluate the derivative in each galaxy of $v(r)$ with its last data points. We conclude that the data is presently compatible with the null hypothesis -no taper off, no decline of $v(r)$ is seen.
In this work we study a robust, $K_s$-band complete, spectroscopically-confirmed sample of 104 unobscured (Type-1) quasars within the COSMOS and XMM-LSS fields of the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey, at 0.60 < $z$ < 3.41. The quasars are selected via $gJK_s$ colour-space and, with 1.3-GHz flux-densities reaching rms ~ 3.0$\mu$Jy beam$^{-1}$, we find a radio-loudness fraction of 5 per cent. Thanks to the deep, multiwavelength datasets that are available over these fields, the properties of radio-loud and radio-quiet quasars can be studied in a statistically-robust way, with the emphasis of this work being on the active-galactic-nuclei (AGN)-related and star-formation-related contributions to the total radio emission. We employ multiple star-formation-rate estimates for the analysis so that our results can be compared more-easily with others in the literature, and find that the fraction of sources that have their radio emission dominated by the AGN crucially depends on the SFR estimate that is derived from the radio luminosity. When redshift dependence is not taken into account, a larger fraction of sources is classed as having their radio emission dominated by the AGN. When redshift dependence $is$ considered, a larger fraction of our sample is tentatively classed as 'starbursts'. We also find that the fraction of (possible) starbursts increases with redshift, and provide multiple suggestions for this trend.
Far-infrared (FIR) and mid-infrared (MIR) fine-structure lines (FSLs) provide key diagnostics of physical conditions in the interstellar medium (ISM). Building on empirical relations established in our previous work, we use photoionization models to systematically investigate the emission from both ionized and neutral gas phases in a coherent structure. By applying power-law fits to model parameters, we quantitatively capture how key FIR FSL ratios scale with physical properties such as density, radiation field strength and hardness, and elemental abundances. Our analysis confirms the primary dependencies behind most observed empirical trends and establishes certain FIR FSL ratios as tracers of physical parameters, while revealing that parameter marginalization-particularly in density and the O/H-$U$-$Q_1/Q_0$ relation-plays a crucial role in shaping tight correlations seen in galaxies. We also identify persistent challenges, including degeneracies between ionization parameter and radiation field hardness, uncertainties in neutral gas density, and difficulties in modeling dusty H II regions. We outline the fundamental observational and theoretical limitations of current FIR FSL diagnostics, and highlight prospects for advancing the field through comprehensive, multi-wavelength studies of diverse galaxy populations.
We report the discovery of a large-scale structure containing multiple overdensities at $z= 4.54 \pm 0.03$ in the COSMOS field, with the most prominent one likely infalling towards the recently identified Taralay protocluster. We use combined wide-band and narrow-band optical photometry to identify Ly$\alpha$ emitters (LAEs) within a 21 cMpc radius from the submm source J1000+0234, at $z= 4.54$, to identify typical star-forming galaxies that may trace an underlying structure. Our approach selects line emitters as narrow-band excess objects and we use the COSMOS2020 photometric redshift catalog to eliminate potential low-redshift interlopers whose line emission (e.g. [OIII] at $z\sim 0.3$) might be responsible for the observed excess in the narrow band. In comparison with the LAE density in the field, our results point to a mean LAE number overdensity of $\bar\delta = 3$ spanning a region of $27 \times 20 \times 36$ cMpc$^3$, probably evolving into a moderate-mass cluster ($3 - 10 \times 10^{14} \, M_\odot$) at $z\sim 0$. This work supports the idea that submm sources, although offset from the major overdensity peaks, serve as traces of moderately massive, potentially infalling structures.
We present the discovery of two intersecting radio shells, likely radio relics, surrounding a compact galaxy group dominated by a massive elliptical galaxy. The shells present as partial, edge-brightened rings with diameters of about 240" (720 kpc) each and resemble a pair of odd radio circles. The central galaxy, WISEA J184105.19-654753.8, which shows signs of interactions, is radio bright, has a stellar mass of 3.1 x 10^11 Msun (for a redshift of zphot = 0.18) and is located in the intersect region. The double radio shell system, which we refer to as ORC J1841-6547 (also known as ORC 6), was detected in 944 MHz radio continuum images obtained with Phased Array Feeds on the Australian Square Kilometre Array Pathfinder (ASKAP). The more prominent, north-western shell may be associated with an X-ray detection, while the weaker, south-eastern shell has no counterpart at non-radio wavelength. We propose outwards moving shocks from galaxy mergers driving into the intragroup medium, re-energising relic radio lobes, as a possible formation scenario for the observed radio shells. We conclude that at least some ORCs are shock-energised relics in the outskirts of galaxy groups, which originate during the merger evolution of the brightest group galaxy.
We calibrate the Tully-Fisher relation (TFR) with data from the DESI Peculiar Velocity (PV) Survey taken during the Survey Validation (SV) period of the DESI galaxy redshift survey. Placing spectroscopic fibers on the centers and major axes of spatially-extended spiral galaxies identified in the 2020 Siena Galaxy Atlas using the DESI Legacy Surveys, we measure the rotational velocities at 0.33R26 for 1163 (1136 + 27 dwarf) spiral galaxies observed during SV. Using 41 spiral galaxies observed in the Coma Cluster, we find a slope for the TFR of -7.96+/-0.13 AB mag in the r-band, with a scatter about the TFR of 1.07+/-0.02 AB mag. We calibrate the zero-point of the TFR using galaxies with independent distances measured using type Ia supernovae via the cosmological distance ladder. From the SN Ia distances, we measure a zero-point of -19.34(+0.30,-0.29) AB mag in the r-band. We produce a public catalog of the distances to these 1136 spiral galaxies observed during DESI SV as part of the DESI PV Survey with our calibrated TFR. This is, to our knowledge, the first catalog of TFR distances produced with velocities measured at a single point in the disk.
We present a view of the stellar halo in the inner-central regions of the Milky Way (R <~ 10 kpc) mapped by RR Lyrae stars. The combined BRAVA-RR/APOGEE RR Lyrae catalog is used to obtain a sample of 281 RR Lyrae stars located in the bulge region of the Galaxy, but with orbits indicating they belong to the inner-central halo. The RR Lyrae stars in the halo are more metal-poor than the bulge RR Lyrae stars and have pulsation properties more consistent with an accreted population. We use the Milky Way-like zoom-in cosmological simulation Auriga to compare the properties of the RR Lyrae stars to those expected from the "Gaia-Enceladus-Sausage" (GES) merger. The integrals of motions and eccentricities of the RR Lyrae stars are consistent with a small fraction of 6-9 +- 2 % of the inner-central halo RR Lyrae population having originated from GES. This fraction, lower than what is seen in the solar neighborhood, is consistent with trends seen in the Auriga simulation, where a GES-like merger would have a decreasing fraction of GES stars at small Galactocentric radii compared to other accreted populations. Very few of the Auriga inner Galaxy GES-18 particles have properties consistent with belonging to a bulge population with (z_max < 1.1 kpc), indicating that no (or very few) RR Lyrae stars with bulge orbits should have originated from GES.
Modern digital sky surveys have been acquiring images of billions of galaxies. While these images often provide sufficient details to analyze the shape of the galaxies, accurate analysis of such high volumes of images requires effective automation. Current solutions often rely on machine learning annotation of the galaxy images based on a set of pre-defined classes. Here we introduce a new approach to galaxy image analysis that is based on generative AI. The method simplifies the galaxy images and automatically converts them into a ``skeletonized" form. The simplified images allow accurate measurements of the galaxy shapes and analysis that is not limited to a certain pre-defined set of classes. We demonstrate the method by applying it to galaxy images acquired by the DESI Legacy Survey. The code and data are publicly available. The method was applied to 125,000 DESI Legacy Survey images, and the catalog of the simplified images is publicly available.
(abridged) We consider a secular orbital evolution of a supermassive binary black hole (SBBH) with unequal masses $M_p$ and $M_s < M_p$ in a central part of a non-spherical nuclear star cluster (NSC). When the mass of NSC inside the orbit is smaller than $M_{s}$ dynamical friction becomes inefficient. The subsequent orbital evolution of SBBH is largely governed by perturbing tidal potential of NSC arising from its non-sphericity. When the perturbing potential is mainly determined by quadrupole harmonics with azimuthal number $|m|=2$ the secular dynamics of the SBBH does not conserve any components of the angular momentum and can lead to the formation of highly eccentric orbits. Such orbits can experience an efficient circularization due to emission of gravitational waves (GW). In this Paper we consider this situation in some detail. We study analytically and numerically the orbital evolution taking into account the important effect of Einstein apsidal precession and estimate the largest possible value of $e$, which can be obtained. We then estimate a possibility of fast orbital circularization through emission of gravitation waves on a highly eccentric orbit, with circularization timescale of the order of the orbital period. We find that our mechanism could result in such events on a time scale of the order of or smaller than a few Gyr. It is stressed that for particular values of the model parameters such events may sometimes be distinguished from the ones expected in more standard scenarios, since in our case the eccentricity may remain substantial all the way down to the final merger. It is also noted that our results can be applied to other astrophysical settings, e.g. to study the orbital evolution of a binary star or a proto planetary system inside a massive deformed gas cloud.
Quasar sightline observations reveal that low-ionization-state gas corotates with the galaxy disk and often at sub-centrifugal velocities, suggesting that the gas is spiraling towards the galaxy disk. However, while observations ubiquitously detect O VI absorption around low-redshift, $\sim L^*$ star-forming galaxies, the relationship between O VI and the galaxy disk, especially the kinematics, is not well-established. This work focuses on the O VI kinematics and its comparison with that of the low ions and galactic disk rotation. We present observations of 18 pairs of quasars and $z\approx0.2$ star-forming galaxies. All quasar sightlines intersect the circumgalactic medium (CGM) within 45$^\circ$ from the galaxy major axes. We show that while individual O VI velocity components do not correlate with disk rotation, the bulk of O VI gas in individual sightlines rarely counter-rotates. We then match O VI velocity components with those of low ions by minimizing the difference of their velocity centroids. The O VI velocity components with successful low-ion matches are typically found at small sightline impact parameters and are more likely to corotate with the disk. We suggest that the low-ion-matched O VI velocity components trace the gas co-spatial with the low ions near the extended disk plane in the inner CGM, whereas those without low-ion matches represent the gas at large 3D radii. While the gas at large radii is theoretically expected to kinematically correlate with the disk angular momentum, this correlation is expected to be weaker due to the higher turbulent to mean rotation velocity ration at large radii, consistent with our results.
We investigate the chemical evolution of N/O using a sample of 45 local star-forming galaxies (SFGs) from the CLASSY survey. This sample spans a wide range of galaxy properties, with robust determinations of nitrogen and oxygen abundances via the direct-$T_{\rm e}$ method. We explore how N/O relates to density structure, stellar mass, star formation rate (SFR), stellar age, compactness, and gas kinematics. In addition, we compare our results with those of galaxies at $z =2-10$ where N/O ratios were derived from optical or UV nitrogen lines, aiming to identify chemical enrichment pathways across cosmic time. Our analysis shows that the N/O-O/H relation in CLASSY galaxies aligns with the trends seen in local galaxies and extragalactic HII regions, and that galaxies at $z = 2-6$ exhibit similar N/O values, indicating no significant redshift evolution in N/O. We identify a significant correlation between electron density $n_{\rm e}$([S II]) and N/O, suggesting that density structure contributes to the scatter in the N/O-O/H relation. The CLASSY galaxies with high SFRs or compact star formation show elevated N/O, though no strong correlation with stellar mass is found. We also find that high-velocity outflows (v$_{out}$ > 350 km/s) and low mass-loading factors are linked to elevated N/O, indicating that feedback plays a significant role. These results highlight the importance of density, star formation, and feedback from young stellar populations in shaping N/O enrichment and provide key insights for interpreting high-$z$ galaxies observed with JWST.
Extra-tidal features around globular clusters (GCs) are tracers of their disruption, stellar stream formation, and their host's gravitational potential. However, these features remain challenging to detect due to their low surface brightness. We conduct a systematic search for such features around 19 GCs in the DECam Local Volume Exploration (DELVE) survey Data Release 2, discovering a new extra-tidal envelope around NGC 5897 and find tentative evidence for an extended envelope surrounding NGC 7492. Through a combination of dynamical modeling and analyzing synthetic stellar populations, we demonstrate these envelopes may have formed through tidal disruption. We use these models to explore the detectability of these features in the upcoming Legacy Survey of Space and Time (LSST), finding that while LSST's deeper photometry will enhance detection significance, additional methods for foreground removal like proper motions or metallicities may be important for robust stream detection. Our results both add to the sample of globular clusters with extra-tidal features and provide insights on interpreting similar features in current and upcoming data.
We present the first results from the COS-EDGES survey, targeting the kinematic connection between the ISM and multi-phase circumgalactic medium (CGM) in nine isolated, edge-on galaxies at z~0.2, each probed along its major axis by a background quasar at impact parameters of 13-38kpc. Using VLT/UVES and HST/COS quasar spectra, we analyse MgI, MgII, HI, CII, CIII, and OVI absorption relative to galaxy rotation curves from Keck/LRIS and Magellan/MagE spectra. We find that at lower $D/R_{vir}$ ($D/R_{vir}\leq 0.2$), over 80% of absorption in all ions lies on the side of systemic velocity matching disk rotation, and the optical-depth-weighted median velocity ($v_{abs}$) is consistent with the peak rotation speed. At higher $D/R_{vir}$ ($D/R_{vir} > 0.2$), the kinematics diverge by ionisation state: For low ionisation gas, the amount of co-rotating absorption remains >80%, yet $v_{abs}$ drops to 60% of the galaxy rotation speed. For high ionisation gas (OVI), only 60% of the absorption is consistent with co-rotation and $v_{abs}$ drops to 20% of the rotation speed. Furthermore, the velocity widths, corresponding to 50% of the total optical depth ($\Delta v_{50}$) for low ionisation gas is 1.8 times larger in the inner halo than at larger radii, while for CIII and OVI $\Delta v_{50}$ remains unchanged with distance. These results suggest a radially dependent CGM kinematic structure: the inner halo hosts cool, dynamically broad gas tightly coupled to disk rotation, whereas beyond 0.2$R_{vir}$, particularly traced by OVI and HI, the CGM shows weaker rotational alignment and lower velocity dispersion. Therefore, low-ionisation gas likely traces extended co-rotating gas, inflows and/or recycled accretion, while high-ionisation gas reflects a mixture of co-rotating, lagging, discrete collisionally ionised structures, indicating a kinematic stratification of the multi-phase CGM. [Abridged]
Atacama Large Millimeter/submillimeter Array (ALMA) observations at 1.3mm have recently revealed surprising complexity in the circumstellar environment of DFK 52, a red supergiant (RSG) located in the Stephenson 2 massive open cluster. We provide an initial characterisation of the star's mass-loss properties by studying its circumstellar emission in continuum, $^{12}$CO, $^{13}$CO, and SiO rotational lines. We find that DFK 52 is surrounded by an extremely large outflow (up to 50,000 au in radius) that shows complex morphologies in both its molecular and dust emission. The size of the circumstellar medium is unprecedented, even when compared with other known extreme RSGs, and its lower luminosity indicates that its mass ejection mechanism may be unique among this population. The molecular emission can be partially reproduced by a two-component model consisting of a fast (27 km/s) detached equatorial component with $M{\sim}0.05$ $M_{\odot}$ and a slow (10 km/s) spherical envelope with $\dot{M}\sim3\times10^{-6}$ $M_{\odot}$ yr$^{-1}$. This suggests that DFK 52 underwent a dramatic mass-loss event $\sim$4000 years ago, but has since transitioned into having a slower more symmetric mass loss. We conservatively estimate a total mass of $0.1-1$ $M_{\odot}$ in the complex extended regions of the outflow. The uncertain nature of the dramatic mass loss warrants extensive follow-up of this likely supernova progenitor.
High-redshift active galactic nuclei (AGN) have long been recognized as key probes of early black hole growth and galaxy evolution. However, modeling this population remains difficult due to the wide range of luminosities and black hole masses involved, and the high computational costs of capturing the hydrodynamic response of gas and evolving radiation fields on-the-fly. In this study, we present a new suite of simulations based on the IllustrisTNG galaxy formation framework, enhanced with on-the-fly radiative transfer, to examine AGN at high redshift (z > 5) in a protocluster environment extracted from the MillenniumTNG simulation. We focus on the co-evolution of black holes and their host galaxies, as well as the radiative impact on surrounding intergalactic gas. The model predicts that black holes form in overdense regions and lie below the local black hole-stellar mass relation, with stellar mass assembly preceding significant black hole accretion. Ionizing photons are primarily produced by stars, which shape the morphology of ionized regions and drive reionization. Given the restrictive black hole growth in the original IllustrisTNG model, we reduce the radiative efficiency from 0.2 to 0.1, resulting in higher accretion rates for massive black holes, more bursty growth, and earlier AGN-driven quenching. However, the resulting AGN remain significantly fainter than observed high-redshift quasars. As such, to incorporate this missing population, we introduce a quasar boosted model, in which we artificially boost the AGN luminosity. This results in strong effects on the surrounding gas, most notably a proximity effect, and large contributions to He ionization.
Emission and absorption line features are important diagnostics for the physics underlying extragalactic astronomy. The interpretation of observed signatures involves comparing against forward modeled spectra from galaxy formation simulations as well as more simplified geometries, while including the complex scattering radiative transfer (RT) of resonant emission lines. Here, we present thor, a modern C++ radiative transfer code focused initially on resonant emission lines. thor is a high-performance, distributed memory MPI-parallel, multi-target code, running on CPUs, GPUs and other accelerators, yielding large $\sim 10-50\rm{x}$ speed-ups compared to previous CPU-only codes. We support multiple grid-based and gridless data structures, enabling comparisons across different hydrodynamical codes as well as toy model geometries. We demonstrate its science capabilities with a number of example use cases across scales: (i) Lyman-alpha RT on simple shell-like gas distributions; (ii) Lyman-alpha RT applied to a high-resolution, high-redshift $z \sim 6$ cosmological hydrodynamical galaxy formation simulation; (iii) Lyman-alpha and Magnesium-II halos, i.e. scattering and emission from the circumgalactic medium of $z \simeq 1-2$ galaxies drawn from cosmological magnetohydrodynamical simulations; (iv) the large-scale cosmic web in gas emission, a $6144^3$ volume-element RT scaling calculation; and (v) synthetic absorption spectra of the Lyman-alpha forest. Extensive verification and benchmarking validates our approach and its computational efficiency.
The chemical abundance patterns of near-pristine objects provide important constraints on the properties of the first generations of stars in the Universe. We present the chemical abundances of five very metal-poor ([M/H]<-2.5) sub damped Lyman alpha systems (subDLAs) covering the redshift range $3.6<z<4.3$, identified with the XQ-100 survey. We find that the subDLAs in our sample show consistent chemical abundance patterns (in particular [C/O], [Al/O], and [Fe/O]) with those of very metal-poor DLAs. Based on Voigt profile fitting, the chemical abundance ratios [C/O], [Al/O], and [Si/O] of individual velocity components in at least three of the subDLAs shows some intrinsic scatter. In order to verify these chemical inhomogeneities in absorption components, we present a novel method for computing ionization corrections (ICs) on a component-by-component basis and show that ICs alone cannot explain the variations in [C/O], [Al/O], and [Si/O] between components of the same absorber at ~2 sigma significance. Comparing the observed abundance ratios to the simulated core-collapse supernovae yields of early stellar populations, we find that all individual components of the subDLAs appears to be enriched by progenitor masses of <30 M_sol. The observed inhomogeneities between components can be reproduced by differences in the progenitor mass or supernova explosion energy. As such, the observed chemical inhomogeneities between components can be explained by poorly mixed gas from different nucleosynthetic events.
We present a detailed kinematic study of a sample of 32 massive ($9.5\leqslant\log(M_*/{\rm M_{\odot}})\leqslant10.9$) main-sequence star-forming galaxies (MS SFGs) at $4<z<6$ from the ALMA-CRISTAL program. The data consist of deep (up to 15hr observing time per target), high-resolution ($\sim1$kpc) ALMA observations of the [CII]158$\mu$m line emission. This data set enables the first systematic kpc-scale characterisation of the kinematics nature of typical massive SFGs at these epochs. We find that $\sim50\%$ of the sample are disk-like, with a number of galaxies located in systems of multiple components. Kinematic modelling reveals these main sequence disks exhibit high-velocity dispersions ($\sigma_0$), with a median disk velocity dispersion of $\sim70{\rm kms^{-1}}$ and $V_{\rm rot}/\sigma_0\sim2$, and consistent with dominant gravity driving. The elevated disk dispersions are in line with the predicted evolution based on Toomre theory and the extrapolated trends from $z\sim0$-$2.5$ MS star-forming disks. The inferred dark matter (DM) mass fraction within the effective radius $f_{\rm DM}(<R_{\rm e})$ for the disk systems decreases with the central baryonic mass surface density, and is consistent with the trend reported by kinematic studies at $z\lesssim3$; roughly half the disks have $f_{\rm DM}(<R_{\rm e})\lesssim30\%$. The CRISTAL sample of massive MS SFGs provides a reference of the kinematics of a representative population and extends the view onto typical galaxies beyond previous kpc-scale studies at $z\lesssim3$.
Interstellar Objects are comets and asteroids that formed around other stars but were ejected before they could accrete into exoplanets. They therefore represent a rare opportunity to compare the building blocks of planets in the Solar System to those in other stellar systems. The third Interstellar Object, 3I/ATLAS, is the newest, brightest, potentially largest, and fastest member of this population. We report observations of 3I/ATLAS taken on 2025 July 3 and 4 with the NASA Infrared Telescope Facility just days after its discovery. In r'-band imaging with 'Opihi, we see no obvious lightcurve variability and derive a g'-i' color of 1.06+/-0.11 which is consistent in spectral slope to recent work by D.Z. Seligman et al. (2025) and C. Opitom (2025). We obtained the first near-infrared (NIR) reflectance spectrum of 3I/ATLAS with SpeX. The visible color and NIR spectrum show a linear, red visible slope, a somewhat less red slope between 0.7 and 1.1 $\mu{m}$, and a neutral or slightly blue slope at longer wavelengths. Challenges in modeling the reflectivity of 3I may indicate that this comet has a complex grain size distribution, grain compositions unlike Solar system comets, or both. Like 2I/Borisov, there are no obvious signatures of water ice in the coma of 3I/ATLAS. Observations closer to perihelion will help elucidate whether 3I has less water than anticipated or whether the Interstellar Objects might retain and release their ices somewhat differently than Solar System comets do.
At this early stage of its passage through our Solar System, 3I/ATLAS, the recently discovered interstellar interloper, has displayed various anomalous characteristics, determined from photometric and astrometric observations. As largely a pedagogical exercise, in this paper we present additional analysis into the astrodynamics of 3I/ATLAS, and hypothesize that this object could be technological, and possibly hostile as would be expected from the 'Dark Forest' resolution to the 'Fermi Paradox'. We show that 3I/ATLAS approaches surprisingly close to Venus, Mars and Jupiter, with a probability of $\lesssim 0.005$\%. Furthermore the low retrograde tilt of 3I/ATLAS's orbital plane to the ecliptic offers various benefits to an Extra-terrestrial Intelligence (ETI), since it allows the object access to our planet with relative impunity. The eclipse by the Sun from Earth of 3I/ATLAS at perihelion, would allow it to conduct a clandestine reverse Solar Oberth Manoeuvre, an optimal high-thrust strategy for interstellar spacecraft to brake and stay bound to the Sun. An optimal intercept of Earth would entail an arrival in late November/early December of 2025, and also, a non-gravitational acceleration of $\sim{5.9} \times 10^{-5}$ au day$^{-2}$, normalized at 1 au from the Sun, would indicate an intent to intercept the planet Jupiter, not far off its path, and a strategy to rendezvous with it after perihelion.
The $\sim 60\,000$ solar-mass (\MSun) star-cluster R136 (NGC~2070) in the Tarantula Nebula in the Large Magellanic Cloud is the host of at least 55 massive stars ($M \apgt 10$\,\MSun) which move away from the cluster at projected velocities $\gteq 27.5$\,km/s \cite{2024Natur.634..809S}. The origin of the high velocities of such runaway stars have been debated since the 1960s, resulting either from dynamical ejections \citep{1961BAN....15..265B,1961BAN....15..291B} or from supernova explosions \citep{1983ApJ...267..322H}. Due to the Gaia satellite's outstanding precision, we can now retrace the most recently ejected binary star, Mel 34, back to the center of R136 and reconstruct the events that 52\,000 years ago let to its removal from R136, i.e., we establish its dynamical interaction and ejection history. We find that this ejection requires the participation of 5 stars in a strong interaction between a triple composed of the tight massive binary Mel~39 orbited by the star VFTS~590, and the binary star Mel~34. The participation of 5 stars is unexpected because runaway stars were not expected to result from triple interactions \cite{2011Sci...334.1380F}. The deterministic nature of the Newtonian dynamics in the scattering enables us to reconstruct the encounter that ejected Mel~34. We then predict that Mel~39 is a binary star with an 80\,\Msun\, companion star that orbits within $\sim 1^\circ$ in the same plane as Mel~34, and escapes the cluster with a velocity of $\sim 64$\,km/s. The five stars will undergo supernova explosions in the coming 5\,Myr at a distance of $\sim 180$\,pc to $\sim 332$\,pc from their birth location (R\,136). The resulting black hole binaries, however, are not expected to merge within a Hubble time.