# Tycho's Supernova: The View from Gaia

Pilar Ruiz-Lapuente, Jonay González Hernández, Roger Mor, Mercè Romero-Gómez, N. Miret-Roig, Francesca Figueras, Luigi Bedin, Ramon Canal, Javier Méndez

► To cite this version:

Pilar Ruiz-Lapuente, Jonay González Hernández, Roger Mor, Mercè Romero-Gómez, N. Miret-Roig, et al.. Tycho's Supernova: The View from Gaia. The Astrophysical Journal, 2019, 870 (2), pp.135. 10.3847/1538-4357/aaf1c1 . hal-02005347

HAL Id: hal-02005347

https://hal.science/hal-02005347v1

Submitted on 29 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Tycho's Supernova: The View from Gaia

Pilar Ruiz-Lapuente<fnref n="1" />,<fnref n="2" /> ID, Jonay I. González Hernández<fnref n="3" />,<fnref n="4" />, Roger Mor<fnref n="2" />, Mercè Romero-Gómez<fnref n="2" />, Núria Miret-Roig<fnref n="2" />,<fnref n="5" />, Francesca Figueras<fnref n="2" />, Luigi R. Bedin<fnref n="6" />, Ramon Canal<fnref n="2" />, and Javier Méndez<fnref n="7" />

Received 2018 July 9; revised 2018 October 26; accepted 2018 November 15; published 2019 January 16

## Abstract

SN 1572 (Tycho Brahe's supernova) clearly belongs to the Ia (thermonuclear) type. It was produced by the explosion of a white dwarf (WD) in a binary system. Its remnant has been the first of this type to be explored in search of a possible surviving companion, the mass donor that brought the WD to the point of explosion. A high peculiar motion with respect to the stars at the same location in the Galaxy, mainly due to the orbital velocity at the time of the explosion, is a basic criterion for the detection of such companions. Radial velocities from the spectra of the stars close to the geometrical center of Tycho's supernova remnant, plus proper motions of the same stars, obtained by astrometry with the *Hubble Space Telescope (HST)*, have been used so far. In addition, a detailed chemical analysis of the atmospheres of a sample of candidate stars had been made. However, the distances to the stars remained uncertain. Now, the Second *Gaia* Data Release (DR2) provides unprecedented accurate distances, and new proper motions for the stars can be compared with those obtained from the *HST*. We consider the Galactic orbits that the candidate stars to the SN companion would have in the future. We do this to explore any kinematic peculiarity. We also locate a representative sample of candidate stars in the Toomre diagram. Using the new data, we re-evaluate the status of the candidates suggested thus far, as well as the larger sample of the stars seen in the central region of the remnant.

**Key words:** supernovae: general – supernovae: individual (SN 1572)

## 1. Introduction

Type Ia supernovae (SNe Ia) are the calibrated standard candles used in the discovery of the accelerated expansion of the universe (Riess et al. 1998; Perlmutter et al. 1999) and they remain a powerful tool in exploring the nature of dark energy. Although a lot of progress has been made in disentangling the nature of the explosions, there are still many points to be addressed concerning the progenitors (see, e.g., the reviews by Wang & Han 2012; Maoz et al. 2014; Ruiz-Lapuente 2014). They appear to be thermonucler explosions of white dwarfs (WDs) made of C+O, and accretion of material by the WD from a companion in a close binary system should be the basic mechanism to induce the explosion, but here the consensus stops. The companion could either be a still thermonuclearly active star in any stage of its evolution (the single-degenerate, SD channel) or another WD (the double-degenerate, DD channel). The explosion could also result from the merging of a WD with the electron-degenerate core of an asymptotic giant branch star. The mode of the accretion could range from steady accretion to violent merger, and the explosion either arise from central ignition of C, when the WD grows close to the Chandrasekhar mass, or be induced by detonation of a He layer near the surface, the mass of the WD being smaller in this case. Observed different SNe Ia may have different origins.

No binary system has ever been discovered in which an SN Ia has later taken place, but some binary systems are considered to be excellent candidates for SN Ia progenitors, such as U Sco, which contains a WD already close to the Chandrasekhar mass. A general prediction for the SD channel is that the companion

star of the WD should survive the explosion and present revealing characteristics.

There are supernova remnants (SNRs) of the explosions of SNe Ia, close and recent enough that their exploration can either detect the presence of a surviving companion or confirm its absence (Ruiz-Lapuente 1997). This has been done for several SNRs of the Ia type, in our own Galaxy and in the Large Magellanic Cloud (Ruiz-Lapuente et al. 2004 (hereafter RL04), 2018; González Hernández et al. 2009 (hereafter GH09), 2012; Kerzendorf et al. 2009, 2012, 2013 (hereafter K13), 2014, 2018; Edwards et al. 2012; Schaefer & Pagnotta 2012; Bedin et al. 2014 (hereafter B14); Pagnotta & Schaefer 2015).

The remnant of SN 1572 (Tycho Brahe's SN) was the first to be explored (RL04), and the findings there were later the subject of several studies (GH09; K13; B14).

Now the *Gaia* Data Release 2 is providing an unprecedented view of the kinematics of the Galactic disk (Brown et al. 2018). It not only gives the 3D location of a very large sample of stars in the Galaxy, but also full velocity information (proper motion and radial velocity) for 7.2 million stars brighter than $G_{RVS}$ = 12 mag, and transverse velocity for an unprecedentedly large number of stars. *Gaia* DR2 provides astrometric parameters (positions, parallaxes, and proper motions) for 1.3 billion sources. The median uncertainty for sources brighter than $G = 14$ mag is 0.03 mas for the parallax and 0.07 mas yr$^{-1}$ for the proper motions. The reference frame is aligned with the International Celestial Reference System (ICRS) and non-rotating with respect to the quasars to within 0.1 mas yr$^{-1}$. The systematics are below 0.1 mas and the

<footnote n="1">¹Instituto de Física Fundamental, Consejo Superior de Investigaciones Científicas, c/. Serrano 121, E-28006, Madrid, Spain</footnote>

<footnote n="2">²Institut de Ciències del Cosmos (UB-IEEC), c/. Martí i Franques 1, E-08028 Barcelon, Spain</footnote>

<footnote n="3">³Instituto de Astrofísica de Canarias, E-38206 La Laguna, Tenerife, Spain</footnote>

<footnote n="4">⁴Universidad de la Laguna, Departamento de Astrofísica. E-38206, La Laguna, Tenerife, Spain</footnote>

<footnote n="5">⁵Laboratoire d'Astrophysique de Bordeaux, Université de Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, F-33615 Pessac, France</footnote>

<footnote n="6">⁶INAF, Osservatorio Astronomico di Padova, Via dell'Osservatorio 3, I-35122 Padova, Italy</footnote>

<footnote n="7">⁷Isaac Newton Group of Telescopes, P.O. Box 321, E-38700 Santa Cruz de la Palma, Spain</footnote>
parallax zeropoint uncertainty is small, about 0.03 mas (Brown et al. 2018).

Previously, the distances to the stars could only be estimated from comparison of the absolute magnitudes deduced from their spectral types and luminosity classes with their photometry, assuming some interstellar extinction in the direction of the SNR. That left considerable uncertainty in many cases (see RL04 and B14). It is here where the Gaia DR2 is most useful.

The situation was better concerning proper motions, where Hubble Space Telescope (HST) astrometry, based on images taken at different epochs, had allowed high precision (see B14). HST proper motions are always relative to a local frame, whereas Gaia DR2 proper motions are absolute, referred to the ICRS. Moreover, Gaia DR2 allows us to calculate the Galactic orbits of the stars. In addition, without a precise knowledge of the distances, the conversion of proper motions into tangential velocities remained uncertain and so was the reconstruction of the total velocities.

The paper is organized as follows. First we describe the characteristics of Tycho's SNR. In Section 3, we examine the distances given by the parallaxes from Gaia, for the surveyed stars, and we compare them with previous estimates. In Section 4, the proper motions from Gaia are compared with those from the HST. Section 5 discusses the position in the Toomre diagram of possible companion stars to SN 1572, as compared with a large sample. In Section 6, we calculate the Galactic orbits of four representative stars and discuss their characteristics. In Section 7 our whole sample is discussed. Section 8 compares the observations with the predictions of models of the evolution of SN Ia companions. Finally, Section 9 gives a summary and conclusions.

## 2. Tycho's SNR

Tycho's SNR lies close to the Galactic plane ($b = 1°4$, which means 59–78 pc above the Galactic plane). The remnant has an angular radius of 4 arcmin. In RL04 a search was performed covering the innermost 0.65 arcmin radius centered on the Chandra X-ray Observatory center of the SN, up to an apparent visual magnitude of 22. Presently we will discuss more stars, roughly doubling the radius of the searched area (see Figure 1). The coordinates of the Chandra geometrical center of the remnant are: R.A. = $00^h25^m19^s9$, decl. = $64°08′18″/2$ (J2000). This is the preferred center, which pratically coincides with that of ROSAT (Hughes 2000), which differs by only 6.5 arcsec. The centroid in radio, from the Very Large Array (Reynoso et al. 1997), is also nearby. The stars closest to the center are A, B, C, D, E, F, and G. These are therefore the preferred candidates.

The distance to SN 1572 has been studied using different methods. The estimated value is converging into a value in the middle of the range from 2 to 4 kpc. Chevalier et al. (1980), using the expansion of the filaments in the remnant and the shock velocity, obtained a distance of $2.3 \pm 0.5$ kpc. A similar distance was obtained by Albinson et al. (1986) through the observation of neutral hydrogen toward the supernova. They place the distance in the range of 1.7–3.7 kpc. Just one year later, Kirshner et al. (1987) revisited the distance through the expansion of the filaments of the remnant and found it to be between 2.0 and 2.8 kpc.

RL04 attempted a different approach. By assembling the records of the historical observations of this supernova in 1572–1574 and evaluating the uncertainties, it was possible to

<figure><img src="image_3.png" /><figcaption>**Figure 1.** B-band image, taken with the 4.2 m William Herschel Telescope, showing all the stars referred to in this paper.</figcaption></figure>

reconstruct the light curve of the SN and the color. After applying the strech factor fitting of light curves of SNe, it was possible to classify this SN within the family of SNe Ia. The derived absolute magnitude was found to be consistent with a distance of $2.8 \pm 0.4$ kpc for the scale of $H_0 \sim 65$ km s$^{-1}$ Mpc$^{-1}$. In this determination, the extinction toward the supernova was derived from the reddening as shown in the color curve of the SN. Given that current estimates of $H_0$ are $67$ km s$^{-1}$ Mpc$^{-1}$, this impacts into a somehow smaller value around 2.7 kpc.

With the acknowledgement of those uncertainties, we take a range of possible distances, in this paper, between 1.7 and $3.7$ kpc ($2.7 \pm 1$ kpc) and study all the stars within this distance range as derived by Gaia as potential candidates. We discuss the distance to the stars in the next section and come back to it when talking about candidate stars. We now see a difference in the distances toward some stars, and compare with those published previously.

## 3. Parallaxes and Gaia Distances

Distances to stars targeted as possible surviving companions of SN 1572 had first been estimated by RL04 (see their Table 1), for 13 of them. Those estimates were made by fitting synthetic spectra (under the assumption of local thermodynamic equilibrium) to the observed ones. The grids of model atmospheres and the atomic data from Kurucz (1993), in combination with the Uppsala Synthetic Spectrum Package (1975), were used in the spectrum synthesis. The atmospheric parameters—effective temperature $T_{eff}$ and surface gravity $g$—were thus determined. Intrinsic colors and absolute visual magnitudes were then deduced from the relationships between spectral type and color and spectral type and absolute magnitude for the different luminosity classes (Schmidt-Kaler 1982). Comparison with BVR photometry obtained with the 2.5 m Isaac Newton Telescope in La Palma yielded the reddening $E(B-V)$, from which the visual extinction $A_V$ and the corrected apparent visual magnitude
$V_0$ were calculated. The high-resolution spectra had been obtained with the UES and ISIS spectrographs at the 4.2 m William Herschel Telescope in La Palma. Low-resolution spectra came, in addition, from the LRIS imaging spectrograph at the 10 m Keck Telescopes in Hawaii. They were compared, after dereddening, with template spectra from Lejeune et al. (1997), and these supplemented the information obtained from the high-resolution spectra.

The detailed characterization of Tycho G (singled out as a likely SN companion in RL04) was done by GH09 using a high-resolution HIRES spectrum obtained at the Keck I telescope. The stellar parameters—effective temperature and gravity—were derived using the excitation and ionization equilibria of Fe together with the fit of the wings of the $H\alpha$ line compared to different synthetic spectra computed. The result pointed to a G2 IV star with metallicity slightly below solar. The individual magnitudes in the different filters were used to estimate a range of possible distances of star G. In addition, low-resolution LIRIS spectra of the stars E, F, G, and D were obtained, which confirmed their spectral types. The best fit for Tycho G gave $T_{\text{eff}} = 5900 \pm 150$ K, $\log g = 3.85 \pm 0.35$ dex, and $[\text{Fe/H}] = -0.05 \pm 0.09$ (see GH09). This result is consistent with that obtained by K13.

K13 recalculated spectrophotometric distances to five of the stars (A, B, C, E, and G). Many of them had large error bars and are compatible with the distance estimate in B14 and the distance values implied by Gaia parallaxes, except for the value given for their star C (which is in fact three stars. Star C1 has a measured Gaia parallax corresponding to a distance of $d = 0.18_{-0.1}^{+0.3}$ kpc. The estimate in B14 is compatible with that value, whereas K13 found a too large distance of $5.5 \pm 3.5$ kpc.) Finally, in B14 there is a list of distances to 23 stars (A to W) (their Table 3), completing the work of RL04. The distances in B14 are in reasonable agreement (within $1\sigma$) with the Gaia parallaxes (see Table 2).

Now the DR2 from Gaia has provided us with precise parallaxes for almost all the stars included in those previous studies. The corresponding distances and their errors are given in column 6 of Table 2, for the stars for which we already had estimates of the distances (column 5), and in column 3 of Table 3 for those for which there were none. Gaia DR2 distances are estimated as the inverse of the parallax.

In Figure 2 are shown the Gaia DR2 distances and their error bars, and they are compared with the estimated distance to Tycho's SNR (blue vertical line) and its error bars (black dashed vertical lines). Proper motions in decl. are on the vertical axis. Solid (blue) error bars mark the stars that, within reasonable uncertainties, might be inside Tycho's SNR. Dashed (black) error bars correspond to the case where, although the error bars of the stars' distances reach the range of distances to the SNR, the errors are so large as to make their association with the SNR very unlikely. Finally, dotted (black) lines correspond to the stars that have incompatible distances with the SNR or parallax errors larger than 100%.

There are 15 stars within the range of possible distances to the SNR ($1.7 < d < 3.7$ kpc). These are stars B, F, G, H, L, N, Q, T, U, Y, AA, AC, AE, AG, and AN. From this list of stars, G and U have significant proper motions in decl. Proper motions of all the targeted stars will be discussed in the next section.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Distances and distance ranges inferred from the parallaxes in the Gaia DR2 and their uncertainties, together with their proper motions in decl. The dashed vertical lines mark the conservative limits of $2.7 \pm 1$ kpc on the distance to Tycho's SNR. Solid (blue) error bars correspond to stars that, within reasonable uncertainties, might be inside the SNR, dashed lines to the case where, although formally the stars' $1\sigma$ error bars reach the distance of the remnant, they are so large as to make it implausible, while dotted lines correspond to the stars that are beyond the limits for the SNR distance or have a parallax relative error higher than 100%.</figcaption></figure>

In Table 1, the Gaia DR2 data on parallaxes, proper motions, and G magnitudes are given as in the Gaia DR2. The proper motions are absolute, referred to in the ICRS (as mentioned above).

In Table 2, the Gaia DR2 distances, as stated, are compared with the distances deduced in B14 from determination of the stellar atmospheric parameters and comparison of the resulting luminosities with the available photometry. We see that there is reasonable agreement in most cases, with a tendency to place the stars at longer distances in B14 as compared with Gaia, which can be attributed to an underestimate, in B14, of the extinction in the direction of Tycho's SNR. Based only on the stars with distance errors $\le 0.5$ kpc in both sets, the underestimate would be by $\Delta\Delta V \simeq 0.5$ mag. Exceptions are stars N, P1, U, V, and W, although in the last two cases the Gaia error bars are so large that the comparison is not really meaningful.

## 4. Proper Motions from Gaia Compared with the HST

High velocities, mainly due to their orbital motion at the time of explosion, must be a salient characteristic of SNe Ia companions. Unless they were mostly moving along the line of sight when the binary system was disrupted, the components of the velocity on the plane of the sky should be observed as high proper motions relative to the stars around the location of the SN. It might happen that the stars were moving at some angle with the line of sight and thus the components of the velocity would be distributed accordingly.

The location of the supernova explosion, in the case of SNe whose remnants still exist, is given in a first, rough approximation, by the centroid of the remnant. High-precision astrometric measurements of the proper motions of the stars within some angular distance from the centroid are the tool
<table><caption>**Table 1**

Gaia IDs, Parallaxes, Proper Motions, and *G* Magnitudes of the Sample of Stars from the Gaia DR2</caption><thead><tr><th>Star</th><th>Gaia ID</th><th>ω<br>(mas)</th><th>μ<sub>α</sub> cos δ<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)</th><th>G<br>(mag)</th></tr><tr><th>(1)</th><th>(2)</th><th>(3)</th><th>(4)</th><th>(5)</th><th>(6)</th></tr></thead><tbody><tr><td>A</td><td>431160565571641856</td><td>1.031 ± 0.052</td><td>-5.321 ± 0.076</td><td>-3.517 ± 0.065</td><td>12.404 ± 0.001</td></tr><tr><td>B</td><td>431160569875463936</td><td>0.491 ± 0.051</td><td>-4.505 ± 0.063</td><td>-0.507 ± 0.049</td><td>15.113 ± 0.001</td></tr><tr><td>C1</td><td>431160359417132800</td><td>5.310 ± 0.483</td><td>-2.415 ± 0.735</td><td>-0.206 ± 0.576</td><td>18.027 ± 0.006</td></tr><tr><td>D</td><td>431160363709280768</td><td>1.623 ± 0.318</td><td>-4.566 ± 0.636</td><td>-2.248 ± 0.376</td><td>19.371 ± 0.003</td></tr><tr><td>E</td><td>431160565573859584</td><td>0.138 ± 0.220</td><td>0.232 ± 0.377</td><td>-0.699 ± 0.265</td><td>18.970 ± 0.002</td></tr><tr><td>F</td><td>431160569875460096</td><td>0.466 ± 0.079</td><td>-5.739 ± 0.130</td><td>-0.292 ± 0.097</td><td>17.036 ± 0.001</td></tr><tr><td>G</td><td>431160359413315328</td><td>0.512 ± 0.021</td><td>-4.417 ± 0.191</td><td>-4.064 ± 0.143</td><td>17.988 ± 0.001</td></tr><tr><td>H</td><td>431160599931508480</td><td>0.620 ± 0.203</td><td>-4.839 ± 0.341</td><td>-0.577 ± 0.248</td><td>18.895 ± 0.002</td></tr><tr><td>I</td><td>431160569867713152</td><td>-0.014 ± 0.566</td><td>-1.479 ± 0.970</td><td>-0.855 ± 0.761</td><td>20.351 ± 0.006</td></tr><tr><td>J</td><td>431160565571749760</td><td>0.134 ± 0.240</td><td>-3.900 ± 0.373</td><td>-1.054 ± 0.292</td><td>18.965 ± 0.002</td></tr><tr><td>K</td><td>431160393780294144</td><td>-0.266 ± 0.290</td><td>-1.735 ± 0.601</td><td>-0.815 ± 0.350</td><td>19.313 ± 0.003</td></tr><tr><td>L</td><td>431160398076768896</td><td>0.689 ± 0.457</td><td>-2.471 ± 0.876</td><td>0.514 ± 0.578</td><td>20.072 ± 0.005</td></tr><tr><td>M</td><td>431160604230502400</td><td>-2.282 ± 1.316</td><td>3.472 ± 1.943</td><td>-1.624 ± 2.070</td><td>20.900 ± 0.011</td></tr><tr><td>N</td><td>431160565571767552</td><td>0.246 ± 0.096</td><td>0.092 ± 0.148</td><td>0.134 ± 0.121</td><td>17.612 ± 0.001</td></tr><tr><td>O</td><td>431160569875457792</td><td>1.169 ± 0.063</td><td>2.607 ± 0.098</td><td>2.108 ± 0.076</td><td>16.542 ± 0.001</td></tr><tr><td>P</td><td>431160565571767424</td><td>0.168 ± 0.092</td><td>-0.889 ± 0.139</td><td>-0.389 ± 0.106</td><td>16.998 ± 0.001</td></tr><tr><td>Q</td><td>431160565575562240</td><td>0.663 ± 0.334</td><td>0.438 ± 0.643</td><td>1.409 ± 0.404</td><td>19.496 ± 0.004</td></tr><tr><td>S</td><td>431160565573859840</td><td>1.235 ± 0.417</td><td>2.091 ± 0.771</td><td>-0.437 ± 0.491</td><td>19.568 ± 0.003</td></tr><tr><td>T</td><td>431159088102994432</td><td>0.565 ± 0.289</td><td>-5.177 ± 0.563</td><td>0.004 ± 0.353</td><td>19.320 ± 0.003</td></tr><tr><td>U</td><td>431159092406721280</td><td>0.504 ± 0.070</td><td>-1.877 ± 0.113</td><td>-5.096 ± 0.083</td><td>17.064 ± 0.001</td></tr><tr><td>V</td><td>431160359413311616</td><td>0.059 ± 1.023</td><td>-2.201 ± 1.184</td><td>1.645 ± 1.279</td><td>20.235 ± 0.007</td></tr><tr><td>W</td><td>431160393773079808</td><td>0.193 ± 0.283</td><td>-2.760 ± 0.600</td><td>0.163 ± 0.343</td><td>19.312 ± 0.003</td></tr><tr><td>X</td><td>431159092398964992</td><td>0.192 ± 0.427</td><td>-1.187 ± 0.812</td><td>-0.836 ± 0.511</td><td>19.812 ± 0.004</td></tr><tr><td>Y</td><td>431159092406717568</td><td>0.631 ± 0.223</td><td>0.144 ± 0.347</td><td>-2.261 ± 0.290</td><td>18.923 ± 0.002</td></tr><tr><td>Z</td><td>431159092398966400</td><td>0.176 ± 0.146</td><td>-1.498 ± 0.233</td><td>-0.294 ± 0.193</td><td>18.082 ± 0.002</td></tr><tr><td>AA</td><td>431159088102995968</td><td>0.957 ± 0.467</td><td>-2.277 ± 0.977</td><td>-1.184 ± 0.595</td><td>19.973 ± 0.005</td></tr><tr><td>AB</td><td>431159088102989056</td><td>-0.090 ± 0.267</td><td>-2.011 ± 0.445</td><td>-1.600 ± 0.316</td><td>19.046 ± 0.002</td></tr><tr><td>AC</td><td>431159088103003520</td><td>0.490 ± 0.160</td><td>-2.376 ± 0.249</td><td>-1.445 ± 0.195</td><td>18.399 ± 0.001</td></tr><tr><td>AE</td><td>431159088102986880</td><td>0.279 ± 0.173</td><td>-0.907 ± 0.268</td><td>-0.241 ± 0.227</td><td>18.559 ± 0.002</td></tr><tr><td>AF</td><td>431158881944551424</td><td>1.323 ± 0.324</td><td>-2.381 ± 0.644</td><td>0.489 ± 0.400</td><td>19.399 ± 0.003</td></tr><tr><td>AG</td><td>431158881944550272</td><td>0.703 ± 0.382</td><td>-2.412 ± 0.779</td><td>0.626 ± 0.453</td><td>19.768 ± 0.004</td></tr><tr><td>AH</td><td>431158881944553216</td><td>0.206 ± 0.087</td><td>-0.704 ± 0.139</td><td>-0.579 ± 0.107</td><td>17.486 ± 0.001</td></tr><tr><td>AII/HP1</td><td>431160359417132928</td><td>2.831 ± 0.273</td><td>71.558 ± 0.530</td><td>-3.030 ± 0.322</td><td>19.159 ± 0.003</td></tr><tr><td>AJ</td><td>431160359413306368</td><td>0.187 ± 0.206</td><td>-1.102 ± 0.333</td><td>0.222 ± 0.249</td><td>18.883 ± 0.002</td></tr><tr><td>AK</td><td>431160393773068032</td><td>-0.476 ± 0.447</td><td>-1.306 ± 0.938</td><td>-0.142 ± 0.518</td><td>20.008 ± 0.005</td></tr><tr><td>AL</td><td>431160398072078592</td><td>0.383 ± 0.618</td><td>-2.827 ± 1.199</td><td>0.727 ± 0.789</td><td>20.461 ± 0.006</td></tr><tr><td>AM</td><td>431160398073281792</td><td>0.752 ± 0.825</td><td>0.303 ± 1.636</td><td>-2.940 ± 1.114</td><td>20.605 ± 0.008</td></tr><tr><td>AN</td><td>431160599931516288</td><td>0.605 ± 0.138</td><td>-4.560 ± 0.221</td><td>-1.330 ± 0.168</td><td>18.284 ± 0.001</td></tr></tbody></table>

needed to detect or discard the presence of possible companions. Until the advent of Gaia, this was only possible with HST astrometry. Radial velocities obtained from high-resolution spectra give complementary information on the velocity along the line of sight.

A first set of measurements of the proper motions of the stars around the centroid of Tycho's SNR was made in RL04. It included 26 stars, labeled from A to W (see their Figure 1). Images from the WFPC2 aboard the HST, taken two months apart (within Cycle 12), were used. It was found that star G was at compatible distance to the supernova (the distance estimate for Tycho G, the most widely named star G, was $3.0_{-0.5}^{+1}$ kpc and its motion was mostly perpendicular to the Galactic plane, with $\mu_b = 6.11 \pm 1.34 \text{ mas yr}^{-1}$ ($\mu_t = -2.6 \pm 1.34 \text{ mas yr}^{-1}$ only). That meant a heliocentric tangential velocity $v_t \sim 94 \text{ km s}^{-1}$ which, combined with a high measured radial velocity of $-87.4 \pm 0.5 \text{ km s}^{-1}$, gave a total heliocentric velocity $v_{tot} \sim 128 \text{ km s}^{-1} \pm 9 \text{ km s}^{-1}$, making star G a likely candidate to have been the companion star of SN 1572.

In B14, proper motions $\mu_\alpha \cos \delta = -2.63 \pm 0.18 \text{ mas yr}^{-1}$ and $\mu_\delta = -3.98 \pm 0.10 \text{ mas yr}^{-1}$ were measured. That, for a distance of $2.83 \pm 0.79 \text{ kpc}$, taken as the SN distance, give a heliocentric tangential velocity $v_t = 64 \pm 11 \text{ km s}^{-1}$. Given the heliocentric radial velocity, $v_r = -87.4 \pm 0.5$, the total velocity is $v_{tot} = 108 \pm 9 \text{ km s}^{-1}$.

Now, from the Gaia DR2, $\mu_\alpha \cos \delta = -4.417 \pm 0.191 \text{ mas yr}^{-1}$ and $\mu_\delta = -4.064 \pm 0.143 \text{ mas yr}^{-1}$. The parallax being $\omega = 0.512 \pm 0.021$, we have $v_\alpha \cos \delta = -40.88 \pm 2.44 \text{ km s}^{-1}$, $v_\delta = -37.61 \pm 2.03 \text{ km s}^{-1}$, and $v_t = 55.55 \pm 2.26 \text{ km s}^{-1}$ (heliocentric).

For the corresponding $v_r$, it results in $v_{tot} = 103.69 \pm 7.52 \text{ km s}^{-1}$. We thus see that, as compared with B14, there is no significant change with the new results. As in B14, they can be interpreted in the framework of a binary model similar to that for U Sco. The excess velocity of star G with respect to the average of the stars at the same location in the Galaxy could come from the orbital velocity it had when the binary was disrupted by the SN explosion.
<table><caption>**Table 2**
BVR Photometry, Distances, and Proper Motions of Stars A–W from B14, Compared with the Distances and Proper Motions from Gaia DR2</caption><thead><tr><th>Star</th><th>B<br>(mag)</th><th>V<br>(mag)</th><th>R<br>(mag)</th><th>d (B14)<br>(kpc)</th><th>d<br>(kpc)</th><th>μ<sub>α</sub> cos δ<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>α</sub> cos δ (B14)<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub> (B14)<br>(mas yr<sup>-1</sup>)</th></tr><tr><th>(1)</th><th>(2)</th><th>(3)</th><th>(4)</th><th>(5)</th><th>(6)</th><th>(7)</th><th>(8)</th><th>(9)</th><th>(10)</th></tr></thead><tbody><tr><td>A</td><td>14.82 ± 0.03</td><td>13.29 ± 0.03</td><td>12.24 ± 0.03</td><td>1.1 ± 0.3</td><td>0.97<sup>+0.05</sup><sub>-0.04</sub></td><td>-3.63 ± 0.08</td><td>...</td><td>-3.06 ± 0.07</td><td>...</td></tr><tr><td>B</td><td>16.35 ± 0.03</td><td>15.41 ± 0.03</td><td>...</td><td>2.6 ± 0.5</td><td>2.03<sup>+0.19</sup><sub>-0.15</sub></td><td>-2.90 ± 0.06</td><td>-1.67 ± 0.06</td><td>-0.09 ± 0.05</td><td>0.59 ± 0.08</td></tr><tr><td>C1</td><td>21.06 ± 0.12</td><td>19.06 ± 0.05</td><td>17.77 ± 0.03</td><td>0.75 ± 0.5</td><td>0.18<sup>+0.03</sup><sub>-0.01</sub></td><td>-0.82 ± 0.73</td><td>-1.98 ± 0.07</td><td>0.39 ± 0.58</td><td>-1.09 ± 0.06</td></tr><tr><td>C2</td><td>22.91 ± 0.20</td><td>20.53 ± 0.15</td><td>...</td><td>~40</td><td>...</td><td>...</td><td>-1.75 ± 0.07</td><td>...</td><td>-1.07 ± 0.07</td></tr><tr><td>C3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.08 ± 0.11</td><td>...</td><td>-0.14 ± 0.10</td></tr><tr><td>D</td><td>22.97 ± 0.28</td><td>20.70 ± 0.10</td><td>19.38 ± 0.06</td><td>0.8 ± 0.2</td><td>0.62<sup>+0.15</sup><sub>-0.11</sub></td><td>-2.97 ± 0.64</td><td>-2.03 ± 0.09</td><td>-1.65 ± 0.38</td><td>-1.28 ± 0.07</td></tr><tr><td>E</td><td>21.24 ± 0.13</td><td>19.79 ± 0.07</td><td>18.84 ± 0.05</td><td>>20</td><td>7.22<sup>+xxy</sup><sub>-4.43</sub></td><td>1.83 ± 0.38</td><td>1.74 ± 0.05</td><td>-0.10 ± 0.26</td><td>0.28 ± 0.05</td></tr><tr><td>F</td><td>19.02 ± 0.05</td><td>17.73 ± 0.03</td><td>16.94 ± 0.03</td><td>1.5 ± 0.5</td><td>2.15<sup>+0.44</sup><sub>-0.32</sub></td><td>-4.14 ± 0.13</td><td>-3.31 ± 0.15</td><td>0.31 ± 0.10</td><td>0.25 ± 0.07</td></tr><tr><td>G</td><td>20.09 ± 0.08</td><td>18.71 ± 0.04</td><td>17.83 ± 0.03</td><td>2.5–5.0</td><td>1.95<sup>+0.60</sup><sub>-0.35</sub></td><td>-2.82 ± 0.19</td><td>-2.63 ± 0.06</td><td>-3.46 ± 0.14</td><td>-3.98 ± 0.04</td></tr><tr><td>H</td><td>21.39 ± 0.14</td><td>19.80 ± 0.07</td><td>18.78 ± 0.05</td><td>~1.8/~24</td><td>1.61<sup>+0.79</sup><sub>-0.40</sub></td><td>-3.24 ± 0.34</td><td>-3.13 ± 0.07</td><td>-0.02 ± 0.25</td><td>-0.84 ± 0.03</td></tr><tr><td>I</td><td>...</td><td>21.75 ± 0.16</td><td>20.36 ± 0.09</td><td>~4</td><td>≥1.81</td><td>0.12 ± 0.97</td><td>0.69 ± 0.06</td><td>-0.25 ± 0.76</td><td>-0.20 ± 0.06</td></tr><tr><td>J</td><td>21.15 ± 0.12</td><td>19.74 ± 0.07</td><td>18.84 ± 0.05</td><td>~9</td><td>7.46<sup>+xxy</sup><sub>-4.77</sub></td><td>-2.30 ± 0.37</td><td>-2.35 ± 0.06</td><td>-0.45 ± 0.29</td><td>-0.28 ± 0.03</td></tr><tr><td>K</td><td>21.64 ± 0.15</td><td>20.11 ± 0.08</td><td>19.15 ± 0.05</td><td>~2.4/~27</td><td>≥41.67</td><td>-0.14 ± 0.60</td><td>0.24 ± 0.12</td><td>-0.21 ± 0.35</td><td>0.03 ± 0.07</td></tr><tr><td>L</td><td>22.77 ± 0.26</td><td>21.08 ± 0.12</td><td>20.00 ± 0.07</td><td>~4</td><td>1.45<sup>+2.87</sup><sub>-0.58</sub></td><td>-0.87 ± 0.88</td><td>0.36 ± 0.12</td><td>1.11 ± 0.58</td><td>-0.08 ± 0.04</td></tr><tr><td>M</td><td>23.49 ± 0.36</td><td>21.82 ± 0.16</td><td>20.72 ± 0.10</td><td>~4</td><td>...</td><td>5.07 ± 1.94</td><td>-0.61 ± 0.12</td><td>-1.02 ± 2.07</td><td>0.44 ± 0.08</td></tr><tr><td>N</td><td>19.59 ± 0.06</td><td>18.29 ± 0.04</td><td>17.47 ± 0.03</td><td>~1.5–2</td><td>4.06<sup>+2.57</sup><sub>-1.14</sub></td><td>1.69 ± 0.15</td><td>2.64 ± 0.13</td><td>0.74 ± 0.12</td><td>0.96 ± 0.04</td></tr><tr><td>O</td><td>18.62 ± 0.04</td><td>17.23 ± 0.03</td><td>16.37 ± 0.03</td><td><1</td><td>0.85<sup>+0.54</sup><sub>-0.22</sub></td><td>4.21 ± 0.10</td><td>5.13 ± 0.20</td><td>2.71 ± 0.08</td><td>2.85 ± 0.14</td></tr><tr><td>P1</td><td>...</td><td>17.61 ± 0.03</td><td>16.78 ± 0.03</td><td>~1</td><td>5.96<sup>+7.23</sup><sub>-2.11</sub></td><td>...</td><td>1.39 ± 0.36</td><td>...</td><td>0.20 ± 0.09</td></tr><tr><td>P2</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>-0.27 ± 0.20</td><td>...</td><td>-1.64 ± 0.21</td></tr><tr><td>Q</td><td>22.35 ± 0.21</td><td>20.59 ± 0.09</td><td>19.41 ± 0.06</td><td>~2</td><td>1.51<sup>+1.53</sup><sub>-0.51</sub></td><td>2.04 ± 0.64</td><td>1.34 ± 0.09</td><td>2.71 ± 0.40</td><td>2.38 ± 0.04</td></tr><tr><td>R</td><td>22.91 ± 0.28</td><td>21.38 ± 0.13</td><td>20.26 ± 0.08</td><td>3.3 ± 0.2</td><td>...</td><td>...</td><td>-0.18 ± 0.10</td><td>...</td><td>0.25 ± 0.05</td></tr><tr><td>S</td><td>...</td><td>21.30 ± 0.13</td><td>19.74 ± 0.07</td><td>1.3 ± 0.1</td><td>0.81<sup>+0.41</sup><sub>-0.20</sub></td><td>3.69 ± 0.77</td><td>3.68 ± 0.09</td><td>0.16 ± 0.49</td><td>0.93 ± 0.05</td></tr><tr><td>T</td><td>21.82 ± 0.17</td><td>20.23 ± 0.08</td><td>19.20 ± 0.05</td><td>~2/~30</td><td>1.77<sup>+1.86</sup><sub>-0.60</sub></td><td>-3.58 ± 0.56</td><td>-2.96 ± 0.04</td><td>0.61 ± 0.35</td><td>-0.53 ± 0.05</td></tr><tr><td>U</td><td>19.03 ± 0.05</td><td>17.73 ± 0.03</td><td>16.95 ± 0.03</td><td>~1</td><td>1.98<sup>+0.32</sup><sub>-0.24</sub></td><td>-0.28 ± 0.11</td><td>0.39 ± 0.10</td><td>-4.49 ± 0.08</td><td>-4.31 ± 0.07</td></tr><tr><td>V</td><td>23.32 ± 0.33</td><td>21.41 ± 0.13</td><td>20.20 ± 0.08</td><td>~3</td><td>16.81<sup>+xxy</sup><sub>-15.89</sub></td><td>-1.16 ± 1.18</td><td>-0.67 ± 0.08</td><td>2.25 ± 1.28</td><td>0.49 ± 0.08</td></tr><tr><td>W</td><td>22.13 ± 0.19</td><td>20.44 ± 0.09</td><td>19.27 ± 0.05</td><td>~2</td><td>5.17<sup>+xxy</sup><sub>-3.07</sub></td><td>-1.16 ± 0.60</td><td>-0.31 ± 0.09</td><td>0.76 ± 0.34</td><td>0.09 ± 0.04</td></tr></tbody></table>

Note. Here the Gaia proper motions have been transformed to the system used in B14; see the text and, for the Gaia original values, Table 1. There are stars that have no upper limit for the Gaia distances. This corresponds to negative parallaxes. They have been marked with xxx. There are other stars (I, K, and AB in Table 3) with negative central value of the parallax, but adding the errors gives positive limits. This corresponds to a lower limit for the distance, therefore they show a $\ge$sign.

Taking as a reference the Besançon model of the Galaxy (Robin et al. 2003), the average heliocentric tangential velocity of disk stars, at the position and distance of star G, is almost negligible, while the average radial velocity is $\langle v_r \rangle \approx -31 \pm 28$ km s$^{-1}$. Then, attributing the excess over average in radial velocity ($\approx -56 \pm 28$ km s$^{-1}$) to orbital motion and the full tangential velocity ($\approx 55 \pm 2$ km s$^{-1}$) to it, we obtain that $v_{\text{orb}} \approx 78 \pm 20$ km s$^{-1}$ (the inclination of the plane of the orbit with respect to the line of sight would thus be $i = 44^{\circ}$).

The evolutionary path giving rise to SN 1572 might have started from a WD with a mass $\sim 0.8 M_{\odot}$ plus a somewhat evolved companion of $\sim 2.0–2.5 M_{\odot}$ filling its Roche lobe (RL04), the system ending up as a WD with Chandrasekhar mass $\sim 1.4 M_{\odot}$, plus a companion of $\sim 1 M_{\odot}$. Using Kepler’s law, that $P^2 = a^3/(M_1+M_2)$ (with $P$ in years, $a$ in astronomical units, $M_1$ and $M_2$ in solar masses), and since $P = 2\pi a/v$, we find a separation $a \approx 25 R_{\odot}$, the period being $P \approx 9$ days. Using Eggleton’s (1983) formula

$$R_L = a \left[ \frac{0.49}{0.6 + q^{-2/3} \ln(1 + q^{1/3})} \right]$$

($q$ being the mass ratio $M_2/M_1$) for the effective Roche lobe radius of the companion just before the explosion, it would thus have been $\approx 9 R_{\odot}$. At present, the radius of the star is only $1–2 R_{\odot}$ (GH09), and would have resulted from the combination

of mass stripping and shock heating by the impact of the SN ejecta, plus subsequent fast cooling of the outer layers up to the present time.

However, an alternative to this hypothesis is discussed later in the paper.

In B14, the proper motions of 872 stars were measured from HST astrometry, using images taken in up to four different epochs and spanning a total of 8 yr. Much higher precision than in RL04 was achieved. The results for 45 of them (all the stars with names in Figure 1) are given in Tables 2 and 3 of B14. The full version was provided as supplementary electronic material.

When comparing the proper motions given by the Gaia DR2 with those obtained by B14 from the astrometry done with the HST, one must take into account that the former are absolute measurements, in the ICRS system, while the latter are relative measurements. This means that the local frame used for the HST astrometry should, in general, move with respect to the ICRS frame. Such a systematic effect is actually seen when we make the comparison. In Tables 2 and 3, the Gaia proper motions have been transformed to the HST frame.

Including only the stars with proper motion errors smaller than 0.25 mas yr$^{-1}$ in B14, we find that, on average, $\mu_{\alpha} \cos \delta$ (Gaia) = $\mu_{\alpha} \cos \delta$ (B14) – 1.599 ± 0.729 mas yr$^{-1}$ and $\mu_{\delta}$ (Gaia) = $\mu_{\delta}$ (B14) – 0.601 ± 0.585 mas yr$^{-1}$.
<table><caption>**Table 3**

G Magnitudes, Distances, and Proper Motions of Stars X–AN from B14, with Comparison of the Proper Motions from B14 and From Gaia DR2</caption><thead><tr><th>Star</th><th>G<br>(mag)</th><th>d<br>(kpc)</th><th>μ<sub>α</sub> cos δ<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>α</sub> cos δ (B14)<br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub> (B14)<br>(mas yr<sup>-1</sup>)</th></tr><tr><th>(1)</th><th>(2)</th><th>(3)</th><th>(4)</th><th>(5)</th><th>(6)</th><th>(7)</th></tr></thead><tbody><tr><td>X</td><td>19.81</td><td>5.20<sup>+xxx</sup><sub>-3.59</sub></td><td>0.41 ± 0.81</td><td>1.32 ± 0.06</td><td>-0.24 ± 0.51</td><td>-0.26 ± 0.06</td></tr><tr><td>Y</td><td>18.92</td><td>1.58<sup>+0.87</sup><sub>-0.41</sub></td><td>1.74 ± 0.35</td><td>2.52 ± 0.03</td><td>-1.66 ± 0.29</td><td>-1.75 ± 0.07</td></tr><tr><td>Z</td><td>18.08</td><td>5.68<sup>+2.72</sup><sub>-2.57</sub></td><td>0.10 ± 0.23</td><td>0.78 ± 0.04</td><td>0.31 ± 0.19</td><td>0.06 ± 0.07</td></tr><tr><td>AA</td><td>19.97</td><td>1.04<sup>+1.00</sup><sub>-0.33</sub></td><td>-0.68 ± 1.00</td><td>-3.19 ± 0.07</td><td>-0.58 ± 0.60</td><td>-1.42 ± 0.04</td></tr><tr><td>AB</td><td>19.05</td><td>≥5.65<sup>+</sup><sub>-</sub></td><td>-0.41 ± 0.45</td><td>-0.37 ± 0.03</td><td>-1.00 ± 0.32</td><td>-1.01 ± 0.06</td></tr><tr><td>AC</td><td>18.40</td><td>2.04<sup>+0.99</sup><sub>-0.50</sub></td><td>-0.78 ± 0.25</td><td>-1.09 ± 0.07</td><td>-0.84 ± 0.20</td><td>-0.87 ± 0.04</td></tr><tr><td>AD1</td><td>17.22</td><td>1.04<sup>+0.84</sup><sub>-0.32</sub></td><td>0.85 ± 0.59</td><td>-1.24 ± 0.14</td><td>...</td><td>-1.58 ± 0.16</td></tr><tr><td>AD2</td><td>...</td><td>...</td><td>...</td><td>-1.12 ± 0.04</td><td>...</td><td>-2.25 ± 0.07</td></tr><tr><td>AE</td><td>19.05</td><td>3.58<sup>+5.84</sup><sub>-1.37</sub></td><td>0.69 ± 0.27</td><td>1.07 ± 0.06</td><td>0.36 ± 0.23</td><td>-0.05 ± 0.06</td></tr><tr><td>AF</td><td>19.40</td><td>0.76<sup>+0.24</sup><sub>-0.15</sub></td><td>-0.78 ± 0.64</td><td>-0.38 ± 0.04</td><td>1.09 ± 0.40</td><td>0.01 ± 0.07</td></tr><tr><td>AG</td><td>19.77</td><td>1.42<sup>+1.69</sup><sub>-0.50</sub></td><td>-0.81 ± 0.78</td><td>-1.11 ± 0.06</td><td>1.23 ± 0.45</td><td>0.93 ± 0.08</td></tr><tr><td>AH</td><td>17.49</td><td>4.85<sup>+3.55</sup><sub>-1.44</sub></td><td>0.89 ± 0.14</td><td>1.26 ± 0.20</td><td>0.02 ± 0.11</td><td>-0.20 ± 0.26</td></tr><tr><td>AI1/HP-1</td><td>19.16</td><td>0.35<sup>+0.04</sup><sub>-0.03</sub></td><td>73.16 ± 0.53</td><td>73.07 ± 0.09</td><td>-2.43 ± 0.32</td><td>-2.82 ± 0.07</td></tr><tr><td>AI2</td><td>...</td><td>...</td><td>...</td><td>1.76 ± 0.28</td><td>...</td><td>0.16 ± 0.21</td></tr><tr><td>AJ</td><td>18.88</td><td>5.35<sup>+xxx</sup><sub>-2.81</sub></td><td>0.50 ± 0.33</td><td>0.18 ± 0.05</td><td>0.82 ± 0.25</td><td>0.73 ± 0.07</td></tr><tr><td>AK</td><td>20.01</td><td>...</td><td>0.29 ± 0.94</td><td>-0.25 ± 0.09</td><td>0.46 ± 0.52</td><td>0.95 ± 0.08</td></tr><tr><td>AL</td><td>20.46</td><td>2.61<sup>+xxx</sup><sub>-1.61</sub></td><td>-1.23 ± 1.20</td><td>-0.14 ± 0.09</td><td>1.33 ± 0.73</td><td>-0.35 ± 0.09</td></tr><tr><td>AM</td><td>20.61</td><td>1.33<sup>+xxx</sup><sub>-0.70</sub></td><td>1.90 ± 1.64</td><td>-1.36 ± 0.10</td><td>-2.34 ± 1.11</td><td>-0.60 ± 0.10</td></tr><tr><td>AN</td><td>18.28</td><td>1.65<sup>+0.49</sup><sub>-0.30</sub></td><td>-2.96 ± 0.22</td><td>-2.83 ± 0.11</td><td>-0.73 ± 0.17</td><td>-0.96 ± 0.05</td></tr></tbody></table>

Note. Here the Gaia proper motions have been transformed to the system used in B14; for the Gaia original values, see Table 1. These stars have not been assigned a distance in B14 or in any other paper. The sign xxx for the upper limit in distance to some stars corresponds to negative upper limit parallax, as in Table 2.

In Table 2 (columns 7 and 9) and Table 3 (columns 4 and 6), we have transformed the Gaia proper motions to the B14 HST frame according to these relations. For our purposes, the local, relative proper motions are most meaningful, since we are interested in the motions of the stars with respect to the average motions of those around their positions. After applying these zero-point shifts, there still are residual differences between the two proper motion sets. On average, $\Delta\mu_{\alpha} \cos\delta = -0.017 \pm 0.788 \text{ mas yr}^{-1}$ and $\Delta\mu_{\delta} = 0.005 \pm 0.630 \text{ mas yr}^{-1}$. The whole set is included here, the dispersion being mainly due to stars which have substantial errors in their Gaia proper motions (see columns 7 and 9 in Table 2 and columns 4 and 6 in Table 3, as well as columns 5 and 6 in Table 4).

## 5. Toomre Diagram

Gaia provides a five-parameter astrometric solution and, for some stars, line-of-sight velocities ($\alpha, \delta, \varpi, \mu_{\alpha}^*$, $\mu_{\delta}, V_r$), together with their associated uncertainties and correlations between the astrometric quantities. For the 13 stars for which we also know their radial velocities, the total space velocities can be derived. It is most useful to see their components in the Galactic coordinate system: $U$ (positive in the direction of the Galactic center), $V$ (positive in the direction of Galactic rotation) and $W$ (positive in the direction of the north Galactic pole) in the Local Standard of Rest (LSR). In Table 4 we give the $U, V$, and $W$ components of the space velocities, as well as the total velocities on the Galactic meridian plane, in the LSR, of these 13 stars, based on the Gaia DR2 parallaxes and proper motions and on the radial velocities from B14 (save for star A, which has a quite precise radial velocity from Gaia). For the transformation of the motions from heliocentric to the LSR, we have adopted, as the peculiar velocity of the Sun with

respect to the LSR, ($U_{\odot}, V_{\odot}, W_{\odot}$) = (11.1, 12.24, 7.25) km s$^{-1}$ (Schönrich et al. 2010).

The Toomre diagram shows the distribution of $U, V, W$ velocities in the LSR. This diagram combines quadratically $U$ and $W$ versus $V$ and allows us to distinguish between stars belonging to different Galactic stellar components (thin disk, transition thin-thick, thick disk, and halo). The Adibekyan et al. (2012) sample, with very high quality spectroscopic data from the HARPS exoplanet search program, is used as a reference in Figure 3. The data have very high precision on radial velocities, stellar parameters, and metallicities. Therefore, the Toomre diagram obtained from this sample combines kinematics information such as orbital motion together with the information derived from the spectroscopic analysis on element abundances. The sample covers a wide range of metallicities [Fe/H] from -1.2 to 0.4.

The sample in the upper left panel of Figure 3 has no imposed boundaries on metallicity, while those in the lower left and right panels include only stars with metallicities equal to or higher than that of star G minus the $1\sigma$ uncertainty, i.e., for [Fe/H] > -0.14. One sees there (with the exception of star J, whose kinematics is very unreliable, with large errors in the Gaia DR2 data), that no other star in our sample moves as fast as star G.

The Gaia DR2 data place star G above the region where most thin disk stars are. The kinematics of star G would locate it among the thick disk stars but its metallicity is that of a thin disk star, while at its location, only 48 pc above the Galactic plane, the density of thick disk stars is very low. Using the Adibekyan et al. (2012) sample, the probability that star G belonged to the thick disk, given its metallicity, is only 2%.

There are some thin disk stars, however, that move fast on the Galactic meridian plane, and thus star G might belong to
<caption>Table 4</caption>

Parallaxes, Heliocentric Radial Velocities, Proper Motions, Galactic *U*, *V*, *W* Velocity Components, and Total Velocities on the Galactic Meridian Plane (Referred to the Local Standard of Rest) of the Stars with Both Radial Velocities from B14 and Parallaxes from Gaia DR2

<table><thead><tr><th>Star</th><th>DR2 number<br>(4311...)<br>(2)</th><th>ω<br>(mas)<br>(3)</th><th>v<sub>r</sub><br>(km s<sup>-1</sup>)<br>(4)</th><th>μ<sub>α</sub> cos δ<br>(mas yr<sup>-1</sup>)<br>(5)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)<br>(6)</th><th>U<br>(km s<sup>-1</sup>)<br>(7)</th><th>V<br>(km s<sup>-1</sup>)<br>(8)</th><th>W<br>(km s<sup>-1</sup>)<br>(9)</th><th>(U<sup>2</sup>+W<sup>2</sup>)<sup>1/2</sup><br>(km s<sup>-1</sup>)<br>(10)</th></tr></thead><tbody><tr><td>A</td><td>60565571641856</td><td>1.03 ± 0.05</td><td>-30 ± 1</td><td>-5.321 ± 0.076</td><td>-3.518 ± 0.065</td><td>48.65 ± 1.12</td><td>-0.74 ± 0.74</td><td>-7.12 ± 0.82</td><td>49.17 ± 1.11</td></tr><tr><td>B</td><td>60569875463936</td><td>0.49 ± 0.04</td><td>-45 ± 8</td><td>-4.505 ± 0.063</td><td>-0.507 ± 0.049</td><td>71.69 ± 5.08</td><td>-5.10 ± 7.15</td><td>5.72 ± 0.53</td><td>71.91 ± 5.07</td></tr><tr><td>C1</td><td>60359417132800</td><td>5.31 ± 0.48</td><td>-40 ± 6</td><td>-2.415 ± 0.735</td><td>-0.206 ± 0.576</td><td>33.22 ± 3.06</td><td>-23.00 ± 5.20</td><td>6.29 ± 0.54</td><td>33.81 ± 3.05</td></tr><tr><td>D</td><td>60363709280768</td><td>1.62 ± 0.32</td><td>-58 ± 0.8</td><td>-4.566 ± 0.636</td><td>-2.248 ± 0.376</td><td>52.16 ± 2.73</td><td>-30.83 ± 1.61</td><td>0.64 ± 1.68</td><td>52.16 ± 2.73</td></tr><tr><td>E</td><td>60565573859584</td><td>0.14 ± 0.22</td><td>-33 ± 18</td><td>0.232 ± 0.377</td><td>-0.699 ± 0.265</td><td>22.80 ± 17.45</td><td>-18.87 ± 17.43</td><td>-18.20 ± 42.45</td><td>29.18 ± 28.82</td></tr><tr><td>F</td><td>60569875460096</td><td>0.47 ± 0.08</td><td>-41 ± 11</td><td>-5.739 ± 0.130</td><td>-0.292 ± 0.097</td><td>82.40 ± 10.16</td><td>5.63 ± 10.72</td><td>9.19 ± 1.02</td><td>82.91 ± 10.10</td></tr><tr><td>G</td><td>60359413315328</td><td>0.51 ± 0.12</td><td>-87 ± 0.5</td><td>-4.417 ± 0.191</td><td>-4.064 ± 0.143</td><td>93.01 ± 8.85</td><td>-40.33 ± 5.37</td><td>-28.20 ± 8.28</td><td>97.19 ± 8.80</td></tr><tr><td>H</td><td>60599931508480</td><td>0.62 ± 0.20</td><td>-78 ± 10</td><td>-4.839 ± 0.341</td><td>-0.577 ± 0.248</td><td>82.61 ± 11.38</td><td>-36.89 ± 10.49</td><td>4.67 ± 2.02</td><td>82.74 ± 11.36</td></tr><tr><td>J</td><td>60565571749760</td><td>0.13 ± 0.24</td><td>-52 ± 6</td><td>-3.900 ± 0.373</td><td>-1.054 ± 0.292</td><td>159.03 ± 214.99</td><td>38.04 ± 125.88</td><td>-17.10 ± 45.78</td><td>159.94 ± 213.81</td></tr><tr><td>N</td><td>60565571767552</td><td>0.25 ± 0.10</td><td>-37 ± 6</td><td>0.092 ± 0.148</td><td>0.134 ± 0.121</td><td>28.12 ± 4.02</td><td>-21.18 ± 5.41</td><td>8.71 ± 2.25</td><td>29.44 ± 3.90</td></tr><tr><td>O</td><td>60569875457792</td><td>1.17 ± 0.05</td><td>-22 ± 7</td><td>2.607 ± 0.098</td><td>2.108 ± 0.076</td><td>12.57 ± 3.57</td><td>-13.00 ± 6.07</td><td>14.12 ± 0.47</td><td>18.90 ± 2.40</td></tr><tr><td>P1</td><td>60565571767424</td><td>0.17 ± 0.09</td><td>-43 ± 10</td><td>-0.889 ± 0.139</td><td>-0.389 ± 0.106</td><td>55.13 ± 13.04</td><td>-11.71 ± 11.20</td><td>-2.17 ± 6.32</td><td>55.17 ± 13.03</td></tr><tr><td>U</td><td>59092406721280</td><td>0.50 ± 0.07</td><td>-45 ± 4</td><td>-1.877 ± 0.113</td><td>-5.096 ± 0.083</td><td>52.70 ± 3.32</td><td>-14.82 ± 3.86</td><td>-39.77 ± 6.62</td><td>66.02 ± 4.78</td></tr></tbody></table>

Note. The proper motions here are in the Gaia system.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Left upper panel: Toomre diagram for a sample of thin disk, thick disk, and transition thin-thick disk stars, covering a wide range of metallicities, with our stars (from Table 4) superimposed (red dots correspond to thin disk stars, green to transition, and blue to thick disk stars). Left lower panel: same as the upper panel, keeping only stars with metallicities equal to or higher than that of star G. Right panel: detail of the lower left panel, leaving out star J. The sample is taken from Adibekyan et al. (2012).</figcaption></figure>

this group although, as we will see, it includes only a small fraction of the thin disk stars.

Quantitatively, in the sample from Adibekyan et al. (2012), of 1111 FGK dwarf stars, there are 601 thin disk stars with metallicities [Fe/H] > −0.14. A total of 446 of them (74.2%) are inside the circle $V^2 + (U^2 + W^2)^{1/2} < 50 \text{ km s}^{-1}$, and 596 are inside the <100 km s⁻¹ circle. Only five (0.8%) have velocities higher than 100 km s⁻¹. That, therefore, is the probability (0.8%), from kinematics alone, that star G is just a fast-moving thin disk star.

In Section 6 the orbits of the stars are discussed and the question of the detailed chemical abundances of star G addressed.

## 6. Stars' Orbits

Using the Gaia DR2 data system, we can calculate the orbits of the stars in the Galaxy. With known distances, proper motions can be translated into tangential velocities. From these and from the radial velocities already obtained, the total velocities of the targeted stars are reconstructed and their orbits as they move across the Galaxy can then be calculated.

The stellar orbits are obtained by integration of the equations of motion. A 3D potential of the Galaxy is required for this. Here we use an axysymmetric potential consisting of a spherical central bulge, a disk, and a massive spherical halo, developed by Allen & Santillán (1991). This is an analytic potential that allows an efficient and accurate numerical computation of the

orbits. In the present case, the total mass of the Galaxy is assumed to be $9 \times 10^{11} M_{\odot}$. We take the Sun to be at 8.5 kpc from the Galactic center and moving circularly with a frequency $\omega_{\odot} = 25.88 \text{ km s}^{-1} \text{ kpc}^{-1}$. We do not consider a Galactic bar or a spiral arms potential.

In Figure 4 we show the orbits of stars B, G, F, and U. We see that only stars G and U reach large distances above and below the Galactic plane, while the other two stars, in contrast, do not appreciably leave it. We also note, in the motion parallel to the Galactic plane, the large eccentricity of the orbit of star G.

The figure is meant to show how far from the Galactic plane the stars with significant proper motions in $\mu_{\alpha}$ and $\mu_{\delta}$ reach. We take four stars with distance compatible with that of the SN. We see that star G will reach up to 500 pc and star U up to almost 700 pc above the Galactic plane within the next 500 Myr. In contrast, B and F (which have an insignificant $\mu_{\delta}$) do not reach 200 pc in their orbits in any lapse of time. We have shown B and F because those have $\mu_{\alpha} \sim 4 \text{ mas yr}^{-1}$ but negligible $\mu_{\delta}$. Their orbits do not look peculiar. These are example of the many stars in a similar situation, which can be seen in our tables. They will be thin disk stars (as seen in the Toomre diagram). The case of star U is unique, in the sense that it has a slightly larger $\mu_{\delta}$ than G. Star U has a negligible $\mu_{\alpha}$, which makes its orbit very circular. Star G has about the same proper motions in $\mu_{\alpha}$ and in $\mu_{\delta}$. This is why it reaches 500 pc above the Galactic plane but, at the same time, unlike star U, its orbit is eccentric.
<figure><img src="image_6.png" /></figure>

**Figure 4.** Orbits of stars B (green), G (red), F (blue), and U (gray), projected on the Galactic meridian plane (left) and on the Galactic plane (right), computed forward in time for the next 500 Myr. The common starting point is marked with a blue square. In the left panel we see that star U reaches the largest distance from the Galactic plane, followed by star G (which corresponds to the respective values of the W component of their velocities in Table 4), while stars B and F scarcely depart from the plane. The behavior of the latter stars is typical of the rest of the sample considered here. In the right panel we see that the orbit of star G on the Galactic plane is highly eccentric (which corresponds to the high value of the U component of its velocity in Table 4), while the other stars (including star U) have orbits close to circular. Also here, the behavior of stars B and F is representative of the whole sample.

<figure><img src="image_7.png" /></figure>

The total velocity of star G is larger than that of star U. This can already be seen from the orbit and more explicitly in the Toomre diagram. When we add the radial velocity vector to obtain the total velocity for star G, we have a $v_r = -87.40 \text{ km s}^{-1}$ (heliocentric), which is larger than that of star U ($-45.40 \text{ km s}^{-1}$). Thus the total velocity for star G is $103.69 \text{ km s}^{-1}$ while for U it is only $68.63 \text{ km s}^{-1}$.

In Figure 5, we have made a histogram of the proper motions in decl. of the stars at distances between $1.7 < d < 3.7$ within $1^\circ$ from the geometrical center of the SN. We see there that the position of star G is anomalous. We also see that around $1^\circ$ in the SN position the proper motions in decl. are small.

In Figure 6 we show the position of the stars compatible with the distance of SN 1572 in proper motions in R.A. and decl. All stars are those from the Gaia DR2 within $1^\circ$ from the SN 1572 center and $1.7 < d < 3.7$ kpc. It is worth noting the logarithmic scale. Basically within $1^\circ$ of the center of SN 1572 and those distances the large majority of stars have very low proper motions.

One might ask what happens with the stars, in the full sample of Tables 2 and 3, that do not have measured heliocentric radial velocities because they were far from the 15% of the radius of the remnant explored in RL04. The Gaia DR2 data show no significant proper motions for any of them.

## 7. Candidate Stars

In order to evaluate the likelihood that a given star is the companion of the SN, we look at the distances provided by Gaia parallaxes and to the proper motions. For some of the stars, we also have radial velocities, obtained from high-resolution spectra. Parallaxes and proper motions allow us to discern whether a given star has received extra momentum from the disruption of the binary system to which it belonged.

<figure><img src="image_8.png" /></figure>

**Figure 5.** Histogram of the distribution in $\mu_\delta$ (in mas yr$^{-1}$) of the stars within $1^\circ$ of the geometrical center of Tycho's SNR in the range of distance compatible with SN 1572 ($1.7 < d < 3.7$ kpc). The data are obtained from Gaia DR2. The red vertical line shows the $\mu_\delta$ of star G.

Such extra momentum (relative to the motion of the center of mass of the binary) mostly comes from the orbital motion of the star when the system was disrupted by the SN explosion. There is also some kick due to the impact of the SN ejecta on the companion, but it is comparatively small and depends on orbital separation and on how compact the companion is. In the hydrodynamic simulations of Marietta et al. (2000), the momentum gained from the kick ranges from 12% to 50% of the momentum the star had before explosion, for main-sequence and subgiant companions, the kick being much
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Proper motions of the candidate stars (see the first column of Table 5) to companion of SN 1572 plotted over the distribution of all the stars with distance compatible with that of the SNR ($1.7 < d < 3.7$ kpc) and included within a radius of 1° of the remnant.</figcaption></figure>

smaller in the case of red giants. Similar values are found by Pan et al. (2014).

We examine first the stars closest to the geometrical center of the remnant.

### Stars A to G.

Star A is the closest to the geometrical center of the SNR. From its stellar atmosphere parameters, this star is at the foreground of the SN. The *Gaia* parallax places it at $0.97_{-0.04}^{+0.05}$ kpc. By comparison of the absolute magnitude corresponding to the stellar parameters with the photometric data, we had derived $d = 1.1 \pm 0.3$ kpc (RL04, B14).

Star B has recently been thought to be a foreground star (Kerzendorf et al. 2018). It is a hot star close to the geometrical center of the remnant. *Gaia* DR2 places this star at a distance compatible with that of the SNR. $d = 2.03_{-0.15}^{+0.19}$ kpc, is now compatible with the SN distance. We already pointed out, in RL04 and B14, that the star likely was at a distance compatible with that of the SNR. We can discard star B, however, on the basis of having no peculiar proper motions nor radial velocity. We have made a reconstruction of the orbit of star B and it moves on the Galactic plane without any disturbance toward upper or lower Galactic latitudes.

Stars C1, C2, and C3 have been observed by the *HST*. *Gaia* could only observe C1, and has determined it to be a very nearby star, at a distance of $0.18_{-0.01}^{+0.03}$ kpc. From the stellar parameters, a distance of $0.75 \pm 0.5$ kpc was estimated. For C2 and C3, RL04, and B14 could not estimate any distance, due to their faint magnitudes. The proper motion values for these three stars obtained from *HST* images in B14 show that they are moving close to the Galactic plane.

Star D is also very nearby. It is at a distance of $0.62_{-0.11}^{+0.15}$ kpc according to *Gaia*. We had calculated a distance of $d = 0.8 \pm 0.2$ kpc. Most of the targeted stars close to the geometrical center of the SNR are at distances below 1.5 kpc.

Star E, though, is at a very large distance. *Gaia* indicates a $d = 7.22_{-4.43}^{+xxx}$ kpc (the upper limit corresponds to a negative parallax). Star E was suggested as the SN companion by Ihara et al. (2007). The authors detected absorption lines in the blue

side, of the spectrum, at 3720 Å, consistent with Fe absorption from Tycho's SNR. They concluded that this might either be due to the Fe I in the SN ejecta or to a peculiarity of the star. In GH09 it was pointed out that it is likely a very distant star. In fact, both its distance and its kinematics exclude it as a possible companion of the SN.

Star F is compatible with the distance to Tycho's SN. *Gaia* measures a distance of $d = 2.15_{-0.32}^{+0.44}$ kpc. We had calculated a distance $d = 1.5 \pm 0.5$ kpc. It is not moving at high radial velocity, nor does it have a high proper motion perpendicular to the Galactic plane. Its orbit does not depart from the Galactic plane.

Star G can also be considered among those close the center of the SNR. *Gaia* has measured a distance $d = 1.95_{-0.35}^{+0.6}$ kpc, which is within the range of distances suggested for the Tycho SN. This star has been a proposed companion in RL04, GH09, and B14. Its kinematics corresponds to that of a thick disk star, but it has a thin disk composition. An enhanced [Ni/Fe] was found in GH09 and questioned in K13. A new calculation was done by B14, which still shows a value [Ni/Fe] larger than solar. We leave this point aside and refer only to the agreed solar metallicity of the star.

All the previous stars have measured $v_r$. They have been placed in the Toomre diagram (Figure 3 and Table 4).

### Proposed stars at the NW of the geometrical center.

Xue & Schaefer (2015) place the site of the explosion of Tycho's SN at the NW of the geometrical center of the SNR. They base their claim, in part, on a reconstruction of the historical center using observations of astronomers that wrote records on SN 1572 in the year of its discovery. Their position is at odds with a previous historically based reconstruction of the location in the sky of SN 1572 from Stephenson & Clark (1977). On the other hand, Xue & Schaefer (2015) use a substitute of a 2D hydrodynamical simulation which, as noted by Williams et al. (2016), would only be valid for perfectly spherical remnants.

Xue & Schaefer (2015) give as the position of the explosion site R.A. = 00°25′15″.36 and decl. = 64°08′40″.2. From this, they suggest that the companion star should be in a small circle around stars O and R. In their Figure 2 they point to stars O, Q, S, and R. From the *Gaia* DR2 data, we see that stars O, Q, and S are at too short distances, incompatible with the distance to Tycho's SN. Star O is at a distance of $0.85_{-0.22}^{+0.54}$ kpc. Star Q is at a distance of $1.51_{-0.51}^{+1.53}$ kpc, only the upper limit being compatible with the distance to the SNR, but it has no high proper motions. Star S is also at a small distance of $0.81_{-0.20}^{+0.41}$ kpc. Star R has no distance measurement or proper motions in the *Gaia* DR2. However, we have proper motions measured with the *HST* (B14) and they are small. Star P is, as seen in the images of the *HST*, two stars P1 and P2. In the *Gaia* DR2 we have only the P1 parallax, though not very well determined but no proper motions are given. In B14, though, we have the proper motions of the stars and they are relatively small. P1 has a $v_r$ of $-43 \pm 10$ and it is possible to place it in the Toomre diagram (Figure 3 and Table 4) where it shows no kinematic peculiarity. Star N is in the same area of the sky as stars O, P1 and P2, Q, R, and S. It has small proper motions and it is possible to place it in the Toomre diagram, where it lies in the region of low kinematical values. In that corner of the sky, as suggested by Xue & Schaefer (2015), there is no star looking like a companion of Tycho's SN in any way.
### The NE proposed center.

In a recent paper by Williams et al. (2016), the expansion center of the remnant is suggested to be at the northeast of the geometrical center. These authors use an extrapolation of the trajectories of different regions of the SNR, but also a 2D hydrodynamical simulation of the expansion of the ejecta in an inhomogeneous medium. They also assume cylindrical symmetry in the initial ejection of the supernova material. However, Krause et al. (2008), from the spectrum of the light echo of SN 1572, suggest that the explosion was aspherical and thus not cylindrically symmetric.

Their suggested center is at R.A. = 00°25′22″6 and decl. = 64°08′32″7. This would be close to stars L and K. Star L is at a distance $d = 1.45_{-0.58}^{+2.87}$ kpc, but with small proper motions. Star K has no distance determined by Gaia. In B14 we suggest it to be around 4 Kpc. The kinematics of the star, with small proper motions, makes it a non-suitable companion of Tycho's SN. In the northeast center there are other stars like W, AK, AL, and AM. These stars do not have accurate parallaxes in the Gaia DR2 release. They have upper limits in distance with negative parallaxes. We have looked at their kinematics and they have moderate and low proper motions.

Therefore, given the various candidates proposed, the best approach is to look for those that are within the range of possible distances to Tycho's SN, show a peculiar kinematics, and are within the region of the sky already explored.

We show in Table 5 the different stars and evaluated their viability as possible companions.

## 8. Luminosities and Models

There are significant differences in the predictions of the characteristics of the surviving companions of the supernova explosion. Podsiadlowski (2003) found that, for a subgiant companion, the object ~400 yr after the explosion might be either significantly overluminous or underluminous, relative to its pre-SN luminosity, depending on the amount of heating and the amount of mass stripped by the impact of the SN ejecta. More recently Shappee et al. (2013) have also followed the evolution of luminosity for years after the impact of the ejecta on a main-sequence companion. The models first rise in temperature and luminosity, peaking at $10^4 L_\odot$ to start cooling and dimming down to $10 L_\odot$ some $10^4$ yr after the explosion. Around 500 days after explosion the companion luminosity would be $10^3 L_\odot$. Pan et al. (2012a, 2012b, 2014) criticize the two preceding approaches for the arbitrariness of the initial models. Starting from their hydrodynamic 3D models, they find lower luminosities for the companions than the previous authors. They find luminosities of the order of only $10 L_\odot$ for the companions, several hundred days after the explosion.

Now, knowing the distances from Gaia, we can derive the luminosities of the stars compatible with being inside the SNR. We take the distance to the SNR from various reliable measurements, which means a value $1.7 < d < 3.7$ kpc. We have 15 candidates compatible with that distance. We already had UV photometry for some of them and now Gaia photometry for all. From these we find that there is no clearly overluminous candidate.

It has been suggested that, within the double-degenerate channel to produce SNe Ia, the explosion can be triggered just at the beginning of the coalescence process of the two WDs, by detonation of a thin helium layer coming from the surface of the less massive one. That would induce a second detonation in

the core of the more massive WD. This hypothetical process has been dubbed the “dynamically driven, double-degenerate, double-detonation scenario” (see Shen et al. 2018 and references therein). In this case, the less massive WD would survive the SN explosion and be ejected at the very high orbital velocity (>1000 km s⁻¹) it had at the moment of the explosion. Those would be seen as “hypervelocity WDs” (Shen et al. 2018). The number of hypervelocity WDs detectable by Gaia depends on the assumed luminosity of these objects. Shen et al. (2018) conclude that, taking into account tidal heating undergone by the WD before the explosion, a typical object would have, after subsequently cooling for ~$10^6$ yr, a luminosity $\ge 0.1 L_\odot$, and thus be detectable by Gaia up to a distance of 1 kpc. Based on that, they predict that ~30 potentially detectable hypervelocity WDs should be found within 1 kpc from the Sun. They have actually found, from Gaia DR2, three objects which, after having been followed up with ground-based telescopes, although not looking like typical WDs, might be the result of heating and bloating of a SN Ia WD companion.

In the case of Tycho's SN, the cooling time of a possible surviving WD companion is only ~450 yr, and thus the luminosity should be significantly higher than the $0.1 L_\odot$ adopted by Shen et al. (2018) for a typical companion having cooled for ~$10^6$ yr.

In order to look for a possible hypervelocity WD companion to Tycho, we must considerably enlarge the search area around the center of the SNR. Taking as an upper limit a velocity perpendicular to the line of sight of 4000 km s⁻¹, the maximum distance traveled in 450 yr, $5.7 \times 10^{13}$ km, translates, at a distance of the SNR, into an angular displacement of 2.1 arcmin (slightly more than 50% of the average radius of the SNR, which is about 4 arcmin).

We have checked that there is no object with unusually high proper motion in the Gaia DR2 data release, within the searched area and up to a G magnitude of 20.7 (V ~ 22). For an extinction $A_V = 2.4$ mag (GH09), and at the distance of Tycho, that means a luminosity $L \sim 0.3 L_\odot$, similar to the lower limit adopted by Shen et al. (2018). That does not take into account the capture of radioactive material by the companion WD predicted by Shen & Schwab (2017). Objects such as the three candidates for hypervelocity former SN Ia companions found by Shen et al. (2018), with G magnitudes ~17–18 mag, would be clearly seen.

## 9. Summary and Conclusions

We have reexamined the distances and proper motions of stars close to the center of Tycho's SNR, using data provided by Gaia DR2. Previously, the distances were only known from determination of the stellar atmosphere parameters and comparison of the corresponding luminosities with the observed apparent magnitudes, with only an approximate knowledge of the extinction and uncertainty about the luminosity classes in a number of cases. More accurate were the proper motions, coming from astrometry made with the HST, but the DR2 has allowed a cross-check here. Besides, only a precise knowledge of the distances allows us to reliably convert proper motions into tangential velocities. Gaia now provides the last word on the distances and kinematics of the previously proposed companions of Tycho's SN.

A good agreement between the distances from Gaia DR2 and those reported in B14 has been found in many cases, but with a general trend to shorter Gaia distances as compared with B14, which can be attributed to an underestimate of
<table><caption>**Table 5**
Criteria Satisfied by the Stars in Table 1</caption><thead><tr><th>Star<br>(1)</th><th>1.7 kpc &le; d &le; 3.7 kpc<br>(2)</th><th>High Velocity<br>(3)</th><th>High v<sub>b</sub><br>(4)</th></tr></thead><tbody><tr><td>A</td><td>No (d = 0.97<sup>+0.05</sup><sub>-0.04</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>B</td><td>Yes (d = 2.03<sup>+0.19</sup><sub>-0.15</sub> kpc)</td><td>No (v = 72 &plusmn; 5 km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -0.5 &plusmn; 0.6 km s<sup>-1</sup>)</td></tr><tr><td>C1</td><td>No (d = 0.18<sup>+0.03</sup><sub>-0.01</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>D</td><td>No (d = 0.62<sup>+0.15</sup><sub>-0.11</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>E</td><td>No (d = 7.22<sup>+xx</sup><sub>-4.43</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>F</td><td>Yes (d = 2.15<sup>+0.44</sup><sub>-0.32</sub> kpc)</td><td>Yes (v = 83 &plusmn; 10 km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = 3 &plusmn; 1 km s<sup>-1</sup>)</td></tr><tr><td>G</td><td>Yes (d = 1.95<sup>+0.60</sup><sub>-0.35</sub> kpc)</td><td>Yes (v = 103 &plusmn; 7 km s<sup>-1</sup>)</td><td>Yes (v<sub>b</sub> = -3.3 &plusmn; 9 km s<sup>-1</sup>)</td></tr><tr><td>H</td><td>Yes (d = 1.61<sup>+0.79</sup><sub>-0.40</sub> kpc)</td><td>Yes (v = 91 &plusmn; 11 km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -0.5 &plusmn; 2.5 km s<sup>-1</sup>)</td></tr><tr><td>I</td><td>...</td><td>...</td><td>...</td></tr><tr><td>J</td><td>No (d = 7.46<sup>+xx</sup><sub>-4.77</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>K</td><td>...</td><td>...</td><td>...</td></tr><tr><td>L</td><td>Yes (d = 1.45<sup>+2.87</sup><sub>-0.58</sub> kpc)</td><td>No (v<sub>t</sub> = 10<sup>+20</sup><sub>-6</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -2 &plusmn; 8 km s<sup>-1</sup>)</td></tr><tr><td>M</td><td>...</td><td>...</td><td>...</td></tr><tr><td>N</td><td>Yes (d = 4.06<sup>+2.57</sup><sub>-1.14</sub> kpc)</td><td>No (v = 36 &plusmn; 4 km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = 2 &plusmn; 3 km s<sup>-1</sup>)</td></tr><tr><td>O</td><td>No (d = 0.85<sup>+0.54</sup><sub>-0.22</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>P1</td><td>No (d = 5.96<sup>+2.23</sup><sub>-2.11</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>Q</td><td>Yes (d = 1.51<sup>+1.53</sup><sub>-0.51</sub> kpc)</td><td>No (v<sub>t</sub> = 24<sup>+24</sup><sub>-8</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = 10 &plusmn; 10 km s<sup>-1</sup>)</td></tr><tr><td>R</td><td>...</td><td>...</td><td>...</td></tr><tr><td>S</td><td>No (d = 0.81<sup>+0.41</sup><sub>-0.20</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>T</td><td>Yes (d = 1.77<sup>+1.86</sup><sub>-0.60</sub> kpc)</td><td>No (v<sub>t</sub> = 30<sup>+32</sup><sub>-10</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = 4 &plusmn; 3 km s<sup>-1</sup>)</td></tr><tr><td>U</td><td>Yes (d = 1.98<sup>+0.32</sup><sub>-0.24</sub> kpc)</td><td>No (v = 68 &plusmn; 5 km s<sup>-1</sup>)</td><td>Yes (v<sub>b</sub> = -46 &plusmn; 8 km s<sup>-1</sup>)</td></tr><tr><td>V</td><td>No (d = 16.81<sup>+xx</sup><sub>-15.89</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>W</td><td>No (d = 5.17<sup>+xx</sup><sub>-3.07</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>X</td><td>No (d = 5.20<sup>+xx</sup><sub>-3.59</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>Y</td><td>Yes (d = 1.58<sup>+0.87</sup><sub>-0.41</sub> kpc)</td><td>No (v<sub>t</sub> = 18<sup>+10</sup><sub>-5</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -17 &plusmn; 9 km s<sup>-1</sup>)</td></tr><tr><td>Z</td><td>No (d = 5.68<sup>+2.72</sup><sub>-2.57</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AA</td><td>Yes (d = 1.04<sup>+1.00</sup><sub>-0.33</sub> kpc)</td><td>No (v<sub>t</sub> = 4<sup>+8</sup><sub>-7</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -5 &plusmn; 6 km s<sup>-1</sup>)</td></tr><tr><td>AB</td><td>...</td><td>...</td><td>...</td></tr><tr><td>AC</td><td>Yes (d = 2.04<sup>+0.99</sup><sub>-0.50</sub> kpc)</td><td>No (v<sub>t</sub> = 11<sup>+6</sup><sub>-10</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -12 &plusmn; 7 km s<sup>-1</sup>)</td></tr><tr><td>AE</td><td>Yes (d = 3.58<sup>+5.84</sup><sub>-1.37</sub> kpc)</td><td>No (v<sub>t</sub> = 16<sup>+13</sup><sub>-5</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = 6<sup>+20</sup><sub>-1</sub> km s<sup>-1</sup>)</td></tr><tr><td>AF</td><td>No (d = 0.76<sup>+0.24</sup><sub>-0.15</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AG</td><td>Yes (d = 1.42<sup>+1.69</sup><sub>-0.50</sub> kpc)</td><td>No (v<sub>t</sub> = 10<sup>+12</sup><sub>-5</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -6 &plusmn; 6 km s<sup>-1</sup>)</td></tr><tr><td>AH</td><td>No (d = 4.85<sup>+3.55</sup><sub>-1.44</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AI1/HP-1</td><td>No (d = 0.35<sup>+0.04</sup><sub>-0.03</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AI2</td><td>...</td><td>...</td><td>...</td></tr><tr><td>AJ</td><td>No (d = 5.35<sup>+xx</sup><sub>-2.81</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AK</td><td>...</td><td>...</td><td>...</td></tr><tr><td>AL</td><td>No (d = 2.61<sup>+xx</sup><sub>-1.61</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AM</td><td>No (d = 1.33<sup>+xx</sup><sub>-0.70</sub> kpc)</td><td>...</td><td>...</td></tr><tr><td>AN</td><td>Yes (d = 1.65<sup>+0.49</sup><sub>-0.30</sub> kpc)</td><td>No (v<sub>t</sub> = 24<sup>+7</sup><sub>-4</sub> km s<sup>-1</sup>)</td><td>No (v<sub>b</sub> = -7 &plusmn; 3 km s<sup>-1</sup>)</td></tr></tbody></table>

Note. *v* refers to the total velocity, *v*<sub>t</sub> to the tangential velocity, and *v*<sub>b</sub> to the velocity perpendicular to the Galactic plane.

extinction in the direction of the remnant in B14. In a few cases, however, the discrepancies are large.

Concerning proper motions, the agreement is very good once due account is made of the systematic effect of the motion of the local frame to which the HST measurements are referred with respect to Gaia's absolute frame.

We find that, within the remaining uncertainties, up to 15 stars are at distances compatible with that of the SNR. The case for Tycho G is that in samples such as the one shown in Figure 3, this star has a thick disk kinematics, but has thin disk metallicity. There is only 0.8% of the stars having similar characteristics. We have inspected the proper motions of all the stars visible up to limit of Gaia DR2 in magnitude, and we have found none with the same peculiar total velocity. We have presented a scenario in which Tycho G could be the companion

of SN 1572. There is, however, the possibility that, after performing several orbits around the Galactic center, and encountering globular clusters and spiral arms, the star orbit becomes eccentric and migrates toward higher Galactic latitudes. This is a suggested explanation for the characteristics of Tycho G. A counterargument is why the other stars that could also have performed several Galactic orbits, at close locations, would not have migrated.

We agree with Kerzendorf et al. (2018) that Tycho B is not a good candidate for a companion of the explosion. We can also exclude, in view of the Gaia DR2 data, that star E could be a companion, since it lies very far away.

In the case that Tycho G were not the companion star, the double-degenerate scenario or the core degenerate scenario are favored, since we have gone well below solar luminosities.
With Gaia DR2, we have also looked for the hypervelocity
stars predicted by some scenarios, but within the magnitudes
reached by Gaia we have found none.

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 DPAC has been provided
by national institutions, in particular the institutions participat-
ing in the Gaia Multilateral Agreement. P.R.-L. is supported
by AYA2015-67854-P from the Ministry of Industry, Science
and Innovation of Spain and the FEDER funds. J.I.G.H.
acknowledges financial support from the Spanish MINECO
(Ministry of Economy of Spain) under the 2013 Ramón y Cajal
program MINECO RyC-2013-14875, also from the MINECO
AYA2014-56359-P and AYA2017-86389-P. This work was
supported as well by the MINECO through grant ESP2016-
80079-C2-1-R (MINECO/FEDER, UE) and ESP2014-55996-
C2-1-R (MINECO/FEDER, UE) and MDM-2014-0369 of
ICCUB (Unidad de Excelencia 'Maria de Maeztu).

ORCID iDs

Pilar Ruiz-Lapuente https://orcid.org/0000-0002-0615-4030

References

Adibekyan, V. Zh., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32

Albinson, J. S., Tuffs, R. J., Swinban, K. E., & Gull, S. F. 1986, MNRAS, 219, 427

Allen, C., & Santillán, A. 1991, RMxA&A, 22, 255

Bedin, L. R., Ruiz-Lapuente, P., González Hernández, J. I., et al. 2014, MNRAS, 439, 354 (B14)

Brown, A. G. A., Vallenari, A., Prusti, T., et al. 2018, A&A, 616, A1

Chevalier, R. A., Kirshner, R. P., & Raymond, J. C. 1980, ApJ, 235, 186

Edwards, Z. I., Pagnotta, A., & Schaefer, B. E. 2012, ApJL, 747, L19

Eggleton, P. P. 1983, ApJ, 268, 368

González Hernández, J. I., Ruiz-Lapuente, P., Filippenko, A. V., et al. 2009, ApJ, 691, 1 (GH09)

González Hernández, J. I., Ruiz-Lapuente, P., Tabernero, H. M., et al. 2012, Natur, 489, 533

Hughes, J. P. 2000, ApJL, 55, L53

Ihara, Y., Ozaki, J., Doi, M., et al. 2007, PASJ, 59, 811

Kerzendorf, W., Childress, M., Scharwächter, J., Do, T., & Schmidt, B. P. 2014, ApJ, 782, 27

Kerzendorf, W. E., Long, K. S., Winkler, P. F., & Do, T. 2018, MNRAS, 479, 5696

Kerzendorf, W. E., Schmidt, B. P., Asplund, M., et al. 2009, ApJ, 701, 1665

Kerzendorf, W. E., Schmidt, B. P., Laird, J. B., Podsiadlowski, P., & Bessell, M. S. 2012, ApJ, 759, 7

Kerzendorf, W. E., Yong, D., Schmidt, B. P., et al. 2013, ApJ, 774, 99 (K13)

Kirshner, R. P., Winkler, P. F., & Chevalier, R. A. 1987, ApJL, 315, L135

Krause, O., Tanaka, M., Usuda, T., et al. 2008, Natur, 456, 617

Kurucz, R. L. 1993, ATLAS9 Stellar Atmospheres Programs, Grids of Model Atmospheres and Line Data [CD-ROM] (Cambridge, MA: Smithsonian Astrophysical Observatory)

Lejeune, Th., Cuisinier, F., & Buser, R. 1997, A&AS, 125, 229

Maoz, D., Mannucci, F., & Nelemans, G. 2014, ARA&A, 52, 107

Marietta, E., Burrows, A., & Fryxell, B. 2000, ApJS, 128, 615

Pagnotta, A., & Schaefer, B. E. 2015, ApJ, 799, 101

Pan, K.-C., Ricker, P. M., & Taam, R. E. 2012a, ApJ, 750, 151

Pan, K.-C., Ricker, P. M., & Taam, R. E. 2012b, ApJ, 760, 21

Pan, K.-C., Ricker, P. M., & Taam, R. E. 2014, ApJ, 792, 71

Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565

Podsiadlowski, P. 2003, arXiv:astro-ph/0303660

Reynoso, E. M., Moffett, D. A., Goss, W. M., et al. 1997, ApJ, 491, 816

Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009

Robin, A. C., Reylé, C., Derrière, S., et al. 2003, A&A, 409, 523

Ruiz-Lapuente, P. 1997, Sci, 276, 1813

Ruiz-Lapuente, P. 2004, ApJ, 612, 357

Ruiz-Lapuente, P. 2014, NewAR, 62, 15

Ruiz-Lapuente, P., Comerón, F., Méndez, J., et al. 2004, Natur, 431, 1069 (RL04)

Ruiz-Lapuente, P., Damiani, F., Bedin, L. R., et al. 2018, ApJ, 862, 124

Schaefer, B. E., & Pagnotta, A. 2012, Natur, 481, 164

Schmidt-Kaler, Th. 1982, in Landolt-Börnstein, New Ser. VI, Vol. 2b, ed. K.-H. Hellwege (Berlin: Springer), 1

Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829

Shappee, B. J., Kochanek, C. S., & Stanek, K. Z. 2013, ApJ, 765, 150

Shen, K. J., Boubert, D., Gänsicke, B. T., et al. 2018, ApJ, 865, 15

Shen, K. J., & Schwab, J. 2017, ApJ, 834, 180

Stephenson, F. R., & Clark, D. H. 1997, QJRAS, 18, 340

Wang, B., & Han, Zh. 2012, NewAR, 56, 122

Williams, B. J., Chomiuk, L., Hewitt, J. W., et al. 2016, ApJL, 823, L32

Xue, Zh., & Schaefer, B. E. 2015, ApJ, 809, 183