Search

We investigate the change in mean stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift ($0.60\leq z\leq0.76$) using the LEGA-C Survey, to low redshift ($0.014\leq z\leq0.10$) using the SAMI Galaxy Survey. We find that, similarly to their low-redshift counterparts, the stellar metallicity of quiescent galaxies at $0.60\leq z\leq 0.76$ closely correlates with $M_*/R_\mathrm{e}$ (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals, by determining the escape velocity required by metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. On the other hand, we find no correlation between stellar age and $M_*/R_\mathrm{e}^2$ (stellar mass surface density $\Sigma$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$\Sigma$ relation in the context of the redshift evolution of the star-forming and quiescent populations in the mass--size plane, and find our results can be explained as a consequence of galaxies forming more compactly at higher redshifts, and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The $z\sim 0$ age--$\Sigma$ relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and so a range of surface densities., Comment: 18 pages, 11 figures, accepted to MNRAS

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment., Comment: 20 pages, 9 figures. Accepted for publication in MNRAS

We present new CO(2–1) observations of a representative sample of 24 local (z, EB and AL acknowledge the support from Comunidad de Madrid through the Atracción de Talento grant 2017-T1/TIC-5213. This research has been partially funded by the Spanish State Research Agency (AEI) Project MDM-2017-0737 Unidad de Excelencia ‘María de Maeztu’– Centro de Astrobiología (INTA-CSIC). MPS acknowledges support from the Comunidad de Madrid through the Atracción de Talento Investigador Grant 2018-T1/TIC-11035 and PID2019-105423GA-I00 (MCIU/AEI FEDER, UE). AAH, MVM and SGB acknowledge support from grant PGC2018-094671-B-I00 (MCIU/AEI/FEDER,UE). SGB acknowledges support from the research project PID2019-106027GA-C44 of the Spanish Ministerio de Ciencia e Innovación. DR acknowledges support from STFC through grant ST/S000488/1. AL, LC and SA acknowledge the support from grant PID2019-106280GB-I00 (MCIU/AEI/FEDER,UE). SC acknowledges financial support from the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). MSG acknowledges support from the Spanish Ministerio de Economía y Competitividad through the grants BES-2016-078922 and ESP2017-83197-P. AP gratefully acknowledges financial support from STFC through grant ST/T000244/1. FV acknowledges support from the Carlsberg Foundation Research Grant CF18-0388 ‘Galaxies: Rise and Death’. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant No. 140. AU acknowledges support from the Spanish grants PGC2018-094671-B-I00, funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe”, and PID2019-108765GB-I00, funded by MCIN/AEI/10.13039/501100011033. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00243.S, ADS/JAO.ALMA#2013.1.00271.S, ADS/JAO.ALMA#2015.1.00714.S, ADS/JAO.ALMA#2017.1.00255.S, ADS/JAO.ALMA#2017.1.00395.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., It has been suggested that gravitational potential can have a significant role in suppressing star formation in nearby galaxies. To establish observational constraints on this scenario, we investigated the connection between the dynamics – taking the circular velocity curves (CVCs) as a proxy for the inner gravitational potential – and star formation quenching in 215 non-active galaxies across the Hubble sequence from the Calar Alto Legacy Integral Field Area (CALIFA) survey. Our results show that galaxies with similar CVCs tend to have a certain star-formation quenching pattern. To explore these findings in more details, we constructed kiloparsec(kpc)-resolved relations of the equivalent width of the Hα (WHα) versus the amplitude (Vc) and shape (β = dlnVc/dlnR) of the circular velocity at given radii. We find that the WHα − Vc is a declining relationship, where the retired regions of the galaxies (the ones with WHα values of below 3 Å) tend to have higher Vc. Concurrently, WHα − β is a bimodal relationship, which is characterised by two peaks: concentration of the star forming regions at a positive β (rising CVC) and a second concentration of the retired regions with a negative β (declining CVC). Our results show that both the amplitude of the CVC – driven by the mass of the galaxies – and its shape – which reflects the internal structure of the galaxies – play an important role in the quenching history of a galaxy. © V. Kalinova et al. 2022., DC acknowledges support by the Deutsche Forschungsgemeinschaft, DFG project number SFB956A. SFS acknowledges Conacyt projects CB-285080 and FC-2016-01 -1916, and PAPIIT-DGAPA IN100519. RGB and RGD acknowledge financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709), and the projects AYA2016-77846-P and PID2019-109067-GB100. ER acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference number RGPIN-2017-03987. KK acknowledges deep gratitude to MPIfR in Bonn for the Guest Researchership 2013–2021. In this study, we made use of the data of the first legacy survey, the Calar Alto Legacy Integral Field Area (CALIFA) survey, based on observations made at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, operated jointly by the Max Planck-Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC). This research made use of the open-source python packages as Astropy (Astropy Collaboration 2018), SciPy (Virtanen et al. 2020), NumPy (Harris et al. 2020), and Matplotlib (Hunter 2007).

We present spatially resolved stellar kinematics for 797 $z=0.6-1$ galaxies selected from the LEGA-C survey and construct axisymmetric Jeans models to quantify their dynamical mass and degree of rotational support. The survey is $K_s$-band selected, irrespective of color or morphological type, and allows for a first assessment of the stellar dynamical structure of the general $L^*$ galaxy population at large lookback time. Using light profiles from Hubble Space Telescope imaging as a tracer, our approach corrects for observational effects (seeing convolution and slit geometry), and uses well-informed priors on inclination, anisotropy and a non-luminous mass component. Tabulated data include total mass estimates in a series of spherical apertures (1, 5, and 10 kpc; 1$\times$ and 2$\times$\re), as well as rotational velocities, velocity dispersions and anisotropy. We show that almost all star-forming galaxies and $\sim$50\% of quiescent galaxies are rotation-dominated, with deprojected $V/\sigma\sim1-2$. Revealing the complexity in galaxy evolution, we find that the most massive star-forming galaxies are among the most rotation-dominated, and the most massive quiescent galaxies among the least rotation-dominated galaxies. These measurements set a new benchmark for studying galaxy evolution, using stellar dynamical structure for galaxies at large lookback time. Together with the additional information on stellar population properties from the LEGA-C spectra, the dynamical mass and $V/\sigma$ measurements presented here create new avenues for studying galaxy evolution at large lookback time., Comment: Accepted for publication in ApJ. Data table will be published with journal article and is now available upon request

We present an orbit-based method of combining stellar and cold gas kinematics to constrain the dark matter profile of early-type galaxies. We apply this method to early-type galaxy NGC 2974, using Pan-STARRS imaging and SAURON stellar kinematics to model the stellar orbits, and introducing HI kinematics from VLA observation as a tracer of the gravitational potential. The introduction of the cold gas kinematics shows a significant effect on the confidence limits of especially the dark halo properties: we exclude more than $95\%$ of models within the $1-\sigma$ confidence level of Schwarzschild modelling with only stellar kinematics, and reduce the relative uncertainty of the dark matter fraction significantly to $10\%$ within $5 R_\mathrm{e}$. Adopting a generalized-NFW dark matter profile, we measure a shallow cuspy inner slope of $0.6^{+0.2}_{-0.3}$ when including the cold gas kinematics in our model. We cannot constrain the inner slope with the stellar kinematics alone., Comment: 13 pages, 11 figures, accepted for publication in MNRAS

We perform spatially resolved stellar population analysis for a sample of 69 early-type galaxies (ETGs) from the CALIFA integral field spectroscopic survey, including 48 ellipticals and 21 S0's. We generate and quantitatively characterize profiles of light-weighted mean stellar age and metallicity within $\lesssim 2R_e$, as a function of radius and stellar-mass surface density $\mu_*$. We study in detail the dependence of profiles on galaxies' global properties, including velocity dispersion $\sigma_e$, stellar mass, morphology. ETGs are universally characterized by strong, negative metallicity gradients ($\sim -0.3\,\text{dex}$ per $R_e$) within $1\,R_e$, which flatten out moving towards larger radii. A quasi-universal local $\mu_*$-metallicity relation emerges, which displays a residual systematic dependence on $\sigma_e$, whereby higher $\sigma_e$ implies higher metallicity at fixed $\mu_*$. Age profiles are typically U-shaped, with minimum around $0.4\,R_e$, asymptotic increase to maximum ages beyond $\sim 1.5\,R_e$, and an increase towards the centre. The depth of the minimum and the central increase anti-correlate with $\sigma_e$. A possible qualitative interpretation of these observations is a two-phase scenario. In the first phase, dissipative collapse occurs in the inner $1\,R_e$, establishing a negative metallicity gradient. The competition between the outside-in quenching due to feedback-driven winds and some form of inside-out quenching, possibly caused by central AGN feedback or dynamical heating, determines the U-shaped age profiles. In the second phase, the accretion of ex-situ stars from quenched and low-metallicity satellites shapes the flatter stellar population profiles in the outer regions., Comment: Accepted version to appear in MNRAS. Improved discussion with respect to original submission and additional tests included

We study the rest-frame optical mass-size relation of bulges and discs from z similar to 2 to z similar to 0 for a complete sample of massive galaxies in the CANDELS fields using two-component Sersic models. Discs and star-forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best-fitting slope of 0.7 for quiescent galaxies against 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc versus bulge dominated) and the star formation activity (star-forming versus quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our results pose a challenge to current cosmological models that predict distinct structural properties for stellar bulges arising from mergers and disc instabilities.

We construct Schwarzschild orbit-based models of NGC 7457, known as a peculiar low-mass lenticular galaxy. Our best-fitting model successfully retrieves most of the unusual kinematics behaviours of this galaxy, in which, the orbital distribution of stars is dominated by warm and hot orbits. The reconstructed surface brightness of the hot component matches fairly well the photometric bulge and the reconstructed LOSVD map of this component shows clear rotation around the major photometric axis of the galaxy. In the absence of a dominant cold component, the outer part of our model is dominated by warm orbits, representing an exponential thick disc. Our orbital analysis also confirms the existence of a counter-rotating orbital substructure in the very centre, reported in previous observational studies. By comparing our model with a variety of simulation studies, and considering the stellar kinematics and populations properties of this galaxy, we suggest that the thick disc is most likely a dynamically heated structure, formed through the interactions and accretion of satellite(s) with near-polar initial inclination. We also suggest a merger-driven process as the most plausible scenario to explain the observed and dynamically modelled properties of the bulge of NGC 7457. We conclude that both the high level of cylindrical rotation and unusually low velocity dispersion reported for the NGC 7457 have most likely external origins. Therefore, NGC 7457 could be considered as a candidate for merger-driven cylindrical rotation in the absence of a strong bar in disc galaxies.

The flat stellar density cores of massive elliptical galaxies form rapidly due to sinking supermassive black holes (SMBHs) in gas-poor galaxy mergers. After the SMBHs form a bound binary, gravitational slingshot interactions with nearby stars drive the core regions towards a tangentially biased stellar velocity distribution. We use collisionless galaxy merger simulations with accurate collisional orbit integration around the central SMBHs to demonstrate that the removal of stars from the centre by slingshot kicks accounts for the entire change in velocity anisotropy. The rate of strong (unbinding) kicks is constant over several hundred Myr at $\sim 3 \ \mathrm{ M}_\odot\, \rm yr^{-1}$ for our most massive SMBH binary (MBH = 1.7 × 1010 M⊙). Using a frequency-based orbit classification scheme (box, x-tube, z-tube, rosette), we demonstrate that slingshot kicks mostly affect box orbits with small pericentre distances, leading to a velocity anisotropy of β ≲ −0.6 within several hundred Myr as observed in massive ellipticals with large cores. We show how different SMBH masses affect the orbital structure of the merger remnants and present a kinematic tomography connecting orbit families to integral field kinematic features. Our direct orbit classification agrees remarkably well with a modern triaxial Schwarzschild analysis applied to simulated mock kinematic maps.

We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new Multi-Gaussian Expansion (MGE) photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a=1.294$\pm$0.039, b= 0.912$\pm$0.025, and zero-point c= 7.067$\pm$0.078. We leverage the high signal-to-noise of SAMI integral field spectroscopy, to determine how structural and stellar-population observables affect the scatter about the FP. The FP residuals correlate most strongly (8$\sigma$ significance) with luminosity-weighted simple-stellar-population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, S\'ersic index and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio $\Upsilon_\star$) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti-)correlate with the residuals of the relation between surface density and $\Upsilon_\star$. This correlation implies that part of the FP scatter is due to the broad age and $\Upsilon_\star$ distribution at any given surface mass density. Using virial mass and $\Upsilon_\star$ we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75%) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane., Comment: 36 pages, 23 figures

We discuss the evolution of the interstellar medium of quiescent galaxies, currently emerging from recent analyses, with the help of a simple model based on well-established empirical relations such as the stellar mass functions and the main sequence of star formation. This model is meant to describe observed quantities without making specific assumptions on the nature of quenching processes, but relying on their observable consequences. We find that the high gas fractions seen or suggested at high redshift in quiescent galaxies, and their apparent mild evolution at early times, can be mostly attributed to a progenitor effect where recently quenched galaxies with ~10% gas fractions dominate the quiescent galaxy population until z~1. In the same context, the much lower gas and dust fractions measured in local early-type galaxies are interpreted as the product of the steady depletion of their interstellar medium on a ~2 Gyr timescale, coupled with a higher fraction of more gas-exhaustive events., 7 pages, 5 figures, 2 tables (including appendices). Accepted for publication in Astronomy & Astrophysics

The morphology of quiescent galaxies has been found to be correlated with the activity of their central super massive black hole. In this work, we use data from the first data release of the LOFAR Two$-$Metre Sky Survey (LoTSS DR1) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) to select more than 15 000 quiescent galaxies at $z10^{21}\rm\,W\,Hz^{-1}$ at fixed stellar mass, velocity dispersion, or surface mass density does not depend on the galaxy projected axis ratio ($q$). However, the high-power ($L_{\rm 150\,MHz}>10^{23}\rm\,W\,Hz^{-1}$) RLAGNs are more likely to be found in massive, round galaxies, while the low- and intermediate-power ($L_{\rm 150\,MHz}\leq10^{23}\rm\,W\,Hz^{-1}$) RLAGNs have similar distributions of $q$ to non-RLAGN galaxies. We argue that our results support the picture that high-power RLAGNs are more easily triggered in galaxies with a merger-rich history, while low-power RLAGNs can be triggered in galaxies growing mainly via secular processes. Our work also supports the idea that the low-luminosity RLAGN may be sufficient for maintenance-mode feedback in low-mass quiescent galaxies with disc-like morphology, which is based on a simple extrapolation from the observed energy balance between cooling and RLAGN-induced cavities in massive clusters. We find no significant difference between the $q$ distributions of RLAGNs likely to be found in clusters and those likely not found in clusters after controlling the radio luminosity and stellar mass of the two samples, indicating that the environment does not significantly influence the morphology--RLAGN correlation., 14 pages, 8 figures, accepted for publication in A&A

We analyze the three-dimensional shapes and kinematics of the young star cluster population forming in a high-resolution GRIFFIN project simulation of a metal-poor dwarf galaxy starburst. The star clusters, which follow a power-law mass distribution, form from the cold ISM phase with an IMF sampled with individual stars down to 4 solar masses at sub-parsec spatial resolution. Massive stars and their important feedback mechanisms are modelled in detail. The simulated clusters follow a surprisingly tight relation between the specific angular momentum and mass with indications of two sub-populations. Massive clusters ($M_\mathrm{cl}\gtrsim 3\times 10^4 M_{\odot})$ have the highest specific angular momenta at low ellipticities ($\epsilon\sim 0.2$) and show alignment between their shapes and rotation. Lower mass clusters have lower specific angular momenta with larger scatter, show a broader range of elongations, and are typically misaligned indicating that they are not shaped by rotation. The most massive clusters $(M \gtrsim 10^5\,M_{\odot})$ accrete gas and proto-clusters from a $ \lesssim 100\,\rm pc$ scale local galactic environment on a $t \lesssim 10\,\rm Myr$ timescale, inheriting the ambient angular momentum properties. Their two-dimensional kinematic maps show ordered rotation at formation, up to $v \sim 8.5\,\rm km s^{-1}$, consistent with observed young massive clusters and old globular clusters, which they might evolve into. The massive clusters have angular momentum parameters $\lambda_R\lesssim 0.5$ and show Gauss-Hermite coefficients $h_3$ that are anti-correlated with the velocity, indicating asymmetric line-of-sight velocity distributions as a signature of a dissipative formation process., Comment: 21 pages, accepted for publication in ApJ

Galaxy environment has been shown to play an important role in transforming late-type, star-forming galaxies to quiescent spheroids. This transformation is expected to be more severe for low-mass galaxies (M < $10^{10}$ \(M_\odot\)) in dense galaxy groups and clusters, mostly due to the influence of their past host halos (also known as pre-processing) and their present-day environments. For the first time, in this study, we investigate a sample of nine early-type dwarf galaxies (dEs) that were accreted as a likely bound group onto the Virgo galaxy cluster about 2-3 Gyr ago. Considering this special condition, these nine dEs may provide a test bed for distinguishing between the influence of the Virgo galaxy cluster and the effects of the previous host halo on their current properties. Specifically, we use VLT/MUSE integral-field unit spectra to derive their kinematics and specific angular momentum ($\lambda_{R}$) profiles. We observe a spread in the $\lambda_{R}$ profiles of our sample dEs, finding that the $\lambda_{R}$ profiles of half of them are as high as those of low-mass field galaxies. The remaining dEs exhibit $\lambda_{R}$ profiles as low as those of Virgo dEs that were likely accreted longer ago. Moreover, we detect nebular emission in one dE with a gas velocity offset suggesting ongoing gas stripping in Virgo. We suggest that the low-$\lambda_{R}$ dEs in our sample were processed by their previous host halo, prior to their infall to Virgo, and that the high-$\lambda_{R}$ dEs may be experiencing ram pressure stripping in Virgo., Comment: 23 pages, 16 figures, Accepted for publication in MNRAS

We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $H{\delta}_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $H{\delta}_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $H{\delta}_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ \r{A}/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ \r{A}$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius., Comment: 18 pages, 9 figures

Observations of neutral hydrogen (HI) and molecular gas show that 50% of all nearby early-type galaxies (ETGs) contain some cold gas. Molecular gas is always found in small gas discs in the central region of the galaxy, while neutral hydrogen is often distributed in a low-column density disc or ring typically extending well beyond the stellar body. Dust is frequently found in ETGs as well. The goal of our study is to understand the link between dust and cold gas in nearby ETGs as a function of HI content. We analyse deep optical $g-r$ images obtained with the MegaCam camera at the Canada-France-Hawaii Telescope for a sample of 21 HI-rich and 41 HI-poor ETGs. We find that all HI-rich galaxies contain dust seen as absorption. Moreover, in 57 percent of these HI-rich galaxies, the dust is distributed in a large-scale spiral pattern. Although the dust detection rate is relatively high in the HI-poor galaxies ($\sim$59 percent), most of these systems exhibit simpler dust morphologies without any evidence of spiral structures. We find that the HI-rich galaxies possess more complex dust morphology extending to almost two times larger radii than HI-poor objects. We measured the dust content of the galaxies from the optical colour excess and find that HI-rich galaxies contain six times more dust (in mass) than HI-poor ones. In order to maintain the dust structures in the galaxies, continuous gas accretion is needed, and the substantial HI gas reservoirs in the outer regions of ETGs can satisfy this need for a long time. We find that there is a good correspondence between the observed masses of the gas and dust, and it is also clear that dust is present in regions further than 3~Reff. Our findings indicate an essential relation between the presence of cold gas and dust in ETGs and offer a way to study the interstellar medium in more detail than what is possible with HI observations., 11 pages, 7 figures, 3 tables, and 4 appendices. Appendices contain dust masses and images of studied galaxies (62). The fits files will be uploaded into the CDS database. Figures in the appendices are of low quality. High-quality images are given in a separated folder. Accepted for publication in A&A

In this letter, we investigate the impact of environment on integrated and spatially-resolved stellar kinematics of a sample of massive, quiescent galaxies at intermediate redshift ($0.6, Comment: 10 pages, 5 figures, submitted to ApJ Letters

Low-power radio sources dominate the radio sky. They tend to be small in size and core-dominated, but the origin of their properties and the evolution of their radio plasma are not well constrained. Interestingly, there is mounting evidence that low-power radio sources can significantly impact their surrounding gaseous medium and, therefore, may be more relevant for galaxy evolution than previously thought. In this paper, we present low radio frequency observations obtained with LOFAR at 147 MHz of the radio source hosted by NGC 3998. This is a rare example of a low-power source which is extremely core-dominated, but which has two large-scale lobes of low surface brightness. We combine the new 147 MHz image with available 1400 MHz data to derive the spectral index over the source. Despite the low surface brightness, reminiscent of remnant structures, the lobes show an optically thin synchrotron spectral index (~ 0.6). We interpret this as being due to rapid decollimation of the jets close to the core, to high turbulence of the plasma flow, and entrainment of thermal gas. This could be the result of intermittent activity of the central AGN, or, more likely, temporary disruption of the jet due to the interaction of the jet with the rich circumnuclear ISM. Both would result in sputtering energy injection from the core which would keep the lobes fed, albeit at a low rate. We discuss these results in connection with the properties of low-power radio sources in general. Our findings show that amorphous, low surface brightness lobes should not be interpreted, by default, as remnant structures. Large, deep surveys (in particular the LOFAR 150 MHz LoTSS and the recently started 1400 MHz Apertif survey) will identify a growing number of objects similar to NGC 3998 where these ideas can be further tested., Accepted for publication in Astronomy & Astrophysics

Galaxy evolution is driven by a variety of physical processes that are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z ∼ 0.6–1 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D n 4000 and Hδ equivalent widths (EWs) distributions over the stellar mass versus environment 2D spaces for the massive galaxy population (M > 1010 M ⊙) at z ∼ 0.6–1.0. D n 4000 and Hδ EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D n 4000 and Hδ EW, when controlling for stellar mass, are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D n 4000 and Hδ EWs, which depend strongly on their stellar mass and are completely independent of the environment at 0.6 < z < 1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point toward massive galaxies forming first in overdensities or the most compact dark matter haloes.

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (~10 higher in its tail than in its stellar body., Accepted for publication in MNRAS. 18 pages, 7 figures, plus 2 appendices

We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies., 21 pages, 16 figures, 1 table. Accepted for publication in MNRAS; minor corrections to the text, and updated Figures 8 and 9 to match published version

We study the rest-frame optical mass-size relation of bulges and discs from z similar to 2 to z similar to 0 for a complete sample of massive galaxies in the CANDELS fields using two-component Sersic models. Discs and star-forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best-fitting slope of 0.7 for quiescent galaxies against 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc versus bulge dominated) and the star formation activity (star-forming versus quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our results pose a challenge to current cosmological models that predict distinct structural properties for stellar bulges arising from mergers and disc instabilities.

We construct Schwarzschild orbit-based models of NGC 7457, known as a peculiar low-mass lenticular galaxy. Our best-fitting model successfully retrieves most of the unusual kinematics behaviours of this galaxy, in which, the orbital distribution of stars is dominated by warm and hot orbits. The reconstructed surface brightness of the hot component matches fairly well the photometric bulge and the reconstructed LOSVD map of this component shows clear rotation around the major photometric axis of the galaxy. In the absence of a dominant cold component, the outer part of our model is dominated by warm orbits, representing an exponential thick disc. Our orbital analysis also confirms the existence of a counter-rotating orbital substructure in the very centre, reported in previous observational studies. By comparing our model with a variety of simulation studies, and considering the stellar kinematics and populations properties of this galaxy, we suggest that the thick disc is most likely a dynamically heated structure, formed through the interactions and accretion of satellite(s) with near-polar initial inclination. We also suggest a merger-driven process as the most plausible scenario to explain the observed and dynamically modelled properties of the bulge of NGC 7457. We conclude that both the high level of cylindrical rotation and unusually low velocity dispersion reported for the NGC 7457 have most likely external origins. Therefore, NGC 7457 could be considered as a candidate for merger-driven cylindrical rotation in the absence of a strong bar in disc galaxies.

We study the dynamical properties of massive quiescent galaxies at $1.4 < z < 2.1$ using deep Hubble Space Telescope WFC3/F160W imaging and a combination of literature stellar velocity dispersion measurements and new near-infrared spectra obtained using KMOS on the ESO VLT. We use these data to show that the typical dynamical-to-stellar mass ratio has increased by $\sim$0.2 dex from $z = 2$ to the present day, and investigate this evolution in the context of possible changes in the stellar initial mass function (IMF) and/or fraction of dark matter contained within the galaxy effective radius, $f_\mathrm{DM}$. Comparing our high-redshift sample to their likely descendants at low-redshift, we find that $f_\mathrm{DM}$ has increased by a factor of more than 4 since $z \approx 1.8$, from $f_\mathrm{DM}$ = $6.6\pm1.0$% to $\sim$24%. The observed increase appears robust to changes in the methods used to estimate dynamical masses or match progenitors and descendants. We quantify possible variation of the stellar IMF through the offset parameter $\alpha$, defined as the ratio of dynamical mass in stars to the stellar mass estimated using a Chabrier IMF. We demonstrate that the correlation between stellar velocity dispersion and $\alpha$ reported among quiescent galaxies at low-redshift is already in place at $z = 2$, and argue that subsequent evolution through (mostly minor) merging should act to preserve this relation while contributing significantly to galaxies overall growth in size and stellar mass., Comment: 40 pages, 21 figures; accepted for publication in ApJ

Unresolved gas and dust observations show a surprising diversity in the amount of interstellar matter in early-type galaxies. Using ALMA observations we resolve the ISM in z$\sim$0.05 early-type galaxies. From a large sample of early-type galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) we selected five of the dustiest cases, with dust masses M$_d\sim$several$\times10^7$M$_\odot$, with the aim of mapping their submillimetre continuum and $^{12}$CO(2-1) line emission distributions. These observations reveal molecular gas disks. There is a lack of associated, extended continuum emission in these ALMA observations, most likely because it is resolved out or surface brightness limited, if the dust distribution is as extended as the CO gas. However, two galaxies have central continuum ALMA detections. An additional, slightly offset, continuum source is revealed in one case, which may have contributed to confusion in the Herschel fluxes. Serendipitous continuum detections further away in the ALMA field are found in another case. Large and massive rotating molecular gas disks are mapped in three of our targets, reaching a few$\times10^{9}$M$_\odot$. One of these shows evidence of kinematic deviations from a pure rotating disc. The fields of our two remaining targets contain only smaller, weak CO sources, slightly offset from the optical galaxy centres. These may be companion galaxies seen in ALMA observations, or background objects. These heterogeneous findings in a small sample of dusty early-type galaxies reveal the need for more such high spatial resolution studies, to understand statistically how dust and gas are related in early-type galaxies., Comment: 15 pages, 10 figures, 3 tables, accepted for publication in MNRAS

We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying substructures within a galaxy, including embedded and counter-rotating discs, that are easily missed in unilateral studies of the photometry alone. Conversely, bars are rarely apparent in the stellar kinematics but are readily identified on images. Consequently, we distinguish the spin-ellipticity radial tracks of seven morpho-dynamical types: elliptical, lenticular, early spiral, late spiral, barred spiral, embedded disc, and 2-sigma galaxies. The importance of probing beyond the inner radii of galaxies is highlighted by the characteristics of galactic features in the spin-ellipticity radial tracks present at larger radii. The density of information presented through spin-ellipticity radial tracks emphasises a clear advantage to representing galaxies as a track, rather than a single point, in spin-ellipticity parameter space., Comment: 23 pages, 15 figures

Post-starburst galaxies are typically considered to be a transition population, en route to the red sequence after a recent quenching event. Despite this, recent observations have shown that these objects typically have large reservoirs of cold molecular gas. In this paper we study the star-forming gas properties of a large sample of post-starburst galaxies selected from the cosmological, hydrodynamical EAGLE simulations. These objects resemble observed high-mass post-starburst galaxies both spectroscopically and in terms of their space density, stellar mass distribution and sizes. We find that the vast majority of simulated post-starburst galaxies have significant gas reservoirs, with star-forming gas masses of ~10$^9$ M$_{\odot}$, in good agreement with those seen in observational samples. The simulation reproduces the observed time evolution of the gas fraction of the post-starburst galaxy population, with the average galaxy losing ~90 per cent of its star-forming interstellar medium in only ~600 Myr. A variety of gas consumption/loss processes are responsible for this rapid evolution, including mergers and environmental effects, while active galactic nuclei play only a secondary role. The fast evolution in the gas fraction of post-starburst galaxies is accompanied by a clear decrease in the efficiency of star formation, due to a decrease in the dense gas fraction. We predict that forthcoming ALMA observations of the gas reservoirs of low-redshift post-starburst galaxies will show that the molecular gas is typically compact and has disturbed kinematics, reflecting the disruptive nature of many of the evolutionary pathways that build up the post-starburst galaxy population., 15 pages, 17 figures, MNRAS accepted

In this second paper of T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) series we study [Mg/Fe] abundance distribution trends of early-type galaxies observed with the SAURON integral field unit, spanning a wide range in mass and local environment densities: 20 low-mass early-types (dEs) of Sybilska et al. (2017) and 258 massive early types (ETGs) of the $ATLAS^{3D}$ project, all homogeneously reduced and analyzed. We show that the [Mg/Fe] ratios scale with velocity dispersion ($\sigma$) at fixed [Fe/H] and that they evolve with [Fe/H] along similar paths for all early-types, grouped in bins of increasing local and global {$\sigma$}, as well as the second velocity moment $V_{rms}$, indicating a common inside-out formation pattern. We then place our dEs on the [Mg/Fe] $vs.$ [Fe/H] diagram of Local Group galaxies and show that dEs occupy the same region and show a similar trend line slope in the diagram as the high-metallicity stars of the Milky Way and the Large Magellanic Cloud. This finding extends the similar trend found for dwarf spheroidal $vs.$ dwarf irregular galaxies and supports the notion that dEs have evolved from late-type galaxies that have lost their gas at a point of their evolution, which likely coincided with them entering denser environments., Comment: Accepted for publication in MNRAS

The distribution of galaxies in the stellar specific angular momentum versus stellar mass plane ($j_{\star}$-$M_{\star}$) provides key insights into their formation mechanisms. In this paper, we determine the location in this plane of a sample of ten field/group unbarred lenticular (S0) galaxies from the CALIFA survey. We performed a bulge-disc decomposition both photometrically and kinematically to study the stellar specific angular momentum of the disc components alone and understand the evolutionary links between S0s and other Hubble types. We found that eight of our S0 discs have a distribution in the $j_{\star}$-$M_{\star}$ plane that is fully compatible with that of spiral discs, while only two have values of $j_{\star}$ lower than the spirals. These two outliers show signs of recent merging. Our results suggest that merger and interaction processes are not the dominant mechanisms in S0 formation in low-density environments. Instead, S0s appear to be the result of secular processes and the fading of spiral galaxies after the shutdown of star formation., 35 pages, 22 figures. Accepted for publication in MNRAS

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population., Accepted for publication in MNRAS

Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time., Comment: 16 pages, 4 main figures main, 7 supplementary figures, 2 supplementary tables; accepted for publication in Nature Astronomy

We present a kinematic analysis of the dust-lane elliptical NGC 5626 based on MUSE observations. These data allow to robustly classify this galaxy as a fast rotator and to infer a virial mass of $10^{11.7} M_\odot$, making it one of the most massive fast rotators known. In addition, the depth and extent of the MUSE data reveal a strong kinematic twist in the stellar velocity field (by up to $45$ degrees beyond $1.5R_e$). A comparison with the ATLAS$^\mathrm{3D}$ sample underlines the rareness of this system, although we show that such a large-scale kinematic twist could have been missed by the ATLAS$^\mathrm{3D}$ data due to the limited spatial sampling of this survey (typically extending to $0.6R_e$ for massive ETGs). MUSE thus has the potential to unveil more examples of this type of galaxies. We discuss the environment and possible formation history of NGC 5626 and finally argue how a merger between the Milky Way and Andromeda could produce a galaxy of the same class as NGC 5626., 5 pages, 3 figures, letter accepted for publication in MNRAS

A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to z~3, suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the local universe. Here we present deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2-1) emission in 8 massive (Mstar ~ 10^11 Msol) galaxies at z~0.7 selected to lie a factor of 3-10 below the star-forming sequence at this redshift, drawn from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. We significantly detect half the sample, finding molecular gas fractions, Comment: 23 pages, 9 figures, accepted to ApJ. Key results in Figures 7, 8, and 9

Galaxies continuously reprocess their interstellar material. One can therefore expect changing dust grain properties in galaxies which have followed different evolutionary pathways. Determining the intrinsic dust grain mix of a galaxy helps in reconstructing its evolutionary history. Early-type galaxies occasionally display regular dust lanes in their central regions. Due to the relatively simple geometry and composition of their stellar bodies, these galaxies are ideal to disentangle dust mix variations from geometric effects. We therefore model the various components of such a galaxy (FCC 167). We reconstruct its recent history, and investigate the possible fate of the dust lane. MUSE and ALMA observations reveal a nested ISM structure. An ionised-gas disk pervades the central regions of FCC 167, including those occupied by the main dust lane. Inward of the dust lane, we also find a disk/ring of cold molecular gas where stars are forming and HII regions contribute to the ionised-gas emission. Further in, the gas ionisation points towards an active galactic nucleus and the fuelling of a central supermassive black hole from its surrounding ionised and molecular reservoir. Observational constraints and radiative transfer models suggest the dust and gas are distributed in a ring-like geometry and the dust mix lacks small grains. The derived dust destruction timescales from sputtering in hot gas are short and we conclude that the dust must be strongly self-shielding and clumpy, or will quickly be eroded and disappear. Our findings show how detailed analysis of individual systems can complement statistical studies of dust-lane ETGs., Comment: 14 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged], Comment: Accepted for publication in A&A, 52 pages, 30 figures

We report a tight linear relation between the HI circular velocity measured at 6 $R_{\rm e}$ and the stellar velocity dispersion measured within 1 $R_{\rm e}$ for a sample of 16 early-type galaxies with stellar mass between $10^{10}$ and $10^{11}$ $\mathrm{M}_\odot$. The key difference from previous studies is that we only use spatially resolved $v_\mathrm{circ}$(HI) measurements obtained at large radius for a sizeable sample of objects. We can therefore link a kinematical tracer of the gravitational potential in the dark-matter dominated outer regions of galaxies with one in the inner regions, where baryons control the distribution of mass. We find that $v_\mathrm{circ}$(HI) = 1.33 $\sigma_\mathrm{e}$ with an observed scatter of just 12 percent. This indicates a strong coupling between luminous and dark matter from the inner- to the outer regions of early-type galaxies, analogous to the situation in spirals and dwarf irregulars. The $v_\mathrm{circ}$(HI)-$\sigma_\mathrm{e}$ relation is shallower than those based on $v_\mathrm{circ}$ measurements obtained from stellar kinematics and modelling at smaller radius, implying that \vcirc\ declines with radius -- as in bulge-dominated spirals. Indeed, the value of $v_\mathrm{circ}$(HI) is typically 25 percent lower than the maximum $v_\mathrm{circ}$ derived at $\sim0.2\ R_\mathrm{e}$ from dynamical models. Under the assumption of power-law total density profiles $\rho \propto r^{-\gamma}$, our data imply an average logarithmic slope $\langle\gamma\rangle=2.18\pm0.03$ across the sample, with a scatter of 0.11 around this value. The average slope and scatter agree with recent results obtained from stellar kinematics alone for a different sample of early-type galaxies., Comment: MNRAS, accepted

We present new evidence for eight early-type galaxies (ETGs) from the CALIFA Survey that show clear rotation around their major photometric axis ("prolate rotation"). These are LSBCF560-04, NGC 0647, NGC 0810, NGC 2484, NGC 4874, NGC 5216, NGC 6173 and NGC 6338. Including NGC 5485, a known case of an ETG with stellar prolate rotation, as well as UGC 10695, a further possible candidate for prolate rotation, we report ten CALIFA galaxies in total that show evidence for such a feature in their stellar kinematics. Prolate rotators correspond to ~9% of the volume-corrected sample of CALIFA ETGs, a fraction much higher than previously reported. We find that prolate rotation is more common among the most massive ETGs. We investigate the implications of these findings by studying N-body merger simulations, and show that a prolate ETG with rotation around its major axis could be the result of a major polar merger, with the amplitude of prolate rotation depending on the initial bulge-to-total stellar mass ratio of its progenitor galaxies. Additionally, we find that prolate ETGs resulting from this formation scenario show a correlation between their stellar line-of-sight velocity and higher order moment h_3, opposite to typical oblate ETGs, as well as a double peak of their stellar velocity dispersion along their minor axis. Finally, we investigate the origin of prolate rotation in polar galaxy merger remnants. Our findings suggest that prolate rotation in massive ETGs might be more common than previously expected, and can help towards a better understanding of their dynamical structure and formation origin., Comment: accepted for publication in A&A

We study the nearby lenticular galaxy NGC 3998. This galaxy is known to host a low-power radio AGN with a kpc-size one-sided jet and a large, nearly polar HI disc. It is therefore a good system to study to understand the relation between the availability of cold-gas and the triggering of AGNs in galaxies. Our new WSRT data reveal two faint, S-shaped radio lobes extending out to $\sim$10 kpc from the galaxy centre. Remarkably, we find that the inner HI disc warps back towards the stellar mid-plane in a way that mirrors the warping of the radio lobes. We suggest that the polar HI disc was accreted through a minor merger, and that the torques causing it to warp in the inner regions are also responsible for feeding the AGN. The "S" shape of the radio lobes would then be due to the radio jets adapting to the changing angular momentum of the accreted gas. The extended radio jets are likely poorly collimated, which would explain their quick fading and, therefore, their rarity in galaxies similar to NGC 3998. The fuelling of the central super-massive black hole is likely occurring via "discrete events", suggested by the observed variability of the radio core and the extremely high core dominance, which we attribute to the formation and ejection of a new jet resulting from a recent fuelling event., 10 pages, 6 figures, accepted for publication in A&A