On the nature of a shell of young stars in the outskirts of the Small Magellanic Cloud

TL;DR

This study characterizes a shell-like stellar over-density in the outskirts of the Small Magellanic Cloud, revealing it as a young, recently formed structure likely triggered by interactions with neighboring galaxies.

Contribution

The paper provides the first detailed analysis of the shell's stellar content, structure, and kinematics, suggesting a recent star formation event rather than a tidal remnant.

Findings

Shell composed of ~150 Myr old stars with no old population.

Rich sub-structure including spiral arm-like features.

Likely formed from recent star formation triggered by galactic interactions.

Abstract

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

Tables

Table 1. Table 1: Young Star Clusters in the SMC shell. RA and DEC are in J2000 system. References: (1) Glatt et al. 2010, (2) Piatti et al. 2014, (3) Alcaino et al. 2003.

Cluster Name	RA	DEC	log(age/yr)	Ref.
B88	14.233	-70.773	8.10	1
HW33	14.346	-70.809	8.10	2
B139	17.617	-71.561	8.30	1
HW64	17.687	-71.338	8.25	1
IC1655	17.971	-71.331	8.30	1
IC1660	18.158	-71.761	8.20	1
L95	18.687	-71.347	8.30	1
NGC458	18.717	-71.550	8.15	3
HW73	19.108	-71.326	8.15	1

Full text

*Aims. *We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at $\sim 1.9\degr$ from the center of the SMC in the north-east direction.

*Methods. *The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2.

Results. The shell is composed of a young population with an age $\sim 150$ Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute $g$-band magnitude of the shell is M${}_{g,shell}=-10.78\pm 0.02$, with a surface brightness of $\mu_{g,shell}=25.81\pm 0.01$ mag arcsec-2

*Conclusions. *We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, $\sim$ 150 Myr ago.

1 Introduction

The Magellanic Clouds are the largest satellites of the Milky Way (MW), and the only irregular galaxies in its immediate surroundings. The existence of one or even two such gas-rich, massive satellites close to a Milky-Way-sized halo has been shown to be quite rare (e.g., Busha et al. 2011; González et al. 2013; Rodríguez-Puebla et al. 2013; Boylan-Kolchin et al. 2011; Patel, Besla & Sohn 2017). There is strong evidence from the motions and distributions of the stellar and gaseous constituents of the Magellanic Clouds that interactions between them have occurred in the past. High-precision proper motion measurements with the Hubble Space Telescope revealed that, contrary to earlier belief, the Magellanic Clouds are probably just completing their first passage around the Milky Way (Besla et al. 2007, 2010; Kallivayalil et al. 2006a, 2006b, 2013; Piatek et al. 2008; Bekki 2011).

Clearly, the past interactions between the Clouds have left their marks. One of the most obvious manifestations of this interaction is the gaseous Magellanic Stream (Mathewson et al. 1974), whose trailing and leading arms have since been traced over more than $200\degr$ across the sky (Putman et al. 1998; Nidever et al. 2010). Both tidal and ram pressure stripping origins have been suggested for the Stream (see D’Onghia & Fox 2016 for a review). The mass of ionized and atomic Magellanic gas found outside of the Clouds exceeds the remaining H i mass in both Clouds combined (Fox et al. 2014) and is of a similar order of magnitude as the total mass of the Small Magellanic Cloud (SMC).

Another prominent feature is the gaseous Magellanic Bridge, which connects the two Clouds. The Bridge also contains stars, specifically an irregularly distributed young stellar population with ages of up to several hundred Myr (e.g., Irwin et al. 1985; Skowron et al. 2014), associations, and star clusters (e.g., Bica et al. 2008). Unequivocal evidence of intermediate-age stars was also found in the Magellanic Bridge area (Noël et al. 2013; Noël et al. 2015) as well as unambiguous evidence of a tidal stripping scenario (Carrera et al. 2017). The density distribution of these inter-cloud stars suggests that they belong to the extended outer regions of the two Clouds (e.g., Skowron et al. 2014; Jacyszyn-Dobrzeniecka et al. 2017; Wagner-Kaiser & Sarajedini 2017). This would support a scenario where the Bridge formed in a close encounter between the Clouds some 250 Myr ago (D’Onghia & Fox 2016 and references therein). On the other hand, Belokurov et al. (2017) report the discovery of a separate, off-set old stellar bridge between the Clouds. However, Jacyszyn-Dobrzeniecka et al. (2019) found the existence of this old bridge as controversial, based on the analysis of the distribution of RR Lyrae stars in the OGLE data.

Moreover, the outer regions of both Clouds show distortions, clumps, arcs, and related overdensities, both in the direction of the Bridge and elsewhere (e.g., Nidever et al. 2013; Casetti-Dinescu et al. 2014; Besla et al. 2016; Belokurov & Koposov 2016; Mackey et al. 2016; Belokurov et al. 2017; Pieres et al. 2017; Subramanian et al. 2017; Carrera et al. 2017; Mackey et al. 2018; Choi et al. 2018a; Choi et al. 2018b). The last close encounter between the two Clouds some 100 to 300 Myr ago is not only consistent with the young ages in the Bridge, but also with a peak in the age distribution of young clusters in both Clouds (Glatt et al. 2010) and with the bimodal ages of Cepheids (Ripepi et al. 2017).

The SMC has been particularly affected by past interactions and shows a highly distorted, amorphous structure in its H i distribution (e.g., Stanimirovic et al. 2004) and in its younger populations, while its old populations are symmetrically and regularly distributed (e.g. Cioni, Habing & Israel 2000; Zaritsky et al. 2000; Haschke et al. 2012; Jacyszyn-Dobrzeniecka et al. 2017; Muraveva et al. 2018). The SMC has a large line-of-sight depth (e.g., Caldwell & Coulson 1986; Mathewson et al. 1988; Crowl et al. 2001; Subramanian & Subramaniam 2012; Haschke et al. 2012; Nidever et al. 2013) and its old stellar population is very extended (Noël & Gallart 2007; Nidever et al. 2011). From their study of the 3-D structure of the SMC using Cepheid stars, Ripepi et al. (2017) found 25-30 kpc (see also Scowcroft et al. 2016; Jacyszyn-Dobrzeniecka et al. 2016). Apart from repeated disruptive encounters between the Clouds, it has been suggested that the complex structure of the SMC may also be due to a dwarf-dwarf merger in the distant past (Bekki & Chiba 2008).

One important approach towards understanding the evolutionary history of the Magellanic Clouds is through deep, multi-colour mapping of the Clouds, especially of their neglected outskirts, which can still contain clues about the times and nature of past interactions. As part of our deep, wide-field imaging survey of faint tidal structures around the Magellanic Clouds and other Milky Way satellites (Besla et al. 2016; Martínez-Delgado et al. in preparation) using telephoto lenses, we have detected a coherent shell-like over-density feature embedded in the diffuse ”outer arm B” (see de Vaucouleurs & Freeman 1972) and previously visible in photographic studies of the Clouds dating back to the 1950s (e.g. see Figure 12b in de Vaucouleurs & Freeman 1972 and Fig. 1 in this work). Albers et al. (1987) also found evidence of shell-type structures on the north-east side of the SMC based on the analysis of star counts from scanned photographic plates, suggesting they could be associated to faint, spiral arm structure in this region of the SMC. The shell-like feature was also previously noted as “a distinct linear feature perpendicular to the main elongation of the SMC” in a stellar density map using only upper main-sequence stars selected from the Magellanic Cloud Photometric Survey in Zaritsky et al. (2000, their Fig. 2). An early color-magnitude diagram of the north-east outer regions of the SMC by Brueck & Marsoglu (1978) concluded the presence of a young population ($\sim$ 60 Myr) and the association of three star clusters in the outer arm B, that were interpreted as the evidence of a recent burst of star formation in this SMC region (Bruck 1980). In this paper, we study the structure, stellar and gas content, kinematics and possible formation scenarios of this over-density by means of its resolved stellar populations as traced by the imaging data obtained in the Survey of the MAgellanic Stellar History (SMASH, Nidever et al. 2017) and the proper motions available from the Gaia Data Release 2 (DR2) (Gaia Collaboration et al. 2018).

2 The Data

2.1 Canon 200 lens data

In our data, the SMC over-density was first detected in deep images taken with an equipment consisting of a Canon EF 200 mm f/2.8L II USM and a SBIG STL-11000M CCD camera during two observing runs in 2009 August and September at the European Southern Observatory (La Silla, Chile). This setup provided a total field of view (FoV) of $10\degr\times 7\degr$ and a pixel scale of 9.27 arcsec pixel*-1*. A set of 58 individual exposures of 300 sec were taken in a Luminance filter (see, e.g., Fig. 1 in Martínez-Delgado et al. 2015), with a total exposure time of 290 minutes. This first image is showed in Fig. 1 (panel B), with the SMC over-density marked with label d.

Standard data reduction procedures for bias subtraction and flat fielding were carried out using the CCDRED package in the Image Reduction and Analysis Facility (IRAF111IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.). A detailed analysis of these observations is presented in Martínez-Delgado et al. (in preparation).

To ensure that this shell-like feature was not an artifact or a reflection, a confirmation wide-field image of the SMC was taken with a different Canon EF 200 mm f/2.8L lens attached to a Canon E0S6D (ISO1600) camera in 2015 October 10 from Hacienda Los Andes (Chile). Figure  2 shows the resulting image, with a total exposure time of 156 minutes (obtained by combining 78 individual images with an exposure time of two minutes). Astrometry of this image was obtained using the astrometric calibration service Astrometry.net (Lang et al. 2010).

2.2 SMASH data

SMASH (Nidever et al. 2017) was carried out using the Dark Energy Camera (DECam, Flaugher et al. 2015) on the 4-m Blanco Telescope at Cerro Tololo Inter-American Observatory, Chile. The survey, covering 480 deg2 of the sky, aims to explore the extended stellar populations of the Clouds and their interaction history. For our current work we use a small subset of the SMASH observations. The region of the observed shell in the Canon images is well covered by SMASH in Fields 9 (central coordinates (J2000): $\alpha=$01:01:27.40, $\delta=$-70:43:05.51), 14 ($\alpha=$01:20:26.78, $\delta=$-71:15:32.76), and 15 ($\alpha=$01:24:33.52, $\delta=$-72:49:30.00), which are part of a ring of fields surrounding the central part of the SMC. Each DECam field covers a FoV of 3 deg2. Each of these fields was observed in the ugriz bands, with deep observations of 999s in $u$, $i$, and $z$, and 801s in $g$ and $r$. In addition, for each field we obtained three shorter exposures (60s) with large offsets to tie all the chips together photometrically and to allow us to cover some of the gaps between the CCDs. Details on the data processing and photometry can be found in Nidever et al. (2017).

The photometry of the DECam images was carried out using the photometry program DAOPHOT (Stetson 1987) as incorporated into IRAF. In order to identify point sources in the photometry we imposed the following cuts in the DAOPHOT parameters: $-1.0<{\rm SHARP}<1.0$ and $\chi^{2}<3.0$. We also required that the Sextractor (Bertin & Arnouts 1996) stellar probability index was ${\rm PROB}>0.8$. Since the Fields 9, 14, and 15 cover a contiguous part of the sky and overlap with each other, we took care to avoid double entries in the final combined catalog. This leaves us with 2,651,378 stars.

3 A shell in the outer region of the SMC

The images of the SMC obtained with two different Canon 200 lenses (Fig.1B and Fig. 2) revealed a clear shell-like feature situated at $1.9\degr$ from the centre of the galaxy, just above the bar of the SMC (Figure 2). The shape of this over-density is similar to those of the tidal features recently reported in the outskirts of the northern LMC (Mackey et al. 2016; Besla et al. 2016) and the SMC (Pieres et al. 2017)222Unfortunately, the position of this shell in the SMC halo (situated at 8 $\deg$ north from the SMC) is outside of the sky area covered in our photographic survey, including the Canon 50 images described in Besla et al. (2016).. In this section, we explore the spatial extent, structure, and stellar populations by analysing the deep photometry obtained in the SMASH survey.

3.1 Spatial distribution as traced with different stellar populations

Figure 3 shows a $g$ versus $(g-i)$ color-magnitude diagram (CMD) of our final catalog. The features of the SMC are clearly seen in this rich diagram, and it shows the wide range of ages of the stellar populations that are typical for this kind of dwarf irregular galaxy. The main sequence (MS) is very extended, tracing the SMC’s young populations. The sub-giant branch (SGB) and a wide red giant branch (RGB) are produced by older populations with typical ages of 1 Gyr or more. The vertical red clump and the red clump (RC) at $19.2\lesssim g\lesssim 20$ are additional obvious features in the diagram, which trace intermediate-age populations. For a detailed exploration of the star formation history of the SMC from wide-field, multi-colour, ground-based imaging we refer the reader to Rubele et al. (2018) and from very deep HST imaging to Cignoni et al. (2012, 2013).

In a first approximation, we selected stars in the main features described above, which are labeled in Figure 3 as boxes 1 (upper MS), 2 (RC) and 3 (lower part of the RGB). Since all of these features are located in the brighter part of the CMD, the spatial coverage is uniform. This is not necessarily true in the lower part of the CMD since the deepest SMASH exposures were not dithered and hence the gaps among CCDs would be noticeable.

To investigate the nature of the potential shell, we constructed density maps of the region based on stars in different parts of the CMD, as shown in Figure 4. In the left panel we show a map containing stars in the RC and RGB of the SMC ($\sim 100,000$ stars from boxes 2 and 3 of Figure 3). The distribution of these stars, which trace the intermediate-age and old populations of the SMC, is quite uniform in the region, showing only a smoothly increasing density toward the centre of the SMC located in the lower right part of the panel, at $\alpha_{J2000}=00$:52:38, $\delta_{J2000}=-72$:48:01. The globular cluster NGC 362 (top right position) is clearly visible in this map since part of its MS population overlaps with our selection box for the SMC RGB stars.

The smooth distribution of the density of stars in this map suggests that there is no significant extinction in this part of the sky that may cause an apparent shell-like feature. This is also supported by the extinction map of the SMC of Haschke et al. (2011, their fig. 4), which shows generally low extinction in the SMC and values of $E(B-V)\simeq 0.05$ mag in the northeastern part of the SMC. The lack of an over-density in the upper panel of Figure 4 is furthermore consistent with Zaritsky et al. (2000) who also found no evidence of the northeastern over-density in their stellar density maps obtained using giants and red clump stars as tracers (see their Fig. 3).

The region has a very different appearance when plotting only young stars belonging to the upper MS of the SMC population, as selected by Box 1 in Figure 3. In this map, made with $\sim 70,000$ stars, the shell is clearly visible (labeled a in the right panel of Fig 4) and a rich structure is associated with it. Besides the main shell, we find a spiral-arm-like feature of young stars attached to the shell (labeled b in Fig 4) and a separated small arc situated $\sim$ 30$\arcmin$ West from the globular cluster NGC 362 (labeled c in the right panel of Fig 4). This last feature has no counterpart in the old population maps, so we conclude that it is not part of a tidal tail from this cluster (see also Carballo-Bello 2019). Instead, it is very interesting that two of the young open clusters discussed in Sec. 3.3 (see Table 1) are embedded in the Southern extreme of this small structure (see Fig.8). The structures a and b (Figure 4, right panel) are also clearly visible in the GALEX image plotted in Fig. 5, which traces mainly hot young stars, showing an excellent agreement with the position and morphology of these features as traced in our stellar density maps.

Our CMD selection also includes a known streamer of young stars into the Bridge at $\delta\sim$ -73.3 $\deg$, which is also evident in the East side of the SMC in our Canon 200 image in Figure 2. This is also the origin of some negative features imprinted in this region in the left panel of Figure 4 at $\delta\sim$ -73.3 $\deg$. This is due to the incompleteness in these more crowded areas, since old red stars are harder to detect near bright blue star formation regions.

3.2 Stellar content

Fig. 4 indicates that the shell contains a rich population of young MS stars with a range of ages, which is much less prominent in the surrounding areas. To study the characteristics of the stellar populations present in the shell, we will use de-reddened $g_{0}$, $(g-i)_{0}$ CMDs and the corresponding colour functions (CF). Given the spatially variable nature of the foreground Galactic reddening, we use the reddening maps of Schlegel, Finkbeiner, & Davis (1998) and combine these with the revised extinction coefficients of Schlafly & Finkbeiner (2011) to de-redden each star individually. We adopt the recommended “mean” distance modulus of 18.96 to the SMC as advocated by de Grijs & Bono (2015).

Given the areal coverage of the SMASH survey in conjunction with different environmental effects in each field across the survey – such as distinct SFHs, differential reddening, crowding effects, depth, etc. – we opt to compare the de-reddened $g_{0}$, $(g-i)_{0}$ CMD and CFs of the shell region against two different “control” fields: ConF1 is located close to the shell (see left panel in Fig. 4), and ConF2 is on the opposite side of the SMC (not shown in Fig. 4). We aim to investigate the characteristics of the stellar content in each field, and identify differences in the stellar populations between them.

In Fig. 6 we plot the CMD of the stars in the shell box (upper left panel) and in ConF1 (upper right panel). Box boxes have the same area. Isochrones from the BaSTI version 5.0.1 library (Pietrinferni et al. 2004333http:/albinone.oa-teramo.inaf.it) of ages of 30 Myr, 150 Myr, 250 Myr, 800 Myr, 2.5 Gyr, 6.0 Gyr (Z=0.002) and 13.5 Gyr (Z=0.001) have been plotted as reference. A value of Z=0.002 is consistent with that derived from H ii regions, young stars, and Cepheids in the SMC (Russell & Dopita 1992; Romaniello et al. 2009; Lemasle et al. 2017). Although the SMASH survey utilizes the 4-m Blanco telescope and DECam imager, who filter system is described in Abbott et al. (2018), the photometric calibration relies on the use of standard stars in the Sloan Digital Sky Survey (SDSS; see Nidever et al. 2017) and as such we adopt the BaSTI isochrones in the SDSS ugriz system for our comparison.

From the two upper panels of Fig. 6 it seems that the shell sample contains a significantly higher number of luminous MS stars than the same region of the CMD in the “control” field ConF1 sample, particularly in the age range 150-250 Myr. It is interesting, however, that the CMD of ConF1 contains a larger quantity of stars younger that 150 Myr, which are very scarce in the shell sample: from the comparison with the isochrones it can be concluded that very little star formation has taken place in the shell during the last 150 Myr. A step in the density of stars in the main sequence can also be observed, both in the shell and in the ConF1 CMD at the approximate position of the 2.5 Gyr isochrone, possibly indicating an enhanced period of star formation at intermediate ages. This enhancement is well-known and corresponds to a common LMC/SMC burst epoch at about 1.5-3 Gyr ago (e.g. see Harris & Zaritsky 2002; Weisz et al. 2013).

The larger number of young, bright objects in the shell results in a $g$-band surface brightness for the shell sample that is more than 0.5 mag arcsec*-2* brighter than the ConF1 control field sample (cf. $\mu_{g,shell}=25.81\pm 0.01$ and $\mu_{g,CF1}=26.68\pm 0.01$ mag arcsec*-2*). We also notice this increased surface brightness in the $i$-band but with a smaller difference (cf. $\mu_{i,shell}=25.55\pm 0.01$ and $\mu_{i,cont}=26.12\pm 0.01$ mag arcsec*-2*). For the “mean” SMC distance modulus of 18.96, we determine that the absolute $g$- and $i$-band magnitudes of the shell sample are M${}_{g,shell}=-10.78\pm 0.02$ and M${}_{i,shell}=-11.05\pm 0.02$ mag, respectively.

The lower left panel of the Fig. 6 shows the CMD of a synthetic population computed using the BasTI on-line Stellar Population Synthesis Program444http://basti.oa-teramo.inaf.it/BASTI/WEB_TOOLS/synth_pop2/, using solar scaled overshooting models. We have assumed a constant star formation rate (SFR) from 13.5 Gyr to 30 Myr ago (the latter is the young age limit for the BasTI library), and a simplified chemical enrichment law, approximately consistent with that obtained by Carrera et al. (2008) from Ca II triplet spectroscopy, that is, [Fe/H]=-0.99 ($\sigma$=0.3) for the second half of the galaxy’s life (age $<$6.75 Gyr), and [Fe/H]=-1.29 ($\sigma$=0.1) for the first half (age$>$6.75 Gyr). Other parameters of the model are a binary fraction of $\beta$=0.4 with a mass ratio q$>$0.5 and a Kroupa et al. (1993) IMF. Different colours for the synthetic stars have been used to highlight the position in the CMD of stars with different ages. The same isochrones as in the observed CMDs were plotted as reference. Finally, the lower right panel of Fig. 6 shows the luminosity function of the shell, control fields, and that of the synthetic CMDs. These luminosity functions allow us to conclude that the observed CMDs are basically complete down to $M_{i}$=3 or $i_{0}$=21.96. The horizontal lines indicate this magnitude limit, which will be used to compute the CFs.

In Fig. 7, we use the CFs to further analyze the stellar population content of the shell and the differences to the SMC field populations at similar galactocentric radius (for an introduction to the use of the CF for stellar population analysis, see Gallart et al. 2005; see also Noël et al. 2007 for an application to study SMC field stellar populations). We compare the shell CF with the CF of the same control field shown in Fig. 4 (ConF1) and control field ConF2 located at the opposite side of the SMC. The upper panel of Fig. 7 shows, in black and different line types, the $(g-i)_{0}$ observed CFs for the shell and the two control fields mentioned above. An absolute magnitude cut of $M_{i}$=3 has been adopted to ensure a high level of completeness in the photometric data used to calculate the CFs (see Fig. 6). The orange solid line shows the CF of the synthetic CMD with a constant star formation rate at all times. The CF of the synthetic CMD was scaled such that the number of stars in a box on its red clump matches that in a box in the same location in the shell CMD.

The lower panel of Figure 7 displays the same observed CFs, together with CFs of the synthetic population in three limited age ranges: 0.03 to 0.8 Gyr, 1.5 to 6 Gyr, and older than 8.0 Gyr. In this case, the normalization of the synthetic CFs is arbitrary, and chosen to approximately match the features present in the observed CF of the shell. The position and width of the blue maximum in the shell CF (approximately centered at $(g-i)_{0}\simeq-0.5$) suggests a very recent star formation event. The second maximum approximately centered at $(g-i)_{0}\simeq 0$ is well reproduced by a population spanning the age range 1.5 to 6 Gyr. The third maximum corresponds to the position of the red clump. These results on the age range that is contributing to each feature of the CF are in good agreement with the conclusions we reached from the comparison of the CMDs with isochrones. The comparison of the observed CFs indicates significant differences in the stellar content of the shell and the two control fields. The very prominent blue maximum in the shell CF has very low counts in both ConF1 and ConF2. The difference in the height of the second, intermediate colour maximum in the shell and the control fields is not striking, but a difference is nonetheless evident. Taking into account the information on the lower panel of Fig. 7 regarding the ages contributing to each feature on the CF (complemented by the isochrone information), we conclude that star formation has been very active in the shell in the last $\simeq$ 1 Gyr and much less active in the control fields. In the intermediate-age range 1.5-6 Gyr ago, star formation was also more active in the shell field than in the control fields. The comparison with the CF of a population with a constant star formation rate (orange line in the upper panel of Fig. 7) discloses that star formation in the shell field has been enhanced (compared to an average constant SFH) in the second half of the galaxy’s life, while star formation in both control fields has been depressed in the last $\simeq$ 1 Gyr compared to a constant star formation rate.

Our analysis of the total stellar content toward the shell field, through the comparison of the CMD with isochrones and the comparison between observed and synthetic CFs, reveals a very complex stellar population, with stars of all ages contributing to the CMD of the shell field. A strong enhancement of the recent star formation rate is, however, indicated by the height of the bluest peak of the CF. Even though a broad age range seems to be necessary to reproduce the width of that peak, the details of this age range and the precise age composition cannot be accurately constrained by this kind of simple analysis. A full SFH derivation, which is beyond the scope of this paper, would be necessary. Additional considerations brought up by the analysis of the young cluster (see Sec.3.3) and Cepheid population (see Sec. 3.4) in that area of the SMC indicate that a conspicuous burst of star formation around 150-250 Myr contributed outstandingly to the young population of the shell and may be the origin of the features that identify the shell morphology and that are visible in the right panel of Figure 4.

3.3 Young Star Clusters

A number of young star clusters lie spatially coincident with the northeastern shell feature, with the most prominent clusters appearing as visible over-densities in the upper MS star count map in Fig. 4. We identify nine young ($<$1 Gyr) star clusters from the Bica et al. (2008) catalog that lie within a $\sim 2.2$ square degree search region. We display their location in Figure 8 (upper panel), and list their properties in Table 1. For this census, we exclude six older (age$>$1 Gyr) clusters from further discussion, and ignore 5-10 diffuse associations and low significance catalog entries. This sample features three relatively massive clusters ($\sim 10^{4}M_{\sun}$; NGC458, IC1655, IC1660), as well as six other less massive systems555 The association of three clusters with the position of this feature was previously mentioned by Brueck & Marsoglou (1978). However, only NGC 458 and L90 actually fall on the shell, while HW62 is clearly offset off the stellar over-density towards the SMC center. The ages assigned to the clusters by these authors, based on photographic plates, have been revisited with our modern observations yielding 160 Myr (instead of 20-50 Myr)..

A striking feature of this young cluster sample is their uniformity in age. As determined from isochrone comparisons to observed CMDs (Alcaino et al. 2003; Glatt, Grebel & Koch 2010; Piatti 2014), the cluster ages pile up around $\sim$ 160 Myr [log(age/yr)=8.2$\pm$0.1] and appear consistent with this single age, within current fitting uncertainties. We demonstrate the agreement with a single age by over-plotting an isochrone with log(age/yr)=8.2 on top of a summed SMASH CMD created for the nine cluster sample (see lower panel in Fig. 8) . This synchronization in cluster ages suggests (or, further reinforces the conclusion) that an important enhancement of the star formation rate took place $\sim$ 160 Myr ago. This epoch broadly agrees with the SMC-wide peak of cluster ages observed by Glatt et al. (2010) as well as the putative age of the most recent LMC-SMC interaction.

3.4 Cepheids

OGLE IV (Udalsky et al. 2015; Soszyński et al. 2016) uncovered almost 5,000 classical Cepheids in the SMC and several tens of them spatially overlap with the shell feature (see Fig. 9). Since Classical Cepheids are young supergiants, their presence in quite large numbers underlines significant star formation in the last few hundred Myr.

Cepheid ages can be computed for individual stars using period-age relations derived from population models (for instance, Bono et al. 2005). Using these relations, the ages vary from $\approx$15 to $\approx$500 Myr666Models that include rotation during the Main Sequence phase lead to Cepheid ages increased by 50 to 100%, depending on the period (Anderson et al. 2016) for SMC Cepheids. The age distribution of SMC Cepheids is known to be bimodal, with two peaks at $\approx$110–130 Myr and $\approx$220–230 Myr separated by a minimum at $\approx$150 Myr (e.g., Inno et al. 2015; Subramanian et al. 2015; Jacyszyn-Dobrzeniecka et al. 2016; Ripepi et al. 2017). The former peak has been associated with star formation triggered by the most recent interaction between the LMC and the SMC. Cepheids lying in the shell region span a relatively wide age range that clearly peaks at 100–130 Myr (see Fig. 10). Such ages are in good agreement with current results for the MS sample and match very well the ages derived for the young clusters in the vicinity of the shell777However, it is important to take into account that the raw age distribution of Cepheids cannot be directly interpreted as an age distribution of the underlying star formation because it is convolved with both the stellar IMF and the lifetime of the star within the instability strip during which it would be identified as a Cepheid variable.. Cepheids therefore support a scenario where the shell population is dominated by young stars formed during a recent star formation event, possibly related to the interaction between both Magellanic Clouds.

3.5 Stellar Kinematics

Using the Gaia DR2 data (Gaia Collaboration et al. 2018), proper motions for thousands of stars in the SMC outer region have recently became available. We use this catalog to investigate the kinematic signature from the shell-like feature. From the DR2 database, we select stars surrounding the SMC and apply a series of astrometric cuts. We start by applying a parallax cut of $\omega<0.2$ mas in order to remove foreground MW stars. Next we use a cut to the renormalized unit weight error (as described in the gaia technical note GAIA-C3-TN-LU-LL-124-01) of 1.40 and a cut for the color excess of the stars (as described in Lindegren et al. 2018 by Equation C.2). As astrometric precision has a strong relationship with the magnitude of the stars, we additionally apply a cut of $\mathrm{G}<18$. Finally, to trim down potential MW contamination, we cut out an area with radius equal to 3 mas yr*-1* around the systemic proper motion (PM) of the SMC, which we will take to be $\mu_{\alpha^{*},0,\mathrm{SMC}},\mu_{\delta^{*},0,\mathrm{SMC}}$ = $0.797\pm 0.030,-1.220\pm 0.030$ mas yr*-1* (Helmi, et al. 2018). This leaves us our final selection of stars (seen in upper Figure 11), where the shell-like over-density can be clearly seen (marked by the purple rectangle).

The CMD of the shell region obtained from the Gaia data (bottom panel of Figure 11, expanded down to G $=20$ for greater context) displays clearly the MS, the red supergiants (RSG), the RC, and the RGB features. Combining the requirement of G $<18$ and the apparent locations of the stellar features, we created masks for each feature. The masks were then applied to the full sample of SMC stars and the resulting spatial distributions were plotted (upper Fig. 12). Similar to the analysis earlier in the paper, the RGB had no apparent correlation with the shell-like region, but both the MS and RSG display over-densities in the location of the shell. Correspondingly, we select these two sequences from the shell region for our kinematic analysis.

For this portion of the analysis, we first subtract the systemic PM from the PM of each star in the shell region. As the feature is noticeably extended on the night sky, we also calculate and remove the viewing perspective for each source, as outlined in van der Marel et al. (2002). These residual proper motions are converted into a Cartesian frame assuming the kinematically-derived center of ($\alpha_{J2000}$, $\delta_{J2000}$) = (16.25∘,$-$72.42∘), using the transformations from the Gaia Collaboration (Helmi et al. 2018). The position angle (PA, which in the context of this analysis is defined as the vector angle of the residual proper motions) is calculated for each source, defined where $0^{\circ}$ points vertically upwards in the spatial plot in Figure 12 and increases in a counter-clockwise direction. We compile all of the PAs for all sources into a histogram for easier interpretation (seen in lower panel of Figure 12). The PAs appear to peak around $70-80^{\circ}$, which points roughly radially outwards from the center of the SMC. The scatter from $0^{\circ}$ to $360^{\circ}$ is expected as the average residual is on the order of $\sim 0.1$ mas yr*-1* with average errors of comparable magnitude, underscoring that the large peak in PA must be a real signal. This coherent radially outward motion of stars in the SMC is consistent with prior studies of the internal kinematics of the SMC (e.g. Zivick et al. 2018).

3.6 H i and H$\alpha$ emission

Using the atomic hydrogen (H i) emission maps from the Galactic All-Sky Survey (GASS) third release888https://www.astro.uni-bonn.de/hisurvey/gass/ (Kalberla & Haud 2015), we scanned the velocity channels available in that survey (from $-495$ km s*-1* out to +495 km s*-1*). The H i gas emission within one square degree centered on $\alpha=1$h10min, $\delta=-71.5{\degr})$ reaches a maximum in the channel $v=198$ km s*-1*. more specifically, the velocities of the gas within that region are close to a Gaussian distribution, with the maximum in $v=198.3$ km s*-1* and a dispersion $\sigma=13.6$ km s*-1*. This velocity is very similar to the mean velocity of the stars measured by Evans & Howarth (2008)(172 km s*-1*) with a dispersion of 30 km s*-1*. The leftmost panel of Fig. 13 shows the emission for this specific velocity channel in a log scale. The second panel of Fig. 13 is the distribution of the mean gas temperature bounded by a larger range of velocities (184-211 km s*-1*) in a linear scale, which is very similar to the previous panel. The rightmost panel shows the gas in the respective velocity channel, in a zoomed region close to the shell discussed in this paper. The young MS stars are over-plotted on the gas density map as white dots. A visual inspection on the last panel suggests a shift between the projected location of the young stars and the gas distribution.

As we do not have constraints to the distance of the H I gas cloud in the Fig. 13, our few arguments are similarities between the velocities and shape of the H I gas cloud and the shell-like feature of young stars in the Fig. 4. Assuming that the ‘Z’ shaped gas cloud was forming stars in a recent past (150-250 Myr ago), the velocities of the young stars and of the gas cloud could provide insights on what is the dominant process (tidal stripping or ram pressure) for decoupling the kinematics of the stars and the gas, at least in that region. However, if the gas cloud formed the shell-like feature, it is not clear what mechanism triggered such an intense star formation episode $\sim$200 Myr ago and why it is no longer acting to form stars. An evidence of the quiescence of the gas in forming stars is the lack of MS stars between the current position of the gas cloud and the position of the shell feature. Nevertheless, we must take into account the possibility of large distances in the line-of-sight between the gas feature and the young stars, being two completely separate substructures.

A comparison of the stellar populations map to the H$\alpha$ map (Fig. 14) from the Magellanic Clouds Emission Line Survey (MCELS; Winkler, Rathore, & Smith 1999) shows a system of filaments that roughly forms a shell with a radius of $0.75\degr$ centered at $\alpha_{J2000}=01$:06, $\delta_{J2000}=-71$:45. This interstellar material is reminiscent of supergiant shells similar to those catalogued in the Magellanic Clouds by Meaburn (1980). Such shells are believed to be driven by the collective action of winds from multiple OB associations and supernovae. When the H$\alpha$ images are compared to the H i ATCA+Parkes observations, there appear to be faint (3-5$\sigma$ significance) H i emission patches coincident with the H$\alpha$ filaments with velocities between $\sim$125-140. We note that this putative shell does not coincide with those originally identified by Stanimirovic et al. (1999) using the same data set.

If this emission corresponds to the limb-brightened edge of a supergiant shell, the lack of detection of a coherent structure including envelopes expanding along the line-of-sight preclude a detailed kinematic comparison with the underlying stellar population. The velocities of the faint H i emission are consistent with the gas (and driving stellar population) being part of the SMC main body and not the material being drawn out by the interaction with the LMC. Based on the size, 750 pc (0.75°at 60 kpc) and assuming this interstellar structure arose from the collective stellar energy feedback from a population originating 150 Myr ago, we expect an expansion velocity of 5 km s*-1* which is low but not unreasonable compared to other Magellanic supergiant shells (15 to 30 km s*-1*, Book et al. 2008).

4 Discussion

Our analysis of the resolved stellar populations of the elongated over-density detected in our Canon 200 images reveals that it is the brighter optical part of a more extended structure mainly traced by blue, young stars distributed in an intricate and complex structure. The structure shows a further, outer arc-like feature observed in the projected proximity of the Galactic globular cluster NGC 362 (see bottom panel of Figure 4). The over-density also contains nine young star clusters with ages tightly clustered at 175 Myr (see Table 1).

We do not detect any counterpart of this over-density in the distribution of the older stars in this area of the SMC. This dominant young age of our SMC over-density and the lack of an old stellar remnant in the stellar density map plotted in Fig. 4 (left panel) suggest that this feature is neither of tidal origin nor the stellar remnant of a tidally disrupted, lower-mass system as observed in some other dwarf galaxies, e.g., in NGC 4449 (Martínez-Delgado et al. 2012) or in Andromeda II (Amorisco et al. 2014). While recently a large number of new, ultra-faint dwarf galaxies were discovered near the Magellanic System (Drlica-Wagner et al. 2015, 2016), these systems are exclusively devoid of gas and contain very old stellar populations. The existence of a former, more gas-rich small dwarf irregular galaxy is not excluded, but we would expect to also see a clear over-density in the old stellar population map if such a system were to merge. With the exception of tidal dwarf galaxies, there are no nearby dwarf galaxies known that do not contain old populations (Grebel & Gallagher 2004). Thus, the nature of this over-density seems to be different from the one discovered at $8\degr$ north of the SMC by Pieres et al. (2017), which is mainly composed of intermediate-age stars and without significant H i gas. Although morphologically similar in appearance to the LMC substructure described by Mackey et al. (2016), the young age of the SMC feature distinguishes it from that older and much larger LMC over-density.

Our stellar density maps of the shell region (Figure 4 and Figure 11) and the GALEX data (Figure 5) confirm with higher resolution the asymmetric structure of young stars originally found by Zaritsky et al. (2000) in the SMC outskirts, which contrasts with the smooth distribution of the older stellar populations. This lack of substructure in the older stellar populations led these authors to conclude that the dominant physical mechanism in determining the current appearance of the SMC must be recent star formation possibly triggered by a hydrodynamic interaction of the SMC with a second gas-rich object (e.g., a hot Milky Way halo or the outer gaseous envelope of the LMC). The contrast between the small-scale spatial structures in the young populations and the smooth distribution of the intermediate-age and old populations are even starker in these new data. This finding strengthens our argument (see also Zaritsky et al. 2000) that a purely tidal origin for the feature, which would have affected stars of all ages similarly, is unlikely. Instead, the origin of the feature is most likely hydrodynamical in nature. Given the coherence of the feature, the absence of an old stellar remnant, and the lack of evidence for a strong shock in the H i distribution, we disfavour models involving a recent accretion event.

The elongated shape, estimated absolute magnitude (Sec. 3.2), and position of the over-density in the outer region of the main body of the SMC could also suggest a tidal dwarf galaxy as another possible formation scenario for the shell-like structure. The stars in such objects are formed from gas stripped in past encounters (Elmegreen, Kaufman & Thomasson 1993) and/or consist of stars that originally formed in the more massive galaxies participating in the interaction (e.g., Duc 2012). Key requirements for tidal dwarfs include that they are made from recycled material, are gravitationally bound, and have decoupled from their former parent (Duc 2012). They may be able to survive for at least 3 Gyr in spite of their lack of dark matter (Ploeckinger et al. 2014). In our case, the over-density appears to be connected with the SMC and there is no evidence to suggest that it is a separate entity. Hence we also discount the possibility of a tidal dwarf galaxy.

A conspicuous feature in the stellar density map is a spiral-arm-like feature emanating from the over-density (labelled b in Figure 4). Because this stellar arm is only seen in young stars, it is unlikely that it was produced via tidal effects. Instead, the arc may consist of young stars formed as a result of a low pitch angle spiral density wave in the outer gas. Such features are common, but challenging to detect, in disks well beyond the optical radius (Ferguson et al. 1998; Herbert-Fort et al. 2012). The presence of a disk population would imply that the young stellar populations of the SMC must have some bulk rotation. However, initial Gaia DR2-based proper motions and rotation measurements of the SMC do not support significant rotation for the young SMC population (van der Marel & Sahlmann 2016), suggesting that, if the density wave scenario is correct, then the observed offset is dominated by the progression of the pattern rather than motion of material through it. In addition, the hypothesis of a density wave also imply an expected offset between current (H$\alpha$) and past star formation. We find not detectable offset between the H$\alpha$ emission and GALEX sources (see Fig. 14) which, given the timescales these tracers are sensitive to (107 vs. 108 yrs), places an upper limit on the pattern speed. Within the new Gaia DR2 catalog, we find among the bright stars the same shell feature in the young stellar populations. Examining this subset of SMC stars reveals coherent motion among the stars but in an outward radial direction. Given the direction of motion, it continues to appear unlikely that the SMC possesses stellar rotation, at least rotation that lays primarily in the plane of the sky as with the neighboring LMC. Though, as the vast majority of these stars do not possess radial velocities, we are unable to fully assess if this motion may be due to rotation with a large inclination. A more holistic analysis of the SMC, including searching for other kinematic substructure, will be required to improve constraints on the degree of stellar rotation.

Our analysis of the stellar content of the shell (see Sec. 3.2) indicates that the main difference between the stellar population of the shell and the surrounding control fields is a recent period of enhanced star formation in the last $\simeq$ 1 Gyr, likely peaking at $\sim$ 150 Myr, as indicated by the clusters and Cepheids age distribution. The comparison of the CFs of the shell and the control fields also provides hints of a comparatively enhanced star formation rate in the total stellar content of this region at intermediate ages. If we compare these hints on the SFH of the shell and the control fields with the SFHs of the fields analyzed by Noël et al. (2009), it can be seen that the SFH of the shell field may present similarities with the SFHs of fields located in the Wing area of the SMC such as qj0112, qj0111, and qj0116, while the SFH of the two control fields may be alike to that of the remaining fields analyzed by Noël et al. (2009), where SFH has been very low in the last $\simeq$ 1 Gyr. The qualitative hints on the SFH of the shell region obtained through the CF analysis are also consistent with the SFH for fields 1 and 4 (located in the bar and the wing area, respectively) analyzed by Cignoni et al. (2012), which present a strong enhancement from $\simeq$ 5 Gyr ago to the present time, and a strong peak at a very recent epoch, $<$ 200 Myr ago, and with the information on detailed SFH maps presented by Rubele et al. (2015), which indicate that the Wing is $<$ 200 Myr old, that metal-poor gas was injected $>$ 1 Gyr ago resulting in the formation of intermediate-age stars, and that the majority of the SMC mass resulted from a star formation episode 5 Gyr ago. The SFH of the Wing field, then, could be then similar to that of the SMC central body and Wing area, but characterized by an enhanced star formation rate in the last 200 Myr with respect to the past average star formation rate. However, since the CMD toward the shell region contains also a component of the field SMC population, it is possible that the intermediate-age populations that make up the CF peak at ($g-i)\sim$ 0 in Fig. 6 are part of that SMC underlying population. Therefore, there are two possible scenarios to explain the nature of the shell: i) a shell composed only of young stars resulting from a huge star forming region active $\sim$ 150 Myr ago; or ii) the shell is a region that had enhanced episodes of star formation at different epochs, with the most striking one happening $\sim$ 150 Myr ago.

It is, thus, interesting to consider the young age of our SMC substructure in the context of the putative recent interaction with the LMC about 100 to 300 Myr ago, which left its signature in the young stellar populations in the Magellanic Bridge (e.g., Skowron et al. 2014) and in the age distribution of young populous star clusters in both Clouds (e.g., Glatt et al. 2010). Recently, Zivick et al. (2018) used mostly Hubble Space Telescope proper motion data to show that the LMC and SMC have had a head-on collision in the recent past, and it is able to constrain both the timescale and the impact parameter of this collision. Based on their measured proper motions and considering the allowed range of masses for the LMC and the Milky Way, they foud that in 97% of all the considered cases, the Clouds experienced a direct collision with each other 147 $\pm$ 33 Myr ago, with a mean impact parameter of 7.5 $\pm$ 2.5 kpc. There is also further evidence of this recent head-on collision in the Gaia DR2 proper motions along the Bridge (Zivick et al. 2019). The age of our feature falls into this age range as well. Moreover, the three-dimensional structure of young populations traced by Cepheids (with ages of about 15–500 Myr) shows a highly asymmetric distribution throughout the SMC: the distribution of classical Cepheids is elongated over 15–20 kpc (e.g., Scowcroft et al. 2016; Jacyszyn-Dobrzeniecka et al. 2016), the NE region being younger, closer to the Sun than the SW region (e.g., Haschke et al. 2012; Subramanian et al. 2015; Ripepi et al. 2017). Haschke et al. (2012) speculate that the displacement and compression of the gas of the SMC through tidal and ram pressure effects caused by the interaction between the Magellanic Clouds (and with the Milky Way) may have locally enhanced star formation, creating some of the irregular, asymmetric features such as the over-density described in the current paper. Among others, Inno et al.(2015) and Jacyszyn-Dobrzeniecka et al.(2016) also propose that the concomitance of the extensive Cepheid formation in the LMC $\approx$140 Myr ago and of the younger episode of Cepheid formation in the SMC may be related to the interaction between the Clouds. Ripepi et al. (2017) elaborate on this scenario and suggest that the relatively young Cepheids ($<$140 Myr) that dominate in the NE region have formed after the dynamical interaction that created the Bridge, from gas already shifted by the interaction.

Unfortunately, there is insufficient resolution in any existing numerical simulations of the LMC-SMC interaction to see the fine structure of our stellar density map in Fig. 4. Our results clearly motivate more detailed studies of the internal structure and star formation induced in the SMC by the LMC-SMC (and the MW) interaction. In particular, if the stellar arc is a spiral structure, these data likely disfavor models where the SMC is originally modeled as a non-rotating spheroid. This result further highlights the discrepant kinematics and spatial distribution of the SMC younger and older stellar populations, the origin of which is currently unknown.

Acknowledgements.

We thank David Hogg for his help with the astrometry solution of the image used in this work. DMD thanks Prof. Ken Freeman for a discussion about the detection of the SMC shell in the photographic plates of the Clouds taken in the 1950s. DMD also thanks the hospitality and fruitful discussion about this work with the European Southern Observatory Garching headquarters staff during his stay as part of the ESO visitor program in September 2017. DMD, EKG, BL and LI acknowledge support by Sonderforschungsbereich (SFB) 881 “The Milky Way System” of the German Research Foundation (DFG), particularly through sub-projects A2, A3 and A5. DMD acknowledges support from the Spanish MINECO grant AYA2016-81065-C2-2. This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) collaboration. M-RC and CB acknowledge support from the European Research Council (ERC) under the European Union*′*s Horizon 2020 research and innovation program (grant agreement No. 682115). BCC acknowledges the support of the Australian Research Council through Discovery project DP150100862. T.d.B. acknowledges support from the European Research Council (ERC StG-335936). This work has been supported by the Spanish Ministry of Economy and Competitiveness (MINECO) under grant AYA2014-56795-P. MS acknowledges support from the ADAP grant NNX14AF81G. R.R.M. acknowledges partial support from project BASAL AFB-$170002$ as well as FONDECYT project N${}^{\circ}1170364$.Y.C. acknowledges support from NSF grant AST 1655677. Based on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. 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 Science and Technology Facilities Council of the United Kingdom, 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, Fundaçã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 Dark Energy Survey. This work has made use of data from the European Space Agency (ESA) mission Gaia (<https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, <https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.