Main Galaxy Research Articles

Angular structure and gravitational imaging

AbstractIn gravitational imaging, the mass model for the main lensing galaxy is one of the main sources of systematic uncertainty. We use subhalo detection models with increasing levels of angular complexity in the lens mass model to analyse 100 HST mock observations. We find that perturbations of just 1% are enough to cause a 20% false positive subhalo detection rate, with order 3 multipoles having the strongest effect. The area in an observation where a substructure can be detected drops by a factor of 10 if multipoles up to 3 per cent amplitude are included in the lens model. The mass of the smallest detectable substructure however is not affected. We find a detection limit of M>108.2M⊙ at 5σ in all models. In order for strong lensing searches for dark matter objects to remain reliable in the future, angular structure beyond the elliptical power-law must be included.

Kinematics of the Local Group gas and galaxies in the hestia simulations

ABSTRACT We investigate the kinematic properties of gas and galaxies in the Local Group (LG) using high-resolution simulations performed by the hestia (High-resolution Environmental Simulations of The Immediate Area) collaboration. Our simulations include the correct cosmography surrounding LG-like regions consisting of two main spiral galaxies of ∼1012 M⊙, their satellites and minor isolated galaxies, all sharing the same large-scale motion within a volume of a few Mpc. We characterize the gas and galaxy kinematics within the simulated LGs, from the perspective of the Sun, to compare with observed trends from recent HST/COS absorption-line observations and LG galaxy data. To analyse the velocity pattern of LG gas and galaxies seen in the observational data, we build sky maps from the local standard of rest, and the Galactic and LG barycentre frames. Our findings show that the establishment of a radial velocity dipole at low/high latitudes, near the preferred barycentre direction, is a natural outcome of simulation kinematics for material outside the Milky Way virial radius after removing Galaxy rotation when the two main LG galaxies are approaching. Our results favour a scenario where gas and galaxies stream towards the LG barycentre producing a velocity dipole resembling observations. While our study shows in a qualitative way the global matter kinematics in the LG as part of its ongoing assembly, quantitative estimates of gas-flow rates and physical conditions of the LG gas have to await a more detailed modelling of the ionization conditions, which will be presented in a follow-up paper.

Cold gas mass measurements for the era of large optical spectroscopic surveys

ABSTRACT Gas plays an important role in many processes in galaxy formation and evolution, but quantifying the importance of gas has been hindered by the challenge to measure gas masses for large samples of galaxies. Data sets of direct atomic and molecular gas measurements are sufficient to establish simple scaling relations, but often not large enough to quantify three-parameter relations, or second-order dependences. As an alternative approach, we derive here indirect cold gas measurements from optical emission lines using photoionization models for galaxies in the Sloan Digital Sky Survey (SDSS) main galaxy sample and the PHANGS-MUSE survey. We calibrate the gas surface density measurements using xCOLD GASS and PHANGS-ALMA molecular gas measurements to ensure that our measurements are reliable. We demonstrate the importance of taking into account the scale dependence of the relation between optical depth (τV) and gas surface density (Σgas) and provide a general prescription to estimate Σgas from τV, metallicity, and the dust-to-metal ratio, at any arbitrary physical resolution. To demonstrate that the indirect cold gas masses are accurate enough to quantify the role of gas in galaxy evolution, we study the mass–metallicity relation of SDSS galaxies and show that as a third parameter gas mass is better than star formation rate at reducing the scatter of the relation, as predicted by models and simulations.

Deciphering the radio–star formation correlation on kpc scales

Given the multiple energy-loss mechanisms of cosmic-ray (CR) electrons in galaxies, the tightness of the infrared (IR)–radio continuum correlation is surprising. As the radio continuum emission at GHz frequencies is optically thin, this offers the opportunity to obtain unbiased star formation rates (SFRs) from radio-continuum flux-density measurements. The calorimeter theory can naturally explain the tightness of the far-infrared (FIR)–radio correlation but makes predictions that do not agree with observations. Noncalorimeter models often have to involve a conspiracy to maintain the tightness of the FIR–radio correlation. We extended a published analytical model of galactic disks by including a simplified prescription for the synchrotron emissivity. The galactic gas disks of local spiral galaxies, low-z starburst galaxies, high-z main sequence star-forming galaxies, and high-z starburst galaxies are treated as turbulent clumpy accretion disks. The magnetic field strength is determined by the equipartition between the turbulent kinetic and the magnetic energy densities. Our fiducial model, which includes neither galactic winds nor CR electron secondaries, reproduces the observed radio continuum spectral energy distributions of most (∼70%) of the galaxies. Except for the local spiral galaxies, fast galactic winds can potentially make the conflicting models agree with observations. The observed IR–radio correlations are reproduced by the model within 2σ of the joint uncertainty of model and data for all datasets. The model agrees with the observed SFR–radio correlations within ∼4σ. Energy equipartition between the CR particles and the magnetic field only approximately holds in our models of main sequence star-forming galaxies. If a CR electron calorimeter is assumed, the slope of the IR–radio correlation flattens significantly. Inverse Compton losses are not dominant in the starburst galaxies because in these galaxies not only the gas density but also the turbulent velocity dispersion is higher than in normal star-forming galaxies. Equipartition between the turbulent kinetic and magnetic field energy densities then leads to very high magnetic field strengths and very short synchrotron timescales. The exponents of our model SFR–radio correlations at 150 MHz and 1.4 GHz are very close to one.

Modeling Redshift-space Clustering with Abundance Matching

We explore the degrees of freedom required to jointly fit projected and redshift-space clustering of galaxies selected in three bins of stellar mass from the Sloan Digital Sky Survey Main Galaxy Sample (SDSS MGS) using a subhalo abundance matching (SHAM) model. We employ emulators for relevant clustering statistics in order to facilitate our analysis, leading to large speed gains with minimal loss of accuracy. We are able to simultaneously fit the projected and redshift-space clustering of the two most massive galaxy samples that we consider with just two free parameters: scatter in stellar mass at fixed SHAM proxy, and the dependence of the SHAM proxy on dark matter halo concentration. We find some evidence for models that include velocity bias, but including orphan galaxies improves our fits to the lower-mass samples significantly. We also model the clustering signals of specific star formation rate (sSFR) selected samples using conditional abundance matching (CAM). We obtain acceptable fits to projected and redshift-space clustering as a function of sSFR and stellar mass using two CAM variants, although the fits are worse than for stellar-mass-selected samples alone. By incorporating nonunity correlations between the CAM proxy and sSFR, we are able to resolve previously identified discrepancies between CAM predictions and SDSS observations of the environmental dependence of quenching for isolated central galaxies.

Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling

Context. Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions (SEDs) into the main luminosity contributors (e.g. stellar populations of distinct age and metallicity or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years, with unknown ramifications to the current understanding of galaxy evolution. Aims. The objective of this work is to estimate the relations between stellar properties (e.g. total mass, mean age, and mean metallicity) of SF galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through the more common purely stellar spectral synthesis approach. Methods. The main galaxy sample from SDSS DR7 was analysed with two distinct population synthesis codes: FADO, which estimates self-consistently both the stellar and nebular contributions to the SED, and the original version of STARLIGHT, as representative of purely stellar population synthesis codes. Results. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as ∼0.06 dex for the overall population of galaxies and ∼0.12 dex for SF galaxies (galaxies with EW(Hα) > 3 Å), with the most prominent difference between both codes in the two populations being in the light-weighted mean stellar age. FADO presents a broader range of mean stellar ages and metallicities for SF galaxies than STARLIGHT, with the latter code preferring metallicity solutions around the solar value (Z⊙ = 0.02). A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(Hα) > 75 Å) reveals that the light contributions of simple stellar populations (SSPs) younger than ≤107 (109) years in STARLIGHT are higher by ∼5.41% (9.11%) compared to FADO. Moreover, FADO presents higher light contributions from SSPs with metallicity ≤Z⊙/200 (Z⊙/50) of around 8.05% (13.51%) when compared with STARLIGHT. This suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to ∼0.17 dex and overestimating the light-weighted metallicity by up to ∼0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(Hα) and decreasing total stellar mass. Moreover, comparing SF results from FADO in a purely stellar mode with the previous results qualitatively suggests that differences between codes are primarily due to mathematical and statistical differences and secondarily due to the impact of the nebular continuum modelling approach (or lack thereof). However, it is challenging to adequately quantify the relative role of each factor since they are likely interconnected. Conclusions. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties (e.g. mean age and mean metallicity) of galaxies with relatively high star formation rates (e.g. dwarf spirals, ‘green peas’, and starburst galaxies). In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity, which are factors worth bearing in mind in light of future high-resolution spectroscopic surveys at higher redshifts (e.g. MOONS and 4MOST-4HS).

An extended [C ii] halo around a massive star-forming galaxy at z = 5.3

ABSTRACT High-redshift observations are often biased towards massive and bright galaxies that are not necessarily representative of the full population. In order to accurately study galaxy evolution and mass assembly at these redshifts, observations of ‘normal’ main sequence galaxies are required. Here we present Atacama Large Millimeter/Submillimeter Array (ALMA) 0.3 arcsec resolution observations of the [C ii] emission line at 158 μm of HZ7, a main sequence galaxy at z = 5.25. Comparing to archival rest-frame UV observations taken by the Hubble Space Telescope (HST), we find strong evidence of the existence of extended [C ii] emission, which we estimate to be twice the size of the rest-frame UV emission, yielding one of the first high-redshift objects where a clear signature of a [C ii] ‘Halo’ has been detected to date. For a matched Sérsic profile with n = 1, we measured a [C ii] effective radius of 0.50 ± 0.04 arcsec (3.07 ± 0.25 kpc) and an average rest-frame UV effective radius of 0.2 ± 0.04 arcsec (1.48 ± 0.16 kpc). The [C ii] morphology and kinematics of the system suggest a merging event resulting in a non-rotating disc system. This event could be responsible for the extended [CII] emission. Alternatively, some potential obscured emission could also explain the [C ii] to UV size ratio. These results contribute to the growing consensus with respect to the existence of extended [C ii] emission around galaxies.

Shaping physical properties of galaxy subtypes in the VIPERS survey: Environment matters

Aims. This study aims to explore the relation between the physical properties of different galaxy sub-classes, from red passive to blue star-forming, and their environment. Our work is based on the analysis of 31 631 galaxies from the VIMOS Public Extragalactic Redshift Survey (VIPERS), observed at 0.5 < z < 0.9. The unprecedented volume of VIPERS and the wealth of auxiliary derived data allow us to associate sub-classes of the main galaxy populations with their possibly different evolutionary paths. This is the first time such a study is performed with such statistical precision. Methods. We use the results of an unsupervised clustering algorithm to distinguish 11 subclasses of VIPERS galaxies, based on the multi-dimensional feature space, defined by rest-frame UV to NIR colours. We investigate the relationship between the properties of these sub-classes of galaxies and their local environment, defined as the galaxy density contrast, δ, derived from the fifth nearest neighbour technique. Results. We confirm that the galaxy population-density relation is already in place at z ∼ 0.9, with the blue galaxy fraction decreasing with density, compensated by an increase in the red fraction. We demonstrate how the properties of red, green, and blue galaxy subclasses are altered as they assemble into denser regions, and we attempt to interpret this in the context of their evolution. On average, red galaxies in the high-density environment are larger by 28% than those in low-density environments. In particular, we find one group of galaxies, sub-class C3, whose increase of size with time can be explained mainly as the result of mergers; for other red subclasses, mergers would not seem to play the major role (subclass C2) or would play a negligible role (sub-class C1). The properties of the green galaxies (sub-classes C4–6) depend on whether their stellar mass is above or below a transition mass, log(Mstar/M⊙)=10.6. Low-mass green (9.5 ≲ log(Mstar/M⊙)≲10.6) galaxies appear to have grown through secular processes, while in high-mass (10.6 ≲ log(Mstar/M⊙)≲11.5) green galaxies, mass assembly appears to be dominated by mergers. When it comes to blue galaxies, the trend of decreasing fraction with denser environments seen for the group as a whole (sub-classes C7–11) is found to be driven mostly by one (the most numerous) group of galaxies; sub-class C10. These are compact low-mass galaxies with high, specific star formation rates, which are preferentially found in low-density environments. However, the remaining blue galaxies (sub-classes C7–9) are larger and appear in denser environments than galaxies within C10.

The Chocolate Chip Cookie Model: Dust Geometry of Milky Way–like Disk Galaxies

We present a new two-component dust geometry model, the Chocolate Chip Cookie model, where the clumpy nebular regions are embedded in a diffuse stellar/interstellar medium disk, like chocolate chips in a cookie. By approximating the binomial distribution of the clumpy nebular regions with a continuous Gaussian distribution and omitting the dust scattering effect, our model solves the dust attenuation process for both the emission lines and stellar continua via analytical approaches. Our Chocolate Chip Cookie model successfully fits the inclination dependence of both the effective dust reddening of the stellar components derived from stellar population synthesis and that of the emission lines characterized by the Balmer decrement for a large sample of Milky Way–like (MW-like) disk galaxies selected from the main galaxy sample of the Sloan Digital Sky Survey. Our model shows that the clumpy nebular disk is about 0.55 times thinner and 1.6 times larger than the stellar disk for MW-like galaxies, whereas each clumpy region has a typical optical depth of τ cl,V ∼ 0.5 in the V band. After considering the aperture effect, our model prediction on the inclination dependence of dust attenuation is also consistent with observations. Not only that, in our model, the dust attenuation curve of the stellar population naturally depends on the inclination, and its median case is consistent with the classical Calzetti law. As the modeling constraints are from the optical wavelengths, our model is unaffected by the optically thick dust component, which however could bias the model’s prediction of the infrared emissions.

Measurements of cosmic expansion and growth rate of structure from voids in the Sloan Digital Sky Survey between redshift 0.07 and 1.0

ABSTRACT We present measurements of the anisotropic cross-correlation of galaxies and cosmic voids in data from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample, Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS luminous red galaxy catalogues from SDSS Data Releases 7, 12, and 16, covering the redshift range 0.07 < z < 1.0. As in our previous work analysing voids in subsets of these data, we use a reconstruction method applied to the galaxy data before void finding in order to remove selection biases when constructing the void samples. We report results of a joint fit to the multipole moments of the measured cross-correlation for the growth rate of structure, fσ8(z), and the ratio DM(z)/DH(z) of the comoving angular diameter distance to the Hubble distance, in six redshift bins. For DM/DH, we are able to achieve a significantly higher precision than that obtained from analyses of the baryon acoustic oscillations and galaxy clustering in the same data sets. Our growth rate measurements are of lower precision but still comparable with galaxy clustering results. For both quantities, the results agree well with the expectations for a Λ cold dark matter model. Assuming a flat Universe, our results correspond to a measurement of the matter density parameter $\Omega _\mathrm{m}=0.337^{+0.026}_{-0.029}$. For more general models, the degeneracy directions obtained are consistent with and complementary to those from other cosmological probes. These results consolidate void–galaxy cross-correlation measurements as a pillar of modern observational cosmology.

Elucidating galaxy assembly bias in SDSS

We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2\sigma$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.

Sub-percentage measure of distances to redshift of 0.1 by a new cosmic ruler

ABSTRACT Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675 ± 0.079 for flat Λ-cold dark matter (ΛCDM) model, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34 per cent, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-$\mathit{ z}$ measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7, and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

Arp 58 and Arp 68: two M 51-type systems

ABSTRACT We study two M 51-type systems Arp 68 and Arp 58, which strongly differ by their stellar masses, gas content, and environment. Long-slit spectral observations obtained at the Russian 6-m telescope were used to trace the distributions of a line-of-sight (LOS) velocity and a gas-phase oxygen abundance along the spectral cuts. Two systems are compared by their observed properties. We found a very strong large-scale non-circular motion of gas in both systems and a kpc-size saw-edged velocity profile along the tidal spiral arm of Arp 68, probably caused by the gas outflow due to the stellar feedback. A deep decrease of LOS velocity is also found in the ‘hinge’ region in Arp 58, where the inner spiral arm transforms into the tidal spiral arm, which was predicted earlier for M 51-type galaxies. Local sites of star formation and the satellites are compared with the evolutionary models at the colour–colour diagrams. Unlike the spiral galaxy Arp 58, the main galaxy in Arp 68 system is experiencing an ongoing burst of star formation. Gas-phase metallicity estimates show that Arp 58 has a higher metal abundance and reveals a shallow negative radial gradient of the gas-phase oxygen abundance. The emission gas in Arp 68 has noticeably lower metallicity than it is expected for a given luminosity of this galaxy, which may be connected with its space position in the local void.

Photometric redshifts from SDSS images with an interpretable deep capsule network

ABSTRACT Studies of cosmology, galaxy evolution, and astronomical transients with current and next-generation wide-field imaging surveys like the Rubin Observatory Legacy Survey of Space and Time are all critically dependent on estimates of photometric redshifts. Capsule networks are a new type of neural network architecture that is better suited for identifying morphological features of the input images than traditional convolutional neural networks. We use a deep capsule network trained on ugriz images, spectroscopic redshifts, and Galaxy Zoo spiral/elliptical classifications of ∼400 000 Sloan Digital Sky Survey galaxies to do photometric redshift estimation. We achieve a photometric redshift prediction accuracy and a fraction of catastrophic outliers that are comparable to or better than current methods for SDSS main galaxy sample-like data sets (r ≤ 17.8 and zspec ≤ 0.4) while requiring less data and fewer trainable parameters. Furthermore, the decision-making of our capsule network is much more easily interpretable as capsules act as a low-dimensional encoding of the image. When the capsules are projected on a two-dimensional manifold, they form a single redshift sequence with the fraction of spirals in a region exhibiting a gradient roughly perpendicular to the redshift sequence. We perturb encodings of real galaxy images in this low-dimensional space to create synthetic galaxy images that demonstrate the image properties (e.g. size, orientation, and surface brightness) encoded by each dimension. We also measure correlations between galaxy properties (e.g. magnitudes, colours, and stellar mass) and each capsule dimension. We publicly release our code, estimated redshifts, and additional catalogues at https://biprateep.github.io/encapZulate-1.

Substructure detection reanalysed: dark perturber shown to be a line-of-sight halo

ABSTRACT Observations of structure at subgalactic scales are crucial for probing the properties of dark matter, which is the dominant source of gravity in the universe. It will become increasingly important for future surveys to distinguish between line-of-sight haloes and subhalos to avoid wrong inferences on the nature of dark matter. We reanalyse a subgalactic structure (in lens JVAS B1938 + 666) that has been previously found using the gravitational imaging technique in galaxy-galaxy lensing systems. This structure has been assumed to be a satellite in the halo of the main lens galaxy. We fit the redshift of the perturber of the system as a free parameter, using the multiplane thin-lens approximation, and find that the redshift of the perturber is $z_\mathrm{int} = 1.42^{+0.10}_{-0.15}$ (with a main lens redshift of z = 0.881). Our analysis indicates that this structure is more massive than the previous result by an order of magnitude. This constitutes the first dark perturber shown to be a line-of-sight halo with a gravitational lensing method.

Spatially resolved gas-phase metallicity in FIRE-2 dwarfs: late-time evolution of metallicity relations in simulations with feedback and mergers

ABSTRACT We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($\rm \mathit{M}_{halo} \approx 10^{11}\, {\rm M}_\odot$, $\rm \mathit{M}_\star \approx 10^{8.8}{-}10^{9.6}\, {\rm M}_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near z ≈ 0, over a period of 1.4 Gyr, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass–metallicity (MZR) and mass–metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95 ${{\ \rm per\ cent}}$ of the gas being within ±0.1 dex between the cold and dense gas (T < 500 K and nH > 1 cm−3), ionized gas (near the H αT ≈ 104 K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy.

The X-shaped Radio Galaxy J0725+5835 is Associated with an AGN Pair

X-shaped radio galaxies (XRGs) are those that exhibit two pairs of unaligned radio lobes (main radio lobes and wings). One of the promising models for the peculiar morphology is jet reorientation. To clarify this, we conducted a 5 GHz observation with the European VLBI Network (EVN) of XRG J0725+5835, which resembles the archetypal binary active galactic nuclei (AGNs) 0402+379 in radio morphology, but it is larger in angular size. In our observation, two milliarcsecond-scale radio components with nonthermal radio emission are detected. Each of them coincides with an optical counterpart with similar photometric redshift and (optical and infrared) magnitude, corresponding to dual active nuclei. Furthermore, with the improved Very Large Array (VLA) images, we find a bridge between the two radio cores and a jet bending in the region surrounding the companion galaxy. This further supports the interplay between the main and companion galaxies. In addition, we also report the discovery of an arcsecond-scale jet in the companion. Given the projected separation of ∼100 kpc between the main and companion galaxies, XRG J0725+5835 is likely associated with a dual jetted-AGN system. In both EVN and VLA observations, we find signatures that the jet is changing its direction, which is likely responsible for the X-shaped morphology. For the origin of jet reorientation, several scenarios are discussed.

An Old, Metal-rich Accreted Stellar Component in the Milky Way Stellar Disk

We study the possibility that the Milky Ways’ cool stellar disk includes mergers with ancient stars. Galaxies are understood to form in a hierarchical manner, where smaller (proto-)galaxies merge into larger ones. Stars in galaxies, like the Milky Way, contain in their motions and elemental abundance tracers of past events and can be used to disentangle merger remnants from stars that formed in the main galaxy. The merger history of the Milky Way is generally understood to be particularly easy to study in the stellar halo. The advent of the ESA astrometric satellite Gaia has enabled the detection of completely new structures in the halo such as the Gaia-Enceladus-Sausage. However, simulations also show that mergers may be important for the build-up of the cool stellar disks. Combining elemental abundances for ∼100 giant branch stars from APOGEE DR17 and astrometric data from Gaia we use elemental abundance ratios to find a hitherto unknown, old stellar component in the cool stellar disk in the Milky Way. We further identify a small sample of RR Lyrae variables with disk kinematics that also show the same chemical signature as the accreted red giant stars in the disk. These stars allow us to date the stars in the accreted component. We find that they are exclusively old.

What is the origin of the stacked radio emission in radio-undetected quasars?

Radio emission in the brightest radio quasars can be attributed to processes inherent to active galactic nuclei (AGN) powered by super massive black holes (SMBHs), while the physical origins of the radio fluxes in quasars without radio detections have not been established with full certainly. Deep radio surveys carried out with the Low Frequency ARray (LOFAR) are at least one order of magnitude more sensitive for objects with typical synchrotron spectra than previous wide-area high-frequency surveys ( > 1.0 GHz). With the enhanced sensitivity that LOFAR offers, we investigate the radio-infrared continuum of LOFAR radio-detected quasars (RDQs) and LOFAR radio-undetected quasars (RUQs) in the 9.3 deg2NOAO Deep Wide-field survey (NDWFS) of the Boötes field; RUQs are quasars that are individually undetected at a level of ≥5σin the LOFAR observations. To probe the nature of the radio and infrared emission, where direct detection is not possible due to the flux density limits, we used a median image stacking procedure. This was done in the radio frequencies of 150 MHz, 325 MHz, 1.4 GHz and 3.0 GHz, and in nine infrared bands between 8 and 500 μm. The stacking analysis allows us to probe the radio-luminosity for quasars that are up to one order of magnitude fainter than the ones detected directly. The radio and infrared photometry allow us to derive the median spectral energy distributions of RDQs and RUQs in four contiguous redshift bins between 0 < z < 6.15. The infrared photometry is used to derive the infrared star-formation rate (SFR) through SED fitting, and is compared with two independent radio-based star-formation (SF) tracers using the far-infrared radio correlation (FIRC) of star-forming galaxies. We find a good agreement between our radio and infrared SFR measurements and the predictions of the FIRC. Moreover, we use the FIRC predictions to establish the level of the contribution due to SMBH accretion to the total radio-luminosity. We show that SMBH accretion can account for ∼5−41% of the total radio-luminosity in median RUQs, while for median RDQs the contribution is ∼50−84%. This implies that vigorous SF activity is coeval with SMBH growth in our median stacked quasars. We find that median RDQs have higher SFRs that agree well with those of massive star-forming main sequence galaxies, while median RUQs present lower SFRs than RDQs. Furthermore, the behavior of the radio-loudness parameter (R = log10(Lrad/LAGN)) is investigated. For quasars withR ≥ −4.5, the radio-emission is consistent with being dominated by SMBH accretion, while for low radio luminosity quasars withR < −4.5 the relative contribution of SF to the radio fluxes increases as the SMBH component becomes weaker. We also find signatures of SF suppression due to negative AGN feedback in the brightest median RDQs at 150 MHz. Finally, taking advantage of our broad spectral coverage, we studied the radio spectra of median RDQs and RUQs. The spectral indices of RUQs and RDQs do not evolve significantly with redshift, but they become flatter towards lower frequencies.

Physics of ULIRGs with MUSE and ALMA: The PUMA project

Context. A classical scenario suggests that ultra-luminous infrared galaxies (ULIRGs) transform colliding spiral galaxies into a spheroid-dominated early-type galaxy. Recent high-resolution simulations have instead shown that, under some circumstances, rotation disks can be preserved during the merging process or rapidly regrown after coalescence. Our goal is to analyse in detail the ionised gas kinematics in a sample of ULIRGs to infer the incidence of gas rotational dynamics in late-stage interacting galaxies and merger remnants. Aims. We analysed integral field spectrograph MUSE data of a sample of 20 nearby (z < 0.165) ULIRGs (with 29 individual nuclei) as part of the Physics of ULIRGs with MUSE and ALMA (PUMA) project. We used multi-Gaussian fitting techniques to identify gaseous disk motions and the 3D-Barolo tool to model them. Methods. We found that 27% (8 out of 29) individual nuclei are associated with kiloparsec-scale disk-like gas motions. The rest of the sample displays a plethora of gas kinematics, dominated by winds and merger-induced flows, which makes the detection of rotation signatures difficult. On the other hand, the incidence of stellar disk-like motions is ∼2 times larger than gaseous disks, as the former are probably less affected by winds and streams. The eight galaxies with a gaseous disk present relatively high intrinsic gas velocity dispersion (σ0 ∈ [30 − 85] km s−1), rotationally supported motions (with gas rotation velocity over velocity dispersion vrot/σ0 ∼ 1 − 8), and dynamical masses in the range (2 − 7)×1010 M⊙. By combining our results with those of local and high-z disk galaxies (up to z ∼ 2) from the literature, we found a significant correlation between σ0 and the offset from the main sequence (δMS), after correcting for their evolutionary trends. Results. Our results confirm the presence of kiloparsec-scale rotating disks in interacting galaxies and merger remnants in the PUMA sample, with an incidence going from 27% (gas) to ≲50% (stars). Their gas σ0 is up to a factor of ∼4 higher than in local normal main sequence galaxies, similar to high-z starbursts as presented in the literature; this suggests that interactions and mergers enhance the star formation rate while simultaneously increasing the velocity dispersion in the interstellar medium.