Search

Galactic bars are important drivers of galactic evolution, and yet how they impact the interstellar medium and correspondingly star formation, remains unclear. We present simulation results for two barred galaxies with different formation mechanisms, bars formed in isolation or via a tidal interaction, to consider the spatially and temporally varying trends of star formation. We focus on the early (< 1Gyr) epoch of bar formation so that the interaction is clearly identifiable. The nearby NGC 4303 (isolated) and NGC 3627 (interaction history) are selected as observational analogues to tailor these simulations. Regardless of formation mechanism, both models show similar internal dynamical features, although the interaction appears to promote bar-arm disconnection in the outer disc velocity structure. Both bars trigger similar boosts in star formation (79%; 66%), while the interaction also triggers an earlier 31% burst. Significant morphological dependence is observed in the relation between surface gas and star formation rate. In both cases, the bar component is notably steepest; the arm is similar to the overall disc average; and the inter-arm clearly the shallowest. A distinguishable feature of the tidal disc is the presence of moderately dense, inefficiently star forming gas mostly confined to tidal debris outside the optical disc. The tidal disc also exhibits a unique trend of radially increasing star formation efficiency and a clear dearth of star formation which persists along the bar between the centre and bar ends. These are potential signatures for identifying a barred system post-interaction., 19 pages, 15 figures, accepted for publication in MNRAS

Clues and traces of the universe's origin and its developmental process are deeply buried in galaxy shapes and formations. Automated galaxies classification from their images is complicated due to the faintness of the galaxy images, conflicting bright background stars, and image noise. For this purpose, the current work proposes a novel logically structured modular algorithm that analyses galaxy morphological raw brightness data to automatically detect galaxy visual center, region, and classification. First, a novel selective brightness threshold is employed to eliminate the effect of bright background stars on detecting galaxy visual centers. Second, a galaxy region detection technique is developed. Finally, a novel technique based on galaxy brightness variation patterns is employed for galaxies classification. The current work has been tested with a run on a collection of 1000 galaxies from the EFIGI catalog. Results demonstrated a success rate of 97.2% in galaxies classification with an average processing time of 0.37 s per galaxy. The high success rates and the low processing time proved the efficiency of the proposed work.

As the closest $L^{\ast}$ galaxy to our own Milky Way, the Andromeda galaxy (M31) is an ideal laboratory for studies of galaxy evolution. The AMIGA project has recently provided observations of the cool ($T\sim10^4$ K) phase of the circumgalactic medium (CGM) of M31, using HST/COS absorption spectra along $\sim40$ background QSO sightlines, located up to and beyond the galaxy virial radius. Based on these data, and by the means of semi-analytic models and Bayesian inference, we provide here a physical description of the origin and dynamics of the cool CGM of M31. We investigate two competing scenarios, in which (i) the cool gas is mostly produced by supernova(SN)-driven galactic outflows or (ii) it mostly originates from infall of gas from the intergalactic medium. In both cases, we take into account the effect of gravity and hydrodynamical interactions with a hot corona, which has a cosmologically motivated angular momentum. We compare the outputs of our models to the observed covering factor, silicon column density and velocity distribution of the AMIGA absorbers. We find that, to explain the observations, the outflow scenario requires an unphysically large (> 100\%) efficiency for SN feedback. Our infall models, on the other hand, can consistently account for the AMIGA observations and the predicted accretion rate, angular momentum and metallicity are consistent with a cosmological infall from the intergalactic medium., Accepted for publication in MNRAS

Chromospherically active binaries’ (CABs) are late-type evolved stars with strong chromosphere/coronae and high flare activity levels. Starspot evolution in these CABs and binaries reveals strong stellar and chromospheric activities. In binaries, the tidal interaction via the internal dynamo motion, during rapid rotation of the component stars generates magnetic fields. The magnetic activity heats the tenuous coronal plasma, and the reconnection of energetic particles leads to stellar flare formation. Cyclic-variation signatures in brightness and emission lines are indicators of coronal/chromospheric activities. The spatial population densities of CABs are large, and their hard energy spectral characteristics are prominent in our galaxy. Photo-ionization/collisional excitation processes in the coronae of CABs generate a strong thermal emission line at 6.7 keV. In the present work, we re-analyzed stellar flare data of two CABs (Algol & GT Muscae) observed with Suzaku. We resolved the 6.7 keV emission spectrum and equivalent width (EW) in each source’s stellar flare data. The resolved energy spectra are remarkably similar to that of the hard energy spectra at 6.7 keV of X-ray point sources observed in different Galactic regions. The EW of the 6.7 keV emission line obtained in this work compares favorably with the EW range (260 eV – 980 eV) of the 6.7 keV energy spectra of Galactic point sources. Hence, cumulative X-ray emission lines of coronal plasma (and their EWs) emitted by numerous point sources in different Galactic plane positions give clues on the Galactic Ridge X-ray Emission (GRXE) mechanism. We are of the view that CABs could account for a large number of fast-transient X-ray sources that emit thermal emission line at 6.7 keV during their quiescent and flared phases, and these underlying population sources could contribute significantly to the total luminosity of 6.7 keV GRXE.

The Lambda-Cold Dark Matter (LCDM) model agrees with most of the cosmological observations, but has some hindrances from observed data at smaller scales such as galaxies. Recently, Berezhiani and Khoury (2015) proposed a new theory involving interacting superfluid dark matter with three model parameters, which explains galactic dynamics with great accuracy. In the present work, we study the cosmological behaviour of this model in the linear regime of cosmological perturbations. In particular, we compute both analytically and numerically the matter linear growth factor and obtain new bounds for the model parameters which are significantly stronger than previously found. These new constraints come from the fact that structures within the superfluid dark matter framework grow quicker than in LCDM, and quite rapidly when the DM-baryon interactions are strong., Comment: 16 pages, 6 figures

Secure confirmation that a gravitational wave (GW) has been gravitationally lensed would bring together these two pillars of General Relativity for the first time. This breakthrough is challenging for many reasons, including: GW sky localization uncertainties dwarf the angular scale of gravitational lensing, the mass and structure of gravitational lenses is diverse, the mass function of stellar remnant compact objects is not yet well constrained, and GW detectors do not operate continuously. We introduce a new approach that is agnostic to the mass and structure of the lenses, compare the efficiency of different methods for lensed GW discovery, and explore detection of lensed kilonova counterparts as a direct method for localising candidates. Our main conclusions are: (1) lensed neutron star mergers (NS-NS) are magnified into the "mass gap" between NS and black holes, therefore selecting candidates from public GW alerts with high mass gap probability is efficient, (2) the rate of detectable lensed NS-NS will approach one per year in the mid-2020s, (3) the arrival time difference between lensed NS-NS images is $1\,\rm sec\lesssim\Delta t\lesssim1\,year$, and thus well-matched to the operations of GW detectors and optical telescopes, (4) lensed kilonova counterparts are faint at peak (e.g.\ $r_{\rm AB}\simeq24-26$ in the mid-2020s), fade quickly ($d, Comment: Accepted by MNRAS. 20 pages, 10 figures

I consider a sample of eight pressure-supported low-surface brightness galaxies in terms of Milgrom's modified Newtonian dynamics (MOND). These objects include seven nearby dwarf spheroidal galaxies -- Sextans, Carina, Leo II, Sculptor, Draco, Leo I, Fornax, and the ultra-diffuse galaxy DF44. The objects are modelled as Milgromian isotropic isothermal spheres characterised by two parameters that are constrained by observations: the constant line-of-sight velocity dispersion and the central surface density. The velocity dispersion determines the total mass, and, with the implied mass-to-light ratio, the central surface brightness. This then specifies the radial run of surface brightness over the entire isothermal sphere. For these objects the predicted radial distribution of surface brightness is shown to be entirely consistent with observations. This constitutes a success for MOND that is independent of the reduced dynamical mass., Comment: 7 pages, 4 figures, minor errors in Table 1 corrected

We describe the possibilities of the “Cluster Cartography” tool which was created for detailed study of the 2D distribution of galaxies in the clusters. The main tasks of the “Cluster Cartography” tool were the detailed study of the morphologyof galaxy clusters using the statistically significant numerical criteria as well as to detect their regular peculiarities. The tool allows to create the 2D map with positions of galaxies in the cluster field and show for each cluster member its shape and orientation as a best-fit ellipse using input catalogue data. The size of symbols for galaxies correspond to input data.It may reflect the galaxy image in arcseconds from catalogue in the map 4000×4000arcsec. Another way connects the size of the symbol with the magnitude of the galaxy. Tool is able to build the map in four modes: the symbols are dots; the symbols are circles with diameters reflected the magnitudes of galaxies; the symbols are ellipses with size reflected the magnitudesand both ellipticities and orientation from the input catalogue; the symbols illustrate the shape of galaxies in projection to the celestial sphere. The “Cluster Cartography” algorithms allow to detect the standard cases in galaxy distribution, suchas the degree of concentration to the cluster center and/or to some line on a statistically significant level using the numerical criteria. Also “Cluster Cartography” allows to detect other features, such as crosses, semi-crosses, complex crosses and short compact chains, as well as to export the list of galaxies forming the peculiarities for the futurestudy. The final version of the “Cluster Cartography” allows to realize the modern scheme for detailed morphological classification of galaxy clusters. The “Cluster Cartography” is powerful and perspective tool for study of features of galaxy clusters.

The existence of massive galaxies lacking a classical bulge has often been proposed as a challenge to $\Lambda$CDM. However, recent simulations propose that a fraction of massive disc galaxies might have had very quiescent merger histories, and also that mergers do not necessarily build classical bulges. We test these ideas with deep MUSE observations of NGC 5746, a massive ($\sim 10^{11}$ M$_\odot$) edge-on disc galaxy with no classical bulge. We analyse its stellar kinematics and stellar populations, and infer that a massive and extended disc formed very early: 80% of the galaxy's stellar mass formed more than 10 Gyr ago. Most of the thick disc and the bar formed during that early phase. The bar drove gas towards the center and triggered the formation of the nuclear disc followed by the growth of a boxy/peanut-shaped bulge. Around $\sim$ 8 Gyr ago, a $\sim$1:10 merger happened, possibly on a low-inclination orbit. The satellite did not cause significant vertical heating, did not contribute to the growth of a classical bulge, and did not destroy the bar and the nuclear disc. It was however an important event for the galaxy: by depositing its stars throughout the whole galaxy it contributed $\sim 30$% of accreted stars to the thick disc. NGC 5746 thus did not completely escape mergers, but the only relatively recent significant merger did not damage the galaxy and did not create a classical bulge. Future observations will reveal if this is representative of the formation histories of massive disc galaxies., Comment: Accepted for publication in MNRAS. 22 pages, 21 figures (including appendix)

The presence of charged dust grains is known to have a profound impact on the physical evolution of the multiphase interstellar medium (ISM). Despite its importance, this process is still poorly explored in numerical simulations due to its complex physics and the tight dependence on the environment. Here we introduce a novel implementation of grain charging in the cosmological radiative transfer code CRASH. We first benchmark the code predictions on a series of idealized dusty HII regions created by a single star, in order to assess the impact of grain properties on the resulting spatial distribution of charges. Second, we perform a realistic radiative transfer simulation of a star forming region extracted from a dusty galaxy evolving in the Epoch of Reionization. We find that $\sim 13$ % of the total dust mass gets negatively charged, mainly silicate and graphite grains of radius micron. A complex spatial distribution of grain charges is also found, primarily depending on the exposure to stellar radiation and strongly varying along different lines of sight, as a result of radiative transfer effects. We finally assess the impact of grain properties (both chemical composition and size) on the resulting charge distribution. The new implementation described here will open up a wide range of possible studies investigating the physical evolution of the dusty ISM, nowadays accessible to observations of high and low redshift galaxies., Comment: MNRAS in press

We present a machine learning model to classify active galactic nuclei (AGNs) and galaxies (AGN-galaxy classifier) and a model to identify type 1 (optically unabsorbed) and type 2 (optically absorbed) AGN (type 1/2 classifier). We test tree-based algorithms, using training samples built from the X-ray Multi-Mirror Mission–Newton (XMM–Newton) catalogue and the Sloan Digital Sky Survey (SDSS), with labels derived from the SDSS survey. The performance was tested making use of simulations and of cross-validation techniques. With a set of features including spectroscopic redshifts and X-ray parameters connected to source properties (e.g. fluxes and extension), as well as features related to X-ray instrumental conditions, the precision and recall for AGN identification are 94 and 93 per cent, while the type 1/2 classifier has a precision of 74 per cent and a recall of 80 per cent for type 2 AGNs. The performance obtained with photometric redshifts is very similar to that achieved with spectroscopic redshifts in both test cases, while there is a decrease in performance when excluding redshifts. Our machine learning model trained on X-ray features can accurately identify AGN in extragalactic surveys. The type 1/2 classifier has a valuable performance for type 2 AGNs, but its ability to generalize without redshifts is hampered by the limited census of absorbed AGN at high redshift., FJC acknowledges financial support from the Spanish Ministry MCIU under project RTI2018-096686-B-C21 (MCIU/AEI/FEDER/UE), cofunded by FEDER funds and from the Agencia Estatal de Investigación, Unidad de Excelencia María de Maeztu, ref. MDM-2017- 0765. Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions.

We present a sample of 16 likely strong gravitational lenses identified in the VST Optical Imaging of the CDFS and ES1 fields (VOICE survey) using Convolutional Neural Networks (CNNs). We train two different CNNs on composite images produced by superimposing simulated gravitational arcs on real Luminous Red Galaxies observed in VOICE. Specifically, the first CNN is trained on single-band images and more easily identifies systems with large Einstein radii, while the second one, trained on composite RGB images, is more accurate in retrieving systems with smaller Einstein radii. We apply both networks to real data from the VOICE survey, taking advantage of the high limiting magnitude (26.1 in the r-band) and low PSF FWHM (0.8" in the r-band) of this deep survey. We analyse $\sim21,200$ images with $mag_r, 15 pages, 12 figures, submitted to MNRAS

We develop a high-performance analytical model of Galactic Chemical Evolution, which accounts for delay time distributions and lock-up of stellar yields in a thermal-phased ISM. The model is capable of searching, for the first time, through the high-dimensional parameter space associated with the r-process enrichment of the Milky Way by its possible sources: Neutron Star Mergers and Collapsar events. Their differing formation mechanisms give these two processes different time dependencies, a property which has frequently been used to argue in favour of collapsars as the dominant r-process source. However, we show that even with large degrees of freedom in the allowed thermal, structural, and chemical properties of the galaxy, large regions of parameter space are in strong tension with the data. In particular, whilst we are able to find models in which neutron star mergers produce the majority of r-process material, the data rule out all models with dominant collapsar yields. With no other identified source, we conclude that Neutron Star Mergers must be the dominant contributors to the modern Milky Way r-process budget., Comment: 20 pages (3 pages appendices), 16 figures

Obtaining accurately calibrated redshift distributions of photometric samples is one of the great challenges in photometric surveys like LSST, Euclid, HSC, KiDS, and DES. We present an inference methodology that combines the redshift information from the galaxy photometry with constraints from two-point functions, utilizing cross-correlations with spatially overlapping spectroscopic samples, and illustrate the approach on CosmoDC2 simulations. Our likelihood framework is designed to integrate directly into a typical large-scale structure and weak lensing analysis based on two-point functions. We discuss efficient and accurate inference techniques that allow us to scale the method to the large samples of galaxies to be expected in LSST. We consider statistical challenges like the parametrization of redshift systematics, discuss and evaluate techniques to regularize the sample redshift distributions, and investigate techniques that can help to detect and calibrate sources of systematic error using posterior predictive checks. We evaluate and forecast photometric redshift performance using data from the CosmoDC2 simulations, within which we mimic a DESI-like spectroscopic calibration sample for cross-correlations. Using a combination of spatial cross-correlations and photometry, we show that we can provide calibration of the mean of the sample redshift distribution to an accuracy of at least 0.002(1+z), consistent with the LSST-Y1 science requirements for weak lensing and large-scale structure probes., Comment: Updated to match the version accepted by the MNRAS, 23 pages, 9 figures, 2 tables

Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3\sigma significance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3\sigma confidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty., Comment: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

Galaxy morphologies, kinematics, and stellar populations are thought to be linked to each other. However, both simulations and observations have pointed out mismatches therein. In this work, we study the nature and origin of the present-day quenched, bulge-dominated, but dynamically cold galaxies within a stellar mass range of $10.3\,\leqslant\,\log\,M_{\ast}/\mathrm{M_{\odot}}\,\leqslant\,11.2$ in the IllustrisTNG-100 Simulation. We compare these galaxies with a population of normal star-forming dynamically cold disc galaxies and a population of normal quenched dynamically hot elliptical galaxies within the same mass range. The populations of the present-day quenched and bulge-dominated galaxies (both being dynamically cold and hot) used to have significantly higher star-formation rates and flatter morphologies at redshift of $z\sim 2$. They have experienced more frequent larger mass-ratio mergers below $z \sim 0.7$ in comparison to their star-forming disc counterparts, which is responsible for the formation of their bulge-dominated morphologies. The dynamically cold populations (both being star-forming and quenched) have experienced more frequent prograde and tangential mergers especially below $z \sim 1$, in contrast to the dynamically hot ellipticals, which have had more retrograde and radial mergers. Such different merging histories can well explain the differences on the cold and hot dynamical status among these galaxies. We point out that the real-world counterparts of these dynamically cold and hot bulge-dominated quenched populations are the fast- and slow-rotating early-type galaxies, respectively, as seen in observations and hence reveal the different evolution paths of these two distinct populations of early-type galaxies., 14 pages, 12 figures; accepted for publication in MNRAS

Different properties of dark matter haloes, including growth rate, concentration, interaction history, and spin, correlate with environment in unique, scale-dependent ways. While these halo properties are not directly observable, galaxies will inherit their host haloes' correlations with environment. In this paper, we show how these characteristic environmental signatures allow using measurements of galaxy environment to constrain which dark matter halo properties are most tightly connected to observable galaxy properties. We show that different halo properties beyond mass imprint distinct scale-dependent signatures in both the galaxy two-point correlation function and the distribution of distances to galaxies' $k$th nearest neighbours, with features strong enough to be accessible even with low-resolution (e.g., grism) spectroscopy at higher redshifts. As an application, we compute observed two-point correlation functions for galaxies binned by half-mass radius at $z=0$ from the Sloan Digital Sky Survey, showing that classic galaxy size models (i.e., galaxy size being proportional to halo spin) as well as other recent proposals show significant tensions with observational data. We show that the agreement with observed clustering can be improved with a simple empirical model in which galaxy size correlates with halo growth., Comment: 24 pages, MNRAS accepted

In the era of precision cosmology and ever-improving cosmological simulations, a better understanding of different galaxy components such as bulges and discs will give us new insight into galactic formation and evolution. Based on the fact that the stellar populations of the constituent components of galaxies differ by their dynamical properties, we develop two simple models for galaxy decomposition using the IllustrisTNG cosmological hydrodynamical simulation. The first model uses a single dynamical parameter and can distinguish 4 components: thin disc, thick disc, counter-rotating disc and bulge. The second model uses one more dynamical parameter, was defined in a probabilistic manner, and distinguishes two components: bulge and disc. The number fraction of disc-dominated galaxies at a given stellar mass obtained by our models agrees well with observations for masses exceeding $ \log_{10}(M_*/M_\odot)=10$. S\'ersic indices and half-mass radii for the bulge components agree well with those for real galaxies. The mode of the distribution of S\'ersic indices for the disc components is at the expected value of $n=1$. However, disc half-mass radii are smaller than those for real galaxies, in accordance with previous findings that the IllustrisTNG simulation produces undersized discs., Comment: 15 pages, 12 figures

A crucial question in galaxy formation is what role new accretion has in star formation. Theoretical models have predicted a wide range of correlation strengths between halo accretion and galaxy star formation. Previously, we presented a technique to observationally constrain this correlation strength for isolated Milky Way-mass galaxies at $z\sim 0.12$, based on the correlation between halo accretion and the density profile of neighbouring galaxies. By applying this technique to both observational data from the Sloan Digital Sky Survey and simulation data from the UniverseMachine, where we can test different correlation strengths, we ruled out positive correlations between dark matter accretion and recent star formation activity. In this work, we expand our analysis by (1) applying our technique separately to red and blue neighbouring galaxies, which trace different infall populations, (2) correlating dark matter accretion rates with $D_{n}4000$ measurements as a longer-term quiescence indicator than instantaneous star-formation rates, and (3) analyzing higher-mass isolated central galaxies with $10^{11.0} < M_*/M_\odot < 10^{11.5}$ out to $z\sim 0.18$. In all cases, our results are consistent with non-positive correlation strengths with $\gtrsim 85$ per cent confidence, suggesting that processes such as gas recycling dominate star formation in massive $z=0$ galaxies., 16 pages; accepted for publication in MNRAS

We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study., 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS

Dual-Active Galactic Nuclei (AGN) are a natural consequence of the hierarchical structure formation scenario, and can provide an important test of various models for black hole growth. However, due to their rarity and difficulty to find at high redshift, very few confirmed dual-AGN are known at the epoch where galaxy formation peaks. Here we report the discovery of a gravitationally lensed dual-AGN system at redshift 2.37 comprising two optical/IR quasars separated by 6.5+/-0.6 kpc, and a third compact (R_eff = 0.45+/-0.02 kpc) red galaxy that is offset from one of the quasars by 1.7+/-0.1 kpc. From Very Large Array imaging at 3 GHz, we detect 600 and 340 pc-scale radio emission that is associated with both quasars. The 1.4 GHz luminosity densities of the radio sources are about 10^24.35 W / Hz, which is consistent with weak jets. However, the low brightness temperature of the emission is also consistent with star-formation at the level of 850 to 1150 M_sun / yr. Although this supports the scenario where the AGN and/or star-formation is being triggered through an ongoing triple-merger, a post-merger scenario where two black holes are recoiling is also possible, given that neither has a detected host galaxy., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We study the relations between the mass of the central black hole (BH) $M_{\rm BH}$, the dark matter halo mass $M_{\rm h}$, and the stellar-to-halo mass fraction $f_\star\propto M_\star/M_{\rm h}$ in a sample of $55$ nearby galaxies with dynamically measured $M_{\rm BH}>10^6\,{\rm M}_\odot$ and $M_{\rm h}>5\times10^{11}\,{\rm M}_\odot$. The main improvement with respect to previous studies is that we consider both early- and late-type systems for which $M_{\rm h}$ is determined either from globular cluster dynamics or from spatially resolved rotation curves. Independently of their structural properties, galaxies in our sample build a well defined sequence in the $M_{\rm BH}$-$M_{\rm h}$-$f_\star$ space. We find that: (i) $M_{\rm h}$ and $M_{\rm BH}$ strongly correlate with each other and anti-correlate with $f_\star$; (ii) there is a break in the slope of the $M_{\rm BH}$-$M_{\rm h}$ relation at $M_{\rm h}$ of $10^{12}\,{\rm M}_\odot$, and in the $f_\star$-$M_{\rm BH}$ relation at $M_{\rm BH}$ of $\sim10^7\!-\!10^8\,{\rm M}_\odot$; (iii) at a fixed $M_{\rm BH}$, galaxies with a larger $f_\star$ tend to occupy lighter halos and to have later morphological types. We show that the observed trends can be reproduced by a simple equilibrium model in the $��$CDM framework where galaxies smoothly accrete dark and baryonic matter at a cosmological rate, having their stellar and black hole build-up regulated both by the cooling of the available gas reservoir and by the negative feedback from star formation and active galactic nuclei (AGN). Feature (ii) arises as the BH population transits from a rapidly accreting phase to a more gentle and self-regulated growth, while scatter in the AGN feedback efficiency can account for feature (iii)., 20 pages, 12 figures, accepted by MNRAS. Key figures are Fig.1 and Fig.6

We present observations and models of the kinematics and distribution of neutral hydrogen (HI) in the superthin galaxy FGC 1440 with an optical axial ratio $a/b = 20.4$. Using the Giant Meterwave Radio telescope (GMRT), we imaged the galaxy with a spectral resolution of 1.7 $\rm kms^{-1}$ and a spatial resolution of $15" \times 13.5"$. We find that FGC 1440 has an asymptotic rotational velocity of 141.8 $\rm kms^{-1}$ . The structure of the HI disc in FGC 1440 is that of a typical thin disc warped along the line of sight, but we can not rule out the presence of a central thick HI disc. We find that the dark matter halo in FGC 1440 could be modeled by a pseudo-isothermal (PIS) profile with $\rm R_{c}/ R_{d}, Comment: Accepted for publication in MNRAS

We apply a spectral stacking technique to Westerbork Synthesis Radio Telescope observations to measure the neutral atomic hydrogen content (HI) of nearby galaxies in and around galaxy groups at $z < 0.11$. Our sample includes 577 optically-selected galaxies (120 isolated galaxies and 457 satellites) covering stellar masses between 10$^{10}$ and 10$^{11.5}$ M$_{\odot}$, cross-matched with Yang's group catalogue, with angular and redshift positions from the Sloan Digital Sky Survey. We find that the satellites in the centres of groups have lower HI masses at fixed stellar mass and morphology (characterised by the inverse concentration index) relative to those at larger radii. These trends persist for satellites in both high-mass ($M_{\rm halo} > 10^{13.5}h^{-1}$M$_{\odot}$) and low-mass ($M_{\rm halo} \leqslant 10^{13.5}h^{-1}$M$_{\odot}$) groups, but disappear if we only consider group members in low local density ($\Sigma, Comment: 12 pages, 13 figures, 4 tables

Relativistic jets can interact with the ambient gas distribution of the host galaxy, before breaking out to larger scales. In the past decade several studies have simulated jet-driven outflows to understand how they affect the nearby environment, and over what spatial and temporal scales such interactions occur. The simulations are able to capture the interaction of the jets with the turbulent clumpy interstellar medium and the resultant energetics of the gas. In this review we summarise the results of such recent studies and discuss their implications on the evolution of the dynamics of the gas distribution and the star formation rate., To appear in Astronomische Nachrichten, proceedings paper from "6th Workshop on Compact Steep Spectrum and GHz-Peaked Spectrum Radio Sources"

We use a sample of $z=0$ galaxies visually classified as slow rotators (SRs) in the EAGLE hydrodynamical simulations to explore the effect of galaxy mergers on their formation, characterise their intrinsic galaxy properties, and study the connection between quenching and kinematic transformation. SRs that have had major or minor mergers (mass ratios $\ge 0.3$ and $0.1-0.3$, respectively) tend to have a higher triaxiality parameter and ex-situ stellar fractions than those that had exclusively very minor mergers or formed in the absence of mergers ("no-merger" SRs). No-merger SRs are more compact, have lower black hole-to-stellar mass ratios and quenched later than other SRs, leaving imprints on their $z=0$ chemical composition. For the vast majority of SRs we find that quenching, driven by active galactic nuclei feedback, precedes kinematic transformation, except for satellite SRs, in which these processes happen in tandem. However, in $\approx 50$\% of these satellites, satellite-satellite mergers are responsible for their SR fate, while environment (i.e. tidal field and interactions with the central) can account for the transformation in the rest. By splitting SRs into kinematic sub-classes, we find that flat SRs prefer major mergers; round SRs prefer minor or very minor mergers; prolate SRs prefer gas-poor mergers. Flat and prolate SRs are more common among satellites hosted by massive halos ($>10^{13.6}\,\rm M_{\odot}$) and centrals of high masses ($M_{\star} > 10^{10.5}\, \rm M_{\odot}$). Although EAGLE galaxies display kinematic properties that broadly agree with observations, there are areas of disagreement, such as inverted stellar age and velocity dispersion profiles. We discuss these and how upcoming simulations can solve them., Accepted for publications in MNRAS. 22 pages and 20 figures. Main difference with previous version is additional panel in Fig. 7 showing Sersic index vs stellar mass for different merger histories

In this work, we use the constrained Gaussian realization technique to study the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations, exploring its relationship with features of the initial density peaks on large scales, ~1 Mpc/h. Our constrained simulations of volume (20 Mpc/h)^3 successfully reconstruct the large-scale structure as well as the black hole growth for the hosts of the rare 10^9 Msun SMBHs found in the BlueTides simulation at z~7. We run a set of simulations with constrained initial conditions by imposing a 5 \sigma_0(R_G) peak on scale of R_G = 1 Mpc/h varying different peak features, such as the shape and compactness as well as the tidal field surrounding the peak. We find that initial density peaks with high compactness and low tidal field induce the most rapid BH growth at early epochs. This is because compact density peaks with a more spherical large scale matter distribution lead to the formation of high density gas clumps in the centers of halos, and thus boost early BH accretion. Moreover, such initially compact density peaks in low tidal field regions also lead to a more compact BH host galaxy morphology. This can explain the tight correlation between BH growth and host galaxy compactness seen in observations., Comment: matches MNRAS accepted version

The physical properties of the accretion flow and of the X-ray emitting plasma, in supermassive black holes accreting at extreme Eddington rates, are still very unclear. Here we present the analysis of simultaneous XMM–Newton and NuSTAR observations of the hyper-Eddington Seyfert 1 galaxy IRAS 04416+1215, carried out in 2020. The main goal of these observations is to investigate the properties of the X-ray corona, as well as the structure of the accretion flow and of the circumnuclear environment, in this regime of extreme accretion. IRAS 04416+1215 has one of the highest Eddington ratio (λEdd ≃ 472) in the local Universe. It shows an interesting spectral shape, very similar to the standard narrow-line Seyfert 1 galaxy’s spectra, with the presence of multiphase absorption structure composed of three phases, whose estimate of the minimum and maximum distances suggests two different interpretations, one consistent with the three X-ray winds being cospatial, and possibly driven by magnetohydrodynamical processes, the other consistent with the multiphase winds being also multiscale. The X-ray spectrum of IRAS 04416+1215 also has a prominent soft excess component and a hard X-ray emission dominated by a reflection component. Moreover, our detailed spectral analysis shows that IRAS 04416+1215 has the lowest coronal temperature measured so far by NuSTAR (kTe = 3–22 keV, depending on the model). This is consistent with a hybrid coronal plasma, in which the primary continuum emission is driven by pair production due to high-energy tail of the energy distribution of non-thermal electrons.

Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules resulting from the union of aromatic rings and related species, which are likely responsible for strong infrared emission features (3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 microns). In this work, using a sample of 50 Seyfert galaxies (DL300-400) as well as neutral species. By subtracting the spectrum of the central source from the total, we compare the PAH emission in the central vs extended region of a small sample of AGN. In contrast to the findings for the central regions of AGN-dominated systems, the PAH ratios measured in the extended regions of both type 1 and type 2 Seyfert galaxies can be explained assuming similar PAH molecular size distribution and ionized fractions of molecules to those seen in central regions of star-forming galaxies (100, Comment: 29 pages, 21 figures, 6 tables, MNRAS, accepted 2021 October 22

Using hydrodynamical $N$-body simulations, we show that the observed structure and kinematics of the extremely metal-poor, dwarf irregular galaxy DDO 68 is compatible with a merger event with at least two smaller satellite galaxies. We were able to obtain a self-consistent model that simultaneously reproduces several of its observed features, including: the very asymmetric and disturbed shape of the stellar component, the overall HI distribution and its velocity field, the arc-like stellar structure to the west, the low-surface brightness stellar stream to the north. The model implies the interaction of the main progenitor of DDO 68 with two systems with dynamical masses $7\times10^8\,M_{\odot}$ and almost $10^8\,M_{\odot}$ -- 1/20 and 1/150 times the dynamical mass of DDO 68, respectively. We show that the merger between DDO 68 and the most massive of its satellites offers a route to explain the large offset of DDO 68 from the mass-metallicity relation. Assuming that the interacting galaxies have metallicities prior to the merger compatible with those of galaxies with similar stellar masses, we provide quantitative evidence that gas mixing alone does not suffice at diluting the gas of the two components; according to our simulations, the HII regions observed along the Cometary Tail trace the low metallicity of the accreted satellite rather than that of DDO 68's main body. In this case, the mass corresponding to the low metallicity is that of the secondary body and DDO 68 becomes consistent with the mass-metallicity relation., Accepted for publication in MNRAS. 17 pages, 10 figures

Many lines of reasoning suggest that external galaxies should host planetary systems, but detecting them by methods typically used in our own Galaxy is not possible. An alternative approach is to study the temporal behaviour of X-rays emitted by bright extragalactic X-ray sources, where an orbiting planet would temporarily block the X-rays and cause a brief eclipse. We report on such a potential event in the X-ray binary M51-ULS-1 in the galaxy M51. We examined a range of explanations for the observed X-ray dip, including a variety of transiting objects and enhancements in the density of gas and dust. The latter are ruled out by the absence of changes in X-ray colours, save any with sharp density gradients that cannot be probed with our data. Instead, the data are well fit by a planet transit model in which the eclipser is most likely to be the size of Saturn. We also find that the locations of possible orbits are consistent with the survival of a planet bound to a mass-transfer binary. A brightness dip in the extragalactic X-ray binary M51-ULS-1 can be well fit by a planet transit model in which the eclipser is most likely Saturn-sized. The locations of possible orbits are consistent with the survival of a planet bound to a mass-transfer binary.

We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

The ratio of the mass of interstellar dust to the total mass of metals (the dust-to-metals/DTM ratio) tends to increase with metallicity. This can be explained by the increasing efficiency of grain growth in the interstellar medium (ISM) at higher metallicities, with a corollary being that the low DTM ratios seen at low metallicities are due to inefficient stellar dust production. This interpretation assumes that the efficiency of dust destruction in the ISM is constant, whereas it might be expected to increase at low metallicity; the decreased cooling efficiency of low-metallicity gas should result in more post-shock dust destruction via thermal sputtering. We show that incorporating a sufficiently strong metallicity dependence into models of galaxy evolution removes the need for low stellar dust yields. The contribution of stellar sources to the overall dust budget may be significantly underestimated, and that of grain growth overestimated, by models assuming a constant destruction efficiency., Comment: 5 pages, 2 figures. MNRAS accepted

According to the standard cosmological scenario, the large galaxies that we observe today have reached their current mass via mergers with smaller galaxy satellites (Moore et al.1999). This hierarchical process is expected to take place on smaller scales for the satellites themselves, that should build-up from the accretion of smaller building blocks (D'Onghia & Lake 2008). The best chance we have to test this prediction is by looking at the most massive satellite of the Milky Way (MW): the Large Magellanic Cloud (LMC). Smaller galaxies have been revealed to orbit around the LMC (Erkal & Belokurov 2020, Patel et al. 2020), but so far the only evidence for mutual interactions is related to the orbital interplay with the nearby Small Magellanic Cloud (SMC), which is the most massive LMC satellite. In this work, we report the likely discovery of a past merger event that the LMC experienced with a galaxy with a low star formation efficiency and likely having a stellar mass similar to those of dwarf spheroidal galaxies. This former LMC satellite has now completely dissolved, depositing the old globular cluster (GC) NGC 2005 as part of its debris. This GC is the only surviving witness of this ancient merger event, recognizable through its peculiar chemical composition. This discovery is the observational evidence that the process of hierarchical assembly has worked also in shaping our closest satellites., 32 pages, 7 figures, 3 table. Published in Nature Astronomy

We analyse the surface brightness profiles of disc-type galaxies in the EAGLE simulations in order to investigate the effects of galaxy mass and environment on galaxy profile types. Following observational works, we classify the simulated galaxies by their disc surface brightness profiles into single exponential (Type I), truncated (Type II) and anti-truncated (Type III) profiles. In agreement with previous observation and theoretical work, we find that Type II discs result from truncated star-forming discs that drive radial gradients in the stellar populations. In contrast, Type III profiles result from galaxy mergers, extended star-forming discs or the late formation of a steeper, inner disc. We find that the EAGLE simulations qualitatively reproduce the observed trends found between profile type frequency and galaxy mass, morphology and environment, such as the fraction of Type III galaxies increasing with galaxy mass, and the the fraction of Type II galaxies increasing with Hubble type. We investigate the lower incidence of Type II galaxies in galaxy clusters, finding, in a striking similarity to observed galaxies, that almost no S0-like galaxies in clusters have Type II profiles. Similarly, the fraction of Type II profiles for disc-dominated galaxies in clusters is significantly decreased relative to field galaxies. This difference between field and cluster galaxies is driven by star formation quenching. Following the cessation of star formation upon entering a galaxy cluster, the young stellar populations of Type II galaxies simply fade, leaving behind Type I galaxies., 11 pages, 10 figures. Accepted for publication in MNRAS

We present a near-infrared spectroscopic survey of narrow-line Seyfert 1 galaxies in the southern hemisphere (using the SOFI instrument on the ESO-NTT telescope), sampled from optical surveys. We examine the kinematics of the broad-line region, probed by the emission line width of hydrogen (Pa$\alpha$ and H$\beta$). We observed 57 objects, of which we could firmly measure Pa$\alpha$ in 49 cases. We find that a single Lorentzian fit (preferred on theoretical grounds) is preferred over multi-component Gaussian fits to the line profiles; a lack of narrow-line region emission, overwhelmed by the pole-on view of the broad line region (BLR) light, supports this. We recompute the catalog black hole (BH) mass estimates, using the values of FWHM and luminosity of H$\beta$, both from catalog values and re-fitted Lorentzian values. We find a relationship slope greater than unity compared to the catalog values. We ascribe this to contamination by galactic light or difficulties with line flux measurements. However, the comparison of masses computed by the fitted Lorentzian and Gaussian measurements show a slope close to unity. Comparing the BH masses estimated from both Pa$\alpha$ and H$\beta$, the line widths and fluxes shows deviations from expected; in general, however, the computed BH masses are comparable. We posit a scenario where an intermixture of dusty and dust-free clouds (or alternately a structured atmosphere) differentially absorbs the line radiation of the BLR, due to dust absorption and hydrogen bound-free absorption., Comment: Accepted for publication by MNRAS. 15 pages, 15 figures

We study the population of active galaxies in void environment in the Sloan Digital Sky Survey. We use optical spectroscopic information to analyse characteristics of the emission lines of galaxies, accomplished by WHAN and BPT diagrams. Also, we study Wide-field Infrared Survey Explorer(WISE) mid-IR colours to assess active galactic nucleus (AGN) activity. We investigate these different AGN classification schemes, both optical and mid-IR, and their dependence on the spatial location with respect to the void centres. To this end, we define three regions: void, the spherical region defined by voidcentric distance relative to void radius (distance/rvoid) smaller than 0.8, comprising overdensities lesser than −0.9, an intermediate/transition shell region (namely, void-wall) 0.8 < distance/rvoid < 1.2, and a region sufficiently distant from voids, the field: distance/rvoid > 2. We find statistical evidence for a larger fraction of AGN and star-forming galaxies in the void region, regardless of the classification scheme addressed (either BPT, WHAN, or WISE). Moreover, we obtain a significantly stronger nuclear activity in voids compared to the field. We find an unusually large fraction of the most massive black holes undergoing strong accretion when their host galaxies reside in voids. Our results suggest a strong influence of the void environment on AGN mechanisms associated with galaxy evolution.

We study the interstellar medium (ISM) properties as a function of the molecular gas size for 77 infrared-selected galaxies at $z \sim 1.3$. Molecular gas sizes are measured on ALMA images that combine CO(2-1), CO(5-4) and underlying continuum observations, and include CO(4-3), CO(7-6)+[CI]($^3 P_2-^3P_1$), [CI]($^3 P_1-^3P_0$) observations for a subset of the sample. The $\gtrsim 46 \%$ of our galaxies have a compact molecular gas reservoir, and lie below the optical disks mass-size relation. Compact galaxies on and above the main sequence have higher CO excitation and star formation efficiency than galaxies with extended molecular gas reservoirs, as traced by CO(5-4)/CO(2-1) and CO(2-1)/$L_{\rm IR, SF}$ ratios. Average CO+[CI] spectral line energy distributions indicate higher excitation in compacts relative to extended sources. Using CO(2-1) and dust masses as molecular gas mass tracers, and conversion factors tailored to their ISM conditions, we measure lower gas fractions in compact main-sequence galaxies compared to extended sources. We suggest that the sub-millimetre compactness, defined as the ratio between the molecular gas and the stellar size, is an unavoidable information to be used with the main sequence offset to describe the ISM properties of galaxies, at least above $M_{\star} \geqslant 10^{10.6}$ M$_{\odot}$, where our observations fully probe the main sequence scatter. Our results are consistent with mergers driving the gas in the nuclear regions, enhancing the CO excitation and star formation efficiency. Compact main-sequence galaxies are consistent with being an early post-starburst population following a merger-driven starburst episode, stressing the important role of mergers in the evolution of massive galaxies., Accepted on MNRAS on October 11th 2021. Added Figure 10, 11, 12 following referee's comments. Abstract slightly modified to fit to arXiv's requirements

In this paper, we investigate the properties of binary neutron stars (BNSs) and their mergers by combining population synthesis models for binary stellar evolution (BSE) with cosmological galaxy formation and evolution models. We obtain constraints on BSE model parameters by using the observed Galactic BNSs and local BNS merger rate density ($R_0$) inferred from Gravitational Wave (GW) observations, and consequently estimate the host galaxy distributions of BNS mergers. We find that the Galactic BNS observations imply efficient energy depletion in the common envelope (CE) phase, a bimodal kick velocity distribution, and low mass ejection during the secondary supernova explosion. However, the inferred $R_0$ does not necessarily require an extremely high CE ejection efficiency and low kick velocities, different from the previous claims, mainly because the latest inferred $R_0$ is narrowed to a lower value ($320_{-240}^{+490}\,{\rm Gpc^{-3}\,yr^{-1}}$). The BNS merger rate density resulting from the preferred model can be described by $R(z)\sim R_0(1+z)^{\zeta}$ at low redshift ($z\lesssim0.5$), with $R_0\sim316$-$784\,{\rm Gpc^{-3}\,yr^{-1}}$ and $\zeta\sim1.34$-$2.03$, respectively. Our results also show that $R_{0}$ and $\zeta$ depend on settings of BSE model parameters, and thus accurate estimates of these parameters by future GW detections will put strong constraints on BSE models. We further estimate that the fractions of BNS mergers hosted in spiral and elliptical galaxies at $z\sim0$ are $\sim81$%-$84$% and $\sim16$%-$19$%, respectively. The BNS merger rate per galaxy can be well determined by the host galaxy stellar mass, star formation rate, and metallicity, which provides a guidance in search for most probable candidates of BNS host galaxies., Comment: 32 pages, 17 figures, accepted for publication in MNRAS

The characteristic mass that sets the peak of the stellar initial mass function (IMF) is closely linked to the thermodynamic behaviour of interstellar gas, which controls how gas fragments as it collapses under gravity. As the Universe has grown in metal abundance over cosmic time, this thermodynamic behaviour has evolved from a primordial regime dominated by the competition between compressional heating and molecular hydrogen cooling to a modern regime where the dominant process in dense gas is protostellar radiation feedback, transmitted to the gas by dust-gas collisions. In this paper we map out the primordial-to-modern transition by constructing a model for the thermodynamics of collapsing, dusty gas clouds at a wide range of metallicities. We show the transition from the primordial regime to the modern regime begins at metallicity $Z\sim 10^{-4} \rm{Z_\odot}$, passes through an intermediate stage where metal line cooling is dominant at $Z \sim 10^{-3}\,\rm{Z_{\odot}}$, and then transitions to the modern dust- and feedback-dominated regime at $Z\sim 10^{-2} \rm{Z_\odot}$. In low pressure environments like the Milky Way, this transition is accompanied by a dramatic change in the characteristic stellar mass, from $\sim 50\,\rm{M_\odot}$ at $Z \sim 10^{-6}\,\rm{Z_{\odot}}$ to $\sim 0.3\,\rm{M_\odot}$ once radiation feedback begins to dominate, which marks the appearance of the modern bottom-heavy Milky Way IMF. In the high pressure environments typical of massive elliptical galaxies, the characteristic mass for the modern, dust-dominated regime falls to $\sim 0.1\,\rm{M_{\odot}}$, thus providing an explanation for the more bottom-heavy IMF observed in these galaxies. We conclude that metallicity is a key driver of variations in the characteristic stellar mass, and by extension, the IMF., 26 pages (main result in figures 7 and 8). MNRAS in press

We place novel constraints on the cosmokinetic parameters by using a joint analysis of the newest Very Large Telescope (VLT)-K-band Multi Object Spectrograph (KMOS) HII galaxies (HIIG) with the Type Ia Supernovae (SNIa) Pantheon sample. We combine the latter data sets in order to reconstruct, in a model-independent way, the Hubble diagram to as high redshifts as possible. Using a Gaussian process, we derive the basic cosmokinetic parameters and compare them with those of Lambda cold dark matter (ΛCDM). In the case of SNIa, we find that the extracted values of the cosmokinetic parameters are in agreement with the predictions of ΛCDM model. Combining SNIa with high-redshift tracers of the Hubble relation, namely HIIG data, we obtain consistent results with those based on ΛCDM as far as the present values of the cosmokinetic parameters are concerned, but find significant deviations in the evolution of the cosmokinetic parameters with respect to the expectations of the concordance ΛCDM model.

We estimate the 3D density profile of the Galactic dark matter (DM) halo within $r \lesssim 30$ kpc from the Galactic centre by using the astrometric data for halo RR Lyrae stars from Gaia DR2. We model both the stellar halo distribution function and the Galactic potential, fully taking into account the survey selection function, the observational errors, and the missing line-of-sight velocity data for RR Lyrae stars. With a Bayesian MCMC method, we infer the model parameters, including the density flattening of the DM halo $q$, which is assumed to be constant as a function of radius. We find that 99\% of the posterior distribution of $q$ is located at $q>0.963$, which strongly disfavours a flattened DM halo. We cannot draw any conclusions as to whether the Galactic DM halo at $r \lesssim 30$ kpc is prolate, because we restrict ourselves to axisymmetric oblate halo models with $q\leq1$. Our result is inconsistent with predictions from cosmological hydrodynamical simulations that advocate more oblate ($\langle{q}\rangle \sim0.8 \pm 0.15$) DM halos within $\sim 15\%$ of the virial radius for Milky-Way-sized galaxies. An alternative possibility, based on our validation tests with a cosmological simulation, is that the true value $q$ of the Galactic halo could be consistent with cosmological simulations but that disequilibrium in the Milky Way potential is inflating our measurement of $q$ by 0.1-0.2. As a by-product of our analysis, our model constrains the DM density in the Solar neighbourhood to be $��_{\mathrm{DM},\odot} = (9.01^{+0.18}_{-0.20})\times10^{-3}M_\odot \mathrm{pc}^{-3} = 0.342^{+0.007}_{-0.007}$ $\;\mathrm{GeV} \mathrm{cm}^{-3}$., MNRAS submitted. 15 pages (plus Appendix 7 pages). 6 figures (plus 10 figures in Appendix). Comments welcome

We measure the projected half-light radii of young star clusters in 31 galaxies from the Legacy Extragalactic UV Survey (LEGUS). We implement a custom pipeline specifically designed to be robust against contamination, which allows us to measure radii for 6097 clusters. This is the largest sample of young star cluster radii currently available. We find that most (but not all) galaxies share a common cluster radius distribution, with the peak at around 3 pc. We find a clear mass-radius relation of the form $R_{\mathrm{eff}} \propto M^{0.24}$. This relation is present at all cluster ages younger than 1 Gyr, but with a shallower slope for clusters younger than 10 Myr. We present simple toy models to interpret these age trends, finding that high-mass clusters are more likely to be not tidally limited and expand. We also find that most clusters in LEGUS are gravitationally bound, especially at older ages or higher masses. Lastly, we present the cluster density and surface density distributions, finding a large scatter that appears to decrease with cluster age. The youngest clusters have a typical surface density of 100 $M_\odot$ pc$^{-2}$., Accepted by MNRAS. Cluster catalogs available at https://gillenbrown.github.io/LEGUS-sizes/

An overwhelming majority of galactic spiral arms trail with respect to the rotation of the galaxy, though a small sample of leading spiral arms has been observed. The formation of these leading spirals is not well understood. Here we show, using collisionless N-body simulations, that a barred disc galaxy in a retrograde dark matter halo can produce long-lived (∼3 Gyr) leading spiral arms. Due to the strong resonant coupling of the disc to the halo, the bar slows rapidly and spiral perturbations are forced ahead of the bar. We predict that such a system, if observed, will also host a dark matter wake oriented perpendicular to the stellar bar. More generally, we propose that any mechanism that rapidly decelerates the stellar bar will allow leading spiral arms to flourish.

With high-quality data from programs like the Hubble Frontier Fields, cluster lensing has reached the point that models are dominated by systematic rather than statistical uncertainties. We introduce a Bayesian framework to quantify systematic effects by determining how different lens modelling choices affect the results. Our framework includes a new two-sample test for quantifying the difference between posterior probability distributions that are sampled by methods like Monte Carlo Markov chains. We use the framework to examine choices related to the selection and treatment of cluster member galaxies in two of the Frontier Field clusters: Abell 2744 and MACS J0416.1–2403. When selecting member galaxies, choices about depth and area affect the models; we find that model results are robust for an I-band magnitude limit of mlim ≥ 22.5 mag and a radial cut of rlim ≥ 90 arcsec (from the centre of the field), although the radial limit likely depends on the spatial extent of lensed images. Mass is typically assigned to galaxies using luminosity/mass scaling relations. We find that the slopes of the scaling relations can have significant effects on lens model parameters but only modest effects on lensing magnifications. Interestingly, scatter in the scaling relations affects the two fields differently. This analysis illustrates how our framework can be used to analyse lens modelling choices and guide future cluster lensing programs.

We investigate a sample of seven edge-on spiral galaxies using Chandra observations. Edge-on spiral galaxies allow us to clearly separate source associated with their star-forming regions versus the outer edges of the system; offering a clear advantage over other systems. We uncover a number of X-ray point sources across these galaxies, and after eliminating contaminating foreground and background sources, we identify 12 candidate ultraluminous X-ray sources. All of these sources are projected onto the central regions, implying that the majority of ULXs in this sample of spiral galaxies are disk/bulge, and thus not halo sources. This also includes two transient ULXs, which may be long-duration transients and low mass X-ray binaries. This finding illustrates the need for further studies of transient ULXs., accepted to MNRAS

Inferring the point spread function (PSF) at galaxy positions is one of the crucial steps of the shear measurement. We introduce a novel method to estimate the PSFs at the galaxy positions by using the galaxy images, which could provide additional constrains for the PSF field variations. We construct the PSF for each star image by using Principal-Components-Analysis (PCA) method, which can capture the most significant characteristics of the data. Our method utilizes the image difference of the same object between multi-exposures to probe the coefficients of the principal components, in which the differences are mainly caused by PSFs. We apply our method to the observed data. The results show that the corresponding PSFs can be properly estimated from multiple images of different exposures. We then use the obtained principal components from the observations to mock multi-exposure images, where the PSFs field of each exposure is constructed by bivariate polynomial on coefficients. We find that our method can reproduce the PSFs consistently with mocked data. Our results show that the multi-exposed galaxy images could provide us additional constraints for the PSF fields in PCA scenario. It offers a promising prospect for combing the information of stars and galaxies together to construct the PSF field when the point sources are sparsely sampled.

H$\alpha$ blobs are off-galaxy emission-line regions with weak or no optical counterparts. They are mostly visible in H$\alpha$ line, appearing as concentrated blobs. Such unusual objects have been rarely observed and studied, and their physical origin is still unclear. We have identified 13 H$\alpha$ blobs in the public data of MaNGA survey, by visually inspecting both the optical images and the spatially resolved maps of H$\alpha$ line for $\sim 4600$ galaxy systems. Among the 13 H$\alpha$ blobs, 2 were reported in previously MaNGA-based studies and 11 are newly discovered. This sample, though still small in size, is by far the largest sample with both deep imaging and integral field spectroscopy. Therefore, for the first time we are able to perform statistical studies to investigate the physical origin of H$\alpha$ blobs. We examine the physical properties of these H$\alpha$ blobs and their associated galaxies, including their morphology, environments, gas-phase metallicity, kinematics of ionized gas, and ionizing sources. We find that the H$\alpha$ blobs in our sample can be broadly divided into two groups. One is associated with interacting/merging galaxy systems, of which the ionization is dominated by shocks or diffuse ionized gas. It is likely that these H$\alpha$ blobs used to be part of their nearby galaxies, but were stripped away at some point due to tidal interactions. The other group is found in gas-rich systems, appearing as low-metallicity star-forming regions that are visually detached from the main galaxy. These H$\alpha$ blobs could be associated with faint disks, spiral arms, or dwarf galaxies., Comment: 26 pages, 18 figures, accepted for publication in MNRAS

We characterize the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one and we also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (rmag ≤ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength r0 and slope γ by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterized and validated with the spectroscopic data of Survey Validation. © 2021 The Author(s)., PZ, OR-M, SC, PN, and CB acknowledge support from the Science Technology Facilities Council through ST/P000541/1 andST/T000244/1. OR-M is supported by the Mexican National Council of Science and Technology (CONACyT) through grant No. 297228/440775 and funding from the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Curie grant agreement No 734374. JM gratefully acknowledges support from the U.S. Department of Energy, Office of Science, Office of High Energy Physics under Award Number DESC0020086 and from the National Science Foundation under grant AST-1616414. Authors want to thank the GAMA collaboration for early access to GAMA DR4 data for this work. Some of the results in this paper have been derived using the healpy and HEALPIX package. We acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr).This work also made extensive use of the NASA Astrophysics Data System and of the astro-ph preprint archive at arXiv.org. This work used the DiRAC@Durham facility managed by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). The equipment was funded by BEIS capital funding via STFC capital grants ST/K00042X/1, ST/P002293/1, and ST/R002371/1, Durham University and STFC operations grant ST/R000832/1. DiRAC is part of the National e-Infrastructure. This research used resources of the National Energy Research Scientific Computing Center (NERSC). NERSC is a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. This research is supported by the Director, Office of Science, Office of High Energy Physics of the U.S. Department of Energy under Contract No. DE-AC02-05CH1123, and by the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility under the same contract; additional support for DESI is provided by the U.S. National Science Foundation, Division of Astronomical Sciences under Contract No. AST-0950945 to the NSF’s National Optical-Infrared Astronomy Research Laboratory; the Science and Technologies Facilities Council of the United Kingdom; the Gordon and Betty Moore Foundation; the Heising-Simons Foundation; the French Alternative Energies and Atomic Energy Commission (CEA); the National Council of Science and Technology of Mexico; the Ministry of Economy of Spain, and by the DESI Member Institutions. The authors are honored to be permitted to conduct astronomical research on Iolkam Du’ag (Kitt Peak), a mountain with particular significance to the Tohono O’odham Nation., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

We report the discovery of [NeV]3426 emission, in addition to HeII4686 emission, in six compact star-forming galaxies. These observations considerably increase the sample of eight such galaxies discovered earlier by our group. For four of the new galaxies, the optical observations are supplemented by near-infrared spectra. All galaxies, but one, have HII regions that are dense, with electron number densities of ~300-700 cm-3. They are all characterised by high Hbeta equivalent widths EW(Hbeta)~190-520A and high O32=[OIII]5007/[OII]3727 ratios of 10-30, indicating young starburst ages and the presence of high ionization radiation. All are low-metallicity objects with 12+logO/H=7.46-7.88. The spectra of all galaxies show a low-intensity broad component of the Halpha line and five out of six objects show Wolf-Rayet features. Comparison with photoionization models shows that pure stellar ionization radiation from massive stars is not hard enough to produce such strong [NeV] and HeII emission in our galaxies. The [NeV]3426/HeII4686 flux ratio of ~1.2 in J1222+3602 is consistent with some contribution of active galactic nucleus ionizing radiation. However, in the remaining five galaxies, this ratio is considerably lower, 20 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society