Search

We present a map of the total intrinsic reddening across ~90 deg$^{2}$ of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222,752 early-type galaxies based on the LEPHARE $\chi^{2}$ minimisation SED-fitting routine. We find excellent agreement between the regions of enhanced intrinsic reddening across the central (4x4 deg$^2$) region of the LMC and the morphology of the low-level pervasive dust emission as traced by far-IR emission. In addition, we are able to distinguish smaller, isolated enhancements that are coincident with known star-forming regions and the clustering of young stars observed in morphology maps. The level of reddening associated with the molecular ridge south of 30 Doradus is, however, smaller than in the literature reddening maps. The reduced number of galaxies detected in this region, due to high extinction and crowding, may bias our results towards lower reddening values. Our map is consistent with maps derived from red clump stars and from the analysis of the star formation history across the LMC. This study represents one of the first large-scale categorisations of extragalactic sources behind the LMC and as such we provide the LEPHARE outputs for our full sample of ~2.5 million sources., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

We present a kinematic analysis of the Small Magellanic Cloud using 3700 spectra extracted from the European Southern Observatory archive. We used data from Gaia and near-infrared photometry to select stellar populations and discard Galactic foreground stars. The sample includes main-sequence, red giant branch and red clump stars, observed with the Fibre Large Array Multi Wavelength Spectrograph. The spectra have a resolving power lambda/Delta(lambda) from 6500 to 38000. We derive radial velocities by employing a full spectrum fitting method using a penalised pixel fitting routine. We obtain a mean radial velocity for the galaxy of 159+/-2 km/s, with a velocity dispersion of 33+/-2 km/s. Our velocities agree with literature estimates for similar (young or old) stellar populations. The radial velocity of stars in the Wing and bar-like structure differ as a consequence of the dynamical interaction with the Large Magellanic Cloud. The higher radial velocity of young main-sequence stars in the bar compared to that of supergiants can be attributed to star formation around 40 Myr ago from gas already influenced by tidal stripping. Similarly, young main-sequence stars in the northern part of the bar, resulting from a prominent episode 25 Myr ago, have a higher radial velocity than stars in the southern part. Radial velocity differences between the northern and southern bar over densities are also traced by giant stars. They are corroborated by studies of the cold gas and proper motion indicating stretching/tidal stripping of the galaxy., Comment: 18 pages, 13 figures, accepted for publication in MNRAS

Multi-mode Cepheids pulsate simultaneously in more than one mode of oscillation. They provide an independent means to test stellar models and pulsation theories. They can also be used to derive metallicities. In recent years, the number of known multi-mode Cepheids has increased dramatically with the discovery of a large number of Galactic double-mode Cepheids. To date, 209 double-mode Cepheids have been detected in the Galactic bulge and disk, mostly based on the Optical Gravitational Lensing Experiment's (OGLE) catalog. In this paper, we conduct a comprehensive search for double-mode Cepheids in the northern sky based on Zwicky Transient Facility Data Release 5. We found 72 such objects in the Milky Way. The periods of the 30 sample objects already included in the OGLE catalog show excellent agreement with the OGLE periods. The period ratios of our new Cepheids are consistent with those of known double-mode Cepheids, as evidenced by their loci in the so-called `Petersen diagram'. Compared with OGLE, the completeness of our double-mode Cepheid sample is around 71\%. The much improved temporal sampling of the Zwicky Transient Facility offers significant scope to find more double-mode Cepheids, especially at the distribution's short-period end., 11 pages, 3 figures, 3 tables; accepted for publication in ApJS

We used data from the near-infrared VISTA survey of the Magellanic Cloud system (VMC) to measure proper motions (PMs) of stars within the Small Magellanic Cloud (SMC). The data analysed in this study comprise 26 VMC tiles, covering a total contiguous area on the sky of ~40 deg$^2$. Using multi-epoch observations in the Ks band over time baselines between 13 and 38 months, we calculated absolute PMs with respect to ~130,000 background galaxies. We selected a sample of ~2,160,000 likely SMC member stars to model the centre-of-mass motion of the galaxy. The results found for three different choices of the SMC centre are in good agreement with recent space-based measurements. Using the systemic motion of the SMC, we constructed spatially resolved residual PM maps and analysed for the first time the internal kinematics of the intermediate-age/old and young stellar populations separately. We found outward motions that point either towards a stretching of the galaxy or stripping of its outer regions. Stellar motions towards the North might be related to the "Counter Bridge" behind the SMC. The young populations show larger PMs in the region of the SMC Wing, towards the young Magellanic Bridge. In the older populations, we further detected a coordinated motion of stars away from the SMC in the direction of the Old Bridge as well as a stream towards the SMC., 20 pages, 20 figures, accepted for publication in MNRAS

Multiple stellar populations (MPs) representing star-to-star light-element abundance variations are common in nearly all ancient Galactic globular clusters. Here we provide the strongest evidence yet that the populous, ~ 1.7 Gyr-old Large Magellanic Cloud cluster NGC 2173 also exhibits light-element abundance variations. Thus, our results suggest that NGC 2173 is the youngest cluster for which MPs have been confirmed to date. Our conclusion is based on the distinct bifurcation at the tip of its red-giant branch in high-quality color--magnitude diagrams generated from Hubble Space Telescope imaging observations. Our results are further supported by a detailed analysis of 'pseudo-$UBI$' maps, which reveal clear evidence of a bimodality in the cluster's red-giant-branch color distribution. Young clusters in the Magellanic Clouds can provide critical insights into galaxy evolution histories. Our discovery of MPs in NGC 2173 suggests that ancient Galactic globular clusters and young massive clusters might share a common formation process., 8 pages, 5 figures

We perform a statistical clustering analysis of upper main-sequence stars in the Large Magellanic Cloud (LMC) using data from the Visible and Infrared Survey Telescope for Astronomy survey of the Magellanic Clouds. We map over 2500 young stellar structures at 15 significance levels across ~120 square degrees centred on the LMC. The structures have sizes ranging from a few parsecs to over 1 kpc. We find that the young structures follow power-law size and mass distributions. From the perimeter-area relation, we derive a perimeter-area dimension of 1.44+-0.20. From the mass-size relation and the size distribution, we derive two-dimensional fractal dimensions of 1.50+-0.10 and 1.61+-0.20, respectively. We find that the surface density distribution is well-represented by a lognormal distribution. We apply the Larson relation to estimate the velocity dispersions and crossing times of these structures. Our results indicate that the fractal nature of the young stellar structures has been inherited from the gas clouds from which they form and that this architecture is generated by supersonic turbulence. Our results also suggest that star formation in the LMC is scale-free from 10 pc to 700 pc., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

We present an extensive catalogue of BY Draconis (BY Dra)-type variables and their stellar parameters. BY Dra are main-sequence FGKM-type stars. They exhibit inhomogeneous starspots and bright faculae in their photospheres. These features are caused by stellar magnetic fields, which are carried along with the stellar disc through rotation and which produce gradual modulations in their light curves (LCs). Our main objective is to characterise the properties of BY Dra variables over a wide range of stellar masses, temperatures and rotation periods. A recent study categorised 84,697 BY Dra variables from Data Release 2 of the Zwicky Transient Facility based on their LCs. We have collected additional photometric data from multiple surveys and performed broad-band spectral energy distribution fits to estimate stellar parameters. We found that more than half of our sample objects are of K spectral type, covering an extensive range of stellar parameters in the low-mass regime (0.1-1.3 M$_{\odot}$ ). Compared with previous studies, most of the sources in our catalogue are rapid rotators, and so most of them must be young stars for which a spin-down has not yet occurred. We subdivided our catalogue based on convection zone depth and found that the photospheric activity index, $S_{\rm ph}$, is lower for higher effective temperatures, i.e., for thinner convective envelopes. We observe a broad range of photospheric magnetic activity for different spectral classes owing to the presence of stellar populations of different ages. We found a higher magnetically active fraction for K- than M-type stars., Comment: 12 pages, 13 figures

We derive the spatially-resolved star formation history (SFH) for a $96$ deg$^2$ area across the main body of the Large Magellanic Cloud (LMC), using the near-infrared photometry from the VISTA survey of the Magellanic Clouds (VMC). The data and analyses are characterised by a great degree of homogeneity and a low sensitivity to the interstellar extinction. 756 subregions of size $0.125$ deg$^2$ -- corresponding to projected sizes of about $296\times322\,\mathrm{pc}^{2}$ in the LMC -- are analysed. The resulting SFH maps, with typical resolution of $0.2$--$0.3$ dex in logarithm of age, reveal main features in the LMC disc at different ages: the patchy star formation at recent ages, the concentration of star formation on three spiral arms and on the Bar up to ages of $\sim\!1.6$ Gyr, and the wider and smoother distribution of older populations. The period of most intense star formation occurred roughly between 4 and 0.5 Gyr ago, at rates of $\sim\!0.3\,\mathrm{M}_{\odot}\mathrm{yr}^{-1}$. We compare young and old star formation rates with the observed numbers of RR Lyrae and Cepheids. We also derive a mean extinction and mean distance for every subregion, and the plane that best describes the spatial distribution of the mean distances. Our results cover an area about 50 per cent larger than the classical SFH maps derived from optical data by Harris & Zaritsky (2009). Main differences with respect to those maps are lower star formation rates at young ages, and a main peak of star formation being identified at ages slightly younger than $1$ Gyr., Accepted for publication in MNRAS

We have collected 2330 Cepheids to establish an intuitive 3D map of the Milky Way’s disk. As regards the warp amplitude, the Cepheid disk agrees well with the gas disk for radii up to 15 kpc. However, the mean line of nodes (LON) of the Cepheid disk deviates from the Galactic Center–Sun direction by 17.5±1.0°. This is a new and different result compared with previous results. The LON is not stable at any given radius, but it twists. The twisted pattern suggests that the formation of the Milky Way’s warp is dominated by the massive inner disk. The kinematic warp defined by the Cepheids is also in concordance with the spatial warp. In the 2020 era, the anticipated increasing number of new Cepheids will provide a key opportunity to view our Milky Way’s disk as a whole, and we expect that our knowledge of the disk’s main structural features will be much improved.

We present a detailed analysis of the projected stellar rotational velocities of the well-separated double main sequence (MS) in the young, $\sim200$Myr-old Milky Way open cluster NGC 2287 and suggest that stellar rotation may drive the split MSs in NGC 2287. We find that the observed distribution of projected stellar rotation velocities could result from a dichotomous distribution of stellar rotation rates. We discuss whether our observations may reflect the effects of tidal locking affecting a fraction of the cluster's member stars in stellar binary systems. The slow rotators are likely stars that initially rotated rapidly but subsequently slowed down through tidal locking induced by low-mass-ratio binary systems. However, the cluster may have a much larger population of short-period binaries than is usually seen in the literature, with relatively low secondary masses., Comment: 4 pages, 1 figure. To appear in the Proceedings of the 351 IAU Symposium "Star Clusters: From the Milky Way to the Early Universe"

The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, e.g. with the SMC about 250 Myr ago. {We calculate proper motions using multi-epoch $K_\mathrm{s}$-band images from the VISTA survey of the Magellanic Clouds system (VMC). Observations span a time baseline of 2$-$5 yr. We combine the VMC data with data from the Gaia early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south east region of the LMC shows a slow rotational speed compared to the overall rotation. $N$-body simulations suggest that this could be caused by a fraction of stripped SMC stars, located in that particular region, that move opposite to the expected rotation., Comment: 21 pages, 21 figures, accepted for publication in A&A

We study the morphology of the stellar periphery of the Magellanic Clouds in search of substructure using near-infrared imaging data from the VISTA Hemisphere Survey (VHS). Based on the selection of different stellar populations using the ($J-K_\mathrm{s}$, $K_\mathrm{s}$) colour-magnitude diagram, we confirm the presence of substructures related to the interaction history of the Clouds and find new substructures on the easter side of the LMC disc which may be owing to the influence of the Milky Way, and on the northern side of the SMC, which is probably associated to the ellipsoidal structure of the galaxy. We also study the luminosity function of red clump stars in the SMC and confirm the presence of a bi-modal distance distribution, in the form of a foreground population. We find that this bi-modality is still detectable in the eastern regions of the galaxy out to a 10 deg distance from its centre. Additionally, a background structure is detected in the North between 7 and 10 deg from the centre which might belong to the Counter Bridge, and a foreground structure is detected in the South between 6 and 8 deg from the centre which might be linked to the Old Bridge., 20 pages, 19 figures, accepted for publication in MNRAS

We independently determine the zero-point offset of the Gaia early Data Release-3 (EDR3) parallaxes based on $\sim 110,000$ W Ursae Majoris (EW)-type eclipsing binary systems. EWs cover almost the entire sky and are characterized by a relatively complete coverage in magnitude and color. They are an excellent proxy for Galactic main-sequence stars. We derive a $W1$-band Period-Luminosity relation with a distance accuracy of $7.4\%$, which we use to anchor the Gaia parallax zero-point. The final, global parallax offsets are $-28.6\pm0.6$ $\mu$as and $-25.4\pm4.0$ $\mu$as (before correction) and $4.2\pm0.5$ $\mu$as and $4.6\pm3.7$ $\mu$as (after correction) for the five- and six-parameter solutions, respectively. The total systematic uncertainty is $1.8$ $\mu$as. The spatial distribution of the parallax offsets shows that the bias in the corrected Gaia EDR3 parallaxes is less than 10 $\mu$as across $40\%$ of the sky. Only $15\%$ of the sky is characterized by a parallax offset greater than 30 $\mu$as. Thus, we have provided independent evidence that the parallax zero-point correction provided by the Gaia team significantly reduces the prevailing bias. Combined with literature data, we find that the overall Gaia EDR3 parallax offsets for Galactic stars are $[-20, -30]$ $\mu$as and 4-10 $\mu$as, respectively, before and after correction. For specific regions, an additional deviation of about 10 $\mu$as is found., Comment: 7 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

The color-magnitude diagrams (CMDs) of young star clusters show that, particularly at ultraviolet wavelengths, their upper main sequences (MSs) bifurcate into a sequence comprising the bulk population and a blue periphery. The spatial distribution of stars is crucial to understand the reasons for these distinct stellar populations. This study uses high-resolution photometric data obtained with the Hubble Space Telescope to study the spatial distributions of the stellar populations in seven Magellanic Cloud star clusters. The cumulative radial number fractions of blue stars within four clusters are strongly anti-correlated with those of the high-mass-ratio binaries in the bifurcated region, with negative Pearson coefficients < -0.7. Those clusters generally are young or in an early dynamical evolutionary stage. In addition, a supporting N-body simulation suggests the increasing percentage of blue-MS stars from the cluster centers to their outskirts may be associated with the dissolution of soft binaries. This study provides a different perspective to explore the MS bimodalities in young clusters and adds extra puzzles. A more comprehensive study combined with detailed simulations is needed in the future., 19 pages, 9 figures, 3 tables. Accepted for publication in ApJ

The impact of stellar rotation on the morphology of star cluster colour-magnitude diagrams is widely acknowledged. However, the physics driving the distribution of the equatorial rotation velocities of main-sequence turn-off (MSTO) stars is as yet poorly understood. Using Gaia Data Release 2 photometry and new Southern African Large Telescope medium-resolution spectroscopy, we analyse the intermediate-age ($\sim1\,$Gyr-old) Galactic open clusters NGC 3960, NGC 6134 and IC 4756 and develop a novel method to derive their stellar rotation distributions based on SYCLIST stellar rotation models. Combined with literature data for the open clusters NGC 5822 and NGC 2818, we find a tight correlation between the number ratio of slow rotators and the clusters' binary fractions. The blue-main-sequence stars in at least two of our clusters are more centrally concentrated than their red-main-sequence counterparts. The origin of the equatorial stellar rotation distribution and its evolution remains as yet unidentified. However, the observed correlation in our open cluster sample suggests a binary-driven formation mechanism., 9 pages, 7 figures, 1 table. Accepted for publication in MNRAS

We have derived high-spatial-resolution metallicity maps covering $\sim$105~deg$^2$ across the Large Magellanic Cloud (LMC) using near-infrared passbands from the VISTA Survey of the Magellanic Clouds. We attempt to understand the metallicity distribution and gradients of the LMC up to a radius of $\sim$ 6~kpc. We identify red giant branch (RGB) stars in spatially distinct $Y, (Y-K_{\rm s})$ colour-magnitude diagrams. In any of our selected subregions, the RGB slope is used as an indicator of the average metallicity, based on calibration to metallicity using spectroscopic data. The mean LMC metallicity is [Fe/H] = $-$0.42~dex ($\sigma$[Fe/H] = 0.04~dex). We find the bar to be mildly metal-rich compared with the outer disc, showing evidence of a shallow gradient in metallicity ($-0.008 \pm 0.001$ dex kpc$^{-1}$) from the galaxy's centre to a radius of 6~kpc. Our results suggest that the LMC's stellar bar is chemically similar to the bars found in large spiral galaxies. The LMC's radial metallicity gradient is asymmetric. It is metal-poor and flatter towards the southwest, in the direction of the Bridge. This hints at mixing and/or distortion of the spatial metallicity distribution, presumably caused by tidal interactions between the Magellanic Clouds., Comment: 13 pages, 15 figures, accepted for publication in the MNRAS

Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed., Comment: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

We derive stellar parameters and abundances (`stellar labels') of 40,034 late-B and A-type main-sequence stars extracted from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope Medium Resolution Survey (LAMOST--MRS). The primary selection of our early-type sample was obtained from LAMOST Data Release 7 based on spectral line indices. We employed the Stellar LAbel Machine (SLAM) to derive their spectroscopic stellar parameters, drawing on Kurucz spectral synthesis models with 6000 K $< T_\mathrm{eff}, Comment: 18 pages, 10 figures, 1 table. Accepted for publication in ApJS

Angular momentum is a key property regulating star formation and evolution. However, the physics driving the distribution of the stellar rotation rates of early-type main-sequence stars is as yet poorly understood. Using our catalog of 40,034 early-type stars with homogeneous $v\sin i$ parameters, we review the statistical properties of their stellar rotation rates. We discuss the importance of possible contaminants, including binaries and chemically peculiar stars. Upon correction for projection effects and rectification of the error distribution, we derive the distributions of our sample's equatorial rotation velocities, which show a clear dependence on stellar mass. Stars with masses less than $2.5\ {M_\odot}$ exhibit a unimodal distribution, with the peak velocity ratio increasing as stellar mass increases. A bimodal rotation distribution, composed of two branches of slowly and rapidly rotating stars, emerges for more massive stars ($M>2.5\ {M_\odot}$). For stars more massive than $3.0\ {M_\odot}$, the gap between the bifurcated branches becomes prominent. For the first time, we find that metal-poor ([M/H] $< -0.2$ dex) stars only exhibit a single branch of slow rotators, while metal-rich ([M/H] $> 0.2$ dex) stars clearly show two branches. The difference could be attributed to unexpectedly high spin-down rates and/or in part strong magnetic fields in the metal-poor subsample., 18 pages, 9 figures. Accepted for publication in ApJ; revised references for section 3.2

The conventional picture of coeval, chemically homogeneous, populous star clusters -- known as `simple stellar populations' (SSPs) -- is a view of the past. Photometric and spectroscopic studies reveal that almost all ancient globular clusters in the Milky Way and our neighbouring galaxies exhibit star-to-star light-element abundance variations, typically known as 'multiple populations' (MPs). Here, we analyse photometric $\it Hubble$ $\it Space$ $\it Telescope$ observations of three young ($, Comment: 11 pages, 5 Figures

The eastern region of the Small Magellanic Cloud (SMC) is found to have a foreground stellar sub-structure, which is identified as a distance bimodality (12 kpc apart) in the previous studies using Red Clump (RC) stars. Interestingly, studies of Red giant branch (RGB) stars in the eastern SMC indicate a bimodal radial velocity (RV) distribution. In this study, we investigate the connection between these two bimodal distributions to better understand the nature and origin of the foreground stellar sub-structure in the eastern SMC. We use the Gaia EDR3 astrometric data and archival RV data of RGB stars for this study. We found a bimodal RV distribution of RGB stars (separated by 35 - 45 km/s) in the eastern and south-western (SW) outer regions. The observed proper motion values of the lower and higher RV RGB components in the eastern regions are similar to those of the foreground and main-body RC stars respectively. This suggests that the two RGB populations in the eastern region are separated by a similar distance as those of the RC stars, and the RGB stars in the lower RV component are part of the foreground sub-structure. Based on the differences in the distance and RV of the two components, we estimated an approximate time of formation of this sub-structure as 307+/-65 Myr ago. This is comparable with the values predicted by simulations for the recent epoch of tidal interaction between the Magellanic Clouds. Comparison of the observed properties of RGB stars, in the outer SW region, with N-body simulations shows that the higher RV component in the SW region is at a farther distance than the main body, indicating the presence of a stellar Counter-Bridge in the SW region of the SMC., Comment: Accepted for publication in MNRAS, 11 pages, 6 figures

GALAH+ is a magnitude-limited survey of high resolution stellar spectra obtained by the HERMES spectrograph at the Australian Astronomical Observatory. Its third data release provides reduced spectra with new derivations of stellar parameters and abundances of 30 chemical elements for 584,015 dwarfs and giants, 88% of them in the Gaia magnitude range 11 < G < 14. Here we use these improved values of stellar parameters to build a library of observed spectra which is useful to study variations of individual spectral lines with stellar parameters. This and other improvements are used to derive radial velocities with uncertainties which are generally within 0.1 km/s or ~25% smaller than in the previous release. Median differences in radial velocities measured here and by the Gaia DR2 or APOGEE DR16 surveys are smaller than 30 m/s, a larger offset is present only for Gaia measurements of giant stars. We identify 4483 stars with intrinsically variable velocities and 225 stars for which the velocity stays constant over >=3 visits spanning more than a year. The combination of radial velocities from GALAH+ with distances and sky plane motions from Gaia enables studies of dynamics within streams and clusters. For example, we estimate that the open cluster M67 has a total mass of ~3300 Msun and its outer parts seem to be expanding, though astrometry with a larger time-span than currently available from Gaia eDR3 is needed to judge if the latter result is real., 15 pages, 4 tables, 17 figures + online data files, accepted for publication in MNRAS

Multiple stellar populations (MPs) with different chemical compositions are not exclusive features of old GCs (older than 10 Gyr). Indeed, recent studies reveal that younger clusters ($\sim$2--6 Gyr-old) in the Magellanic Clouds also exhibit star-to-star chemical variations among evolved stars. However, whether MPs are present among less evolved dwarfs of these intermediate-age clusters is still unclear. In this work, we search for chemical variations among GK-type dwarfs in the $\sim$2 Gyr-old cluster NGC 1978, which is the youngest cluster with MPs. We exploit deep ultraviolet and visual observations from the Hubble Space Telescope to constrain the nitrogen (N) and oxygen (O) variations among MS stars. To do this, we compare appropriate photometric diagrams that are sensitive to N and O with synthetic diagrams of simple stellar populations and MPs. We conclude that the G- and K-type MS stars in NGC\,1978 host MPs. Our statistical analysis shows that the fraction of N-rich stars ranges from $\sim$40\% to $\sim$80\%, depending on the detailed distributions of nitrogen and oxygen., 16 pages, 10 figures, ApJ accepted

Recent observational studies identified a foreground stellar sub-structure traced by red clump (RC) stars (~ 12 kpc in front of the main body) in the eastern regions of the Small Magellanic Cloud (SMC) and suggested that it formed during the formation of the Magellanic Bridge (MB), due to the tidal interaction of the Magellanic Clouds. Previous studies investigated this feature only up to 4.0 deg from the centre of the SMC due to the limited spatial coverage of the data and hence could not find a physical connection with the MB. To determine the spatial extent and properties of this foreground population, we analysed data from the Gaia data release 2 (DR2) of a ~ 314 sq. deg region centred on the SMC, which cover the entire SMC and a significant portion of the MB. We find that the foreground population is present only between 2.5 deg to ~ 5-6 deg from the centre of the SMC in the eastern regions, towards the MB and hence does not fully overlap with the MB in the plane of the sky. The foreground stellar population is found to be kinematically distinct from the stellar population of the main body with ~ 35 km/s slower tangential velocity and moving to the North-West relative to the main body. Though the observed properties are not fully consistent with the simulations, a comparison indicates that the foreground stellar structure is most likely a tidally stripped counterpart of the gaseous MB and might have formed from the inner disc (dominated by stars) of the SMC. A chemical and 3D kinematic study of the RC stars along with improved simulations, including both tidal and hydro-dynamical effects, are required to understand the offset between the foreground structure and MB., Accepted for publication in MNRAS

Following Columbus' voyages to the Americas, Castilian (Spanish) and Portuguese rulers engaged in heated geopolitical competition, which was eventually reconciled through a number of treaties that divided the world into two unequal hemispheres. However, the early-sixteenth-century papal demarcation line was poorly defined. Expressed in degrees with respect to a vague reference location, determination of longitude at sea became crucial in the nations' quest for expanding spheres of influence. In Spain, King Philip II and his son, Philip III, announced generous rewards for those whose solutions to the longitude problem performed well in sea trials and which were suitable for practical implementation. The potential reward generated significant interest from scientist-scholars and opportunists alike. The solutions proposed and the subset taken to sea provided important physical insights that still resonate today. None of the numerous approaches based on compass readings ("magnetic declination") passed the exacting sea trials, but the brightest sixteenth-century minds already anticipated that lunar distances and the use of marine timepieces would eventually enable more precise navigation. With most emphasis in the English-language literature focused on longitude solutions developed in Britain, France and the Low Countries, the earlier yet groundbreaking Spanish efforts have, undeservedly, largely been forgotten. Yet, they provided a firm basis for the development of an enormous "cottage industry" that lasted until the end of the eighteenth century., 32 pages. Journal of Astronomical History and Heritage, vol. 23, in press (December 2020 issue); author's version

The Magellanic Clouds are a nearby pair of interacting dwarf galaxies and satellites of the Milky Way. Studying their kinematic properties is essential to understanding their origin and dynamical evolution. They have prominent tidal features and the kinematics of these features can give hints about the formation of tidal dwarfs, galaxy merging and the stripping of gas. In addition they are an example of dwarf galaxies that are in the process of merging with a massive galaxy. The goal of this study is to investigate the kinematics of the Magellanic Bridge, a tidal feature connecting the Magellanic Clouds, using stellar proper motions to understand their most recent interaction. We calculated proper motions based on multi-epoch $K_{s}$-band aperture photometry, which were obtained with the Visible and Infrared Survey Telescope for Astronomy (VISTA), spanning a time of 1-3 yr, and we compared them with $Gaia$ Data Release 2 (DR2) proper motions. We tested two methods for removing Milky Way foreground stars using $Gaia$~DR2 parallaxes in combination with VISTA photometry or using distances based on Bayesian inference. We obtained proper motions for a total of 576,411 unique sources over an area of $23$ deg$^{2}$ covering the Magellanic Bridge including mainly Milky Way foreground stars, background galaxies, and a small population of possible Magellanic Bridge stars ($, Comment: 17 pages, 12 figures, accepted for publication in A&A

We present new near-infrared ($JHK_s$) time-series observations of RR Lyrae variables in the Messier 3 (NGC 5272) globular cluster using the WIRCam instrument at the 3.6-m Canada France Hawaii Telescope. Our observations cover a sky area of $\sim 21'\times 21'$ around the cluster center and provide an average of twenty epochs of homogeneous $JHK_s$-band photometry. New homogeneous photometry is used to estimate robust mean magnitudes for 175 fundamental-mode (RRab), 47 overtone-mode (RRc), and 11 mixed-mode (RRd) variables. Our sample of 233 RR Lyrae variables is the largest thus far obtained in a single cluster with time-resolved, multi-band near-infrared photometry. Near-infrared to optical amplitude ratios for RR Lyrae in Messier 3 exhibit a systematic increase moving from RRc to short-period ($P < 0.6$~days) and long-period ($P \gtrsim 0.6$~days) RRab variables. We derive $JHK_s$-band Period--Luminosity relations for RRab, RRc, and the combined sample of variables. Absolute calibrations based on the theoretically predicted Period--Luminosity--Metallicity relations for RR Lyrae stars yield a distance modulus, $��= 15.041 \pm 0.017~(\textrm{statistical}) \pm 0.036~(\textrm{systematic})$~mag, to Messier 3. When anchored to trigonometric parallaxes for nearby RR Lyrae stars from the {\it Hubble Space Telescope} and the {\it Gaia} mission, our distance estimates are consistent with those resulting from the theoretical calibrations, albeit with relatively larger systematic uncertainties., 22 pages, 14 figures, Accepted for publication in The Astronomical Journal

Longitude determination at sea gained increasing commercial importance in the late Middle Ages, spawned by a commensurate increase in long-distance merchant shipping activity. Prior to the successful development of an accurate marine timepiece in the late-eighteenth century, marine navigators relied predominantly on the Moon for their time and longitude determinations. Lunar eclipses had been used for relative position determinations since Antiquity, but their rare occurrences precludes their routine use as reliable way markers. Measuring lunar distances, using the projected positions on the sky of the Moon and bright reference objects--the Sun or one or more bright stars--became the method of choice. It gained in profile and importance through the British Board of Longitude's endorsement in 1765 of the establishment of a Nautical Almanac. Numerous 'projectors' jumped onto the bandwagon, leading to a proliferation of lunar ephemeris tables. Chronometers became both more affordable and more commonplace by the mid-nineteenth century, signaling the beginning of the end for the lunar distance method as a means to determine one's longitude at sea., Review article, 31 pages; Journal of Astronomical History and Heritage, in press (December 2020 issue); author's version

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars., Comment: accepted to MNRAS

Studies of young stellar objects (YSOs) in the Galaxy have found that a significant fraction exhibit photometric variability. However, no systematic investigation has been conducted on the variability of extragalactic YSOs. Here we present the first variability study of massive YSOs in a $\sim 1.5\,\mathrm{deg^2}$ region of the Large Magellanic Cloud (LMC). The aim is to investigate whether the different environmental conditions in the metal-poor LMC ($\sim$ 0.4-0.5 Z_sun) have an impact on the variability characteristics. Multi-epoch near-infrared (NIR) photometry was obtained from the VISTA Survey of the Magellanic Clouds (VMC) and our own monitoring campaign using the VISTA telescope. By applying a reduced $\chi^2$-analysis, stellar variability was identified. We found 3062 candidate variable stars from a population of 362 425 stars detected. Based on several Spitzer studies, we compiled a sample of high-reliability massive YSOs: a total of 173 massive YSOs have NIR counterparts ($K_{\mathrm{s}}\sim 18.5\,$mag) in the VMC catalogue, of which 39 display significant ($>3\sigma$) variability. They have been classified as eruptive, fader, dipper, short-term variable and long-period variable YSOs based mostly on the appearance of their $K_{\mathrm{s}}$ band light curves. The majority of YSOs are aperiodic, only five YSOs exhibit periodic lightcurves. The observed amplitudes are comparable or smaller than those for Galactic YSOs (only two Magellanic YSOs exhibit $\Delta K_{\mathrm{s}}>1\,$mag), not what would have been expected from the typically larger mass accretion rates observed in the Magellanic Clouds., Comment: Accepted for publication in MNRAS, 32 pages, 24 figures

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69., Comment: 21 pages, 10 figures, accepted by MNRAS for publication

We consider the published distance moduli to the Fornax and Coma galaxy clusters, with emphasis on the period since 1990. We have carefully homogenized our catalogs of distance moduli onto the distance scale established in the previous papers in this series. We assessed systematic differences associated with the use of specific tracers, and discarded results based on application of the Tully--Fisher relation and of globular cluster and planetary nebula luminosity functions. We recommend `best' weighted relative distance moduli for the Fornax and Coma clusters with respect to the Virgo cluster of $��(m-M)_0^{\rm Fornax - Virgo} = 0.18 \pm 0.28 $ mag and $��(m-M)_0^{\rm Coma - Virgo} = 3.75 \pm 0.23$ mag. The set of weighted mean distance moduli (distances) we derived as most representative of the clusters' distances is, \begin{eqnarray} (m-M)_0^{\rm Fornax} &=& 31.41 \pm 0.15 \mbox{ mag } (D = 19.1^{+1.4}_{-1.2} \mbox{ Mpc) and} \nonumber &=& 31.21 \pm 0.28 \mbox{ mag } (D = 17.5^{+2.4}_{-2.2} \mbox{ Mpc)}; \nonumber \\ (m-M)_0^{\rm Coma} &=& 34.99 \pm 0.38 \mbox{ mag } (D = 99.5^{+19.0}_{-15.9} \mbox{ Mpc) and} \nonumber &=& 34.78 \pm 0.27 \mbox{ mag } (D = 90.4^{+11.9}_{-10.6} \mbox{ Mpc)}, \nonumber \end{eqnarray} where the first value for each cluster is the result of our analysis of the direct distance moduli, while the second modulus is based on distance moduli relative to the Virgo cluster. The absolute and relative distance moduli for both clusters are mutually consistent within the uncertainties; the relative distance moduli yield shorter distances by $\sim$1$��$. Lingering uncertainties in the underlying absolute distance scale appear to have given rise to a systematic uncertainty on the order of 0.20 mag., 15 pages, 3 figures; accepted for publication in ApJS

We present the physical parameters of 2335 late-type contact binary (CB) systems extracted from the Catalina Sky Survey (CSS). Our sample was selected from the CSS Data Release 1 by strictly limiting the prevailing temperature uncertainties and light-curve fitting residuals, allowing us to almost eliminate any possible contaminants. We developed an automatic Wilson--Devinney-type code to derive the relative properties of CBs based on their light-curve morphology. By adopting the distances derived from CB (orbital) period--luminosity relations (PLRs), combined with the well-defined mass--luminosity relation for the systems' primary stars and assuming solar metallicity, we calculated the objects' masses, radii, and luminosities. Our sample of fully eclipsing CBs contains 1530 W-, 710 A-, and 95 B-type CBs. A comparison with literature data and with the results from different surveys confirms the accuracy and coherence of our measurements. The period distributions of the various CB subtypes are different, hinting at a possible evolutionary sequence. W-type CBs are clearly located in a strip in the total mass versus mass ratio plane, while A-type CBs may exhibit a slightly different dependence. There are no significant differences among the PLRs of A- and W-type CBs, but the PLR zero points are affected by their mass ratios and fill-out factors. Determination of zero-point differences for different types of CBs may help us improve the accuracy of the resulting PLRs. We demonstrate that automated approaches to deriving CB properties could be a powerful tool for application to the much larger CB samples expected to result from future surveys., Comment: 2o pages, 10 figures, 2 tables. Accepted for publication in ApJS

We present a map of the total intrinsic reddening across similar or equal to 34 deg(2) of the Small Magellanic Cloud (SMC) derived using optical (ugriz) and near-infrared (IR; YJK(s)) spectral energy distributions (SEDs) of background galaxies. The reddening map is created using a subsample of 29 274 galaxies with low levels of intrinsic reddening based on the LEPHARE chi(2) minimization SED-fitting routine. We find statistically significant enhanced levels of reddening associated with the main body of the SMC compared with regions in the outskirts [Delta E(B - V) similar or equal to 0.3 mag]. A comparison with literature reddening maps of the SMC shows that, after correcting for differences in the volume of the SMC sampled, there is good agreement between our results and maps created using young stars. In contrast, we find significant discrepancies between our results and maps created using old stars or based on longer wavelength far-IR dust emission that could stem from biased samples in the former and uncertainties in the far-IR emissivity and the optical properties of the dust grains in the latter. This study represents one of the first large-scale categorizations of extragalactic sources behind the SMC and as such we provide the LEPHARE outputs for our full sample of similar to 500 000 sources. © 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement no. 682115). SR acknowledges support from the ERC consolidator grant project STARKEY(grant agreement no. 615604). YC acknowledges support from NSF grant AST 1655677. DMDacknowledges financial support from the Spanish Ministry for Science, Innovation and Universities and FEDER funds through grant AYA2016-81065-C2-2, the State Agency for Research of the Spanish MCIU through the 'Centre of Excellence Severo Ochoa' award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709) and from grant PGC2018-095049B-C21. RRM acknowledges partial support from project BASAL AFB-170002 as well as FONDECYT project no. 1170364. SS acknowledges support from the Science and Engineering Research Board, India through a Ramanujan Fellowship. The authors would like to thank the Cambridge Astronomy Survey Unit (CASU) and theWide Field Astronomy Unit (WFAU) in Edinburgh for providing the necessary data products under the support of the Science and Technology Facilities Council (STFC) in theU.K. The authorswould also like to thank K. Gordon and J. Roman-Duval for discussions related to the use of the HERITAGE dust maps. The authors would like to extend their gratitude to the referee, Geoff Clayton, who provided several comments that improved the manuscript. This study was based on observations made with VISTA at the La Silla Paranal Observatory under programme ID 179.B-2003. This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) collaboration. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the STFC, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo 'a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia, Tecnologia e Inovacao, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi NationalAccelerator Laboratory, the University of Illinois atUrbanaChampaign, the Institut de Ciencies de l'Espai (IEEC/CSIC), the Institut de Fisica d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, the Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. Based on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory (NOAO Prop. ID: 2013A-0411 and 2013B-0440; PI: Nidever), which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. Finally, this project has made extensive use of the Tool for OPerations on Catalogues And Tables (TOPCAT) software package (Taylor 2005) as well as the following open-source PYTHON packages: Astropy (The Astropy Collaboration 2018), MATPLOTLIB (Hunter 2007), NUMPY (Oliphant 2015), PANDAS (McKinney 2010), and SCIPY (Virtanen et al. 2020).