Search

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

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

We report here on the discovery with XMM–Newton of pulsations at 22 ms from the central compact source associated with IKT 16, a supernova remnant in the Small Magellanic Cloud (SMC). The measured spin period and spin period derivative correspond to 21.7661076(2) ms and 2.9(3) × 10−14 s s−1, respectively. Assuming standard spin-down by magnetic dipole radiation, the spin-down power corresponds to 1.1 × 1038 erg s−1 implying a Crab-like pulsar. This makes it the most energetic pulsar discovered in the SMC so far and a close analogue of PSR J0537–6910, a Crab-like pulsar in the Large Magellanic Cloud. The characteristic age of the pulsar is 12 kyr. Having for the first time a period measure for this source, we also searched for the signal in archival data collected in radio with the Parkes telescope and in γ-rays with the Fermi/LAT, but no evidence for pulsation was found in these energy bands.

Warps are vertical distortions of the stellar or gaseous disks of galaxies. One of the proposed scenarios for the formation of warps involves tidal interactions among galaxies. A recent study identified a stellar warp in the outer regions of the south-western (SW) disk of the Large Magellanic Cloud (LMC) and suggested that it might have originated due to the tidal interaction between the LMC and the Small Magellanic Cloud (SMC). Due to the limited spatial coverage of the data, the authors could not investigate the counterpart of this warp in the north-eastern (NE) region, which is essential to understanding the global shape, nature, and origin of the outer LMC warp. In this work, we study the structure of the LMC disk using data on red clump stars from the Gaia Early Data Release 3 (EDR3), which cover the entire Magellanic system. We detected a warp in the NE outer LMC disk which is deviated from the disk plane in the same direction as that of the SW outer warp, but with a lower amplitude. This suggests that the outer LMC disk has an asymmetric stellar warp, which is likely to be a U-shaped warp. Our result provides an observational constraint to the theoretical models of the Magellanic system aimed at improving the understanding the LMC-SMC interaction history., Comment: 9 pages, 6 figures. Accepted for publication in A&A and in press

We show that Smoothed Particle Hydrodynamics (SPH) simulations of dwarf galaxies interacting with a Milky Way-like disk produce moving groups in the simulated stellar disk. We analyze three different simulations: one that includes dwarf galaxies that mimic the Large Magellanic Cloud, Small Magellanic Cloud and the Sagittarius dwarf spheroidal; another with a dwarf galaxy that orbits nearly in the plane of the Milky Way disk; and a null case that does not include a dwarf galaxy interaction. We present a new algorithm to find large moving groups in the $V_R, V_\phi$ plane in an automated fashion that allows us to compare velocity sub-structure in different simulations, at different locations, and at different times. We find that there are significantly more moving groups formed in the interacting simulations than in the isolated simulation. A number of dwarf galaxies are known to orbit the Milky Way, with at least one known to have had a close pericenter approach. Our analysis of simulations here indicates that dwarf galaxies like those orbiting our Galaxy produce large moving groups in the disk. Our analysis also suggests that some of the moving groups in the Milky Way may have formed due to dynamical interactions with perturbing dwarf satellites. The groups identified in the simulations by our algorithm have similar properties to those found in the Milky Way, including similar fractions of the total stellar population included in the groups, as well as similar average velocities and velocity dispersions., Comment: 15 pages, 12 figures, Accepted for publication in MNRAS. Major revisions in the accepted version of the paper

We report the discovery of a new high mass X-ray binary pulsar, XMMU J050722.1-684758, possibly associated with the supernova remnant MCSNR J0507-6847 in the Large Magellanic Cloud, using XMM-Newton X-ray observations. Pulsations with a periodicity of 570 s are discovered from the Be X-ray binary XMMU J050722.1-684758 confirming its nature as a HMXB pulsar. The HMXB is located near the geometric centre of the supernova remnant MCSNR J0507-6847 (0.9 arcmin from the centre) which supports the XRB-SNR association. The estimated age of the supernova remnant is 43-63 kyr which points to a middle aged to old supernova remnant. The large diameter of the supernova remnant combined with the lack of distinctive shell counterparts in optical and radio indicates that the SNR is expanding into the tenous environment of the superbubble N103. The estimated magnetic field strength of the neutron star is $B\gtrsim10^{14}$ G assuming a spin equilibrium condition which is expected from the estimated age of the parent remnant and assuming that the measured mass-accretion rate remained constant throughout., 12 pages, 13 figures, accepted for publication in MNRAS

We present the analysis of the colour-magnitude diagram (CMD) morphology of the ~ 800 Myr old star cluster NGC1831 in the Large Magellanic Cloud, exploiting deep, high-resolution photometry obtained using the Wide Field Camera 3 onboard the Hubble Space Telescope. We perform a simultaneous analysis of the wide upper main sequence and main sequence turn-off observed in the cluster, to verify whether these features are due to an extended star formation or a range of stellar rotation rates, or a combination of these two effects. Comparing the observed CMD with Monte Carlo simulations of synthetic stellar populations, we derive that the morphology of NGC1831 can be fully explained in the context of the rotation velocity scenario, under the assumption of a bimodal distribution for the rotating stars, with ~40% of stars being slow-rotators ($\Omega$ / $\Omega_{crit}$ < 0.5) and the remaining ~ 60% being fast rotators ($\Omega$ / $\Omega_{crit}$ > 0.9). We derive the dynamical properties of the cluster, calculating the present cluster mass and escape velocity, and predicting their past evolution starting at an age of 10 Myr. We find that NGC1831 has an escape velocity $v_{esc}$ = 18.4 km/s, at an age of 10 Myr, above the previously suggested threshold of 15 km/s, below which the cluster cannot retain the material needed to create second-generation stars. These results, combined with those obtained from the CMD morphology analysis, indicate that for the clusters whose morphology cannot be easily explained only in the context of the rotation velocity scenario, the threshold limit should be at least ~ 20 km/s., Comment: 11 pages, 10, figures, 2 tables. Accepted for publication in MNRAS

Modelling dust formation in single stars evolving through the carbon-star stage of the asymptotic giant branch (AGB) reproduces well the mid-infrared colours and magnitudes of most of the C-rich sources in the Large Magellanic Cloud (LMC), apart from a small subset of extremely red objects (EROs). An analysis of the spectral energy distributions of EROs suggests the presence of large quantities of dust, which demand gas densities in the outflow significantly higher than expected from theoretical modelling. We propose that binary interaction mechanisms that involve common envelope (CE) evolution could be a possible explanation for these peculiar stars; the CE phase is favoured by the rapid growth of the stellar radius occurring after C/O overcomes unity. Our modelling of the dust provides results consistent with the observations for mass-loss rates $\dot{M} \sim 5\times 10^{-4}\,{\rm M}_{\odot }$ yr−1, a lower limit to the rapid loss of the envelope experienced in the CE phase. We propose that EROs could possibly hide binaries with orbital periods of about days and are likely to be responsible for a large fraction of the dust production rate in galaxies.

Massive stars are usually called cosmic engines due to their well-known importance at different scales, from shaping their surroundings to even cosmological implications. Binarity is one of the main ingredients in the evolution of massive stars as it has been demonstrated by the high binary fraction of OB-type stars found by several studiesin the Galaxy and the Magellanic Clouds. Binarity can also play a crucial role onother physical properties of massive stars such as rotation and chemical evolution.Furthermore, massive binaries might end their lives as double-compact objects, which are progenitors of gravitational waves. To understand the connection between binarity and the evolution and final fate of massive stars, it is imperative to characterise large and homogeneous samples of massive binaries. In particular, the distributions of their orbital parameters are important inputsfor evolutionary and population-synthesis models used to investigate stellar populations and the chemical evolution of galaxies. Among massive stars, B-type stars are the most numerous representatives and the spectral range between B0 and B3 dominates the rate of core-collapse supernovae. Here I present a detailed study of different aspects of B-type stars in this important spectral range. Firstly, I describe my work for the Galactic O-type Stars Spectroscopic Survey (GOSSS) in Chapter 2. In line with previous instalments of the GOSSS series, a new atlas of galactic B-type giants and supergiants and a new grid of standard stars for spectral classification are presented. Taking one step further, I have developed a quantitative approach to spectral classification of B-type stars, based on the equivalent widths ofspectral lines that have been used historically for this purpose. I propose a classification procedure that uses in first instance the quantitative criteria, and can then use thestandard stars defined in this work for complicated cases. In Chapter 3, I present the B-type Binaries Characterisation (BBC) programme, the largest multi-epoch, spectroscopic study of early B-type binaries in the 30 Doradus region of the Large Magellanic Cloud. In the context of this programme, I havedetermined orbital solutions for a sample of 88 B-type Binaries and the distributions oftheir orbital parameters that are contrasted to other samples of massive stars. Interesting similarities are found when comparing the period distribution of the BBC stars witha sample of O-type stars from the same region of 30 Doradus. This evidence, added to the similar multiplicity fractions of B- and O-type stars, argues in favour of similar formations mechanisms for the complete range of core-collapse-supernova progenitors, but the intrinsic distributions need further study to draw more definite conclusions. The fourth Chapter of this thesis delves into the search of compact companions to B-type stars. Taking advantage of the numerous BBC sample and their derived orbitalproperties, I used the K-P diagram (semi-amplitude velocity vs. orbital period) in Chapter 3 to identify stars with large velocities that could be caused by the presence of massive companions that remained hidden in the spectra of the system. The K-P diagram is explored further in Chapter 4, and strategies are discussed to limit the inclusion of false-positives among the candidates to binaries with compact companions. I have used spectral classifications to estimate primary masses and derive the mass of the secondaries. The best candidates (those with high mass ratios) are discussed further. Finally, Chapter 5 presents a summary of the results presented in this thesis, and describes the possible lines of research that could follow the work developed during my PhD.

Despite mounting evidence that dark matter (DM) exists in the Universe, its fundamental nature remains unknown. We present sensitivity estimates to detect DM particles with a future very-high-energy ($\gtrsim$ TeV) wide field-of-view gamma-ray observatory in the Southern Hemisphere, currently in its research and development phase under the name Southern Wide field-of-view Gamma-ray Observatory (SWGO). This observatory would search for gamma rays from the annihilation or decay of DM particles in many key targets in the Southern sky, such as the Galactic halo, several dwarf galaxies, including the promising Reticulum II, and the Large Magellanic Cloud. With a wide field of view and long exposures, such observatory will have unprecedented sensitivity to DM in the mass range of $\sim$100 GeV to a few PeV from observations of a large fraction of the Galactic halo around the Galactic Center and from Galactic subhalos targets. These results, combined with those from other present and future gamma-ray observatories, will likely probe the thermal relic annihilation cross section of Weakly Interacting Massive Particles for all masses from $\sim$80 TeV down to the GeV range in most annihilation channels.

We present the first satellite system of the Large Binocular Telescope Satellites Of Nearby Galaxies Survey (LBT-SONG), a survey to characterize the close satellite populations of Large Magellanic Cloud to Milky Way-mass, star-forming galaxies in the Local Volume. In this paper, we describe our unresolved diffuse satellite finding and completeness measurement methodology and apply this framework to NGC 628, an isolated galaxy with $\sim1/4$ the stellar mass of the Milky Way. We present two new dwarf satellite galaxy candidates: NGC 628 dwA, and dwB with $\text{M}_{\text{V}}$ = $-12.2$ and $-7.7$, respectively. NGC 628 dwA is a classical dwarf while NGC 628 dwB is a low-luminosity galaxy that appears to have been quenched after reionization. Completeness corrections indicate that the presence of these two satellites is consistent with CDM predictions. The satellite colors indicate that the galaxies are neither actively star-forming nor do they have the purely ancient stellar populations characteristic of ultrafaint dwarfs. Instead, and consistent with our previous work on the NGC 4214 system, they show signs of recent quenching, further indicating that environmental quenching can play a role in modifying satellite populations even for hosts smaller than the Milky Way., 17 pages, 11 figures

We present observations in X-ray and optical emission of the supernova remnant (SNR) 0520–69.4 in the Large Magellanic Cloud. Using XMM–Newton observatory data, we produced images of the diffuse X-ray emission and spectra to obtain the X-ray parameters, such as luminosity and temperature, of hot plasma in the SNR. Diffuse X-ray emission with filled-centre morphology goes beyond the Hα region, suggesting that the hot gas escapes through the pores of the Hα shell. We fitted a model that has a plasma temperature of 1.1 × 107 K for an X-ray thermal luminosity of 3.3 × 1035 erg s−1. However, from Hα and [O iii] Fabry–Perot observations obtained with the Marseille Hα Survey of the Magellanic Clouds and the Milky Way at La Silla, European Southern Observatory, we are able to obtain physical parameters such as the velocity of the shock induced in the cloudlets emitting at optical wavelengths and the electron density of this gas. With the parameters described above, we test the model proposed by White & Long (1991, ApJ, 373, 543) for explaining the mixed-morphology observed.

Recent measurements suggest that the Large Magellanic Cloud (LMC) may weigh as much as 25\% of the Milky Way. In this work we explore how such a large satellite affects mass estimates of the Milky Way based on equilibrium modelling of the stellar halo or other tracers. In particular, we show that if the LMC is ignored, the Milky Way mass is overestimated by as much as 50\%. This bias is due to the bulk motion in the outskirts of the Galaxy's halo and can be, at least in part, accounted for with a simple modification to the equilibrium modelling. Finally, we show that the LMC has a substantial effect on the orbit Leo I which acts to increase its present day speed relative to the Milky Way. We estimate that accounting for a $1.5\times10^{11} M_\odot$ LMC would lower the inferred Milky Way mass to $\sim10^{12} M_\odot$., 7 pages, 6 figures. Submitted to MNRAS. Comments welcome

Nearly all of the well studied ancient globular clusters (GCs), in the Milky Way and in nearby galaxies, show star-to-star variations in specific elements (e.g., He, C, N, O, Na, Al), known as "multiple populations" (MPs). However, MPs are not restricted to ancient clusters, with massive clusters down to $\sim2$ Gyr showing signs of chemical variations. This suggests that young and old clusters share the same formation mechanism but most of the work to date on younger clusters has focused on N-variations. Initial studies even suggested that younger clusters may not host spreads in other elements beyond N (e.g., Na), calling into question whether these abundance variations share the same origin as in the older GCs. In this work, we combine HST photometry with VLT/MUSE spectroscopy of a large sample of RGB stars (338) in the Large Magellanic Cloud cluster NGC 1978, the youngest globular to date with reported MPs in the form of N-spreads. By combining the spectra of individual RGB stars into N-normal and N-enhanced samples, based on the "chromosome map" derived from HST, we search for mean abundance variations. Based on the NaD line, we find a Na-difference of $\Delta$[Na/Fe]$=0.07\pm0.01$ between the populations. While this difference is smaller than typically found in ancient GCs (which may suggest a correlation with age), this result further confirms that the MP phenomenon is the same, regardless of cluster age and host galaxy. As such, these young clusters offer some of the strictest tests for theories on the origin of MPs., Comment: 10 pages, 9 figures, 1 table. Accepted for publication in MNRAS

The dominant gaseous structure in the Galactic halo is the Magellanic Stream. This extended network of neutral and ionized filaments surrounds the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), the two most massive satellite galaxies of the Milky Way1–4. Recent observations indicate that the LMC and SMC are on their first passage around the Galaxy5, that the Magellanic Stream is made up of gas stripped from both clouds2,6,7 and that the majority of this gas is ionized8,9. Although it has long been suspected that tidal forces10,11 and ram-pressure stripping12,13 contributed to the formation of the Magellanic Stream, models have not been able to provide a full understanding of its origins3. Several recent developments—including the discovery of dwarf galaxies associated with the Magellanic group14–16, determination of the high mass of the LMC17, detection of highly ionized gas near stars in the LMC18,19 and predictions of cosmological simulations20,21—support the existence of a halo of warm (roughly 500,000 kelvin) ionized gas around the LMC (the ‘Magellanic Corona’). Here we report that, by including this Magellanic Corona in hydrodynamic simulations of the Magellanic Clouds falling onto the Milky Way, we can reproduce the Magellanic Stream and its leading arm. Our simulations explain the filamentary structure, spatial extent, radial-velocity gradient and total ionized-gas mass of the Magellanic Stream. We predict that the Magellanic Corona will be unambiguously observable via high-ionization absorption lines in the ultraviolet spectra of background quasars lying near the LMC. Embedding the Magellanic Clouds in a corona of ionized gas allows the gaseous Magellanic Stream to be modelled accurately.

Stars slingshotted by the supermassive black hole at the Galactic Centre escape from the Milky Way so quickly that their trajectories are almost straight lines. Previous works have shown how these ‘hypervelocity stars’ (stars moving faster than the local Galactic escape speed) are subsequently deflected by the gravitational field of the Milky Way and the Large Magellanic Cloud (LMC), but have neglected to account for the reflex motion of the Milky Way in response to the fly-by of the LMC. A consequence of this motion is that the hypervelocity stars we see in the outskirts of the Milky Way today were ejected from where the Milky Way centre was hundreds of millions of years ago. This change in perspective causes large apparent deflections of several degrees in the trajectories of the hypervelocity stars. We quantify these deflections by simulating the ejection of hypervelocity stars from an isolated Milky Way (with a spherical or flattened dark matter halo), from a fixed-in-place Milky Way with a passing LMC, and from a Milky Way that responds to the passage of the LMC, finding that LMC passage causes larger deflections than can be caused by a flattened Galactic dark matter halo in ΛCDM. The $10\, \mu \mathrm{as}\, \mathrm{yr}^{-1}$ proper motion precision necessary to measure these deflections will be possible with the combination of Gaia with the proposed GaiaNIR successor mission, and these measurements will directly probe the shape of the Milky Way, the mass of the LMC, and the dance of these two galaxies.

Full author list: A Pieres, L Girardi, E Balbinot, B Santiago, L N da Costa, A Carnero Rosell, A B Pace, K Bechtol, M A T Groenewegen, A Drlica-Wagner, T S Li, M A G Maia, R L C Ogando, M dal Ponte, H T Diehl, A Amara, S Avila, E Bertin, D Brooks, D L Burke, M Carrasco Kind, J Carretero, J De Vicente, S Desai, T F Eifler, B Flaugher, P Fosalba, J Frieman, J García-Bellido, E Gaztanaga, D W Gerdes, D Gruen, R A Gruendl, J Gschwend, G Gutierrez, D L Hollowood, K Honscheid, D J James, K Kuehn, N Kuropatkin, J L Marshall, R Miquel, A A Plazas, E Sanchez, S Serrano, I Sevilla-Noarbe, E Sheldon, M Smith, M Soares-Santos, F Sobreira, E Suchyta, M E C Swanson, G Tarle, D Thomas, V Vikram, A R Walker, We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from TRILEGAL models to real data. We apply this technique, which we call MWFITTING, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ~2300 deg2. We further explore a smaller subset (~1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of MWFITTING in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid., The DESDM is supported by the National Science Foundation under grants AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2015-71825, ESP2015-66861, FPA2015-68048, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020, and the Brazilian Instituto Nacional de Ci?ncia e Tecnologia (INCT) do e-Universe (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. 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 Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Centre for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Centre for Cosmology and AstroParticle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundac¸ão Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovac¸ão, the Deutsche Forschungsgemeinschaft, and the Collaborating Institutions in the DES.

Galaxy interactions distort the distribution of baryonic matter and can affect star formation. The nearby Magellanic Clouds are a prime example of an ongoing galaxy interaction process. Here we use the intermediate-age ($\sim1$-$10$ Gyr) red clump stars to map the three-dimensional structure of the Small Magellanic Cloud (SMC) and interpret it within the context of its history of interaction with the Large Magellanic Cloud (LMC) and the Milky Way. Red clump stars are selected from near-infrared colour-magnitude diagrams based on data from the VISTA survey of the Magellanic Clouds. Interstellar reddening is measured and removed, and the corrected brightness is converted to a distance, on a star-by-star basis. A flat plane fitted to the spatial distribution of red clump stars has an inclination $i=35\deg$-$48\deg$ and position angle PA$=170\deg$-$186\deg$. However, significant deviations from this plane are seen, especially in the periphery and on the eastern side of the SMC. In the latter part, two distinct populations are present, separated in distance by as much as 10 kpc. Distant red clump stars are seen in the North of the SMC, and possibly also in the far West; these might be associated with the predicted `Counter-Bridge'. We also present a dust reddening map, which shows that dust generally traces stellar mass. The structure of the intermediate-age stellar component of the SMC bears the imprints of strong interaction with the LMC a few Gyr ago, which cannot be purely tidal but must have involved ram pressure stripping., Comment: 16 pages, 11 figures; accepted for publication in MNRAS

We use Gaia DR2 systemic proper motions of 45 satellite galaxies to constrain the mass of the Milky Way using the scale free mass estimator of Watkins et al. (2010). We first determine the anisotropy parameter $\beta$, and the tracer satellites' radial density index $\gamma$ to be $\beta$=$-0.67^{+0.45}_{-0.62}$ and $\gamma=2.11\pm0.23$. When we exclude possible former satellites of the Large Magellanic Cloud, the anisotropy changes to $\beta$=$-0.21^{+0.37}_{-0.51}$. We find that the index of the Milky Way's gravitational potential $\alpha$, which is dependent on the mass itself, is the parameter with the largest impact on the mass determination. Via comparison with cosmological simulations of Milky Way-like galaxies, we carried out a detailed analysis of the estimation of the observational uncertainties and their impact on the mass estimator. We found that the mass estimator is biased when applied naively to the satellites of simulated Milky Way halos. Correcting for this bias, we obtain for our Galaxy a mass of $0.58^{+0.15}_{-0.14}\times10^{12}$M$_\odot$ within 64 kpc, as computed from the inner half of our observational sample, and $1.43^{+0.35}_{-0.32}\times10^{12}$M$_\odot$ within 273 kpc, from the full sample; this latter value extrapolates to a virial mass of $M_\mathrm{vir\,\Delta=97}$=$1.51^{+0.45}_{-0.40} \times 10^{12}M_{\odot}$ corresponding to a virial radius of R$_\mathrm{vir}$=$308\pm29$ kpc. This value of the Milky Way mass lies in-between other mass estimates reported in the literature, from various different methods., Comment: 17 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

Context. It is well known that massive O stars are frequently (if not always) found in binary or higher-order multiple systems, but this fact has been less robustly investigated for the lower mass range of the massive stars, represented by B-type stars. Obtaining the binary fraction and orbital parameter distributions of B-type stars is crucial to understand the impact of multiplicity on the archetypal progenitor of core-collapse supernovae as well as to properly investigate formation channels for gravitational wave progenitors. Aims. This work aims to characterise the multiplicity of the B star population of the young open cluster NGC 6231 through multi-epoch optical spectroscopy of 80 B-type stars. Methods. We analyse 31 FLAMES/GIRAFFE observations of 80 B-type stars, monitoring their radial velocities (RVs) and performing a least-squares spectral analysis (Lomb-Scargle) to search for periodicity in those stars with statistically significant variability in their RVs. Results. We constrained an observed spectroscopic binary fraction of 33 ± 5% for the B-type stars of NGC 6231, with a first order bias correction giving a true spectroscopic binary fraction of 52 ± 8%. Out of 27 B-type binary candidates, we obtained orbital solutions for 20 systems: 15 single-lined (SB1) and five double-lined spectroscopic binaries (SB2s). We present these orbital solutions and the orbital parameter distributions associated with them. Conclusions. Our results indicate that Galactic B-type stars are less frequently found in binary systems than their more massive O-type counterparts, but their orbital properties generally resemble those of B- and O-type stars in both the Galaxy and Large Magellanic Cloud.

The highly-substructured outskirts of the Magellanic Clouds provide ideal locations for studying the complex interaction history between both Clouds and the Milky Way (MW). In this paper, we investigate the origin of a >20$^\circ$ long arm-like feature in the northern outskirts of the Large Magellanic Cloud (LMC) using data from the Magellanic Edges Survey (MagES) and Gaia EDR3. We find that the arm has a similar geometry and metallicity to the nearby outer LMC disk, indicating that it is comprised of perturbed disk material. Whilst the azimuthal velocity and velocity dispersions along the arm are consistent with those in the outer LMC, the in-plane radial velocity and out-of-plane vertical velocity are significantly perturbed from equilibrium disk kinematics. We compare these observations to a new suite of dynamical models of the Magellanic/MW system, which describe the LMC as a collection of tracer particles within a rigid potential, and the SMC as a rigid Hernquist potential. Our models indicate the tidal force of the MW during the LMC's infall is likely responsible for the observed increasing out-of-plane velocity along the arm. Our models also suggest close LMC/SMC interactions within the past Gyr, particularly the SMC's pericentric passage ~150 Myr ago and a possible SMC crossing of the LMC disk plane ~400 Myr ago, likely do not perturb stars that today comprise the arm. Historical interactions with the SMC prior to ~1 Gyr ago may be required to explain some of the observed kinematic properties of the arm, in particular its strongly negative in-plane radial velocity., Accepted by MNRAS

Measurements of [Fe/H] and radial velocity are presented for 89 RR Lyrae (RRL) candidates within 6 kpc of the Sun. After removal of two suspected non-RRL, these stars were added to an existing database, which yielded 464 RRL with [Fe/H] on a homogeneous scale. Using data from the Gaia satellite (Data Release 2), we calculated the positions and space velocities for this sample. These data confirm the existence of a thin-disc of RRL with [alpha/Fe] ~ solar. The majority of the halo RRL with large total energies have near zero angular momenta about the Z axis. Kinematically these stars closely resemble the Gaia-Sausage/Gaia-Enceladus stars that others have proposed are debris from the merger of a large galaxy with the Milky Way. The metallicity and period distributions of the RRL and their positions in the period-amplitude diagram suggest that this disrupted galaxy was as massive as the Large Magellanic Cloud and possibly greater., Comment: 18 pages, 13 figures

Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way. A local wake is predicted to trail the orbit of the LMC, while a large-scale over-density is predicted to exist across a large area of the northern Galactic hemisphere. Here we present the detection of both the local wake and Northern over-density (hereafter the "collective response") in an all-sky star map of the Galaxy based on 1301 stars at 60, Comment: To appear in the 22 April 2021 issue of Nature

Because the galaxies of the Local Group have such large angular sizes, much of their diffuse, large-angular-scale emission is filtered out by the Herschel data reduction process. In this work, we restore this previously missed dust in Herschel observations of the Large Magellanic Cloud, Small Magellanic Cloud, M31, and M33. We do this by combining Herschel data (including new reductions for the Magellanic Clouds), in Fourier space, with lower-resolution data from all-sky surveys (Planck, IRAS, and COBE) that did not miss the extended emission. With these new maps, we find that a significant amount of emission was missing from uncorrected Herschel data of these galaxies; over 20% in some bands. Our new photometry also resolves the disagreement between fluxes reported from older HERITAGE Magellanic Cloud Herschel reductions, and fluxes reported from other telescopes. More emission is restored in shorter wavelength bands, especially in the galaxies' peripheries, making these regions 20-40% bluer than before. We also find that the Herschel-PACS instrument response conflicts with the all-sky data, over the 20-90' angular scales to which they are both sensitive, by up to 31%. By binning our new data based on hydrogen column density, we are able to detect emission from dust at low ISM densities (at $\Sigma_{\rm H} < 1\,{\rm M_{\odot} pc^{-2}}$ in some cases), and are able to detect emission at much lower densities (a factor of 2.2 lower on average, and more than a factor of 7 lower in several cases) than was possible with uncorrected data., Comment: Accepted for publication in ApJ

Most of the observed and simulated galactic bars are symmetric in the face-on view. However, there are indeed cases of bars that are off-center with respect to the disk or have an asymmetric shape. The only well-known example showing both these features is the Large Magellanic Cloud. We report on the identification of several lopsided galactic bars in the Illustris TNG100 simulation found among a sample of elongated bar-like galaxies studied in the past. The bars show a clear asymmetry in the face-on view, which is in the shape of a footprint. We measured the evolution of the different parameters of the bars' shape and asymmetry as a function of time and find that the asymmetry is preserved for a few Gyr. It can grow together with a bar or appear later, after bar formation. We considered two scenarios leading to the formation of lopsided bars using controlled simulations. In the first, a Milky Way-like galaxy interacts with a massive companion placed on a radial orbit in the plane of the disk and perpendicular to the orientation of the bar at the time of the first passage. In the second, the galaxy initially has an off-center disk and the growth of the bar and its asymmetry is more similar to the one found in IllustrisTNG galaxies, where it is also preceded by the presence of an asymmetric disk. It is thus possible that lopsided bars are born in lopsided disks, although in some cases, the time difference between the occurrence of the asymmetry in the two components is quite large., 12 pages, 11 figures, accepted for publication in A&A

We develop an automatic bubble-recognition routine based on Minkowski functionals (MF) and tensors (MT) to detect bubble-like interstellar structures in optical emission line images. Minkowski functionals and MT are powerful mathematical tools for parameterizing the shapes of bodies. Using the papaya2-library, we created maps of the desired MF or MT of structures at a given window size. We used maps of the irreducible MT $\psi_2$, which is sensitive to elongation, to find filamentary regions in H$\alpha$, [SII], and [OIII] images of the Magellanic Cloud Emission Line Survey (MCELS). Using the phase of $\psi_2$, we were able to draw lines perpendicular to each filament and thus obtain line-density maps. This allowed us to find the center of a bubble-like structure and to detect structures at different window sizes. The detected bubbles in all bands are spatially correlated to the distribution of massive stars, showing that we indeed detect interstellar bubbles without large spatial bias. Eighteen out of 59 supernova remnants in the Large Magellanic Cloud (LMC) and 13 out of 20 superbubbles are detected in at least one wavelength. The lack of detection is mostly due to surrounding emission that disturbs the detection, a too small size, or the lack of a (circular) counterpart in our emission line images. In line-density maps at larger scales, maxima can be found in regions with high star formation in the past, often inside supergiant shells (SGS). In SGS LMC 2, there is a maximum west of the shell where a collision of large gas clouds is thought to have occurred. In the Small Magellanic Cloud (SMC), bubble detection is impaired by the more complex projected structure of the galaxy. Line maps at large scales show large filaments in the SMC in a north-south direction, especially in the [SII] image. The origin of these filaments is unknown., Comment: 19 pages, 27 figures. Accepted for publication in Astronomy & Astrophysics

We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galaxies, we find quenching times older than 11.5 Gyr ago, compatible with the scenario in which reionization suppresses the star formation of small dark matter halos. However, our analysis also reveals some differences in the SFHs of candidate Magellanic Cloud satellites, i.e., galaxies that are likely satellites of the Large Magellanic Cloud and that entered the Milky Way potential only recently. Indeed, Magellanic satellites show quenching times about 600 Myr more recent with respect to those of other Milky Way satellites, on average, even though the respective timings are still compatible within the errors. This finding is consistent with theoretical models that suggest that satellites' SFHs may depend on their host environment at early times, although we caution that within the error bars all galaxies in our sample are consistent with being quenched at a single epoch., 7 pages, 3 figures, 2 tables. Accepted for publication in ApJL

We present Period-Luminosity Relations (PLRs) for 55 Cepheids in M31 with periods ranging from 4 to 78 days observed with the Hubble Space Telescope (HST) using the same three-band photometric system recently used to calibrate their luminosities. Images were taken with the Wide Field Camera 3 in two optical filters (F555W and F814W) and one near-infrared filter (F160W) using the Drift And SHift (DASH) mode of operation to significantly reduce overheads and observe widely-separated Cepheids in a single orbit. We include additional F160W epochs for each Cepheid from the Panchromatic Hubble Andromeda Treasury (PHAT) and use light curves from the Panoramic Survey Telescope and Rapid Response System of the Andromeda galaxy (PAndromeda) project to determine mean magnitudes. Combined with a 1.28$\%$ absolute calibration of Cepheid PLRs in the Large Magellanic Cloud from Riess et al. (2019a) in the same three filters, we find a distance modulus to M31 of $\mu_0$ = 24.407 $\pm$ 0.032, corresponding to 761 $\pm$ 11 kpc and 1.49$\%$ uncertainty including all error sources, the most precise determination of its distance to date. We compare our results to past measurements using Cepheids and the Tip of the Red Giant Branch (TRGB). This study also provides the groundwork for turning M31 into a precision anchor galaxy in the cosmic distance ladder to measure the Hubble constant together with efforts to measure a fully geometric distance to M31., Comment: Accepted by ApJ, 19 pages, 5 figures, 6 tables

We comprehensively study the variability of Miras in the Large Magellanic Cloud (LMC) by simultaneous analysing light curves in 14 bands in the range of 0.5$-$24 microns. We model over 20-years-long, high cadence $I$-band light curves collected by The Optical Gravitational Lensing Experiment (OGLE) and fit them to light curves collected in the remaining optical/near-infrared/mid-infrared bands to derive both the variability amplitude ratio and phase-lag as a function of wavelength. We show that the variability amplitude ratio declines with the increasing wavelength for both oxygen-rich (O-rich) and carbon-rich (C-rich) Miras, while the variability phase-lag increases slightly with the increasing wavelength. In a significant number of Miras, mostly the C-rich ones, the spectral energy distributions (SEDs) require a presence of a cool component (dust) in order to match the mid-IR data. Based on SED fits for a golden sample of 140 Miras, we calculated synthetic period-luminosity relations (PLRs) in 42 bands for the existing and future sky surveys that include OGLE, The VISTA Near-Infrared $YJK_\mathrm{s}$ Survey of the Magellanic Clouds System (VMC), Legacy Survey of Space and Time (LSST), Gaia, Spitzer, The Wide-field Infrared Survey Explorer (WISE), The James Webb Space Telescope (JWST), The Nancy Grace Roman Space Telescope (formerly WFIRST), and The Hubble Space Telescope (HST). We show that the synthetic PLR slope decreases with increasing wavelength for both the O-rich and C-rich Miras in the range of 0.1$-$40 microns. Finally, we show the location and motions of Miras on the color-magnitude (CMD) and color-color (CCD) diagrams., Accepted for publication in The Astrophysical Journal Supplement Series. 27 pages, 16 figures, 7 tables. The full machine readable tables will be available in the online journal just after publication

Within $��$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultra faint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association., 19 pages, 9 figures, 4 tables. Accepted for publication in MNRAS

A standard prediction of galaxy formation theory is that the ionizing background suppresses galaxy formation in haloes with peak circular velocities smaller than Vpeak ~ 20 km/s, rendering the majority of haloes below this scale completely dark. We use a suite of cosmological zoom simulations of Milky Way-like haloes that include central Milky Way disk galaxy potentials to investigate the relationship between subhaloes and ultrafaint galaxies. We find that there are far too few subhaloes within 50 kpc of the Milky Way that had Vpeak > 20 km/s to account for the number of ultrafaint galaxies already known within that volume today. In order to match the observed count, we must populate subhaloes down to Vpeak ~ 6 km/s with ultrafaint dwarfs. The required haloes have peak virial temperatures as low as 1,500 K, well below the atomic hydrogen cooling limit of 10^4 K. Allowing for the possibility that the Large Magellanic Cloud contributes several of the satellites within 50 kpc could potentially raise this threshold to 10 km/s (4,000 K), still below the atomic cooling limit and far below the nominal reionization threshold., 10 pages, 5 figures, 1 table. Submitted to MNRAS

The properties of carbon stars in the Magellanic Clouds (MCs) and their total dust production rates are predicted by fitting their spectral energy distributions (SED) over pre-computed grids of spectra reprocessed by dust. The grids are calculated as a function of the stellar parameters by consistently following the growth for several dust species in their circumstellar envelopes, coupled with a stationary wind. Dust radiative transfer is computed taking as input the results of the dust growth calculations. The optical constants for amorphous carbon are selected in order to reproduce different observations in the infrared and optical bands of \textit{Gaia} Data Release 2. We find a tail of extreme mass-losing carbon stars in the Large Magellanic Cloud (LMC) with low gas-to-dust ratios that is not present in the Small Magellanic Cloud (SMC). Typical gas-to-dust ratios are around $700$ for the extreme stars, but they can be down to $\sim160$--$200$ and $\sim100$ for a few sources in the SMC and in the LMC, respectively. The total dust production rate for the carbon star population is $\sim 1.77\pm 0.45\times10^{-5}$~M$_\odot$~yr$^{-1}$, for the LMC, and $\sim 2.52\pm 0.96 \times 10^{-6}$~M$_\odot$~yr$^{-1}$, for the SMC. The extreme carbon stars observed with the Atacama Large Millimeter Array and their wind speed are studied in detail. For the most dust-obscured star in this sample the estimated mass-loss rate is $\sim 6.3 \times 10^{-5}$~M$_\odot$~yr$^{-1}$. The grids of spectra are available at: https://ambrananni085.wixsite.com/ambrananni/online-data-1 and included in the SED-fitting python package for fitting evolved stars https://github.com/s-goldman/Dusty-Evolved-Star-Kit ., 22 pages, 16 figures, accepted for publication in MNRAS

Three decades have passed since the supernova SN 1987A was observed in the Large Magellanic Cloud, inside which the product is most likely a neutron star (NS) formed in the core collapse explosion.Although lots of observations with sensitive radio telescopes have taken place, astronomers have not yet detected any evidence for a radio pulsar around the remnant of 1987A. To investigate pulsars inside the SN remnants, we calculate the cut-off oscillation frequency of the plasma around the presumed NS inside SN1987A, as shown to be about 33 GHz at present (2018 CE), which is much higher than the favorite “searching window” (e.g. L-band ∼ 1.4 GHz) of radio pulsar surveys that have been commonly exploited by astronomers. Since radio waves with frequencies lower than the plasma cut-off frequency cannot penetrate the SN remnant media, we suggest that astronomers use higher frequency bands to search for a pulsar in SN 1987A.Furthermore, with the expansion of SN remnant media, we find that the plasma cut-off frequency can decay to the L-band (1.4 GHz) in the future. The strategy of finding a pulsar of SN 1987A is that either the high frequency bands of radio telescopes, or the high energy detections at Gamma-ray and X-ray bands by space satellites are applied.

Observational work have shown that the two brightest satellite galaxies of the Milky Way (MW), the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), are rare amongst MW analogues. It is then interesting to know whether the presence of massive satellite has any effect on the whole satellite population in MW analogues. In this article, we investigate this problem using a semi-analytical model combined with the Millennium-II Simulation. MW-analogous galaxies are defined to have similar stellar mass or dark matter halo mass to the MW. We find that, in the first case, the halo mass is larger and there are, on average, twice as many satellites in Milky Way analogs if there is a massive satellite galaxy in the system. This is mainly from the halo formation bias. The difference is smaller if MW analogues are selected using halo mass. We also find that the satellites distribution is slightly asymmetric, being more concentrated on the line connecting the central galaxy and the massive satellite and that, on average, LMC have brought in 14.7 satellite galaxies with $M_{r}, 9 pages, 10 figures, published in MNRAS

Variable stars with well-calibrated period-luminosity relationships provide accurate distance measurements to nearby galaxies and are therefore a vital tool for cosmology and astrophysics. While these measurements typically rely on samples of Cepheid and RR-Lyrae stars, abundant populations of luminous variable stars with longer periods of $10 - 1000$ days remain largely unused. We apply machine learning to derive a mapping between lightcurve features of these variable stars and their magnitude to extend the traditional period-luminosity (PL) relation commonly used for Cepheid samples. Using photometric data for long period variable stars in the Large Magellanic cloud (LMC), we demonstrate that our predictions produce residual errors comparable to those obtained on the corresponding Cepheid population. We show that our model generalizes well to other samples by performing a blind test on photometric data from the Small Magellanic Cloud (SMC). Our predictions on the SMC again show small residual errors and biases, comparable to results that employ PL relations fitted on Cepheid samples. The residual biases are complementary between the long period variable and Cepheid fits, which provides exciting prospects to better control sources of systematic error in cosmological distance measurements. We finally show that the proposed methodology can be used to optimize samples of variable stars as standard candles independent of any prior variable star classification., Comment: 14 pages, 10 figures, 1 table, updated to match the version accepted by the MNRAS

We carried out a spectroscopic follow-up program of the four new stellar stream candidates detected by Belokurov & Koposov (2016) in the outskirts of the Large Magellanic Cloud (LMC) using FORS2 (VLT). The medium-resolution spectra were used to measure the line-of-sight velocities, estimate stellar metallicities and to classify stars into Blue Horizontal Branch (BHB) and Blue Straggler (BS) stars. Using the 4-D phase-space information, we attribute approximately one half of our sample to the Magellanic Clouds, while the rest is part of the Galactic foreground. Only two of the four stream candidates are confirmed kinematically. While it is impossible to estimate the exact levels of MW contamination, the phase-space distribution of the entire sample of our Magellanic stars matches the expected velocity gradient for the LMC halo and extends as far as 33 deg (angular separation) or 29 kpc from the LMC center. Our detections reinforce the idea that the halo of the LMC seems to be larger than previously expected, and its debris can be spread in the sky out to very large separations from the LMC center. Finally, we provide some kinematic evidence that many of the stars analysed here have likely come from the Small Magellanic Cloud., 15 pages, 13 figures, accepted for publication in MNRAS

PSR J0540$-$6919 is the second-most energetic radio pulsar known and resides in the Large Magellanic Cloud. Like the Crab pulsar it is observed to emit giant radio pulses (GPs). We used the newly-commissioned PTUSE instrument on the MeerKAT radio telescope to search for GPs across three observations. In a total integration time of 5.7 hrs we detected 865 pulses above our 7$\sigma$ threshold. With full polarisation information for a subset of the data, we estimated the Faraday rotation measure, $\rm{RM}=-245.8 \pm 1.0$ rad m$^{-2}$ toward the pulsar. The brightest of these pulses is $\sim$ 60% linearly polarised but the pulse-to-pulse variability in the polarisation fraction is significant. We find that the cumulative GP flux distribution follows a power law distribution with index $-2.75 \pm 0.02$. Although the detected GPs make up only $\sim$ 10% of the mean flux, their average pulse shape is indistinguishable from the integrated pulse profile, and we postulate that there is no underlying emission. The pulses are scattered at L-band frequencies with the brightest pulse exhibiting a scattering time-scale of $\tau = 0.92 \pm 0.02$ ms at 1.2 GHz. We find several of the giants display very narrow-band "flux knots" similar to those seen in many Fast Radio Bursts, which we assert cannot be due to scintillation or plasma lensing. The GP time-of-arrival distribution is found to be Poissonian on all but the shortest time-scales where we find four GPs in six rotations, which if GPs are statistically independent is expected to occur in only 1 of 7000 observations equivalent to our data., Comment: 15 pages, 10 figures

In a failed supernova, partial ejection of the progenitor's outer envelope can occur due to weakening of the core's gravity by neutrino emission in the protoneutron star phase. We consider emission when this ejecta sweeps up the circumstellar material, analogous to supernova remnants (SNRs). We focus on failed explosions of blue supergiants, and find that the emission can be bright in soft X-rays. Due to its soft emission, we find that sources in the Large Magellanic Cloud (LMC) are more promising to detect than those in the Galactic disk. These remnants are characteristic in smallness ($\lesssim 10$ pc) and slowness (100s of ${\rm km\ s^{-1}}$) compared to typical SNRs. Although the expected number of detectable sources is small (up to a few by eROSITA 4-year all-sky survey), prospects are better for deeper surveys targeting the LMC. Detection of these failed SNRs will realize observational studies of mass ejection upon black hole formation., 6 pages, 4 figures (content identical to v1). Accepted for publication in PASJ

Massive star evolution at low metallicity is closely connected to many fields in high-redshift astrophysics, but poorly understood. The Small Magellanic Cloud (SMC) is a unique laboratory to study this because of its metallicity of 0.2 Zsol, its proximity, and because it is currently forming stars. We used a spectral type catalog in combination with GAIA magnitudes to calculate temperatures and luminosities of bright SMC stars. By comparing these with literature studies, we tested the validity of our method, and using GAIA data, we estimated the completeness of stars in the catalog as a function of luminosity. This allowed us to obtain a nearly complete view of the most luminous stars in the SMC. When then compared with stellar evolution predictions. We also calculated the extinction distribution, the ionizing photon production rate, and the star formation rate. Our results imply that the SMS hosts only 30 very luminous main-sequence stars (M > 40 Msol; L > 10^5 Lsol), which are far fewer than expected from the number of stars in the luminosity range 3*10^4 < L/Lsol < 3*10^5 and from the typically quoted star formation rate in the SMC. Even more striking, we find that for masses above M > 20 Msol, stars in the first half of their hydrogen-burning phase are almost absent. This mirrors a qualitatively similar peculiarity that is known for the Milky Way and Large Magellanic Cloud. This amounts to a lack of hydrogen-burning counterparts of helium-burning stars, which is more pronounced for higher luminosities. We argue that a declining star formation rate or a steep initial mass function are unlikely to be the sole explanations for the dearth of young bright stars. Instead, many of these stars might be embedded in their birth clouds, although observational evidence for this is weak. We discuss implications for cosmic reionization and the top end of the initial mass function., 26 pages, 23 figures. Astronomy and Astrophysics, in press

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($\sigma_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $\sigma_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries., Comment: 7 pages. 5 figures. Accepted for publication in Astronomy & Astrophysics

We perform a novel test of general relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids' periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler's third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be ${G}_{\mathrm{LMC}}/{G}_{\mathrm{SS}}=0.9{3}_{\ensuremath{-}0.04}^{+0.05}$. We discuss the implications of this 5% measurement of Newton's constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers ${G}_{\mathrm{LMC}}l{G}_{\mathrm{SS}}$ at $2.6\ensuremath{\sigma}$. We also obtain new bounds on critical parameters that appear in semianalytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $\ensuremath{\alpha}=0.9{0}_{\ensuremath{-}0.26}^{+0.36}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

This paper presents results from photometric and statistical-parallax analysis of a sample of 850 field RR Lyrae (RRL) variables. The photometric and spectroscopic data for sample RRLs are obtained from (1) our new spectroscopic observations (for 448 RRLs) carried out with the Southern African Large Telescope (SALT); (2) our photometric observations using the 1.0-m telescope of the South African Astronomical Observatory (SAAO), and (3) literature. These are combined with accurate proper motion data from the second release of \textit{Gaia} mission (DR2). This study primarily determines the velocity distribution of solar neighborhood RRLs, and it also calibrates the zero points of the RRLs visual V-band luminosity-metallicity (LZ or $M_V-$[\text{Fe/H}]) relation and their period-luminosity-metallicity (PLZ) relations in the \textit{WISE} $W_1-$ and \textit{2MASS} $Ks-$band. The calibrated PLZ and LZ relations are used to estimate the Galactic Center distance and the distance modulus of the Large Magellanic Cloud (LMC), which are found to be 7.99$\pm$0.49\,kpc and 18.46$\pm$0.09 \,mag, respectively. All our results are in excellent agreement with available literature based on statistical parallax analysis, but are considerably more accurate and precise. Moreover, the zero-points of our calibrated PLZ and LZ relations are quite consistent with current results found by other techniques and yield the LMC distance modulus that is within 0.04\,mag of the current most precise estimate., 20 pages,21 figueres

We report on the extinction properties in the fields around the clusters NGC 1854, NGC 1856, and NGC 1858 in the bar of the Large Magellanic Cloud. The colour-magnitude diagrams of the stars in all these regions show an elongated red giant clump that reveals a variable amount of extinction across these fields, ranging from Av~0.2 to Av~1.9, including Galactic foreground extinction. The extinction properties nonetheless are remarkably uniform. The slope of the reddening vectors measured in the (V-I,V) and (B-I,B) colour-magnitude planes is fully in line with the Av/E(B-V)~5.5 value found in the outskirts of 30 Dor. This indicates the presence of an additional grey extinction component in the optical requiring big grains to be about twice as abundant as in the diffuse Galactic interstellar medium (ISM). Areas of higher extinction appear to be systematically associated with regions of more intense star formation, as measured by the larger number of stars more massive than 8 Msun, thus making injection of big grains into the ISM by SNII explosion the likely mechanism at the origin of the observed grey extinction component., Comment: 12 pages, 7 figures, accepted for publication in The Astrophysical Journal

We present time-series photometry of 21 nearby Type II Cepheids in the near-infrared J, H and Ks passbands. We use this photometry, together with the Third Gaia Early Data Release parallaxes, to determine for the first time period-luminosity relations (PLRs) for Type II Cepheids from field representatives of these old pulsating stars in the near-infrared regime. We found PLRs to be very narrow for BL Herculis stars, which makes them candidates for precision distance indicators. We then use archival photometry and the most accurate distance obtained from eclipsing binaries to recalibrate PLRs for Type II Cepheids in the Large Magellanic Cloud (LMC). Slopes of our PLRs in the Milky Way and in the LMC differ by slightly more than 2{\sigma} and are in a good agreement with previous studies of the LMC, Galactic Bulge and Galactic Globular Clusters Type II Cepheids samples. We use PLRs of Milky Way Type II Cepheids to measure the distance to the LMC and we obtain a distance modulus of 18.540$\pm$0.026(stat.)$\pm$0.034(syst.)mag in the WJKs Wesenheit index. We also investigate the metallicity effect within our Milky Way sample and we find rather significant value of about -0.2mag/dex in each band meaning that more metal-rich Type II Cepheids are intrinsically brighter than their more metal-poor counterparts, in agreement with the value obtained from Type II Cepheids in Galactic Globular Clusters. The main source of systematic error on our Milky Way PLRs calibration and the LMC distance is the current uncertainty of the Gaia parallax zero point., Comment: Accepted for ApJ

In the first months after the launch in July 2019, eROSITA onboard Spektr-RG (SRG) performed long-exposure observations in the regions around SN 1987A and SNR N132D in the Large Magellanic Cloud (LMC). We analyse the distribution and the spectrum of the diffuse X-ray emission in the observed fields to determine the physical properties of the hot phase of the interstellar medium (ISM). The eROSITA data are complemented by newly derived column density maps for the Milky Way and the LMC, 888 MHz radio continuum map from the Australian Square Kilometer Array Pathfinder (ASKAP), and optical images of the Magellanic Cloud Emission Line Survey (MCELS). We detect significant emission from thermal plasma with kT=0.2 keV in all the regions. There is also an additional higher-temperature emission component from a plasma with kT = 0.7 keV. In addition, non-thermal X-ray emission is significantly detected in the superbubble 30 Dor C. The absorbing column density NH in the LMC derived from the analysis of the X-ray spectra taken with eROSITA is consistent with the NH obtained from the emission of the cold medium over the entire area. Neon abundance is enhanced in the regions in and around 30 Dor and SN 1987A, indicating that the ISM has been chemically enriched by the young stellar population. Emission from the stellar cluster RMC 136 and the Wolf-Rayet stars RMC 139 and RMC 140 is best modelled with a high-temperature (kT>1 keV) non-equilibrium ionisation plasma emission and a non-thermal component with a photon index of {\Gamma} =1.3. In addition, the optical SNR candidate J0529-7004 is also detected with eROSITA and we thus confirm the source as an SNR., Comment: Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We report on the detection of a black hole (NGC 1850 BH1) in the ∼100-Myr-old massive cluster NGC 1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off star (4.9 ± 0.4 M·), which is starting to fill its Roche lobe and is becoming distorted. Using 17 epochs of Very Large Telescope/Multi-Unit Spectroscopic Explorer observations, we detected radial velocity variations exceeding 300 km s-1 associated with the target star, linked to the ellipsoidal variations measured by the fourth phase of the Optical Gravitational Lensing Experiment in the optical bands. Under the assumption of a semidetached system, the simultaneous modelling of radial velocity and light curves constrains the orbital inclination of the binary to 38° ± 2°, resulting in a true mass of the unseen companion of 11.1 -2.4+2.1\M⊙. This represents the first direct dynamical detection of a black hole in a young massive cluster, opening up the possibility of studying the initial mass function and the early dynamical evolution of such compact objects in high-density environments., SS, NB and ICZ acknowledge financial support from the European Research Council (ERC-CoG-646928, Multi-Pop). SK acknowledges funding from the UKRI in the form of a Future Leaders Fellowship (grant no. MR/T022868/1). CU acknowledges the support of the Swedish Research Council, Vetenskapsrådet. MG acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). VHB acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990.

We report the possible discovery of a new stellar system (YMCA-1), identified during a search for small scale overdensities in the photometric data of the YMCA survey. The object's projected position lies on the periphery of the Large Magellanic Cloud about $13^\circ$ apart from its center. The most likely interpretation of its color-magnitude diagram, as well as of its integrated properties, is that YMCA-1 may be an old and remote star cluster of the Milky Way at a distance of 100 kpc from the Galactic center. If this scenario could be confirmed, then the cluster would be significantly fainter and more compact than most of the known star clusters residing in the extreme outskirts of the Galactic halo, but quite similar to Laevens~3. However, much deeper photometry is needed to firmly establish the actual nature of the cluster and the distance to the system., 5 Pages, 3 figures. Published on Research Notes of the AAS

We assemble a catalogue of Magellanic Cloud red giants from Data Release 2 of the $Gaia$ mission and, utilising machine learning methods, obtain photometric metallicity estimates for them. In doing so, we are able to chemically map the entirety of the Magellanic System at once. Our high resolution maps reveal a plethora of substructure, with the Large Magellanic Cloud (LMC) bar and spiral arm being readily apparent. We uncover a curious spiral-like feature in the southern portion of the LMC disc, hosting relatively metal-rich giants and likely a by-product of historic encounter with the Small Magellanic Cloud (SMC). Modelling the LMC as an inclined thin disc, we find a shallow metallicity gradient of $-0.048 \pm 0.001$ dex/kpc out to $\sim 12^{\circ}$ from the centre of the dwarf. We see evidence that the Small Magellanic Cloud is disrupting, with its outer iso-density contours displaying the S-shape symptomatic of tidal stripping. On studying the proper motions of the SMC giants, we observe a population of them being violently dragged towards the larger Cloud. The perturbed stars predominately lie in front of the SMC, and we interpret that they exist as a tidal tail of the dwarf, trailing in its motion and undergoing severe disruption from the LMC. We find the metallicity structure in the Magellanic Bridge region to be complex, with evidence for a composite nature in this stellar population, consisting of both LMC and SMC debris., MNRAS, submitted

The Fermi-LAT observational data of the diffuse $\gamma$ ray emission from the Large Magellanic Cloud (LMC), were examined to test for the existence of underlying long range correlations. A statistical test applied to the data indicated that the probability that data are random is $\sim 10^{-99}$. Thus we proceeded and have used the counts-number data to compute 2D spatial autocorrelation, power spectrum, and structure function. The most important result of the present study is a clear indication for large scale spatial underlying correlations. This is evident in {\bf all} the functions mentioned above. The 2D power spectrum has a logarithmic slope of -3 on large spatial scales and a logarithmic slope of -4 on small spatial scales. The structure function has logarithmic slopes equaling 1 and 2 for the large and small scales respectively. The logarithmic slopes of the structure function and the power spectrum are consistent. A plausible interpretation of these results is the existence of a large scale {\it compressible turbulence} with a 3D logarithmic slope of -4 extending over the entire extent of the LMC. This may reflect the fact that the $\gamma$ Ray emission is in star forming regions, where jets and shocks are abundant. Both the power spectrum and structure function exhibit steeper logarithmic slopes for smaller spatial scales. This is interpreted as an indication that the turbulent region has an effective depth of about 1.5 kpc., Comment: 8 pages, 6 figures, submitted to AAS journals

CONTEXT. Gaia EDR3 has produced parallaxes for 1.468x10^9 sources but there are calibration issues that require corrections to the published values and uncertainties. AIMS. We want to characterize the behavior of the uncertainties of the Gaia EDR3 parallaxes. We also aim to provide a procedure for the calculation of distances to stars and stellar clusters. METHODS. We reanalyze some of the data in the calibration papers for QSO and LMC parallaxes and combine those results with measurements for six bright GCs. We calculate the angular covariance of EDR3 parallaxes at small separations based on the LMC results and combine it with the results for larger angles using QSOs to obtain an analytical formula for the angular covariance over the whole sky. The results for the six GCs are used to validate the parallax bias correction as a function of magnitude, color, and ecliptic latitude and to determine the constant used to convert internal uncertainties to external ones. RESULTS. The angular covariance at zero separation is 106.2 muas^2, yielding a minimum uncertainty for EDR3 parallaxes of 10.3 muas for individual stars. That value can be only slightly reduced for GCs after considering the behavior of the angular covariance of the parallaxes for small separations. The Lindegren et al. parallax bias correction works quite well, except for the brighter magnitudes, suggesting improvements may be possible there. The value of k is 1.1-1.7 and depends on G. Stars with moderately large values of RUWE can still provide useful parallaxes albeit with larger values of k. We give accurate and precise Gaia EDR3 distances to the six GCs and for the specific case of 47 Tuc we are able to beat the angular covariance limit and derive a high-precision distance of 4.53+-0.06 kpc. Finally, a recipe for the derivation of distances to stars and stellar clusters using Gaia EDR3 parallaxes is given. [ABRIDGED], Version accepted in A&A