Search

Star formation is a complex process that typically occurs in dense regions of molecular clouds mainly regulated by magnetic fields, magnetohydrodynamic (MHD) turbulence, and self-gravity. However, it remains a challenging endeavor to trace the magnetic field and determine regions of gravitational collapse where the star is forming. Based on the anisotropic properties of MHD turbulence, a new technique termed Velocity Gradient Technique (VGT) has been proposed to address these challenges. In this study, we apply the VGT to two regions of the giant California Molecular Cloud (CMC), namely, L1478 and L1482, and analyze the difference in their physical properties. We use the $^{12}$CO (J = 2 - 1), $^{13}$CO (J = 2 - 1), and C$^{18}$O (J = 2 - 1) emission lines observed with the Heinrich Hertz Submillimeter Telescope. We compare VGT results calculated in the resolutions of $3.3'$ and $10'$ to Planck polarization at 353 GHz and $10'$ to determine areas of MHD turbulence dominance and self-gravity dominance. We show that the resolution difference can introduce misalignment between the two measurements. We find the VGT-measured magnetic fields globally agree with that from Planck in L1478 suggesting self-gravity's effect is insignificant. The best agreement appears in VGT-$^{12}$CO. As for L1482, the VGT measurements are statistically perpendicular to the Planck polarization indicating the dominance of self-gravity. This perpendicular alignment is more significant in VGT-$^{13}$CO and VGT-C$^{18}$O., Comment: 14 pages, 16 figures, accepted for publication in MNRAS

Classification of galaxy morphology is a hot issue in astronomical research. Although significant progress has been made in the last decade in classifying galaxy morphology using deep learning technology, there are still some deficiencies in spatial feature representation and classification accuracy. In this study, we present a multiscale convolutional capsule network (MSCCN) model for the classification of galaxy morphology. First, this model improves the convolutional layers using a multibranch structure to extract the multiscale hidden features of galaxy images. In order to further explore the hidden information in the features, the multiscale features are encapsulated and fed into the capsule layer. Second, we use a sigmoid function to replace the softmax function in dynamic routing, which can enhance the robustness of MSCCN. Finally, the classification model achieves 97 per cent accuracy, 96 per cent precision, 98 per cent recall, and 97 per cent F1-score under macroscopic averaging. In addition, a more comprehensive model evaluation was accomplished in this study. We visualized the morphological features for the part of sample set that used the t-distributed stochastic neighbour embedding (t-SNE) algorithm. The results show that the model has a better generalization ability and robustness, and it can be effectively used in the galaxy morphological classification.

In this paper we analyse a sample of 46 barred galaxies of MaNGA. Our goal is to investigate the stellar kinematics of these galaxies and obtain their rotation curves. Additionally, we aim to derive the total stellar and dynamical masses, as well as the maximum rotation velocity, in order to examine their distributions and scaling relations. Using the Pipe3D dataproducts publicly available we obtained the rotation curves, which were fitted considering two components of an axisymmetric Miyamoto–Nagai gravitational potential. We found a wide range of the maximum rotation velocities (117–340 ${\rm km\, s^{-1}}$), with a mean value of 200 ${\rm km\, s^{-1}}$. In addition we found that the total stellar and dynamical masses are in the range of log(Mstar/M⊙) = 10.1−11.5, with a mean value of log(Mstar/M⊙) = 10.8, and log(Mdyn/M⊙) = 10.4−12.0, with a mean value of log(Mdyn/M⊙) = 11.1, respectively. We found a strong correlation between dynamical mass and maximum velocity, between maximum velocity and magnitude, and between stellar mass and maximum velocity. According to these results, barred galaxies exhibit similar behaviour to that of normal spiral galaxies with respect to these relations, as well as in terms of the distribution of their dynamical mass and maximum rotation velocity. However, we found that the distribution of stellar masses of barred galaxies is statistically different from other samples including non-barred galaxies. Finally, analysing the galaxies that show nuclear activity, we find no difference with the rotation curves of normal galaxies.

Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require high-resolution multi-wavelength data. The new sensitive and high-resolution MIGHTEE survey of the XMM-LSS field along with data from the Low Frequency Array (LOFAR) and the Giant Metrewave Radio Telescope (GMRT) provide data taken at different frequencies with similar resolution, which enables us to determine the spectral age distribution for radio loud AGN in the survey field. In this study we present a sample of 28 radio galaxies with their best fitting spectral age distribution analyzed using the Jaffe-Perola (JP) model on a pixel-by-pixel basis. Fits are generally good and objects in our sample show maximum ages within the range of 2.8 Myr to 115 Myr with a median of 8.71 Myr. High-resolution maps over a range of frequencies are required to observe detailed age distributions for small sources and high-sensitivity maps will be needed in order to observe fainter extended emission. We do not observe any correlation between the total physical size of the sources and their age and we speculate both dynamical models and the approach to spectral age analysis may need some modification to account for our observations., 20 pages, 9 figures

We present a study of rotation across 30 square degrees of the Orion Star-forming Complex, following a $\sim$200 d photometric monitoring campaign by the Next Generation Transit Survey (NGTS). From 5749 light curves of Orion members, we report periodic signatures for 2268 objects and analyse rotation period distributions as a function of colour for 1789 stars with spectral types F0$-$M5. We select candidate members of Orion using $\textit{Gaia}$ data and assign our targets to kinematic sub-groups. We correct for interstellar extinction on a star-by-star basis and determine stellar and cluster ages using magnetic and non-magnetic stellar evolutionary models. Rotation periods generally lie in the range 1$-$10 d, with only 1.5 per cent of classical T Tauri stars or Class I/II young stellar objects rotating with periods shorter than 1.8 d, compared with 14 per cent of weak-line T Tauri stars or Class III objects. In period$-$colour space, the rotation period distribution moves towards shorter periods among low-mass (>M2) stars of age 3$-$6 Myr, compared with those at 1$-$3 Myr, with no periods longer than 10 d for stars later than M3.5. This could reflect a mass-dependence for the dispersal of circumstellar discs. Finally, we suggest that the turnover (from increasing to decreasing periods) in the period$-$colour distributions may occur at lower mass for the older-aged population: $\sim$K5 spectral type at 1$-$3 Myr shifting to $\sim$M1 at 3$-$6 Myr., This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. Published by Oxford University Press on behalf of the Royal Astronomical Society. 20 pages. 21 figures

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core-collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper-motion data from Gaia EDR3 and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions ($\beta\approx -0.4$$\pm$$0.1$) in the outskirts. We also robustly measure a dark central mass of roughly $800\pm300 \,$M$_{\odot}$, but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), or a dark population of stellar remnants of extent $\approx 0.016\,\rm pc \simeq 3300\,AU$. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended ($\sim 0.034\, \rm{pc} \approx 7000\,\rm AU$). We use Monte Carlo $N$-body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e., IMBH vs. remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core., Comment: 19 page, 15 figures, 3 tables. Accepted for publication in MNRAS

We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the PHATTER catalogue of 1214 clusters with ages measured via colour-magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution ALMA $^{12}$CO(2-1) survey. In M31, we compare the PHAT catalogue of 1249 clusters to 251 GMCs measured from a CARMA $^{12}$CO(1-0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are $\leq$10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than $\approx 18$ Myr are uncorrelated with the molecular ISM. Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend $\approx$4-6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of $\approx 11$-15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5-10 km s$^{-1}$ velocity dispersion relative to the molecular ISM., Comment: 14 pages, 13 figures, 1 tables, accepted for publication at MNRAS

The Granger causality is an econometric test for determining whether one time series is useful for forecasting another one with a certain Granger lag. Here, the light curves in the 0.3-1 keV (reflection dominated, soft) and 1.2-5 keV (continuum dominated, hard) bands of Active Galactic Nuclei (AGNs) are produced, and the Granger lags are estimated and compared to the traditional lag-frequency spectra. We find that the light curves in the hard band Granger-cause (lead) those in the soft band, whereas the obtained Granger lags could be interpreted as the intrinsic reverberation lags. Then, we extract the Granger-lag profiles from 14 XMM-Newton observations of IRAS 13224-3809, and find that the lags are significant in 12 observations. The majority of the obtained Granger (intrinsic) lags are ~ 200-500 s. With the IRAS 13224-3809 mass of 2 $\times 10^6$ $M_{\odot}$, these lags correspond to the true light-travel distance of ~ 20-50 $r_g$. Hence, by assuming a lamp-post geometry and a face-on disc, this places the corona at ~ 10-25 $r_{g}$ above the central black hole. Moreover, multiple Granger lags consisting of the small and large lags of < 500 s and > 1000 s are detected in 4 observations. This means that the corona height can significantly evolve from ~ 10-25 $r_{g}$ to ~ 55 $r_{g}$, or vice versa, towards the end of the observation. Granger lags are a promising way to measure the intrinsic lags, and provide evidence of coronal height variability within each individual observation., 13 pages, 10 figures, 1 Table, accepted for publication in MNRAS

We present an interstellar medium and stellar population analysis of three spectroscopically confirmed $z>7$ galaxies in the ERO JWST NIRCam and JWST NIRSpec data of the SMACS J0723.3-7327 cluster. We use the Bayesian spectral energy distribution (SED) fitting code \texttt{Prospector} with a flexible star-formation history (SFH), a variable dust attenuation law, and a self-consistent model of nebular emission (continuum and emission lines). Importantly, we self-consistently fit both the emission line fluxes from JWST NIRSpec and the broad-band photometry from JWST NIRCam, taking into account slit-loss effects. We find that these three $z=7.6-8.5$ galaxies ($M_{\star}\approx10^{8}~M_{\odot}$) are young with rising SFHs and mass-weighted ages of $3-4$ Myr, though we find indications for underlying older stellar populations. The inferred gas-phase metallicities broadly agree with the direct metallicity estimates from the auroral lines. The galaxy with the lowest gas-phase metallicity ($\mathrm{Z}_{\rm gas}=0.06~\mathrm{Z}_{\odot}$) has a steeply rising SFH, is very compact ($, 14 pages, 9 figures; accepted for publication in MNRAS

We measure the black hole mass and investigate the accretion flow around the local ($z=0.0502$) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H$\beta$ lags and linewidths from which we estimate a black hole mass of $\log \left(M_{\bullet}/\mathrm{M}_{\odot} \right) = 7.0$, uncertain by $\sim0.4$ dex. High cadence photometric monitoring over two years with the Las Cumbres Observatory provides lightcurves in 7 optical bands suitable for intensive continuum reverberation mapping. We identify variability on two timescales. Slower variations on a 100-day timescale exhibit excess flux and increased lag in the $u'$ band and are thus attributable to diffuse bound-free continuum emission from the broad line region. Faster variations that we attribute to accretion disc reprocessing lack a $u'$-band excess and have flux and delay spectra consistent with either $\tau \propto \lambda^{4/3}$, as expected for a temperature structure of $T(R) \propto R^{-3/4}$ for a thin accretion disc, or $\tau \propto \lambda^{2}$ expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_{\nu} \sim \rm{const}$. Modelling the SED predicts an Eddington ratio of $\lambda_{\rm Edd} > 1$, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc which requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear, however a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations., Comment: 23 pages, 16 figures, 10 tables. Accepted for publication in MNRAS

We analyse the quenched fractions, gas content, and star formation histories of ~1200 satellite galaxies with $M_* \geq 5 \times 10^6~{\rm M}_\odot$ around 198 Milky Way- (MW) and Andromeda-like (M31) hosts in TNG50, the highest-resolution simulation of IllustrisTNG. Satellite quenched fractions are larger for smaller masses, for smaller distances to their host galaxy, and in the more massive M31-like compared to MW-like hosts. As satellites cross their host's virial radius, their gas content drops: most satellites within 300 kpc lack detectable gas reservoirs at $z=0$, unless they are massive like the Magellanic Clouds and M32. Nevertheless, their stellar assembly exhibits a large degree of diversity. On average, the cumulative star formation histories are more extended for brighter, more massive satellites with a later infall, and for those in less massive hosts. Based on these relationships, we can even infer infall periods for observedMWand M31 dwarfs: e.g. 0-4 Gyr ago for the Magellanic Clouds and Leo I, 4-8 and 0-2 Gyr ago for M32 and IC 10, respectively. Ram pressure stripping (in combination with tidal stripping) deprives TNG50 satellites of their gas reservoirs and ultimately quenches their star formation, even though only a few per cent of the present-day satellites around the 198 TNG50 MW/M31-like hosts appear as jellyfish. The typical time since quenching for currently quenched TNG50 satellites is $6.9^{+2.5}_{-3.3}~{\rm Gyr}$ ago. The TNG50 results are consistent with the quenched fractions and stellar assembly of observed MW and M31 satellites, however, satellites of the SAGA survey with $M_* \sim 10^{8-9}~{\rm M}_\odot$ exhibit lower quenched fractions than TNG50 and other, observed analogues., Comment: Accepted for publication in MNRAS

We present a novel technique called "photometric IGM tomography" to map the intergalactic medium (IGM) at $z\simeq4.9$ in the COSMOS field. It utilizes deep narrow-band (NB) imaging to photometrically detect faint Ly$\alpha$ forest transmission in background galaxies across the Subaru/Hyper-Suprime Cam (HSC)'s $1.8\rm\,sq.\,deg$ field of view and locate Ly$\alpha$ emitters (LAEs) in the same cosmic volume. Using ultra-deep HSC images and Bayesian spectral energy distribution fitting, we measure the Ly$\alpha$ forest transmission at $z\simeq4.9$ along a large number ($140$) of background galaxies selected from the DEIMOS10k spectroscopic catalogue at $4.98, Comment: 27 pages, 21 figures, accepted by MNRAS, full tables available as online supplementary material

We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range ($1 \lesssim \mu \lesssim 5$) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, $\mu=5$, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMCs ($\mu=1$): a first burst of star formation with duration, intensity and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remaining trapped in the GMC by the large-scale B-field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters., Comment: 14 Pages, 10 Figures, Submitted to MNRAS

Galaxy pairs constitute the initial building blocks of galaxy evolution, which is driven through merger events and interactions. Thus, the analysis of these systems can be valuable in understanding galaxy evolution and studying structure formation. In this work, we present a new publicly available catalogue of close galaxy pairs identified using photometric redshifts provided by the Physics of the Accelerating Universe Survey (PAUS). To efficiently detect them we take advantage of the high-precision photo$-z$ ($\sigma_{68} < 0.02$) and apply an identification algorithm previously tested using simulated data. This algorithm considers the projected distance between the galaxies ($r_p < 50$ kpc), the projected velocity difference ($\Delta V < 3500$ km/s) and an isolation criterion to obtain the pair sample. We applied this technique to the total sample of galaxies provided by PAUS and to a subset with high-quality redshift estimates. Finally, the most relevant result we achieved was determining the mean mass for several subsets of galaxy pairs selected according to their total luminosity, colour and redshift, using galaxy-galaxy lensing estimates. For pairs selected from the total sample of PAUS with a mean $r-$band luminosity $10^{10.6} h^{-2} L_\odot$, we obtain a mean mass of $M_{200} = 10^{12.2} h^{-1} M_\odot$, compatible with the mass-luminosity ratio derived for elliptical galaxies. We also study the mass-to-light ratio $M/L$ as a function of the luminosity $L$ and find a lower $M/L$ (or steeper slope with $L$) for pairs than the one extrapolated from the measurements in groups and galaxy clusters., Comment: 13 pages, 10 figures, accepted by MNRAS

ASCC 123 is a little-studied young and dispersed open cluster. Recently, we conducted the first research devoted to it. In this paper, we complement our previous work with TESS photometry for the 55 likely members of the cluster. We pay special attention to seven of these high-probability members, all with FGK spectral types, for which we have high-resolution spectra from our preceding work. By studying the TESS light curves of the cluster members we determine the rotational period and the amplitude of the rotational modulation for 29 objects. The analysis of the distribution of the periods allows us to estimate a gyrochronogical age for ASCC 123 similar to that of the Pleiades, confirming the value obtained in our previous investigation. A young cluster age is also suggested by the distribution of variation amplitudes. In addition, for those stars with spectroscopic data we calculate the inclination of their rotation axis. These values appear to follow a random distribution, as already observed in young clusters, with no indication of spin alignment. However, our sample is too small to confirm this on more solid statistical grounds. Finally, for these seven stars we study the level of magnetic activity from the H$\alpha$ and CaII H&K lines. Despite the small number of data points, we find a correlation of the H$\alpha$ and CaII flux with Rossby number. The position of these stars in flux--flux diagrams follows the general trends observed in other active late-type stars., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 14 pages, 17 figures

Our ability to interpret observations of galaxies and trace their stellar, gas, and dust content over cosmic time critically relies on our understanding of how the dust abundance and properties vary with environment. Here, we compute the dust surface density across cosmic times to put novel constraints on simulations of the build-up of dust. We provide observational estimates of the dust surface density consistently measured through depletion methods across a wide range of environments, going from the Milky Way up to $z=5.5$ galaxies. These conservative measurements provide complementary estimates to extinction-based observations. In addition, we introduce the dust surface density distribution function -- in analogy with the cold gas column density distribution functions. We fit a power law of the form: $\log f( \Sigma_{\rm Dust})=-1.92 \times \log \Sigma_{\rm Dust} - 3.65$ which proves slightly steeper than for neutral gas and metal absorbers. This observed relation, which can be computed by simulations predicting resolved dust mass functions through 2D projection, provides new constraints on modern dust models., Comment: 11 pages. Accepted for publication in MNRAS

We report a detailed analysis of the XMM-Newton spectra of six Narrow-Line Seyfert 1 (NLS1) galaxies at redshift z = 0.35--0.92. Compared with the NLS1s at lower redshift in the previously most-studied sample, these NLS1s have larger black hole (BH) masses ($\log\,M_\text{BH}>7.5$) with similar or even lower Eddington ratios. Our extended XMM-Newton sample of NLS1s shows strong soft X-ray excess emission below 2 keV. The quantified soft excess strength does not show an obvious discrepancy from previous studies of the lower-redshift NLS1s. The systematic effect in the measurement of the Eddington ratio mainly lies in the bolometric correction factor. We also tentatively fit the spectra assuming two more physical models for the soft excess: warm Comptonization and relativistic reflection from the inner accretion disk. In the first scenario, we confirm the ubiquity of a warm and optically thick corona. The behavior of a single source can be better explained by relativistic reflection, although we cannot distinguish which model is a more favorable explanation for the soft excess based on the best-fit statistics., 13 pages, 10 figures

We present a CO dynamical estimate of the mass of the super-massive black hole (SMBH) in three nearby early-type galaxies: NGC0612, NGC1574 and NGC4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3-6 observations of the $^{12}$CO(2-1) emission line with spatial resolutions of $14-58$ pc ($0.01''-0.26''$). We detect disc-like CO distributions on scales from $\lesssim200$ pc (NGC1574 and NGC4261) to $\approx10$ kpc (NGC0612). In NGC0612 and NGC1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC1574, allowing us to obtain the first measure of its mass, $M_{\rm BH}=(1.0\pm0.2)\times10^{8}$ M$_{\odot}$ (1$\sigma$ uncertainty). In NGC4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of $(1.62{\pm 0.04})\times10^{9}$ M$_{\odot}$ ($1\sigma$ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC4261. In NGC0612, we are only able to obtain a (conservative) upper limit of $M_{\rm BH}\lesssim3.2\times10^{9}$ M$_{\odot}$. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the $M_{\rm BH}-\sigma_{\star}$ relation., Comment: Main text: 20 pages, 14 Figures; Appendix: 7 pages, 6 Figures. Accepted for publication in MNRAS on 2023 March 28

Ultra-faint dwarf galaxies are among the oldest and most metal-poor galaxies in the cosmos, observed to contain no gas and a high dark matter mass fraction. Understanding the chemical abundance dispersion in such extreme environments could shed light on the very first generations of stars. We present a novel inhomogeneous chemical evolution model, {\tt i-GEtool}, that we apply to two ultra-faint dwarf galaxies, Carina II and Reticulum II, both satellites of the Large Magellanic Cloud. Our model is based on the Monte Carlo sampling of the initial mass function as star formation proceeds in different gas cells of the galaxy volume. We account for the chemical enrichment of Supernova bubbles as they spread in the interstellar medium, causing dispersion in the elemental abundances. We recreate the abundance patterns of $\alpha$- and odd-$\textit{Z}$ elements, predicting two sequences in [C/Fe] and [N/Fe] at all metallicities. Our models underestimate [C/Fe] and [Ti/Fe] because of the large uncertainty in the adopted stellar nucleosynthesis yields. We discuss that the observed C and N abundances had likely been affected by internal mixing processes, which changed the initial surface abundances in the red giants. Our Supernova feedback scheme is responsible for driving galactic outflows, which quench the star formation activity at early times. We predict an average outflow mass-loading factor $\approx 10^{3}$, which extrapolates towards very low galaxy stellar masses the trend observed at high masses. Finally, by combining our model with the MIST isochrone database, we compare our synthetic colour-magnitude diagrams to observations., Comment: 19 Pages, 12 Figures, 1 Table, Accepted to MNRAS

We simulate the formation and collapse of prestellar cores at few-AU resolution in a set of radiation-magneto-hydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming-in onto individual massive prestellar cores within the GMC. We identify two distinct modes of fragmentation: "quasi-spherical" and "filamentary". In both modes the fragments eventually become embedded in a quasi-steady accretion disk or toroid with radii ~ 500-5000 AU and opening angles $H/R \sim 0.5-1$. The disks/toroids are Toomre stable but the accreted pre-existing fragments are found orbiting the outer disk, appearing as disk fragmentation. Each core converts nearly 100 percent of the gas mass into a few massive stars forming near the disk center. Large and massive disks around high-mass stars are supported by magnetic pressure in the outer disk, at radii >200-1000 AU, and turbulent pressure in the inner disk. The most massive core accretes several times more mass than its initial mass, forming a (proto)star cluster of 8 massive stars enshrouded by a toroid, suggesting a competitive accretion scenario for ultra-high-mass star formation. We also find that the HII regions produced by a single massive star remain trapped in the dense circumstellar disks for a few hundred kiloyears, while the dynamic motions of massive stars in wide binaries or multiple systems displace the stars from the densest parts of the disk, allowing UV radiation to escape producing steady or pulsating bipolar HII regions., 20 pages, 13 figures, submitted to MNRAS. Comments welcome. Simulation movie available at https://chongchonghe.github.io/data/#cores

We present an empirical model of age-dependent photospheric lithium depletion, calibrated using a large, homogeneously-analysed sample of 6200 stars in 52 open clusters, with ages from 2--6000 Myr and $-0.3, Accepted for publication in Monthly Notices of the Royal Astronomical Society. Code available at https://github.com/robdjeff/eagles

This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically-aligned JAM$_{\rm cyl}$ or spherically-aligned JAM$_{\rm sph}$. We also make four assumptions for the models' dark vs. luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically-constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results., Comment: 30 pages, 15 figures, 4 tables, accepted by MNRAS, the catalogue is available here https://manga-dynpop.github.io/

We report the detection of ammonia masers in the non-metastable (6, 3), (7, 5) and (6, 5) transitions, the latter is the first unambiguous maser detection of that transition ever made. Our observations include the first VLBI detection of ammonia maser emission, which allowed effective constrain of the (6, 5) maser brightness temperature. The masers were detected towards G358.931-0.030, a site of 6.7-GHz class~II methanol maser emission that was recently reported to be undergoing a period of flaring activity. These ammonia masers appear to be flaring contemporaneously with the class~II methanol masers during the accretion burst event of G358.931-0.030. This newly detected site of ammonia maser emission is only the twelfth such site discovered in the Milky Way. We also report the results of an investigation into the maser pumping conditions, for all three detected masing transitions, through radiative transfer calculations constrained by our observational data. These calculations support the hypothesis that the ammonia (6, 5) maser transition is excited through high colour temperature infrared emission, with the (6, 5) and (7, 5) transition line-ratio implying dust temperatures >400K. Additionally, we detect significant linearly polarised emission from the ammonia (6, 3) maser line. Alongside our observational and radiative transfer calculation results, we also report newly derived rest frequencies for the ammonia (6, 3) and (6, 5) transitions., Accepted into MNRAS 2023 April 24. 13 pages, 5 figures, 3 tables

Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $ii_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution., 9 pages, 7 figures, 1 table

Double white dwarfs (DWDs) will be the most numerous gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). Most of the Galactic DWDs will be unresolved and will superpose to form a confusion noise foreground, the dominant LISA noise source around $\sim 0.5\mathrm{-}3\,\mathrm{mHz}$. A small fraction of these sources will stand out from the background and be individually detectable. Uniquely among GW sources, a handful of these binaries will be known in advance from electromagnetic (EM) observations and will be guaranteed sources of detectable GWs in the LISA band; these are known as verification binaries (VBs). High-cadence photometric surveys are continuously discovering new VB systems, and their number will continue to grow ahead of the launch of LISA. We analyse, in a fully Bayesian framework, all the currently known VB candidates with the latest design requirements for the LISA mission and find that 25 of the considered sources can be detected within a $4\,\mathrm{yr}$ observation time. We explore what can be expected from GW observations, both alone and in combination with EM observations, and estimate the VB's time to detection in the early months of LISA operations. We also show how VBs can be analysed in the case where their GW signals compete with many other unknown binary signals (both resolved and unresolved) from a realistic Galactic population of DWDs., Comment: 13 pages, 8 figures, plus appendices (accepted in Monthly Notices of the Royal Astronomical Society on 26 April 2023)

Galaxy populations are known to exhibit a strong colour bimodality, corresponding to blue star-forming and red quiescent subpopulations. The relative abundance of the two populations has been found to vary with stellar mass and environment. In this paper, we explore the effect of environment considering different types of measurements. We choose a sample of $49, 911$ galaxies with $0.05 < z < 0.18$ from the Galaxy And Mass Assembly survey. We study the dependence of the fraction of red galaxies on different measures of the local environment as well as the large-scale "geometric" environment defined by density gradients in the surround- ing cosmic web. We find that the red galaxy fraction varies with the environment at fixed stellar mass. The red fraction depends more strongly on local environmental measures than on large-scale geometric environment measures. By comparing the different environmental densities, we show that no density measurement fully explains the observed environmental red fraction variation, suggesting the different measures of environmental density contain different information. We test whether the local environmental measures, when combined together, can explain all the observed environmental red fraction variation. The geometric environment has a small residual effect, and this effect is larger for voids than any other type of geometric environment. This could provide a test of the physics applied to cosmological-scale galaxy evolution simulations as it combines large-scale effects with local environmental impact., Accepted for publication in MNRAS; 16 pages; 10 figures; 2 tables

We model the effect of grain size distribution in a galaxy on the evolution of CO and H$_2$ abundances. The formation and dissociation of CO and H$_2$ in typical dense clouds are modelled in a manner consistent with the grain size distribution. The evolution of grain size distribution is calculated based on our previous model, which treats the galaxy as a one-zone object but includes various dust processing mechanisms in the interstellar medium (ISM). We find that typical dense clouds become fully molecular (H$_2$) when the dust surface area increases by shattering while an increase of dust abundance by dust growth in the ISM is necessary for a significant rise of the CO abundance. Accordingly, the metallicity dependence of the CO-to-H$_2$ conversion factor, $X_\mathrm{CO}$, is predominantly driven by dust growth. We also examine the effect of grain size distribution in the galaxy by changing the dense gas fraction, which controls the balance between coagulation and shattering, clarifying that the difference in the grain size distribution significantly affects $X_\mathrm{CO}$ even if the dust-to-gas ratio is the same. The star formation time-scale, which controls the speed of metal enrichment also affects the metallicity at which the CO abundance rapidly increases (or $X_\mathrm{CO}$ drops). We also propose dust-based formulae for $X_\mathrm{CO}$, which need further tests for establishing their usefulness., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

Measuring the proper motion of the emission component in radio-quiet quasars (RQQs) could help to distinguish between the origins of the radio emission and to understand whether the jet production mechanism is the same in radio-loud quasars (RLQs) and RQQs. PG 1351+640 is one of the few RQQs suitable for proper motion studies: it has two compact components on milli-arcsecond scales, a flat-spectrum core and a steep-spectrum jet; both components are >2 mJy at 5 GHz and are well suited for Very Long Baseline Array (VLBA) observations. We compare recent VLBA observations with that made seventeen years ago and find no significant change in the core-jet separation between 2005 and 2015 (a proper motion of 0.003 mas yr-1). However, the core-jet separation increased significantly between 2015 and 2022, inferring a jet proper motion velocity of 0.063 mas yr-1, which corresponds to an apparent transverse velocity of 0.37c. The result suggests that the jet of the RQQ PG 1351+640 is mildly relativistic and oriented at a relatively small viewing angle., The article has been published by Oxford University Press: https://academic.oup.com/mnrasl/advance-article-abstract/doi/10.1093/mnrasl/slad051/7146829?utm_source=advanceaccess&utm_campaign=mnrasl&utm_medium=email

A complete demographic of active galactic nuclei (AGN) is essential to understand the evolution of the Universe. Optical surveys estimate the population of AGN in the local Universe to be of $\sim$ 4%. However, these results could be biased towards bright sources, not affected by the host galaxy attenuation. An alternative method for detecting these objects is through the X-ray emission. In this work, we aim to complement the AGN population of the optical CALIFA survey (941 sources), by using X-ray data from Chandra, which provides the best spatial resolution to date, essential to isolate the nuclear emission from the host galaxy. We study a total of 138 sources with available data. We find 34 new bonafide AGN and 23 AGN candidates, which could increase the AGN population to 7-10\% among the CALIFA survey. X-rays are particularly useful for low-luminosity AGN since they are excluded by the criterion of large equivalent width of the $H\alpha$ emission line when applied to optical selections. Indeed, placing such a restrictive criteria might cause a loss of up to 70% of AGN sources. X-ray detected sources are preferentially located in the right side of the [$OIII$]/$H\beta$ versus [$NII$]/$H\alpha$ diagram, suggesting that this diagram might be the most reliable at classifying AGN sources. Our results support the idea that multi-wavelength studies are the best way to obtain a complete AGN population., Comment: 18 pages, 8 figures, 7 tables, 5 pages of supplementary online material available. Accepted for publication in MNRAS

One of the most crucial yet poorly constrained parameters in modelling the ionizing emissivity is the escape fraction of photons from star-forming galaxies. Several theoretical and observational studies have been conducted over the past few years, but consensus regarding its redshift evolution has yet to be achieved. We present here the first non-parametric reconstruction of this parameter as a function of redshift from a data-driven reionization model using a Gaussian Process Regression method. Our finding suggests a mild redshift evolution of escape fraction with a mean value of $4\%,7\%,\sim10\%$ at $z=2,6,12$. However, a constant escape fraction of $6-10\%$ at $z\gtrsim 6$ is still allowed by current data and also matches other reionization-related observations. With the detection of fainter high redshift galaxies from upcoming observations of JWST, the approach presented here will be a robust tool to put the most stringent constraint on escape fraction as well as reionization histories., Comment: accepted in MNRAS Letter

We estimated the angular momentum, $J$, of $720$ galaxies in the Local Volume with distances $r < 12$ Mpc. The distribution of the average angular momentum along the Hubble sequence has a maximum at the morphological type $T= 4$ (Sbc), while the dispersion of the $J$-values for galaxies is minimal. Among the Local Volume population, 27 elite spiral galaxies stand out, with an angular momentum greater than 0.15 of the Milky Way, $J > 0.15 J_{MW}$, making the main contribution ($ > 90\%$) to the total angular momentum of galaxies in the considered volume. Using observational data on the kinematics and structure of these galaxies, we determined the direction of their spins. We present the first map of the distribution of the spins of 27 nearby massive spiral galaxies in the sky and note that their pattern does not exhibit significant alignment with respect to the Local Sheet plane. The relationship between the magnitude of the angular momentum and stellar mass of the local galaxies is well represented by a power law with an exponent of ($5/3$) over an interval of $6$ orders of magnitude of the mass of galaxies., 8 pages, 9 figures, accepted in MNRAS

Constraints on the formation and evolution of the Milky Way Galaxy require multidimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large data set and often require careful analysis on a star-by-star basis in asteroseismology. Because spectra are easier to obtain for large samples, being able to determine precise ages from spectra allows for large age samples to be constructed, but spectroscopic ages are often imprecise and contaminated by abundance correlations. Here we present an application of a variational encoder–decoder on cross-domain astronomical data to solve these issues. The model is trained on pairs of observations from APOGEE and Kepler of the same star in order to reduce the dimensionality of the APOGEE spectra in a latent space while removing abundance information. The low dimensional latent representation of these spectra can then be trained to predict age with just ∼1000 precise seismic ages. We demonstrate that this model produces more precise spectroscopic ages ($\sim 22~{{\ \rm per\ cent}}$ overall, $\sim 11~{{\ \rm per\ cent}}$ for red-clump stars) than previous data-driven spectroscopic ages while being less contaminated by abundance information (in particular, our ages do not depend on [α/M]). We create a public age catalogue for the APOGEE DR17 data set and use it to map the age distribution and the age-[Fe/H]-[α/M] distribution across the radial range of the Galactic disc.

We investigate the effects of an oscillating gas filament on the dynamics of its embedded stellar clusters. Motivated by recent observational constraints, we model the host gas filament as a cylindrically symmetrical potential, and the star cluster as a Plummer sphere. In the model, the motion of the filament will produce star ejections from the cluster, leaving star cluster remnants that can be classified into four categories: a) Filament Associated clusters, which retain most of their particles (stars) inside the cluster and inside the filament; b) destroyed clusters, where almost no stars are left inside the filament, and there is no surviving bound cluster; c) ejected clusters, that leave almost no particles in the filament, since the cluster leaves the gas filament; and d) transition clusters, corresponding to those clusters that remain in the filament, but that lose a significant fraction of particles due to ejections induced by filament oscillation. Our numerical investigation predicts that the Orion Nebula Cluster is in the process of being ejected, after which it will most likely disperse into the field. This scenario is consistent with observations which indicate that the Orion Nebula Cluster is expanding, and somewhat displaced from the Integral Shaped Filament ridgeline., Comment: 13 pages, 17 figures

We utilize high-resolution cosmological simulations to reveal that high-redshift galaxies tend to undergo a robust `wet compaction' event when near a `golden' stellar mass of $\sim 10^{10} M_{\odot}$. This is a gaseous shrinkage to a compact star-forming phase, a `blue nugget' (BN), followed by central quenching of star formation to a compact passive stellar bulge, a `red nugget' (RN), and a buildup of an extended gaseous disc and ring. Such nuggets are observed at cosmic noon and seed today's early-type galaxies. The compaction is triggered by a drastic loss of angular momentum due to, e.g., wet mergers, counter-rotating cold streams, or violent disc instability. The BN phase marks drastic transitions in the galaxy structural, compositional and kinematic properties. The transitions are from star-forming to quenched inside-out, from diffuse to compact with an extended disc-ring and a stellar envelope, from dark matter to baryon central dominance, from prolate to oblate stellar shape, from pressure to rotation support, from low to high metallicity, and from supernova to AGN feedback. The central black hole growth, first suppressed by supernova feedback when below the golden mass, is boosted by the compaction, and the black hole keeps growing once the halo is massive enough to lock in the supernova ejecta., 33 pages, 26 figures in the main body (49 pages, 45 figures including appendix)

Despite its crucial role in galaxy evolution, the complex circumgalactic medium (CGM) remains underexplored. Although it is known to be multi-phase, the importance of the molecular gas phase to the total CGM mass budget is, to date, unconstrained. We present the first constraints on the molecular gas covering fraction in the CGM of low-redshift galaxies, using measurements of CO column densities along sightlines towards mm-bright background quasars with intervening galaxies. We do not detect molecular absorption against the background quasars. For the individual, low-redshift, 'normal' galaxy haloes probed here, we can therefore rule out the presence of an extremely molecular gas-rich CGM, as recently reported in high-redshift protoclusters and around luminous active galactic nuclei. We also set statistical limits on the volume filling factor of molecular material in the CGM as a whole, and as a function of radius. ISM-like molecular clouds of ~30 pc in radius with column densities of N(CO) >~ 10^16 cm^-2 have volume filling factors of less than 0.2 per cent. Large-scale smooth gas reservoirs are ruled out much more stringently. The development of this technique in the future will allow deeper constraining limits to be set on the importance (or unimportance) of molecular gas in the CGM., 5 pages, 5 figures, accepted for publication in MNRAS

Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement errors and data sparsity into the stream path uncertainty. The technique is Galactic-model independent, non-parametric, and works on phase-wrapped streams. With this technique, we can uniformly analyze and compare data with simulations, enabling both comparison of simulation techniques and ensemble analysis with stream tracks of many stellar streams. Our method is implemented in the public Python package TrackStream, available at https://github.com/nstarman/trackstream., 16 pages, 11 figures, preprint

We present the splashback radius measurements around the SRG/eROSITA eFEDS X-ray selected galaxy clusters by cross-correlating them with HSC S19A photometric galaxies. The X-ray selection is expected to be less affected by systematics related to projection that affects optical cluster finder algorithms. We use a nearly volume-limited sample of 109 galaxy clusters selected in 0.5-2.0 keV band having luminosity $L_X > 10^{43.5}\,{\rm erg s^{-1} h^{-2}}$ within the redshift $z, Comment: 15 pages, 10 figures

Throughout the passage within the Galactic halo, high-velocity clouds (HVCs) sweep up ambient magnetic fields and form stretched and draped configurations of magnetic fields around them. Many earlier numerical studies adopt spherically symmetric uniform-density clouds as initial conditions for simplicity. However, observations demonstrate that HVCs are clumpy and turbulent. In this paper, we perform 3D magnetohydrodynamic simulations to study the evolution of clouds with initial density distributions described by power-law spatial power spectra. We systematically study the role of (i) the initial density structure, (ii) halo magnetic fields, and (iii) radiative cooling efficiency upon infalling HVCs. We find that (i) the clouds' density structure regulates mixing and mass growth. Uniform clouds grow from the onset of the simulations while clumpy clouds initially lose gas and then grow at later times. Along the same lines, the growth curve of clumpy clouds depends on the slope of the initial density power spectra. (ii) Magnetic fields suppress hydrodynamic instabilities and the growth of small-scale structures. As a result, magnetized clouds develop long filaments extended along the streaming direction whereas non-magnetized clouds are fragmented into many small clumps. (iii) Efficient cooling keeps the main cloud body more compact and produces decelerated dense clumps condensed from the halo gas. This work potentially helps us understand and predict the observed properties of HVCs such as the detectability of magnetized clouds, the presence of decelerated HI structures associated with HVC complexes and small-scale features, and a possible link between the origin and the fate of HVCs., Comment: 21 pages, 13 figures, Accepted to MNRAS

The evolution of star-forming galaxies at high redshifts is very sensitive to the strength and nature of stellar feedback. Using two sets of cosmological, zoom-in simulations from the VELA suite, we compare the effects of two different models of feedback: with and without kinetic feedback from the expansion of supernovae shells and stellar winds. At a fixed halo mass and redshift, the stellar mass is reduced by a factor of 1-3 in the models with stronger feedback, so the stellar-mass-halo-mass relation is in better agreement with abundance matching results. On the other hand, the three-dimensional shape of low-mass galaxies is elongated along a major axis in both models. At a fixed stellar mass, Ms 300 Myr), quenched, stellar (or gas-poor) clumps is absent in the model with strong feedback. On the other hand, giant star-forming clumps with intermediate ages (age_c= 100 - 300 Myr) can survive for several disc dynamical times, independently of feedback strength. The evolution through compaction followed by quenching in the plane of central surface density and specific star-formation rate is similar under the two feedback models., 14 pages, 12 figures, accepted version at MNRAS

Financial support from the Spanish Ministerio de Economía y Competitividad under grant PID2019-107408GB-C44, from Junta de Andalucia under project P18-FR-2664, and also from the grant CEX2021-001131-S funded by MCIN/AEI/ 10.13039/501100011033. P P thanks Fundaçao para a Ciencia e a Tecnologia (FCT) for managing research funds graciously provided to Portugal by the EU. This work was supported through FCT grants UID/FIS/04434/2019, UIDB/04434/2020, UIDP/04434/2020 and the project ‘Identifying the Earliest Supermassive Black Holes with ALMA (IdEaS with ALMA) (PTDC/FIS-AST/29245/2017)’ SDP is grateful to the Fonds de Recherche du Québec - Nature et Technologies. SDP also acknowledges financial support from Juan de la Cierva Formación Fellowship (FJC2021-047523-I) financed by MCIN/AEI/10.130395011000110333 and by the European Union "NextGenerationEU"/PRTR. This study is based on observations collected at the Centro As- tronómico Hispano en Andalucía (CAHA) at Calar Alto, Spain, operated jointly by the Instituto de Astrofísica de Andalucía (CSIC) and Junta de Andalucía. The CAHA Archive is part of the Spanish Virtual Observatory project funded by MCIN/AEI/10.13039/501100011033 through grant PID2020-112949GB-I00 CAB (INTA-CSIC). Part based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Fa-cility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). This research has made use of the NASA/IPAC Extragalactic Data base (NED), which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology., Using the PMAS Integral Field Unit on the Calar Alto 3.5 m telescope, we observed the southern component (Markarian 59) of the ‘cometary’ starburst galaxy NGC 4861. Mrk 59 is centred on a giant nebula and concentration of stars 1 kpc in diameter. Strong H α emission points to a star-formation rate (SFR) at least 0.47 M⊙ yr−1. Mrk 59 has a very high [O III] λ5007/H β ratio, reaching 7.35 in the central nebula, with a second peak at a star-forming hotspot further north. Fast outflows are not detected but nebular motion and galaxy rotation produce relative velocities up to 40 km s−1. Spectral analysis of different regions with ‘Fitting Analysis using Differential evolution Optimization’ (FADO) finds that the stars in the central and ‘spur’ nebulae are very young, ≤ 125 Myr with a large < 10 Myr contribution. Older stars (∼ 1 Gyr) make up the northern disk component, while the other regions show mixtures of 1 Gyr age with very young stars. This and the high specific SFR ∼ 3.5 Gyr−1 imply a bimodal star formation history, with Mrk 59 formed in ongoing starbursts fuelled by a huge gas inflow, turning the galaxy into an asymmetric ‘green pea’ or blue compact dwarf. We map the He II λ4686 emission, and identify a broad component from the central nebula, consistent with the emission of ∼300 Wolf–Rayet stars. About a third of the He II λ4686 flux is a narrow line emitted from a more extended area covering the central and spur nebulae, and may have a different origin., Centro Astronómico Hispano en Andalucía, IPAC, National Aeronautics and Space Administration NASA, European Commission PTDC/FIS-AST/29245/2017, UID/FIS/04434/2019 EC, European Space Agency ESA, Ministerio de Economía y Competitividad PID2019-107408GB-C44 MINECO, Ministerio de Ciencia e Innovación MICINN, Junta de Andalucía CEX2021-001131-S, P18-FR-2664, Agencia Estatal de Investigación AEI, Instituto de Astrofísica de Andalucía PID2020-112949GB-I00 CAB (INTA-CSIC IAA

Low mass star formation regions are unlikely to fully populate their initial mass functions, leading to a deficit of massive stars. In binary stellar populations, the full range of binary separations and mass ratios will also be underpopulated. To explore the effects of stochastic sampling in the integrated light of stellar clusters, we calculate models at a broad range of cluster masses, from 10^2 to 10^7 M_sun, using a binary stellar population synthesis code. For clusters with stellar masses less than 10^5 M_sun, observable quantities show substantial scatter and their mean properties reflect the expected deficit of massive stars. In common with previous work, we find that purely stochastic sampling of the initial mass function appears to underestimate the mass of the most massive star in known clusters. However, even with this constraint, the majority of clusters likely inject sufficient kinetic energy to clear their birth clusters of gas. For quantities which directly measure the impact of the most massive stars, such as N_{ion}, xi_{ion} and beta_{UV}, uncertainties due to stochastic sampling dominate over those from the IMF shape or distribution of binary parameters, while stochastic sampling has a negligible effect on the stellar continuum luminosity density., Comment: 15 pages and 8 page appendix. Accepted for publication in MNRAS

Characterizing prestellar cores in star-forming regions is an important step towards the validation of theoretical models of star formation. Thanks to their sub-arcsecond resolution, ALMA observations can potentially provide samples of prestellar cores up to distances of a few kpc, where regions of massive star formation can be targeted. However, the extraction of real cores from dust-continuum observations of turbulent star-forming clouds is affected by complex projection effects. In this work, we study the problem of core extraction both in the idealized case of column-density maps and in the more realistic case of synthetic 1.3\,mm ALMA observations. The analysis is carried out on 12 regions of high column density from our 250 pc simulation. We find that derived core masses are highly unreliable, with only {\em a weak correlation between the masses of cores selected in the synthetic ALMA maps and those of the corresponding three-dimensional cores}. The fraction of real three-dimensional cores detected in the synthetic maps increases monotonically with mass and remains always below 50\%. Above $\sim 1\,M_{\odot}$, the core mass function derived from the column-density maps is steeper than that of the three-dimensional cores, while the core mass function from the synthetic ALMA maps has a slope closer to that of the real three-dimensional cores. Because of the mass uncertainties, proper guidance from realistic simulations is essential if ALMA observations of protoclusters at kpc distances are to be used to test star-formation models., Comment: 20 pages, 21 figures, MNRAS, in press

We present polarization observations of the young supernova remnant (SNR) Cas A using the High-resolution Airborne Wideband Camera-Plus (HAWC+) instrument onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization map at 154 microns reveals dust grains with strong polarization fractions (5 - 30 percent), supporting previous measurements made over a smaller region of the remnant at 850 microns. The 154 microns emission and the polarization signal is coincident with a region of cold dust observed in the southeastern shell and in the unshocked central ejecta. The highly polarized far-IR emission implies the grains are large (greater than 0.14 microns) and silicate-dominated. The polarization level varies across the SNR, with an inverse correlation between the polarization degree and the intensity and smaller polarization angle dispersion for brighter SNR emission. Stronger polarization is detected between the bright structures. This may result from a higher collision rate between the gas and dust producing a lower grain alignment efficiency where the gas density is higher. We use the dust emission to provide an estimate of the magnetic field strength in Cas A using the Davis-Chandrasekhar-Fermi method. The high polarization level is direct evidence that grains are highly elongated and strongly aligned with the magnetic field of the SNR. The dust mass from the polarized region is 0.14+-0.04 Msun, a lower limit of the amount of dust present within the ejecta of Cas A. This result strengthens the hypothesis that core-collapse SNe are an important contributor to the dust mass in high redshift galaxies., Comment: MNRAS, accepted (18 pages with 14 figures)

In this work, we report the presence of rapid intra-night optical variations in both -- flux and polarization of the blazar BL Lacertae during its unprecedented 2020--2021 high state of brightness. The object showed significant flux variability and some color changes, but no firmly detectable time delays between the optical bands. The linear polarization was also highly variable in both -- polarization degree and angle (EVPA). The object was observed from several observatories throughout the world, covering in a total of almost 300 hours during 66 nights. Based on our results, we suggest, that the changing Doppler factor of an ensemble of independent emitting regions, travelling along a curved jet that at some point happens to be closely aligned with the line of sight can successfully reproduce our observations during this outburst. This is one of the most extensive variability studies of the optical polarization of a blazar on intra-night timescales., 23 pages,7 figures, 5 Tables (2 as appendix). Accepted for publication in MNRAS

We report on the case of optical periodic modulation discovered in two Active Galactic Nuclei (AGN) and one candidate AGN. Analyzing the archival optical data obtained from large transient surveys, namely the Catalina Real-Transient Survey (CRTS) and the Zwicky Transient Facility (ZTF), we find periodicities of 2169.7, 2103.1, and 1462.6\,day in sources J0122+1032, J1007+1248 (or PG~1004+1248), and J2131$-$1127, respectively. The optical spectra of the first two indicate that the first is likely a blazar and the second a type 1 Seyfert galaxy, and while no spectroscopic information is available for the third one, its overall properties suggest that it is likely an AGN. In addition, mid-infrared (MIR) light curve data of the three sources, taken by the Wide-field InfraredSurvey Explorer (WISE), are also analyzed. The light curves show significant variations, but not appearing related to the optical periodicities. Based on the widely-discussed super-massive black-hole binary (SMBHB) scenario, we discuss the origin of the optical modulation. Two possible interesting features, an additional 162-day short optical periodicity in J2131$-$1127 and the consistency the X-ray flux variations of J1007+1248 with its optical periodicity, are also discussed within the SMBHB scenario., Comment: 10 pages, 6 figures, 5 tables, accepted for publication in MNRAS

Recent data from the James Webb Space Telescope suggest that there are realistic prospects for detecting the earliest generation of stars at redshift ~20. These metal-poor, gaseous Population III (Pop III) stars are likely in the mass range 10-1000 solar masses. We develop a framework for calculating the abundances of Pop III stars as well as the distribution of the most massive Pop III stars based on an application of extreme-value statistics. Our calculations use the star formation rate density from a recent simulation to calibrate the star-formation efficiency from which the Pop III stellar abundances are derived. Our extreme-value modelling suggests that the most massive Pop III stars at redshifts 10-20 are likely to be $\gtrsim10^3-10^4\,{\rm M}_\odot$. Such extreme Pop III stars were sufficiently numerous to be the seeds of supermassive black holes at high redshifts and possibly source detectable gravitational waves. We conclude that the extreme-value formalism provides an effective way to constrain the stellar initial mass function., Comment: MNRAS accepted

We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of AGN feedback at z~6, and a kinetic feedback mode acts in each galaxy by z~4.7 to complete the quenching process, which occurs on a time-scale of ~0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z=0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z=0, these descendants are the centres of galaxy clusters with average stellar masses of 8 x 10^{11} M_sun. We make predictions for the first quenched galaxies to be located by the James Webb Space Telescope (JWST)., 6 pages, 4 figures

We have found evidence of bulk velocities following active galactic nucleus (AGN) bubbles in the Virgo cluster and galaxy motions in the Centaurus cluster. In order to increase the sample and improve our understanding of the intracluster medium (ICM), we present the results of a detailed mapping of the Ophiuchus cluster with {\it XMM-Newton} to measure bulk flows through very accurate Fe~K measurements. To measure the gas velocities we use a novel EPIC-pn energy scale calibration, which uses the Cu K$\alpha$ instrumental line as reference for the line emission. We created 2D spectral maps for the velocity, metallicity, temperature, density, entropy and pressure with a spatial resolution of 0.25$'$ ($\sim 26$~kpc). The ICM velocities in the central regions where AGN feedback is most important are similar to the velocity of the brightest cluster galaxy (BCG). We have found a large interface region where the velocity changes abruptly from blueshifted to redshifted gas which follows a sharp surface brightness discontinuity. We also found that the metallicities and temperatures do not change as we move outwards from the giant radio fossil previously identified in radio observations of the cluster. Finally, we have found a contribution from the kinetic component of $, Comment: 15 pages, 20 figures. This paper is part of a series on the ICM velocity structure using XMM-Newton observations. Related series papers: arXiv:2109.06213, arXiv:2203.12635

I describe the selection and initial characterisation of 20 eclipsing binary stars that are suitable for calibration and testing of stellar models and data analysis algorithms used by the PLATO mission and spectroscopic surveys. The binary stars selected are F-/G-type dwarf stars with M-type dwarf companions that contribute less than 2% of the flux at optical wavelengths. The light curves typically show well-defined total eclipses with very little variability between the eclipses. I have used near-infrared spectra obtained by the APOGEE survey to measure the spectroscopic orbit for both stars in HD22064. Combined with an analysis of the TESS light curve, I derive the following masses and radii: $M_1 = 1.35 \pm 0.03 M_{\odot}$, $M_2 = 0.58 \pm 0.01 M_{\odot}$, $R_1 = 1.554 \pm 0.014 R_{\odot}$, $R_2 = 0.595 \pm 0.008 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3, I find that the primary star's angular diameter is $\theta = 0.1035 \pm 0.0009 $ mas. The apparent bolometric flux of the primary star is ${\mathcal F}_{\oplus,0} = (7.51\pm 0.09)\times10^{-9}$ erg cm$^{-2}$ s$^{-1}$. Hence, this F2V star has an effective temperature $T_{\rm eff,1} = 6763{\rm\,K} \pm 39{\rm \,K}$. HD22064 is an ideal benchmark star that can be used for "end-to-end" tests of the stellar parameters measured by large-scale spectroscopic surveys, or stellar parameters derived from asteroseismology with PLATO. The techniques described here for HD22064 can be applied to the other eclipsing binaries in the sample in order to create an all-sky network of such benchmark stars., Comment: 13 pages, 9 figures. Accepted for publication in MNRAS. Revised and re-submitted version. Includes supplementary information in the source files. arXiv admin note: text overlap with arXiv:2205.01466

Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies (H IIs and absorption line galaxies, ALGs). Dividing our sample into active and inactive, they define distinct (radio core luminosity)$-$(bulge mass), L_R,core-M_*,bulge, relations, with a mass turnover at M_*, bulge ~ 10^(9.8 +- 0.3) M_sun (supermassive black hole mass M_BH ~ 10^(6.8 +- 0.3) M_sun), which marks the transition from AGN-dominated nuclear radio emission in more massive bulges to that mainly driven by stellar processes in low-mass bulges. None of our 10/173 bulgeless galaxies host an AGN. The AGN fraction increases with increasing M_*, bulge such that f_optical_AGN $\propto$ M_*,bulge^(0.24 +- 0.06) and f_radio_AGN $\propto$ M_*,bulge^(0.24 +- 0.05). Between M_*,bulge ~ 10^8.5 and 10^11.3 M_sun, f_optical_AGN steadily rises from 15 +- 4 to 80 +- 5 per cent. We find that at fixed bulge mass, the radio loudness, nuclear radio activity and the (optical and radio) AGN fraction exhibit no dependence on environment. Radio-loud hosts preferentially possess an early-type morphology than radio-quiet hosts, the two types are however indistinguishable in terms of bulge S\'ersic index and ellipticity, while results on the bulge inner logarithmic profile slope are inconclusive. We finally discuss the importance of bulge mass in determining the AGN triggering processes, including potential implications for the nuclear radio emission in nearby galaxies., Comment: 27 pages, 15 figures, 4 tables, accepted for publication in MNRAS