Search

The Gaia HR diagram shows the presence of apparently young stars at high tangential velocities. Using a simple analytical model I show that these stars are likely to be blue stragglers. Once normalized to red giant stars, the fraction of nearby halo blue stragglers is of order 20 percent, and remarkably close to that measured in dwarf galaxies. Motivated by this similarity, I apply to field blue stragglers scaling relations inferred from blue stragglers in dwarf galaxies. Doing this for the Milky Way halo returns an average stellar density of 3.4 x 10^-5 Msun/pc^3 and a dark matter density of ~0.006 Msun/pc^3 ~ 0.22 GeV/cm^3 within 2 kpc from the Sun. These values compare favourably to other determinations available in the literature, but are based on an independent set of assumptions. A few considerations of this methodology are discussed, most notably that the correlation between the dark matter halo core-density and stellar mass seen in dwarf galaxies seems to hold also for the nearby Milky Way halo., Comment: ApJ accepted. Video of Figure 1 at https://www.dropbox.com/s/3z183bjti0p7vi4/video.mp4?dl=0

We present results on the extra-tidal features of the Milky Way globular cluster NGC 7099, using deep gr photometry obtained with the Dark Energy Camera (DECam). We reached nearly 6 mag below the cluster Main Sequence (MS) turnoff, so that we dealt with the most suitable candidates to trace any stellar structure located beyond the cluster tidal radius. From star-by-star reddening corrected color-magnitude diagrams (CMDs) we defined four adjacent strips along the MS, for which we built the respective stellar density maps, once the contamination by field stars was properly removed. The resulting field star cleaned stellar density maps show a short tidal tail and some scattered debris. Such extra-tidal features are hardly detected when much shallower Gaia DR2 data sets are used and the same CMD field star cleaning procedure is applied. Indeed, by using 2.5 magnitudes below the cluster MS turnoff as the faintest limit (G < 20.5 mag), cluster members turned out to be distributed within the cluster's tidal radius, and some hints for field star density variations are found across a circle of radius 3.5deg centered on the cluster and with similar CMD features as cluster stars. The proper motion distribution of these stars is distinguishable from that of the cluster, with some superposition, which resembles that of stars located beyond 3.5deg from the cluster center., 7 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

We study the mass fallback rate of tidally disrupted stars on marginally bound and unbound orbits around a supermassive black hole (SMBH) by performing three-dimensional smoothed particle hydrodynamic (SPH) simulations with three key parameters: the star is modeled by a polytrope with two different indexes ($n=1.5$ and $3$). The stellar orbital properties are characterized by five orbital eccentricities ranging from $e=0.98$ to $1.02$ and five different penetration factors ranging from $\beta=1$ to $3$, where $\beta$ represents the ratio of the tidal disruption to pericenter distance radii. We derive analytic formulae for the mass fallback rate as a function of the stellar density profile, orbital eccentricity, and penetration factor. Moreover, two critical eccentricities to classify tidal disruption events (TDEs) into five different types: eccentric ($e, Comment: 42 pages 11 figures, accepted for publication in ApJ

The stellar birth environment can significantly shorten protoplanetary disc (PPD) lifetimes due to the influence of stellar feedback mechanisms. The degree to which these mechanisms suppress the time and mass available for planet formation is dependent on the local far-ultraviolet (FUV) field strength, stellar density, and ISM properties. In this work, we present the first theoretical framework quantifying the distribution of PPD dispersal time-scales as a function of parameters that describe the galactic environment. We calculate the probability density function for FUV flux and stellar density in the solar neighbourhood. In agreement with previous studies, we find that external photoevaporation is the dominant environment-related factor influencing local stellar populations after the embedded phase. Applying our general prescription to the Central Molecular Zone of the Milky Way (i.e. the central ~250 pc), we predict that 90% of PPDs in the region are destroyed within 1 Myr of the dispersal of the parent molecular cloud. Even in such dense environments, we find that external photoevaporation is the dominant disc depletion mechanism over dynamical encounters between stars. PPDs around low-mass stars are particularly sensitive to FUV-induced mass loss, due to a shallower gravitational potential. For stars of mass ~1 $M_\odot$, the solar neighbourhood lies at approximately the highest gas surface density for which PPD dispersal is still relatively unaffected by external FUV photons, with a median PPD dispersal timescale of ~4 Myr. We highlight the key questions to be addressed to further contextualise the significance of the local galactic environment for planet formation., Comment: 22 pages, 12 figures; accepted for publication in MNRAS

We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes., Comment: 22 pages, 13 figures, accepted for publication in MNRAS

We present a large suite of MHD simulations of turbulent, star-forming giant molecular clouds(GMCs) with stellar feedback, extending previous work by simulating 10 different random realizations for each point in the parameter space of cloud mass and size. It is found that oncethe clouds disperse due to stellar feedback, both self-gravitating star clusters and unbound stars generally remain, which arise from the same underlying continuum of substructured stellar density, ie. the hierarchical cluster formation scenario. The fraction of stars that are born within gravitationally-bound star clusters is related to the overall cloud star formation efficiency set by stellar feedback, but has significant scatter due to stochastic variations in the small-scale details of the star-forming gas flow. We use our numerical results to calibrate a model for mapping the bulk properties (mass, size, and metallicity) of self-gravitating GMCs onto the star cluster populations they form, expressed statistically in terms of cloud-level distributions. Synthesizing cluster catalogues from an observed GMC catalogue in M83, we find that this model predicts initial star cluster masses and sizes that are in good agreement with observations, using only standard IMF and stellar evolution models as inputs for feedback. Within our model, the ratio of the strength of gravity to stellar feedback is the key parameter setting the masses of star clusters, and of the various feedback channels direct stellar radiation(photon momentum and photoionization) is the most important on GMC scales., Comment: Submitted to MNRAS. 20 pages, 14 figures. Feedback welcome :)

Most stars form in a clustered environment. Both single and binary stars will sometimes encounter planetary systems in such crowded environments. Encounter rates for binaries may be larger than for single stars, even for binary fractions as low as 10-20 per cent. In this work, we investigate scatterings between a Sun-Jupiter pair and both binary and single stars as in young clusters. We first perform a set of simulations of encounters involving wide ranges of binaries and single stars, finding that wider binaries have larger cross sections for the planet's ejection. Secondly, we consider such scatterings in a realistic population, drawing parameters for the binaries and single stars from the observed population. The scattering outcomes are diverse, including ejection, capture/exchange and collision. The binaries are more effective than single stars by a factor of several or more in causing the planet's ejection and collision. Hence, in a cluster, as long as the binary fraction is larger than about 10 per cent, the binaries will dominate the scatterings in terms of these two outcomes. For an open cluster of a stellar density 50 pc$^{-3}$, a lifetime 100 Myr and a binary fraction 0.5, we estimate that of the order of 1 per cent of the Jupiters are ejected, 0.1 per cent collide with a star, 0.1 per cent change ownership and 10 per cent of the Sun-Jupiter pairs acquire a stellar companion during scatterings. These companions are typically 1000s of au distant and in half of the cases (so 5 per cent of all Sun-Jupiter pairs), they can excite the planet's orbit through Kozai--Lidov mechanism before stripped by later encounters. Our result suggests that the Solar System may have once had a companion in its birth cluster., Comment: to appear in MNRAS, 9 figures

Statistics of low-mass pre-main sequence binaries in the Orion OB1 association with separations ranging from 0.6" to 20" (220 to 7400 au at 370 pc) are studied using images from the VISTA Orion mini-survey and astrometry from Gaia. The input sample based on the CVSO catalog contains 1137 stars of K and M spectral types (masses between 0.3 and 0.9 Msun), 1021 of which are considered to be association members. There are 135 physical binary companions to these stars with mass ratios above ~0.13. The average companion fraction is 0.09+-0.01 over 1.2 decades in separation, slightly less than, but still consistent with, the field. We found a difference between the Ori OB1a and OB1b groups, the latter being richer in binaries by a factor 1.6+-0.3. No overall dependence of the wide-binary frequency on the observed underlying stellar density is found, although in the Ori OB1a off-cloud population these binaries seem to avoid dense clusters. The multiplicity rates in Ori OB1 and in sparse regions like Taurus differ significantly, hinting that binaries in the field may originate from a mixture of diverse populations., Comment: Accepted for publicaion in The Astronomical Journal. 19 pages, 22 figures, 7 tables

We study the atomic gas content and the central stellar mass concentration for a sample of almost 1500 nearby galaxies to investigate the nature of starbursts and the influence of galaxy-galaxy interactions on star formation (SF). We used a sample of catalogued interacting and non-interacting galaxies in the S$^4$G survey - along with archival HI gas masses, stellar masses ($M_{\ast}$), and SF rates from $IRAS$ far-infrared fluxes - and calculate depletion times and gas fractions. We traced the central stellar mass concentration from the inner slope of the stellar component of the rotation curves. Starbursts are defined as galaxies with a factor $>4$ enhanced SF rate relative to a control sample of non-interacting galaxies which are $\pm 0.2$ dex in stellar mass and $\pm 1$ in $T$-type. Starbursts are mainly early-type ($T\lesssim 5$), massive spiral galaxies ($M_{\ast}\gtrsim 10^{10}M_{\odot}$) that are not necessarily interacting. For a given stellar mass bin, starbursts are characterised by lower gas depletion times, similar gas fractions, and larger central stellar mass concentrations than non-starburst galaxies. The global distributions of gas fraction and gas depletion time are not statistically different for interacting and non-interacting galaxies. However, in the case of currently merging galaxies, the median gas depletion time is a factor of $0.4 \pm 0.2$ that of control sample galaxies, and their SF rates are a factor of $1.9 \pm 0.5$ enhanced, even though the median gas fraction is similar. Starbursts present long-lasting SF in the circumnuclear regions that causes an enhancement of the central stellar density at $z\approx0$ in both interacting and non-interacting systems. Starbursts have low gas depletion timescales, yet similar gas fractions as normal main-sequence galaxies. Galaxy mergers cause a moderate enhancement of the SF efficiency (Abridged)., Accepted for publication in A&A (February 21, 2020). 6 pages, 5 figures

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$-$340$ nBHBs to be observed by Gaia, although the numbers may drop to $\sim10$-$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For LAMOST we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy., Comment: 23 pages, 15 figures, 8 tables

We present a new estimate of the mass of the Milky Way, inferred via a Bayesian approach by making use of tracers of the circular velocity in the disk plane and stars in the stellar halo, as from the publicly available {\tt galkin} compilation. We use the rotation curve method to determine the dark matter distribution and total mass under different assumptions for the dark matter profile, while the total stellar mass is constrained by surface stellar density and microlensing measurements. We also include uncertainties on the baryonic morphology via Bayesian model averaging, thus converting a potential source of systematic error into a more manageable statistical uncertainty. We evaluate the robustness of our result against various possible systematics, including rotation curve data selection, uncertainty on the Sun's velocity $V_0$, dependence on the dark matter profile assumptions, and choice of priors. We find the Milky Way's dark matter virial mass to be $\log_{10}M_{200}^{\rm DM}/ {\rm M_\odot} = 11.92^{+0.06}_{-0.05}{\rm(stat)}\pm{0.28}\pm0.27{\rm(syst)}$ ($M_{200}^{\rm DM}=8.3^{+1.2}_{-0.9}{\rm(stat)}\times10^{11}\,{\rm M_\odot}$). We also apply our framework to Gaia DR2 rotation curve data and find good statistical agreement with the above results., 32 pages, 10 figures, 7 tables, matches published version

Knowledge about the large-scale distribution of galaxies is far from complete in the Zone of Avoidance, which is mostly due to high interstellar extinction and to source confusion at lower Galactic latitudes. Past near-infrared (NIR) surveys, such as the Two Micron All Sky Survey (2MASS), have shown the power of probing large-scale structure at these latitudes. Our aim is to map the galaxy distribution across the Southern Galactic plane using the VISTA Variables in the V\'ia L\'actea Survey (VVV), which reach 2 to 4 magnitudes deeper than 2MASS. We used SExtractor + PSFEx to identify extended objects and to measure their sizes, the light concentration index, magnitudes, and colours. Morphological and colour constraints and visual inspection were used to confirm galaxies. We present the resulting VVV NIR Galaxy Catalogue of 5563 visually confirmed galaxies, of which only 45 were previously known. This is the largest catalogue of galaxies towards the Galactic plane, with 99% of these galaxies being new discoveries. We found that the galaxy density distribution closely resembled the distribution of low interstellar extinction of the existing NIR maps. We also present a description of the 185 2MASS extended sources observed in the region, of which 16% of these objects had no previous description, which we have now classified. We conclude that interstellar extinction and stellar density are the main limitations for the detection of background galaxies in the Zone of Avoidance. The VVV NIR Galaxy Catalogue is a new data set providing information for extragalactic studies in the Galactic plane., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal 21 pages

We describe the circumstances that led to the discovery of Kepler-36b, and the subsequent characterization of its host planetary system. The Kepler-36 system is remarkable for its physical properties: the close separation of the planets, the contrasting densities of the planets despite their proximity, and the short chaotic timescale. Its discovery and characterization was also remarkable for the novelty of the detection technique and for the precise characterization due to the large transit-timing variations caused by the close proximity of the planets, as well as the precise stellar parameters due to asteroseismology. This was the first multi-planet system whose transit data was processed using a fully consistent photometric-dynamical model, using population Markov Chain Monte Carlo techniques to precisely constrain system parameters. Amongst those parameters, the stellar density was found to be consistent with a complementary, concurrent asteroseismic analysis. In a first, the 3D orientation of the planets was constrained from the lack of transit-duration variations. The system yielded insights into the composition and evolution of short-period planet systems. The denser planet appears to have an Earth-like composition, with uncertainties comparable to the highest precision rocky exoplanet measurements, and the planet densities foreshadowed the rocky/gaseous boundary. The formation of this system remains a mystery, but should yield insights into the migration and evolution of compact exoplanet systems., 31 pages, 9 figures. In press with New Astronomy Reviews special issue on key discoveries with Kepler

We use astrometry and broad-band photometry from Data Release 2 of the ESA's Gaia mission to map out low surface-brightness features in the stellar density distribution around the Large and Small Magellanic Clouds. The LMC appears to have grown two thin and long stellar streams in its Northern and Southern regions, highly reminiscent of spiral arms. We use computer simulations of the Magellanic Clouds' in-fall to demonstrate that these arms were likely pulled out of the LMC's disc due to the combined influence of the SMC's most recent fly-by and the tidal field of the Milky Way., Comment: Submitted to MNRAS Letters. Comments welcome!

We calculate the mean density profiles for luminous and dark matter on distance scales $D \sim(1 - 100)$ Mpc around us using recent all-sky catalogs of galaxy groups. Within the Local Volume $( D < 11 ~\rm Mpc)$ we derived the mean stellar density $\Omega_*= 0.44\%$ in the critical density units and the mean total matter density $\Omega_m = 0.17$. In the sphere with a radius of 40 Mpc these quantities drop to $\Omega_* = 0.24 - 0.32\%$ and $\Omega_m = 0.09 - 0.14$. In a larger volume within $D\sim135$ Mpc the discussed densities become more uncertain: $\Omega_* = 0.20 - 0.24\%$ and $\Omega_m = 0.05 - 0.16$. We summarize that the major part of the cosmic dark matter locates outside the virial and collapsing zones of groups and clusters., Comment: 11 pages, 8 figures, 4 tables, accepted by AN

We investigate the effect of microlensing on the standardisation of strongly lensed Type Ia supernovae (GLSNe Ia). We present predictions for the amount of scatter induced by microlensing across a range of plausible strong lens macromodels. We find that lensed images in regions of low convergence, shear and stellar density are standardisable, where the microlensing scatter is < 0.15 magnitudes, comparable to the intrinsic dispersion of for a typical SN Ia. These standardisable configurations correspond to asymmetric lenses with an image located far outside the Einstein radius of the lens. Symmetric and small Einstein radius lenses (< 0.5 arcsec) are not standardisable. We apply our model to the recently discovered GLSN Ia iPTF16geu and find that the large discrepancy between the observed flux and the macromodel predictions from More et al. (2017) cannot be explained by microlensing alone. Using the mock GLSNe Ia catalogue of Goldstein et al. (2017), we predict that ~ 22% of GLSNe Ia discovered by LSST will be standardisable, with a median Einstein radius of 0.9 arcseconds and a median time-delay of 41 days. By breaking the mass-sheet degeneracy the full LSST GLSNe Ia sample will be able to detect systematics in H0 at the 0.5% level., 11 pages, 8 Figures. Accepted by MNRAS May 17 2018

We measure and model the wavelength dependence of the PSF in the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey. We find that PSF chromaticity is present in that redder stars appear smaller than bluer stars in the $g, r,$ and $i$-bands at the 1-2 per cent level and in the $z$ and $y$-bands at the 0.1-0.2 per cent level. From the color dependence of the PSF, we fit a model between the monochromatic PSF trace radius, $R$, and wavelength of the form $R(\lambda)\propto \lambda^{b}$. We find values of $b$ between -0.2 and -0.5, depending on the epoch and filter. This is consistent with the expectations of a turbulent atmosphere with an outer scale length of $\sim 10-100$ m, indicating that the atmosphere is dominating the chromaticity. We find evidence in the best seeing data that the optical system and detector also contribute some wavelength dependence. Meyers and Burchat (2015) showed that $b$ must be measured to an accuracy of $\sim 0.02$ not to dominate the systematic error budget of the Large Synoptic Survey Telescope (LSST) weak lensing (WL) survey. Using simple image simulations, we find that $b$ can be inferred with this accuracy in the $r$ and $i$-bands for all positions in the LSST field of view, assuming a stellar density of 1 star arcmin$^{-2}$ and that the optical PSF can be accurately modeled. Therefore, it is possible to correct for most, if not all, of the bias that the wavelength-dependent PSF will introduce into an LSST-like WL survey., Comment: 14 pages, 10 figures. Submitted to MNRAS. Comments welcome

We present {complexity} of the Galactic outer disc {by fitting the stellar volume densities of the red giant branch stars with a two-disc component model}. {The discs are confirmed to} extend to $R\sim19$\,kpc. The radial density profile of the discs shows {two breaks at $R\sim11$ and $\sim14$\,kpc, respectively, which} separate the radial profile into three segments with different scale lengths of $2.12\pm0.26$, $1.18\pm0.08$, and $2.72$\,kpc at $R14$\,kpc, respectively. The first {break} is likely due to the sudden drop in the radial profile of the thin disc, which may be an evidence of the radial migration. {Beyond $14$\,kpc, the thick disc becomes prominent and the transition from thin to thick disc leads to the second break.} This implies that the geometrically defined thick disc is more {radially} extended than the thin disc. This is also supported by the larger scale length of the thick disc than that of the thin disc. Meanwhile, {the scale height of the thicker component increases from $0.637_{-0.036}^{+0.056}$ at $R=8$ to $1.284_{-0.079}^{+0.086}$\,kpc at $R=19$\,kpc, showing an intensive flared disc}. Moreover, rich substructures are displayed in the residuals of the stellar density. Among them, the substructures $D14+2.0$ and $O14-1.5$ show a north-south asymmetry, which can be essentially explained by southward shifting of the thick disc. However, no significant overdensity is found for the Monoceros ring. Finally, the thick disc shows a ripple-like feature with unclear origin at $9, Comment: 15 pages, 12 figures, accepted by MNRAS

The Galactic centre Nuclear Star Cluster is one of the densest stellar clusters in the Galaxy. The stars in its inner portions orbit the supermassive black hole associated with compact radio source Sgr~A* at the orbital speeds of several thousand km/s. The B-type star S2 is currently the best case to test the general relativity as well as other theories of gravity based on its stellar orbit. Yet its orbital period of $\sim 16\,{\rm yr}$ and the eccentricity of $\sim 0.88$ yields the relativistic pericentre shift of $\sim 11'$, which is observationally still difficult to reliably measure due to possible Newtonian perturbations as well as reference-frame uncertainties. A naive way to solve this problem is to find stars with smaller pericentre distances, $r_{\rm p}\lesssim 1529$ Schwarzschild radii ($120\,{\rm AU}$), and thus more prominent relativistic effects. In this contribution, we show that to detect stars on relativistic orbits is progressively less likely given the volume shrinkage and the expected stellar density distributions. Finally, one arrives to a sparse region, where the total number of bright stars is expected to fall below one. One can, however, still potentially detect stars crossing this region. In this contribution, we provide a simple formula for the detection probability of a star crossing a sparse region. We also examine an approximate time-scale on which the star reappears in the sparse region, i.e. a "waiting" time-scale for observers., 7 pages, 6 figures, 2 tables, accepted by the Astronomical Notes as a contribution to the proceedings of the IBWS 2018 workshop (Carlsbad, https://www.ibws.cz/), added references

We present for the first time extended stellar density and/or surface brightness radial profiles for almost all the known Large Magellanic Cloud (LMC) old globular clusters (GCs). These were built from DECam images and reach out to ~ 4 times the GCs' tidal radii. The background subtracted radial profiles reveal that the GCs located closer than ~ 5 kpc from the LMC centre contain an excess of stars in their outermost regions with respect to the stellar density expected from a King profile. Such a residual amount of stars - not seen in GCs located farther than ~ 5 kpc from the LMC centre-, as well as the GCs' dimensions, show a clear dependence with the GCs' positions in the galaxy, in the sense that, the farther the GC from the centre of the LMC, the larger both the excess of stars in its outskirts and size. Although the masses of GCs located inside and outside ~ 5 kpc are commensurate, the outermost regions of GCs located closer than ~ 5 kpc from the LMC centre appear to have dynamically evolved more quickly. These outcomes can be fully interpreted in the light of the known GC radial velocity disc-like kinematics, from which GCs have been somehow mostly experiencing the influence of the LMC gravitational field at their respective mean distances from the LMC centre., 13 pages, 16 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society. arXiv admin note: substantial text overlap with arXiv:1712.05198

Context. The measurement of apsidal motion rates in eccentric eclipsing binaries is a unique way to gain insight into the internal structure of stars through the internal density concentration parameter, k2. High-quality physical parameters of the stellar components, together with precise measurements of the advance of the periastron, are needed for the comparison with values derived from models. Aims. As a product of the Transiting Exoplanet Survey Satellite (TESS) mission, high-precision light curves of a large number of eclipsing binaries are now available. Using a selection of well-studied, double-lined eccentric eclipsing binary systems, we aim to determine their apsidal motion rates and place constraints on the internal density concentration and compare with the predictions from state-of-the-art theoretical models. Methods. We computed times of minimum light using the TESS light curves of 34 eclipsing binaries with precise absolute parameters. We used the changing difference over time between primary and secondary eclipse timings to determine the apsidal motion rate. To extend the time baseline, we combined the high-precision TESS timings with reliable archival data. On the other hand, for each component of our sample of double-lined eclipsing binaries, we computed grids of evolutionary stellar models for the observed stellar mass exploring ranges of values of the overshooting parameter fov, the mixing-length parameter, and the metallicity. To find the best solution for the two components we adopted a χ2 statistic to infer the optimal values of the overshooting parameter and the mixing-length parameter. The theoretical internal structure constants to be compared with the observed values were calculated by integrating the differential equations of Radau for each stellar model. Results. We have determined the apsidal motion rate of 27 double-lined eclipsing binaries with precise physical parameters. The obtained values, corrected for their relativistic contribution, yield precise empirical parameters of the internal stellar density concentration. The comparison of these results with the predictions based on new theoretical models shows very good agreement. Small deviations are identified but remain within the observational uncertainties and the path for a refined comparison is indicated. © ESO 2021., The Spanish MEC (ESP2017-87676-C5-2-R, PID2019-107061GB-C64, and PID2019-109522GB-C52) is gratefully acknowledged for its support during the development of this work. A.C. also acknowledges financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). D.B., I.R., and J.C.M. acknowledge support from the Spanish Ministry of Science and Innovation and the European Regional Development Fund through grant PGC2018-098153-B-C33, from the Generalitat de Catalunya/CERCA programme, and from the Agancia de Gestio d'Ajuts Universitaris i de Recerca of the Generalitat de Catalunya. This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France, and of NASA's Astrophysics Data System Abstract Service.

We present the results of our search for variable stars using the long-term Las Cumbres Observatory (LCO) monitoring of white dwarf ZTF J0139+5245 with the two 1.0-m telescope nodes located at McDonald Observatory using the Sinistro imaging instrument. In this search, we find 38 variable sources, of which 27 are newly discovered or newly classified (71%) based on comparisons with previously published catalogs, thereby increasing the number of detections in the field-of-view under consideration by a factor of $\approx$ 2.5. We find that the improved photometric precision per-exposure due to longer exposure time for LCO images combined with the greater time-sampling of LCO photometry enables us to increase the total number of detections in this field-of-view. Each LCO image covers a field-of-view of $26' \times 26'$ and observes a region close to the Galactic plane ($b = -9.4^\circ$) abundant in stars with an average stellar density of $\approx 8$ arcmin$^{-2}$. We perform aperture photometry and Fourier analysis on over 2000 stars across 1560 LCO images spanning 537 days to find 28 candidate BY Draconis variables, 3 candidate eclipsing binaries of type EA, and 7 candidate eclipsing binaries of type EW. In assigning preliminary classifications to our detections, we demonstrate the applicability of the Gaia color-magnitude diagram (CMD) as a powerful classification tool for variable star studies., Accepted for publication in The Astronomical Journal, 14 pages, 7 figures, 1 table

Context.In previous work, we showed how the planar and vertical eccentricities of disc stars,eande′, could be used as indicators of the stars’ kinematic populations. For a local stellar sample drawn from theGaiaDR2 catalogue, these populations were represented geometrically in the eccentricity diagram,e′2vs.e2, approximately separated by straight lines.Aims.In the current work, we propose a new relationship between the star’s perpendicular velocity and its vertical eccentricity, allowing for a reevaluation of the critical vertical eccentricity and maximum height,zmax, specific to each population component.Methods.We approximated the local potential function to be consistent with the actual shape of the curve that relates the maximum vertical speed of a star and its maximum height. The curve corresponds to a non-linear restoring vertical force, where the stiffness decreases with an increase in the maximum height. The constants involved in this fitting, together with the population velocity dispersions, determine the specific region for each population in the eccentricity diagram.Results.The new classification determines 88% of the sample is made up of thin disc stars and 9% of thick disc stars, whereby 3% of the stars have been relabelled, by providing thinner thin and thick discs. Nested thin disc subsamples allow us to estimate Strömberg’s asymmetric drift equation, leading to a heliocentric velocity of the circular orbit ofVc ≈ −12.9 km s−1, an absolute rotation velocity of Θc ≈ 227 km s−1, and a rotation component of the Galactocentric velocity of the Sun at Θ⊙ ≈ 240 km s−1.Conclusions.The thin disc stars of our local sample are characterised based on values 0 ≤ e ≤ 0.32, 0 ≤ e′ ≤ 0.09, andzmax = 0.7 kpc. Disc stars satisfy 0 ≤ e ≤ 0.44, 0 ≤ e′ ≤ 0.18,zmax = 1.5 kpc. The maximum vertical peculiar velocity for disc stars is found to bew0 = 115 km s−1. The assumed potential provides a stellar density of the disc vanishing atz0 = 1.8 kpc. The approximate behaviour in the local disc is that a small decrease in the stiffness is associated with a relative decrease in the limiting velocity, which produces a thinner disc and a loss of stars in the local cylinder, both in a similar proportion to the limiting velocity.

In Stage-IV imaging surveys, a significant amount of the cosmologically useful information is due to sources whose images overlap with those of other sources on the sky. The cosmic shear signal is primarily encoded in the estimated shapes of observed galaxies and thus directly impacted by overlaps. We introduce a framework based on the Fisher formalism to analyze effects of overlapping sources (blending) on the estimation of cosmic shear. For the Rubin Observatory Legacy Survey of Space and Time (LSST), we present the expected loss in statistical sensitivity for the ten-year survey due to blending. We find that for approximately 62% of galaxies that are likely to be detected in full-depth LSST images, at least 1% of the flux in their pixels is from overlapping sources. We also find that the statistical correlations between measures of overlapping galaxies and, to a much lesser extent the higher shot noise level due to their presence, decrease the effective number density of galaxies, $N_{eff}$, by $\sim$18%. We calculate an upper limit on $N_{eff}$ of 39.4 galaxies per arcmin$^2$ in $r$ band. We study the impact of varying stellar density on $N_{eff}$ and illustrate the diminishing returns of extending the survey into lower Galactic latitudes. We extend the Fisher formalism to predict the increase in pixel-noise bias due to blending for maximum-likelihood (ML) shape estimators. We find that noise bias is sensitive to the particular shape estimator and measure of ensemble-average shape that is used, and properties of the galaxy that include redshift-dependent quantities such as size and luminosity., Comment: Accepted for publication in JCAP. 45 pages, 19 figures

We report on observational evidence of an extra-tidal clumpy structure around NGC 288 from an homogeneous coverage of a large area with the Pan-STARRS PS1 database. The extra-tidal star population has been disentangled from that of the Milky Way field by using a cleaning technique that successfully reproduced the stellar density, luminosity function and colour distributions of MW field stars. We have produced the cluster stellar density radial profile and a stellar density map from independent approaches, from which we found results in excellent agreement : the feature extends up to 3.5 times the cluster tidal radius. Previous works based on shallower photometric data sets have speculated on the existence of several long tidal tails, similar to that found in Pal 5. The present outcome shows that NGC 288 could hardly have such tails, but favours the notion that interactions with the MW tidal field has been a relatively inefficient process for stripping stars off the cluster. These results point to the need of a renewed overall study of the external regions of Galactic globular clusters (GGCs) in order to reliably characterise them. Hence, it will be possible to investigate whether there is any connection between detected tidal tails, extra-tidal stellar populations, extent diffuse halo-like structures with the GGCs' dynamical histories in the Galaxy., Comment: 7 pages, 3 figures. Accepted for publication in MNRAS

We characterize the stellar and gas volume density, potential, and gravitational field profiles in the central $\sim$ 0.5 pc of the Orion Nebula Cluster (ONC), the nearest embedded star cluster (or rather, proto-cluster) hosting massive star formation available for detailed observational scrutiny. We find that the stellar volume density is well characterized by a Plummer profile $\rho_{stars}(r) = 5755\,{\rm M}_{\odot}\,{\rm pc}^{-3}\,(1+(r/a)^2)^{-5/2}$, where $a = 0.36$ pc. The gas density follows a cylindrical power law $\rho_{gas}(R) = 25.9\,{\rm M}_{\odot}\,{\rm pc}^{-3}\,(R/{\rm pc})^{-1.775}$. The stellar density profile dominates over the gas density profile inside $r\,\sim\,1$ pc. The gravitational field is gas-dominated at all radii, but the contribution to the total field by the stars is nearly equal to that of the gas at $r\,\sim\,a$. This fact alone demonstrates that the proto-cluster cannot be considered a gas-free system or a virialized system dominated by its own gravity. The stellar proto-cluster core is dynamically young, with an age of $\sim$ 2-3 Myr, a 1D velocity dispersion of $\sigma_{\rm obs} = 2.6$ km s$^{-1}$, and a crossing time of $\sim$ 0.55 Myr. This timescale is almost identical to the gas filament oscillation timescale estimated recently by Stutz & Gould (2016). This provides strong evidence that the proto-cluster structure is regulated by the gas filament. The proto-cluster structure may be set by tidal forces due to the oscillating filamentary gas potential. Such forces could naturally suppress low density stellar structures on scales $\gtrsim\,a$. The analysis presented here leads to a new suggestion that clusters form by an analog of the "slingshot mechanism" previously proposed for stars., Comment: 11 pages, 11 figures, MNRAS accepted

We present extra-tidal features of spatial configuration of stars around three metal-poor globular clusters (NGC 6266, NGC 6273, NGC 6681) located in the Galactic bulge. The wide-field photometric data were obtained in BVI bands with the MOSAIC II camera at CTIO 4 m Blanco telescope. The derived color-magnitude diagrams (CMDs) contain stars in a total 71’×71’area including a cluster and its surrounding field outside of the tidal radius of the cluster. Applying statistical filtering technique, we minimized the field star contaminations on the obtained cluster CMDs and extracted the cluster members. On the spatial stellar density maps around the target clusters, we found overdensity features beyond the tidal radii of the clusters. We also found that the radial density profiles of the clusters show departures from the best-fit King model for their outer regions which support the overdensity patterns.

The astrophysical parameters of four poorly studied open star clusters namely Teutsch 126, Teutsch 54, Teutsch 61 and Czernik 3, have been estimated using the Two Micron All Sky Survey (2MASS) database. The stellar density distributions and color-magnitude diagrams are used to determine their structural parameters (cluster center, cluster radius, core radius, tidal radius, Galactocenteric coordinates and the distance from the Galactic plane). We have also derived age, color excesses, total mass, relaxation time, luminosity and mass function for each clusters. The mass function slopes for these clusters are derived as 1.59 ± 0.62, 1.31 ± 0.60, 1.22 ± 0.75 and 1.62 ± 0.56 for Teutsch 126, Teutsch 54, Teutsch 61 and Czernik 3 respectively. These values are very close with the Salpeter value ( x = 1.35 ) within the errors. The effect of mass-segregation are observed in the clusters Teutsch 126 and Teutsch 61. Estimated values of dynamical relaxation time are less than age of the clusters under study. This concludes that these objects are dynamically relaxed. The possible reason for relaxation may be due to the dynamical evolution or imprint of star formation or both.

Although still under debate, observations generally suggest that dwarf spheroidal (dSph) galaxies exhibit large constant-density cores in the centers, which can hardly be explained by dissipationless cold dark matter simulations without baryonic feedback. Wave dark matter (${\psi {\rm DM}}$), characterized by a single parameter, the dark matter particle mass $m_{\psi}$, predicts a central soliton core in every galaxy arising from quantum pressure against gravity. Here we apply Jeans analysis assuming a soliton core profile to the kinematic data of eight classical dSphs so as to constrain $m_{\psi}$, and obtain $m_{\psi}=1.18_{-0.24}^{+0.28}\times10^{-22}{\,\rm eV}$ and $m_{\psi}=1.79_{-0.33}^{+0.35}\times10^{-22}{\,\rm eV}~(2\sigma)$ using the observational data sets of Walker et al. (2007) and Walker et al. (2009b), respectively. We show that the estimate of $m_{\psi}$ is sensitive to the dSphs kinematic data sets and is robust to various models of stellar density profile. We also consider multiple stellar subpopulations in dSphs and find consistent results. This mass range of $m_{\psi}$ is in good agreement with other independent estimates, such as the high-redshift luminosity functions, the reionization history, and the Thomson optical depth to the cosmic microwave background., Comment: 12 pages, 11 figures, 2 tables, add appendix B about the joint analysis, submitted to MNRAS

We explore the possibility of the formation of globular clusters under ultraviolet (UV) background radiation. One-dimensional spherical symmetric radiation hydrodynamics (RHD) simulations by Hasegawa et al. have demonstrated that the collapse of low-mass (10^6-10^7 solar masses) gas clouds exposed to intense UV radiation can lead to the formation of compact star clusters like globular clusters (GCs) if gas clouds contract with supersonic infall velocities. However, three-dimensional effects, such as the anisotropy of background radiation and the inhomogeneity in gas clouds, have not been studied so far. In this paper, we perform three-dimensional RHD simulations in a semi-cosmological context, and reconsider the formation of compact star clusters in strong UV radiation fields. As a result, we find that although anisotropic radiation fields bring an elongated shadow of neutral gas, almost spherical compact star clusters can be procreated from a "supersonic infall" cloud, since photo-dissociating radiation suppresses the formation of hydrogen molecules in the shadowed regions and the regions are compressed by UV heated ambient gas. The properties of resultant star clusters match those of GCs. On the other hand, in weak UV radiation fields, dark matter-dominated star clusters with low stellar density form due to the self-shielding effect as well as the positive feedback by ionizing photons. Thus, we conclude that the "supersonic infall" under a strong UV background is a potential mechanism to form GCs., 16 pages, 12 figures, 3 tables, accepted for publication in MNRAS

We investigate the morphology of the colour-magnitude diagram (CMD) of the open cluster NGC 2509 in comparison with other Galactic open clusters of similar age using Gaia photometry. At $\sim900$ Myr Galactic open clusters in our sample all show an extended main sequence turn off (eMSTO) with the exception of NGC 2509, which presents an exceptionally narrow CMD. Our analysis of the Gaia data rules out differential extinction, stellar density, and binaries as a cause for the singular MSTO morphology in this cluster. We interpret this feature as a consequence of the stellar rotation distribution within the cluster and present the analysis with MIST stellar evolution models that include the effect of stellar rotation on which we based our conclusion. In particular, these models point to an unusually narrow range of stellar rotation rates ($\Omega/\Omega_{\rm{crit,ZAMS}} = [0.4, 0.6]$) within the cluster as the cause of this singular feature in the CMD of NGC 2509. Interestingly, models that do not include rotation are not as good at reproducing the morphology of the observed CMD in this cluster., Comment: 9 pages, 7 figures

Planet-forming circumstellar disks are a fundamental part of the star formation process. Since stars form in a hierarchical fashion in groups of up to hundreds or thousands, the UV radiation environment that these disks are exposed to can vary in strength by at least six orders of magnitude. This radiation can limit the masses and sizes of the disks. Diversity in star forming environments can have long lasting effects in disk evolution and in the resulting planetary populations. We perform simulations to explore the evolution of circumstellar disks in young star clusters. We include viscous evolution, as well as the impact of dynamical encounters and external photoevaporation. We find that photoevaporation is an important process in destroying circumstellar disks: in regions of stellar density $\rho \sim 100 \mathrm{\ M}_\odot \mathrm{\ pc}^{-3}\mathrm{\ }$ around 80% of disks are destroyed before 2 Myr of cluster evolution. Our findings are in agreement with observed disk fractions in young star forming regions and support previous estimations that planet formation must start in timescales < 0.1 - 1 Myr., Comment: 13 pages, 14 figures. Accepted for publication in MNRAS

We apply a flexible parametric model, a combination of generalized Plummer profiles, to infer the shapes of the stellar density profiles of the Milky Way's satellite dwarf spheroidal galaxies (dSphs). We apply this model to 40 dSphs using star counts from the Sloan Digital Sky Survey, PanStarrs-1 Survey, Dark Energy Survey, and Dark Energy Camera Legacy Survey. Using mock data, we examine systematic errors associated with modelling assumptions and identify conditions under which our model can identify `non-standard' stellar density profiles that have central cusps and/or steepened outer slopes. Applying our model to real dwarf spheroidals, we do not find evidence for centrally cusped density profiles among the fifteen Milky Way satellites for which our tests with mock data indicate there would be sufficient detectability. We do detect steepened (with respect to a standard Plummer model) outer profiles in several dSphs--Fornax, Leo I, Leo II, and Reticulum II--which may point to distinct evolutionary pathways for these objects. However, the outer slope of the stellar density profile does not yet obviously correlate with other observed galaxy properties., 25 pages, 17 figures. Published in ApJ

Author(s): Titus, R; Ney, EM; Zegers, RGT; Bazin, D; Belarge, J; Bender, PC; Brown, BA; Campbell, CM; Elman, B; Engel, J; Gade, A; Gao, B; Kwan, E; Lipschutz, S; Longfellow, B; Lunderberg, E; Mijatovic, T; Noji, S; Pereira, J; Schmitt, J; Sullivan, C; Weisshaar, D; Zamora, JC | Abstract: Background: In the late stages of stellar core collapse just prior to core bounce, electron captures on medium-heavy nuclei drive deleptonization. Therefore, simulations require the use of accurate reaction rates. Nuclei with neutron number near N=50 above atomic number Z=28 play an important role. Rates presently used in astrophysical simulations rely primarily on a relatively simple single-state approximation. In order to improve the accuracy of the astrophysical simulations, experimental data are needed to test the electron-capture rates and to guide the development of better theoretical models and astrophysical simulations. Purpose: The purpose of the present work was to measure the Gamow-Teller transition strength from Kr86 to Br86, to derive the stellar electron-capture rates based on the extracted strengths, and to compare the derived rates with rates based on shell-model and quasiparticle random-phase approximation (QRPA) Gamow-Teller strengths calculations, as well as the single-state approximation. An additional purpose was to test the impact of using improved electron-capture rates on the late evolution of core-collapse supernovae. Method: The Gamow-Teller strengths from Kr86 were extracted from the Kr86(t,He3+γ) charge-exchange reaction at 115MeV/u. The electron-capture rates were calculated as a function of stellar density and temperature. Besides the case of Kr86, the electron-capture rates based on the QRPA calculations were calculated for 78 additional isotopes near N=50 above Z=28. The impact of using these rates instead of those based on the single-state approximation is studied in a spherically symmetrical simulation of core collapse just prior to bounce. Results: The derived electron-capture rates on Kr86 from the experimental Gamow-Teller strength distribution are much smaller than the rates estimated based on the single-state approximation. Rates based on Gamow-Teller strengths estimated in shell-model and QRPA calculations are more accurate. The core-collapse supernova simulation with electron-capture rates based on the QRPA calculations indicate a significant reduction in the deleptonization during the collapse phase. Conclusions: It is important to utilize microscopic theoretical models that are tested by experimental data to constrain and estimate Gamow-Teller strengths and derived electron-capture rates for nuclei near N=50 that are inputs for astrophysical simulations of core-collapse supernovae and their multimessenger signals, such as the emission of neutrinos and gravitational waves.

We present the results of a binary population study in the Orion Nebula Cluster (ONC) using archival Hubble Space Telescope (HST) data obtained with the Advanced Camera for Surveys (ACS) in Johnson V filter (HST Proposal 10246, PI M. Robberto). Young clusters and associations hold clues to the origin and properties of multiple star systems. Binaries with separations $< 100 $ AU are useful as tracers of the initial binary population since they are not as likely to be destroyed through dynamical interactions. Low mass, low stellar density star-forming regions such as Taurus-Auriga, reveal an excess of multiples compared to the Galactic Field. Studying the binary population of higher mass, higher stellar density star-forming regions like the ONC provides useful information concerning the origin of the Galactic Field star population. In this survey, we characterize the previously unexplored (and incomplete) separation parameter space of binaries in the ONC (15 - 160 AU) by fitting a double-PSF model built from empirical PSFs. We identified 14 candidate binaries (11 new detections) and find that 8$_{-2\%}^{+4\%}$ of our observed sample are in binary systems, complete over mass ratios and separations of 0.6 $< $ q $< $ 1.0 and 30 $< $ a $< $ 160 AU. This is consistent with the Galactic Field M-dwarf population over the same parameter ranges, 6.5\% $\pm$ 3\%. Therefore, high mass star forming regions like the ONC would not require further dynamical evolution for their binary population to resemble the Galactic Field, as some models have hypothesized for young clusters., 14 pages, 3 figures

The harvest of exoplanet discoveries has opened the area of exoplanet characterisation. But this cannot be achieved without a careful analysis of the host star parameters. The system of HD219134 hosts two transiting exoplanets and at least two additional non-transiting exoplanets. We used the VEGA/CHARA interferometer to measure the angular diameter of HD219134, leading to a stellar radius of $R_{\star}=0.726\pm0.014 R_{\odot}$. We also derived the stellar density from the transits light curves ($\rho_{\star}=1.82\pm0.19 \rho_{\odot}$), which finally gives a direct estimate of the mass ($M_{\star}=0.696\pm0.078 M_{\odot}$) with a correlation of 0.46 between $R_{\star}$ and $M_{\star}$. This new mass is smaller than that derived from the C2kSMO stellar evolutionary model, which provides a mass range of 0.755$-$0.810 ($\pm 0.040$) $M_{\odot}$. This allows us to infer the mass, radius and density of the two transiting exoplanets of the system. We then use an inference model to obtain the internal parameters of these two transiting exoplanets. Moreover, we find that planet $b$ and $c$ have smaller radii than previously estimated ($1.500\pm0.057$ and $1.415\pm0.049 R_{\oplus}$, respectively); this clearly puts these planets out of the gap in the exoplanetary radii distribution and validates their super-Earth nature. Planet $b$ is more massive than planet $c$, but possibly less dense. We investigate whether this could be caused by partial melting of the mantle and find that tidal heating due to non-zero eccentricity of planet $b$ may be powerful enough. The system of HD219134 constitutes a very valuable benchmark for both stellar physics and exoplanetary science. The direct determination of the stellar density, radius and mass should be more extensively applied to provide accurate exoplanets properties and calibrate stellar models., Comment: 12 pages, 6 figures, 7 tables; published in A$\&$A. (This version includes language editing corrections.)

Prompted by some inconsistencies in the photometry of the VISTA Variables in the V\'ia L\'actea (VVV) survey, we conduct a revision of the standard calibration procedure of VISTA data in the $J$, $H$, and $K_S$ passbands. Two independent sources of bias in the photometric zero-points are identified: First, high sky backgrounds severely affect the $H$-band measurements, but this can mostly be minimized by strict data vetting. Secondly, during the zero-point calibration, stars serving as standards are taken from the 2MASS catalog, which can suffer from high degrees of blending in regions of high stellar density, affecting both the absolute photometric calibration, as well as the scatter of repeated observations. The former affects studies that rely on an absolute magnitude scale, while the latter can also affect the shapes and amount of scatter in the VVV light curves, thus potentially hampering their proper classification. We show that these errors can be effectively eliminated by relatively simple modifications of the standard calibration procedure, and demonstrate the effect of the recalibration on the VVV survey's data quality. We give recommendations for future improvements of the pipeline calibration of VISTA photometry, while also providing preliminary corrections to the VVV $JHK_\mathrm{S}$ observations as a temporary measure., Comment: Accepted for publication in Experimental Astronomy; Table 1 is available from the first author upon request

We present a calibration of the Tip of the Red Giant Branch (TRGB) in the Large Magellanic Cloud (LMC) on the HST/ACS F814W system. We use archival HST observations to derive blending corrections and photometric transformations for two ground-based wide-area imaging surveys of the Magellanic Clouds. We show that these surveys are biased bright by up to ~0.1 mag in the optical due to blending, and that the bias is a function of local stellar density. We correct the LMC TRGB magnitudes from Jang & Lee (2017) and use the geometric distance from Pietrzynski et al. (2019) to obtain an absolute TRGB magnitude of M_F814W=-3.97+/-0.046 mag. Applying this calibration to the TRGB magnitudes from Freedman et al. (2019) in SN Ia hosts yields a value for the Hubble constant of H_0=72.4+/-2.0 km/s/Mpc for their TRGB+SNe Ia distance ladder. The difference in the TRGB calibration and the value of H_0 derived here and by Freedman et al. (2019) primarily results from their overestimate of the LMC extinction, caused by inconsistencies in their different sources of TRGB photometry for the Magellanic Clouds. Using the same source of photometry (OGLE) for both Clouds and applying the aforementioned corrections yields a value for the LMC I-band TRGB extinction that is lower by 0.06 mag, consistent with independent OGLE reddening maps used by us and by Jang & Lee (2017) to calibrate TRGB and determine H_0., Accepted for publication in ApJ. This version incorporates all changes introduced during the refereeing process

We study the evolution of star clusters in the Galactic tidal field starting from their birth in molecular clumps. Our model clusters form according to the local-density-driven cluster formation model in which the stellar density profile is steeper than that of gas. As a result, clusters resist the gas expulsion better than predicted by earlier models. We vary the impact of the Galactic tidal field {\lambda}, considering different Galactocentric distances (3-18 kpc), as well as different cluster sizes. Our model clusters survive the gas expulsion independent of {\lambda}. We investigated the relation between the cluster mass at the onset of secular evolution and their dissolution time. The model clusters formed with a high star-formation efficiency (SFE) follow a tight mass-dependent dissolution relation, in agreement with previous theoretical studies. However, the low-SFE models present a shallower mass-dependent relation than high-SFE clusters, and most dissolve before reaching 1 Gyr (cluster teenage mortality)., Comment: 4 pages, 2 figures. Submitted for publication in the Proceedings of IAU Symposium 351

We study the effect of dynamical friction on globular clusters and on the stars evaporated from the globular clusters (stellar streams) moving in a galactic halo. Due to dynamical friction, the position of a globular cluster (GC) as a stream progenitor starts to shift with respect to its original position in the reference frame of initial GC orbit. Therefore the stars that have evaporated at different times have different mean position with respect to the GC position. This shifting results in a certain asymmetry in stellar density distribution between the leading and trailing arms of the stream. The degree of the asymmetry depends on the characteristics of the environment in which the GC and the stream stars move. As GCs are located mainly in outer parts of a galaxy, this makes dynamical friction a unique probe to constrain the underlying dark matter spatial density and velocity distributions. For a GC NGC 3201 we compared our theoretical shift estimates with available observations. Due to large uncertainties in current observation data, we can only conclude that the derived estimates have the same order of magnitude., 11 pages, 5 figures, accepted to MNRAS

Establishing precisely how stars and interstellar medium distribute within the central 100 pc area around an AGN, down to the pc scale, is key for understanding how the very final transfer of matter from kpc scale to the sub-parsec size of the accretion disc is achieved. Using AO-assisted (SPHERE-VLT) near-IR images in H and Ks and narrow-band of the Seyfert 2 galaxy NGC1068 we analyse the radial distribution of brightness in the central r < 100 pc area down to the pc scale. The median-averaged radial profiles are adjusted by a cusp (power-law) plus a central point-source. A simple radiative transfer model is used to interpret the data. We find that the fit of profiles beyond 10pc is done quite precisely at Ks by a cusp of exponent -2.0 plus a central point-source and by a cusp of exponent -1.2 at H. The difference between H and Ks can be explained by differential extinction, provided that the distribution of dust is itself cuspy as r^-1. However, the required stellar density follows a r^-4 cusp, much steeper than any other cusp theoretically predicted and the mass of the cluster is unreasonable, even introducing a segregation in the stellar population with an excess of giant stars inward. A much more acceptable solution is found with a K profile dominated by warm dust emission, while the H profile corresponds to a stellar cusp. NGC1068 is shown to satisfy a relationship between half-light radius, cusp luminosity and exponent which suggests that the cusp is the remnant of a recent starbust. We identify the central point-like source with the very hot dust at the internal wall of the putative torus and derive an intrinsic luminosity that requires an overall extinction AK ~ 8, a value consistent with predictions by several models., 20 pages. Revised version after referee's comments. Recommended for publication in A&A, after revision (this version)

We present results on the stellar density radial profile of the outer regions of NGC6779, a Milky Way globular cluster recently proposed as a candidate member of the Gaia Sausage structure, a merger remnant of a massive dwarf galaxy with the Milky Way. Taking advantage of the Pan-STARRS PS1 public astrometric and photometric catalogue, we built the radial profile for the outermost cluster regions using horizontal branch and main sequence stars, separately, in order to probe for different profile trends because of difference stellar masses. Owing to its relatively close location to the Galactic plane, we have carefully treated the chosen colour-magnitude regions properly correcting them by the amount of interstellar extinction measured along the line-of-side of each star, as well as cleaned them from the variable field star contamination observed across the cluster field. In the region spanning from the tidal to the Jacobi radii the resulting radial profiles show a diffuse extended halo, with an average power law slope of -1. While analysing the relationships between the Galactocentric distance, the half-mass density, the half-light radius, the slope of the radial profile of the outermost regions, the internal dynamical evolutionary stage, among others, we found that NGC6779 shows structural properties similar to those of the remaining Gaia Sausage candidate globular clusters, namely, they are massive clusters (>10^5Mo) in a moderately early dynamical evolutionary stage, with observed extra-tidal structures., Comment: 8 pages, 6 figures. Accepted for publication in Monthly Noticies of the Royal Astronomical Society

The Galactic center (GC) is the densest region of the Milky Way. Variability surveys towards the GC potentially provide the largest number of variable stars per square degree within the Galaxy. However, high stellar density is also a drawback due to blending. Moreover, the GC is affected by extreme reddening, therefore near infrared observations are needed. We plan to detect new variable stars towards the GC, focusing on type II Cepheids (T2Cs) which have the advantage of being brighter than RR Lyrae stars. We perform parallel Lomb-Scargle and Generalized Lomb-Scargle periodogram analysis of the $K_s$-band time series of the VISTA variables in the Via Lactea survey, to detect periodicities. We employ statistical parameters to clean our sample. We take account of periods, light amplitudes, distances, and proper motions to provide a classification of the candidate variables. We detected 1,019 periodic variable stars, of which 164 are T2Cs, 210 are Miras and 3 are classical Cepheids. We also found the first anomalous Cepheid in this region. We compare their photometric properties with overlapping catalogs and discuss their properties on the color-magnitude and Bailey diagrams. We present the most extensive catalog of T2Cs in the GC region to date. Offsets in E($J-K_s$) and in the reddening law cause very large ($\sim$1-2 kpc) uncertainties on distances in this region. We provide a catalog which will be the starting point for future spectroscopic surveys in the innermost regions of the Galaxy., A&A, accepted

Understanding how the birthplace of stars affects planet-forming discs is important for a comprehensive theory of planet formation. Most stars are born in dense star-forming regions where the external influence of other stars, particularly the most massive stars, will affect the survival and enrichment of their planet-forming discs. Simulations suggest that stellar dynamics play a central role in regulating how external feedback affects discs, but comparing models to observations requires an estimate of the initial stellar density in star-forming regions. Structural analyses constrain the amount of dynamical evolution a star-forming region has experienced; regions that maintain substructure and do not show mass segregation are likely dynamically young, and therefore close to their birth density. In this paper, we present a structural analysis of two clusters in the Carina Nebula, Tr14 and Tr16. We show that neither cluster shows evidence for mass segregation or a centrally concentrated morphology, suggesting that both regions are dynamically young. This allows us to compare to simulations from Nicholson et al. (2019) who predict disc survival rates in star-forming regions of different initial densities. The surviving disc fractions in Tr14 and Tr16 are consistent with their predictions (both are $\sim 10$%), supporting a growing body of evidence that the star-forming environment plays an important role in the survival and enrichment of protoplanetary discs., MNRAS, accepted

The aim of this work is to investigate the characteristics of the young stellar population in the spiral galaxy NGC 247. In particular, we focused our attention in its hierarchical clustering distributions and the properties of the smallest groups. We used multiband Hubble Space Telescope (HST) data from three fields covering more than half of NGC 247 to select the young population. Then, through the path linkage criterion (PLC), we found compact young star groups, and estimated their fundamental parameters, such as their stellar densities, sizes, amount of members and luminosity function (LF) slopes. We also performed a fractal analysis to determinate the clustering properties of this population. We build a stellar density map and the corresponding dendrograms corresponding to the galactic young population to detect large structures and draw their main characteristics. We detected 339 young star groups, for which we computed a mean radius of $\sim$ 60 pc and a maximum in the size distribution between 30 and 70 pc. We also obtained LF slopes with a bimodal distribution showing peaks at $\sim$ 0.1 and $\sim$ 0.2. We identified several candidates to HII regions which follow an excellent spatial correlation with the young groups found by the PLC. We observed that the young population are hierarchically organized, where the smaller and denser structures are within larger and less dense ones. We noticed that all these groups presented a fractal subclustering, following the hierarchical distribution observed in the corresponding stellar density map. For the large young structures observed in this map, we obtained a fractal dimension of $\sim$ 1.6-1.8 using the perimeter-area relation and the cumulative size distribution. These values are consistent with a scenario of hierarchical star formation., Accepted for publication in A&A

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs., Accepted for publication in MNRAS. 22 pages (+22 pages of appendix), 10 (+21 in appendix) figures, 5 (+1 in appendix) tables

In the local universe, black holes of $10^{5-6}$ M$_{\odot}$ are hosted in galaxies displaying a variety of stellar profiles and morphologies. These black holes are the anticipated targets of LISA, the Laser Interferometer Space Antenna that will detect the low-frequency gravitational-wave signal emitted by binary black holes in this mass interval. In this paper, we infer upper limits on the lifetime of binary black holes of $10^{5-6}$ M$_{\odot}$ and up to $10^8$ M$_{\odot}$, forming in galaxy mergers, exploring two underlying stellar density profiles, by Dehnen and by Prugniel & Simien, and by exploiting local scaling relations between the mass of the black holes and several quantities of their hosts. We focus on the phase of the dynamical evolution when the binary is transitioning from the hardening phase ruled by the interaction with single stars to the phase driven by the emission of gravitational waves. We find that different stellar profiles predict very distinct trends with binary mass, with lifetimes ranging between fractions of a Gyr to more than 10 Gyr, and with a spread of about one order of magnitude, given by the uncertainties in the observed correlations, which are larger in the low-mass tail of the observed black hole population., Comment: 11 pages, 6 figures; published on MNRAS

Recent observational programmes are providing a global view of the Milky Way bulge that serves as template for detailed comparison with models and extragalactic bulges. A number of surveys (i.e. VVV, GIBS, GES, ARGOS, BRAVA, APOGEE) are producing comprehensive and detailed extinction, metallicity, kinematics and stellar density maps of the Galactic bulge with unprecedented accuracy. However, the still missing key ingredient is the distribution of stellar ages across the bulge. To overcome this limitation, we aim to age-date the stellar population in several bulge fields with the ultimate goal of deriving an age map of the Bulge. This paper presents the methodology and the first results obtained for a field along the Bulge minor axis, at $b=-6^\circ$. We use a new PSF-fitting photometry of the VISTA Variables in the V\'{i}a L\'{a}ctea (VVV) survey data to construct deep color-magnitude diagrams of the bulge stellar population down to $\sim$ 2 mag below the Main Sequence turnoff. We find the bulk of the bulge stellar population in the observed field along the minor axis to be at least older than $\sim$ 7.5 Gyr. In particular, when the metallicity distribution function spectroscopically derived by GIBS is used, the best fit to the data is obtained with a combination of synthetic populations with ages in between $\sim$ 7.5 Gyr and 11 Gyr. However, the fraction of stars younger than $\sim$ 10 Gyr strongly depends upon the number of Blue Straggler Stars present in the bulge. Simulations show that the observed color-magnitude diagram of the bulge in the field along the minor axis is incompatible with the presence of a conspicuous population of intermediate-age/young (i.e. $\lesssim 5$ Gyr) stars., Comment: 14 pages, 7 figues, Accepted for publication on A&A

Robust measurements of cosmological parameters from galaxy surveys rely on our understanding of systematic effects that impact the observed galaxy density field. In this paper we present, validate, and implement the idea of adopting the systematics mitigation method of Artificial Neural Networks for modeling the relationship between the target galaxy density field and various observational realities including but not limited to Galactic extinction, seeing, and stellar density. Our method by construction allows a wide class of models and alleviates over-training by performing k-fold cross-validation and dimensionality reduction via backward feature elimination. By permuting the choice of the training, validation, and test sets, we construct a selection mask for the entire footprint. We apply our method on the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Emission Line Galaxies (ELGs) selection from the Dark Energy Camera Legacy Survey (DECaLS) Data Release 7 and show that the spurious large-scale contamination due to imaging systematics can be significantly reduced by up-weighting the observed galaxy density using the selection mask from the neural network and that our method is more effective than the conventional linear and quadratic polynomial functions. We perform extensive analyses on simulated mock datasets with and without systematic effects. Our analyses indicate that our methodology is more robust to overfitting compared to the conventional methods. This method can be utilized in the catalog generation of future spectroscopic galaxy surveys such as eBOSS and Dark Energy Spectroscopic Instrument (DESI) to better mitigate observational systematics., Comment: 31 pages, 25 figures, accepted for publication in MNRAS. Moderate revision throughout the paper. The new version includes a quantitative evaluation of the remaining systematic effects in the DECaLS DR7 data in Summary and Discussion and in Conclusion. Pipeline available at https://github.com/mehdirezaie/SYSNet

We present a homogeneous 2MASS bright galaxy catalogue at low Galactic latitudes ($|b| \le 10.0^{\rm o}$, called Zone of Avoidance) which is complete to a Galactic extinction-corrected magnitude of $K^o_s \le 11.25^{\rm m}$. It also includes galaxies in regions of high foreground extinctions ($E(B-V) > 0.95^{\rm m}$) situated at higher latitudes. This catalogue forms the basis of studies of large-scale structures, flow fields and extinction across the ZoA and complements the ongoing 2MASS Redshift and Tully-Fisher surveys. It comprises 3763 galaxies, 70% of which have at least one radial velocity measurement in the literature. The catalogue is complete up to star density levels of $\log N_*/{\rm deg}^2, 29 pages, 16 figures; accepted for publication in MNRAS; full catalogue available as ancillary file in sourcefile package