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

HAL
open science

# The impact of improved estimates of radio star astrometric models on the alignment of the Gaia bright reference frame to ICRF3

Susanne Lunz, James M. Anderson, Ming H. Xu, Robert Heinkelmann, Oleg Titov, Jean-François Lestrade, Megan C. Johnson, Fengchun Shu, Wen Chen, Jamie Mccallum, et al.

► To cite this version:

Susanne Lunz, James M. Anderson, Ming H. Xu, Robert Heinkelmann, Oleg Titov, et al. The impact of improved estimates of radio star astrometric models on the alignment of the Gaia bright reference frame to ICRF3. Astronomy & Astrophysics - A&A, 2024, 689, 10.1051/0004-6361/202142081. insu-04822443

HAL Id: insu-04822443

https://insu.hal.science/insu-04822443v1

Submitted on 6 Dec 2024

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The impact of improved estimates of radio star astrometric models on the alignment of the Gaia bright reference frame to ICRF3

Susanne Lunz<sup>1,2 D</sup>, James M. Anderson<sup>1,2,* D</sup>, Ming H. Xu<sup>2,3,4 D</sup>, Robert Heinkelmann<sup>2</sup>, Oleg Titov<sup>5 D</sup>, Jean-François Lestrade<sup>6</sup>, Megan C. Johnson<sup>7</sup>, Fengchun Shu<sup>8</sup>, Wen Chen<sup>9</sup>, Alexey Melnikov<sup>10</sup>, Andrei Mikhailov<sup>10</sup>, Jamie McCallum<sup>10</sup>, Yulia Lopez<sup>10</sup>, Pablo de Vicente Abad<sup>11 D</sup>, and Harald Schuh<sup>1,2 D</sup>

Received 24 August 2021 / Accepted 20 June 2024

## ABSTRACT

**Aims.** We investigated the residual orientation offset and spin between the bright ($G \le 13$ mag) frame of the Gaia Early Data Release 3 (Gaia EDR3) and the third realization of the International Celestial Reference Frame (ICRF3). For this purpose, six rotation parameters (orientation offset and its time derivative, the spin), as well as corrections to the Gaia astrometric model for each star involved, are fitted to the differences in the astrometric models derived from very long baseline interferometry (VLBI) and Gaia. This study aims to find reliable estimates for the rotation parameters between the two frames.

**Methods.** We reprocessed our previous analyses while taking into account the effect of Galactocentric acceleration on the VLBI observations. Furthermore, we replaced VLBI data for 12 stars by improved estimates of models of stellar motion from combining historical data with the new positions, rather than including the new observations directly as single-epoch positions in the analysis of the rotation parameters. Additionally, we replaced the model positions by positions obtained without correcting the calibrator data for source structure whenever possible to better reference the star position to ICRF3. In the same fashion, the VLBI proper motion and parallax were included for two of the stars for the first time, and data for five new stars were added.

**Results.** The iterative solutions for the spin parameters show less scatter in the X component when the new models of stellar motion from VLBI are applied. The mean formal errors of the spin parameters decrease by about 15%, whereas those of the orientation offsets increase by about 15%. Small additional improvements in the mean formal error were achieved by including new VLBI data and by excluding stars that produce offsets in the iterative rotation parameter estimates from the beginning. The orientation offset [$\epsilon_X(T)$, $\epsilon_Y(T)$, $\epsilon_Z(T)$] and the spin [$\omega_X$, $\omega_Y$, $\omega_Z$] of the final baseline solution of this work were found to be (+0.322, +0.228, +0.163)±(0.203, 0.251, 0.155) mas and (+0.034, +0.072, -0.026)±(0.023, 0.025, 0.023) mas yr$^{-1}$. As a consequence, no significant orientation offset of Gaia EDR3 toward ICRF3 is detected; however, the spin $\omega_Y$ is statistically significant at the $3\sigma$ level. The rotation parameters between the Gaia and VLBI frames in the Y direction remain the least well determined in terms of formal errors. The impact of Galactocentric acceleration on the rotation parameter analysis was found to be negligible with the currently available VLBI data. As a result, it was found that the uncorrected bright Gaia frame exhibits a closer alignment with ICRF3 compared to the corrected frame.

**Key words.** techniques: interferometric – astrometry – reference systems – radio continuum: stars

## 1. Introduction

In December 2020 the third data release (Early Data Release 3; EDR3; Gaia Collaboration 2021a) of the Gaia spacecraft (Gaia Collaboration 2016), operated by the European Space Agency

(ESA), was published online. It contains astrometric and photometric data of about 1.8 billion objects in the Milky Way and Local Group and therefore is the largest catalog of this kind to date. Since 01 January 2022, the celestial reference frame of Gaia EDR3, Gaia-CRF3 (Gaia Collaboration 2022), has become the fundamental realization of the International Celestial Reference System (ICRS; Arias et al. 1995) in the optical domain (IAU resolution B3, 2021). The counterpart in the radio domain is the third realization of the International Celestial

<footnote n="1">¹ Technische Universität Berlin, Chair of Satellite Geodesy, Kaiserin-Augusta-Allee 104-106, 10553 Berlin, Germany</footnote>

<footnote n="2">² GFZ German Research Centre for Geosciences, 14473 Potsdam, Germany
e-mail: susanne.lunz@gfz-potsdam.de</footnote>

<footnote n="3">³ Aalto University Metsähovi Radio Observatory, Metsähovintie 114, 02540 Kylmälä, Finland</footnote>

<footnote n="4">⁴ Aalto University Department of Electronics and Nanoengineering, PL15500, 00076 Aalto, Finland</footnote>

<footnote n="5">⁵ Geoscience Australia, PO Box 378, Canberra 2601, Australia</footnote>

<footnote n="6">⁶ LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 75014 Paris, France</footnote>

<footnote n="7">⁷ United States Naval Observatory (USNO), 3450 Massachusetts Ave NW, Washington, DC 20392, USA</footnote>

<footnote n="8">⁸ Shanghai Astronomical Observatory, 80 Nandan Road, Shanghai, 200030, China</footnote>

<footnote n="9">⁹ Yunnan Observatories, Yangfangwang 396th, Kunming 650011, Yunnan, China</footnote>

<footnote n="10">¹⁰ University of Tasmania, Private Bag 37, Hobart, Tasmania, 7001, Australia</footnote>

<footnote n="11">¹¹ Observatorio de Yebes (IGN), Cerro de la Palera S/N, 19141, Guadalajara, Spain</footnote>

<footnote>* Current address: Leibniz Institute for Agricultural Engineering and Bioeconomy, 14469 Potsdam, Germany.</footnote>
Reference Frame (ICRF3; Charlot et al. 2020), which has been
in use since 01 January 2019 (IAU resolution B2, 2018). The
ICRF3 was produced from very long baseline interferometry
(VLBI) observations of extragalactic radio sources at radio
wavelengths (S/X, K, and X/Ka). The Gaia-CRF3 is aligned to
the ICRF3 with an uncertainty of 10 µas at epoch T = 2016.0
and it is nonrotating with respect to the ICRF3 and the All-
WISE (mid-infrared) active galactic nucleus catalog (Secrest
et al. 2015, 2016) with a spin of less than ±10 µas yr⁻<fnref n="1" /> (Gaia
Collaboration 2021a). The alignment of the Gaia catalog to
ICRF3 is known to be G magnitude dependent due to the Gaia
internal calibration strategies, and thus, the above numbers are
only valid for the optical magnitude G = 19 since the coun-
terparts used for calculations lie mostly around this magnitude
(Lindegren et al. 2021b). None of the ICRF3 sources that were
detected by Gaia have optical magnitudes G ≤ 13 mag. As a
result, the spin correction for the Gaia EDR3 bright frame (G ≤
13 mag) was derived by comparing proper motions of bright
sources calculated from the position differences between Gaia
Data Release 2 (DR2; Gaia Collaboration 2018) and HIPPARCOS
(van Leeuwen 2007; ESA 1997), divided by the epoch differ-
ence of 24.25 yr, and from Gaia EDR3. It was then applied to
the Gaia bright frame in EDR3. Thus no residual spin should
be present between the Gaia bright frame and ICRS after the
correction, and therefore also no residual spin should be present
between the Gaia bright frame and the Gaia faint frame. Using
this approach, the orientation offset (to the HIPPARCOS reference
frame) was assumed to be zero (Lindegren et al. 2021b). How-
ever, HIPPARCOS positions and the ICRS at epoch J1991.25 are
aligned to no better than 0.6 mas in each axis (Kovalevsky et al.
1997), which introduces a systematic uncertainty of 24 µas yr⁻<fnref n="1" />
for Gaia EDR3. This is above the expected uncertainties for
bright objects with G ≤ 13 mag, which are predicted to be less
than 6 µas for positions and less than 2 µas yr⁻<fnref n="1" /> for proper
motions in the final Gaia data release<fnref n="1" />. Because the error budget
with this approach cannot be reduced by enhanced Gaia data and
obviously no improvements of the HIPPARCOS positions are to
be expected, the alignment to ICRS needs to be tested by other
methods for the final Gaia data release.

Such an alternative method for the verification of the align-
ment is by using VLBI observations of optically bright radio
stars, as introduced by Lindegren (2020a). Positions, proper
motions, and parallaxes from VLBI are compared to those from
Gaia. From the differences in these astrometric parameters,
residual rotations (orientation offsets εₓ, εᵧ, and ε_z and spins
ωₓ, ωᵧ, and ω_z) about the three orthogonal axes X, Y, and
Z and adjustments to the Gaia 5-parameter astrometric models
for the individual stars are determined in a combined fit. With
this method, not only do the proper motions contribute to the
spin parameter adjustment, but also position offsets at a com-
mon epoch divided by the epoch difference between VLBI and
Gaia positions. In Lunz et al. (2023) this method was adopted
and the dataset provided by Lindegren (2020a) was homogenized
to consistently reference historical VLBI positions of radio stars
to ICRF3 at epoch 2015.0. Furthermore, a realistic error budget
was applied to the positions. The main concern for increasing
the uncertainties was the difference in the observational method
for radio-bright extragalactic objects used to create the ICRF3
(absolute geodetic VLBI based on dual-frequency group delays)
and the observational method for radio stars. The latter radio

stars are bright at optical frequencies but faint at radio frequen-
cies, thus, only the technique of phase referencing, which allows
for the determination of relative star positions with respect to
radio-bright calibrators in the ICRF3 based on phase delays,
can be used. New single-epoch positions for 32 stars, observed
in January 2020 in project UL005 with the Very Long Base-
line Array (VLBA) at X and C bands, were added in Lunz
et al. (2023) to enhance the fit of orientation and spin parame-
ters. However, the analysis showed that the realistic single-epoch
position uncertainties were too large for the new positions to have
a significant effect on the spin determination.

In this study, we replace the data for 12 stars that were
detected in January 2020 with improved models of stellar
motion. These models were derived from a fit of the five astro-
metric parameters to a combined position time series consisting
of data found in literature and the January 2020 position for each
star. In addition, data from five stars recently observed with other
VLBI networks have been incorporated into the analysis.

In ICRF3 the effect of Galactocentric acceleration was
taken into account for the first time in ICRF history (Charlot
et al. 2020). In Gaia EDR3 the Galactocentric acceleration has
also been detected due to reduced internal systematics (Gaia
Collaboration 2021b). Thus, the impact of this effect on the mod-
els of stellar motion, and therefore also on the rotation parameter
analysis, must be addressed.

In Sect. 2 we consider the effect of Galactocentric accelera-
tion and investigate its impact on the phase referencing results.
Section 3 describes the algorithm used to fit models of stellar
motions to VLBI position time series, the reprocessing of the
phase referencing observations from January 2020, and the new
model estimates including comparisons to literature for each of
the respective stars. The choice of the absolute VLBI positions
to be used for comparison with the Gaia positions is also dis-
cussed. Section 4 presents rotation parameter results for various
scenarios that test different aspects of the analysis, including the
impact of Galactocentric acceleration, the impact of single stars
on the analysis, the impact of the new estimates for models of
stellar motion, and the impact of the additional data of the five
new stars. Finally, in Sect. 5 we evaluate the results and in Sect. 6
conclusions are drawn.

2. Impact of Galactocentric acceleration on phase referencing results

To account for the global systematic effect of Galactocentric
acceleration in ICRF3, a time-dependent correction in the shape
of a dipole pattern of magnitude |D| = 5.8 µas yr⁻<fnref n="1" /> in the Galac-
tocentric direction (α_D = 266.4°, δ_D = −29.0°) was added for
the first time to the definition of the ICRF positions, which
are given for the epoch 2015.0 (Charlot et al. 2020). While
for comparison of the VLBI results to Gaia DR2 it was suffi-
cient to neglect this effect because internal systematics of the
Gaia DR2 dataset were well above its magnitude, it is appro-
priate to investigate the impact when comparing to Gaia EDR3.
Using Gaia EDR3, the effect was detected to have a magnitude
|D| = 5.05 ± 0.35 µas yr⁻<fnref n="1" /> in direction α_D = 269.1° ± 5.1°, δ_D =
−31.6° ± 4.1° solely from the proper motions of extragalac-
tic compact objects (Gaia Collaboration 2021b). These results
match the VLBI-based result within the error bounds.

In Lunz et al. (2023) all star positions were homogeneously
referred to the calibrator positions in ICRF3 at epoch 2015.0.
To investigate the effect of Galactocentric acceleration on the
phase referencing results, we distinguish between the correction

<footnote n="1">¹ https://www.cosmos.esa.int/web/gaia/science-performance</footnote>
for absolute positions and for proper motions: The absolute position of a star is determined relative to the catalog position of its phase calibrator in ICRF3. Therefore, if this position accounts for the effect of Galactocentric acceleration, the star position also includes this effect. Considering this relation, the homogenized star positions in ICRF3 at epoch 2015.0 can be corrected for the effect of Galactocentric acceleration (see e.g., Titov et al. 2011; MacMillan et al. 2019) by

$$ 
\begin{aligned}
\Delta\alpha &= \Delta t \cdot (-D_1 \sin \alpha + D_2 \cos \alpha) / \cos \delta, \\
\Delta\delta &= \Delta t \cdot (-D_1 \cos \alpha \sin \delta - D_2 \sin \alpha \sin \delta + D_3 \cos \delta),
\end{aligned}
\quad (1) $$

where $D_1 = |D| \cdot \cos \delta_D \cos \alpha_D$, $D_2 = |D| \cdot \cos \delta_D \sin \alpha_D$, $D_3 = |D| \cdot \sin \delta_D$ and $\Delta t = t_B - 2015.0$, where $t_B$ is the epoch of the star position, and $\alpha$ and $\delta$ are the calibrator coordinates. The correction needs to be added to the star position.

At the same time, the proper motions of the stars from position time series are also systematically biased by not accounting for the effect of Galactocentric acceleration in the calibrator positions. Calibrator positions in phase referencing analysis have never been corrected for the effect of Galactocentric acceleration before, since the catalogs used for deriving positions for the phase calibrators never accounted for this effect prior to ICRF3. Therefore, the historical absolute positions of the stars also do not account for the effect of Galactocentric acceleration. The effect needs to be considered for VLBI results because also the Gaia proper motions include it. The corrections

$$ 
\begin{aligned}
\Delta\mu_{\alpha} &= (-D_1 \sin \alpha + D_2 \cos \alpha) / \cos \delta, \\
\Delta\mu_{\delta} &= -D_1 \cos \alpha \sin \delta - D_2 \sin \alpha \sin \delta + D_3 \cos \delta,
\end{aligned}
\quad (2) $$

need to be added to the star proper motions $\mu_\alpha$ and $\mu_\delta$ to correct for the time-varying calibrator coordinates that should have been used during the data processing with respect to ICRF3. In this work, all $\alpha \cdot \cos(\delta)$ terms are abbreviated by $\alpha^*$.

## 3. Fitting models of stellar motion to multi-epoch star positions

In this section, we first describe the fitting process of the main astrometric parameters (Sect. 3.1). Since the structure corrected absolute positions of the stars to their calibrators are needed for the model-fitting, the corresponding data processing for the January 2020 positions is also outlined (Sect. 3.2). Then, for each star, a detailed description of input positions and estimated model parameters is given along with a comparison to astrometric models from the literature and Gaia EDR3 (Sect. 3.3). Finally, only the proper motions and parallaxes are selected from the newly derived models as input for the computation of the rotation parameters in the following section. To be consistent with ICRF3, the absolute positions as determined in Lunz et al. (2023) are used instead of the fitted model positions. These single-epoch absolute positions are corrected for the parallax effect.

### 3.1. Functional models and fitting procedure

The new estimates for models of stellar motion for each star were determined by fitting its single-epoch coordinates $\alpha(t)$ and $\delta(t)$ at observation epoch $t$ to the functional model as, for example, described in Loinard et al. (2007):

$$ 
\begin{aligned}
\alpha(t) &= \alpha_0 + \mu_{\alpha*}t + \bar{\omega}f_{\alpha}(t), \\
\delta(t) &= \delta_0 + \mu_{\delta}t + \bar{\omega}f_{\delta}(t).
\end{aligned}
\quad (3) $$

The model consists of positions $\alpha_0$ and $\delta_0$ at a reference epoch $t_0$, parallax $\bar{\omega}$, as well as linear proper motion terms $\mu_{\alpha*}$ and $\mu_{\delta}$. The terms $f_\alpha(t)$ and $f_\delta(t)$ are the projections of the parallactic ellipse over $\alpha$ and $\delta$ (Seidelmann 1992), calculated as

$$ 
\begin{aligned}
f_{\alpha}(t) &= [X_E(t) \sin \alpha_1 - Y_E(t) \cos \alpha_1] / \cos \delta_1, \\
f_{\delta}(t) &= X_E(t) \cos \alpha_1 \sin \delta_1 + Y_E(t) \sin \alpha_1 \sin \delta_1 - Z_E(t) \cos \delta_1,
\end{aligned}
\quad (4) $$

where $\alpha_1 = \alpha(t) - \bar{\omega}f_\alpha(t)$, $\delta_1 = \delta(t) - \bar{\omega}f_\delta(t)$, and $(X_E(t), Y_E(t), Z_E(t))$ is the vector between the solar system barycenter and Earth in units of AU. The equations therefore depend on a prior knowledge of $\bar{\omega}$ and need to be solved iteratively. Using Eq. (3) a linear least squares fit was performed to the data. According to Loinard et al. (2007), it gives identical results as if using singular value decomposition. Additional variances $E_{\alpha*}^2$ and $E_{\delta*}^2$ are added to the variances of the individual positions for weighting so that the reduced $\chi^2$ values per coordinate direction approximate unity in the fit, yielding realistic uncertainty estimates for the parameters.

### 3.2. Analysis of phase referencing observations

We used the single-epoch observations of 32 stars detected in phase referencing mode at X and C bands (8.11225 GHz and 4.61175 GHz) with the VLBA in project UL005. The observation setup, the data calibration, and the analysis are described in Lunz et al. (2023). In that work, absolute positions for stars were derived from a fringe fit of the calibrator with a point-like calibrator model for its structure. This scheme is consistent with that used for ICRF3 where source structure was not considered. However, because of systematic errors, including the unknown difference between group-delay and phase delay positions, the error budget of the absolute star position was inflated.

For the present analysis, which aims at fitting models of stellar motions to the position time series, the calibrator fringe fits were repeated by applying self-calibrated images of the phase calibrators, instead of a point source model. The self-calibrated images were derived from the same data. This reanalysis minimizes the potential star position jitter coming from the calibrator structure and yields better consistency with positions from other studies, where the same procedure was applied. The positions of the stars $\alpha(t)$ and $\delta(t)$ at the observation epoch $t$, which is the mean of the first and last scans involved, were determined using the task modelfit in the Caltech Difmap imaging package (Shepherd 1997). Their uncertainties $\sigma_{\alpha*,random}$ and $\sigma_{\delta*,random}$ are due to thermal noise and are derived from the beam shape and the root mean square (RMS) noise of the image based on formulas for elliptical Gaussians in Condon (1997). The new positions are listed in Table 1.

In the January 2020 data, the RMS of the differences in the positions of the 32 stars using one or the other method for fringe-fitting (i.e., based on a point-source model or a self-calibrated image for the calibrator structure) were only on the order of 19 µas in $\alpha*$ and 35 µas in $\delta$. The deviation in coordinates is at most 112 µas in $\delta$ for the star HD 283572. All other coordinate deviations are within ±75 µas. Therefore, the impact of the calibrator structure on the results remains small. The differences in the dynamic range are negligible.

### 3.3. New estimates for models of stellar motion

For all stars in Lunz et al. (2023), we searched the literature for single-epoch positions and added to this list the positions
<table><caption>**Table 1.** Absolute positions in ICRF3 for 14 stars based on observations conducted in January 2020 and analysis including the modeling of the phase calibrator structure.</caption><thead><tr><th>Star</th><th>Calibrator</th><th>Δ<sub>c</sub><br>[°]</th><th>Epoch t<br>(Julian year)</th><th>α (t)<br>[°]</th><th>σ<sub>α*,random</sub><br>(mas)</th><th>δ (t)<br>[°]</th><th>σ<sub>δ,random</sub><br>(mas)</th></tr></thead><tbody><tr><td>HD 283572</td><td>J0429+2724</td><td>2.0</td><td>2020.01547</td><td>65.495256413</td><td>0.033</td><td>28.301662253</td><td>0.072</td></tr><tr><td>V410 Tau</td><td>J0429+2724</td><td>2.7</td><td>2020.01546</td><td>64.629671673</td><td>0.059</td><td>28.454348213</td><td>0.126</td></tr><tr><td>SS Cyg</td><td>J2136+4301</td><td>1.3</td><td>2020.01451</td><td>325.679208442</td><td>0.071</td><td>43.586257890</td><td>0.112</td></tr><tr><td>Brun 334</td><td>J0529-0519</td><td>1.2</td><td>2020.01560</td><td>83.665661454</td><td>0.055</td><td>-5.407113858</td><td>0.134</td></tr><tr><td>TYC 5346-538-1</td><td>J0542-0913</td><td>1.1</td><td>2020.01561</td><td>85.640320721</td><td>0.101</td><td>-8.120884576</td><td>0.248</td></tr><tr><td>Haro 1-6</td><td>J1633-2557</td><td>2.3</td><td>2020.01429</td><td>246.512495679</td><td>0.055</td><td>-24.393536459</td><td>0.148</td></tr><tr><td>CoKu HP Tau G2</td><td>J0438+2153</td><td>1.2</td><td>2020.01549</td><td>68.975719883</td><td>0.061</td><td>22.903662815</td><td>0.134</td></tr><tr><td>BH CVn</td><td>J1324+3622</td><td>2.1</td><td>2020.01370</td><td>203.699807328</td><td>0.020</td><td>37.18235859</td><td>0.050</td></tr><tr><td>BH CVn</td><td>J1340+3754</td><td>1.3</td><td>2020.01370</td><td>203.699807118</td><td>0.019</td><td>37.18235822</td><td>0.046</td></tr><tr><td>σ<sup>2</sup> CrB</td><td>J1613+3412</td><td>0.4</td><td>2020.01416</td><td>243.668436110</td><td>0.018</td><td>33.85812083</td><td>0.031</td></tr><tr><td>HD 199178</td><td>J2102+4702</td><td>3.0</td><td>2020.01452</td><td>313.473754649</td><td>0.019</td><td>44.38640432</td><td>0.037</td></tr><tr><td>AR Lac</td><td>J2153+4322</td><td>3.5</td><td>2020.01460</td><td>332.169653348</td><td>0.029</td><td>45.74250851</td><td>0.046</td></tr><tr><td>IM Peg</td><td>J2253+1608</td><td>0.7</td><td>2020.01469</td><td>343.259318422</td><td>0.019</td><td>16.84104056</td><td>0.041</td></tr><tr><td>HD 22468</td><td>J0339-0146</td><td>2.5</td><td>2020.01545</td><td>54.196857495</td><td>0.016</td><td>0.58685785</td><td>0.038</td></tr><tr><td>del Lib</td><td>J1456-0617</td><td>2.1</td><td>2020.01405</td><td>225.242763738</td><td>0.037</td><td>-8.51898195</td><td>0.082</td></tr></tbody></table>

**Notes.** Star coordinates at epoch *t* are denoted *α*(*t*) and *δ*(*t*), and their uncertainties *σ*<sub>α*,random</sub> and *σ*<sub>δ,random</sub>. The spatial distance between star and phase calibrator is Δ<sub>c</sub>.

we derived in January 2020 (Sect. 3.2). New estimates for models of stellar motion, obtained through the scheme presented in Sect. 3.1, were then computed for each star when the data allowed. This computation was possible for 15 stars.

Each star position was corrected to be referenced to the calibrator position in ICRF3 by accounting for the difference between the calibrator position in ICRF3 and the calibrator position used in the literature. Except where noted, all positions are free from systematic errors due to the calibrator structure, which was corrected as described in Sect. 3.2. Variances based on thermal noise σ<sub>random</sub> were used to weight the individual positions. The new estimates for models of stellar motion obtained in this study are listed in Table A.1 along with estimates found in the literature. The corresponding sky motions are shown in Fig. A.1. The estimates for the models of the five astrometric parameters α₀, δ₀, ϖ, μ<sub>α*</sub>, and μ<sub>δ</sub> derived using only data from the literature are denoted M<sub>old</sub>, while the estimates for the models containing the 2020 positions are labeled M<sub>new</sub>. The table also provides information on the additional uncertainties E<sub>α*</sub> and E<sub>δ</sub>. The RMS of the post-fit residuals *v* is denoted by RMS<sub>v</sub>, specifically RMS<sub>v,α*</sub> and RMS<sub>v,δ</sub> for the two coordinate directions. The epoch of α₀ and δ₀ was chosen as the mean epoch t₀ of all observation epochs *t* to minimize the correlations between the parameters. For the projection of the parallax ellipse in Eq. (4), the barycentric coordinates of the Earth's center of mass at the observation epoch were derived from the DE 421 ephemeris (Folkner et al. 2009).

Because of the presence of systematic errors, special attention should be paid to the combination of star positions referenced to different phase calibrators. For the VLBA observations in 2020, position differences with an RMS of 0.28 mas in α* and 0.69 mas in δ were determined from referencing seven stars to two different phase calibrators (Lunz et al. 2023). To account for this effect, all observations not referenced to a primary phase calibrator have been down-weighted by adding the respective additional uncertainties to the star position uncertainties through variance propagation. A mean value of 0.5 mas in both directions was adopted and assumed to be not too pessimistic in the

presence of different uv-coverage, source structure, differences between group- and phase delay positions, and residual delay model errors in the historical observations. The resulting models are denoted “w0.5”. Details of the data and analysis for each star are discussed in the following subsections.

### 3.3.1. HD 283572

The star HD 283572 (HDE 283572) was observed relative to J0429+2724 at 6 epochs from September 2004 to December 2005 at 8.4 GHz (Torres et al. 2007). These data were reproduced in Galli et al. (2018) using the same scheme as that described in Sect. 3.2. Observations in the same band and using the same primary calibrator were carried out in project UL005. The shifts of -0.1333 mas in α and 0.1413 mas in δ (as given in Table E.1 of Lunz et al. 2023) were subtracted from the single-epoch star positions from Galli et al. (2018) to reference them to ICRF3.

The work of Galli et al. (2018) does not report a model position, whereas Lindegren (2020a) determined a full five-parameter model for the star using the single-epoch positions of Galli et al. (2018). The solution M<sub>old</sub> is not significantly different from that of Galli et al. (2018) – however, the formal errors of M<sub>old</sub> are larger. Adding the 2020 position to the analysis (M<sub>new</sub>) reduces the formal errors by 10% to 20% for α₀, δ₀, and ϖ, and by about 90% for μ<sub>α*</sub> and μ<sub>δ</sub> compared to M<sub>old</sub>. The latter is due to the much longer time span involved for the estimation. The additional uncertainties E<sub>α*</sub> and E<sub>δ</sub> decrease by 10% to 15%, while the values of RMS<sub>v</sub> do not change significantly. The differences in model estimates between M<sub>new</sub> and Galli et al. (2018) are not significant.

### 3.3.2. V410 Tau

The 2020 position of V410 Tau (HD 283518) was added to the 9 epochs listed in Galli et al. (2018). These epochs span from November 2013 to September 2017 and used the same antenna network as the January 2020 observations. The first 5
epochs were observed at 8.4 GHz, whereas the last 4 epochs were observed at 5.0 GHz, consistent with the January 2020 observation at 4.6 GHz. The same phase calibrator, J0429+2724, was used. The shifts of 0.0167 mas in α and 0.1813 mas in δ (as given in Table E.1 of Lunz et al. 2023) were subtracted from the star positions from Galli et al. (2018) to reference them to ICRF3.

The formal errors are not improved between $M_{old}$ and $M_{new}$. The additional observation uncertainties increase from $E_{α*} = 116$ µas and $E_{δ} = 79$ µas to $E_{α*} = 199$ µas and $E_{δ} = 93$ µas when the January 2020 position is included in the analysis, and the scatter of the residuals increases. The residuals at the X band and C band do not differ. Since V410 Tau is a multiple star, this difference and the increase in the residuals may be due to the signature of orbital motion, even though only one component is detected at radio frequencies (Harris et al. 2012). However, we defer the testing of this hypothesis to when further observations (at more epochs) are acquired.

### 3.3.3. SS Cyg

SS Cyg was observed by Miller-Jones et al. (2013) from April 2010 to October 2012 in 7 epochs with the VLBA at 8.4 GHz and in 2 epochs at 5 GHz with the European VLBI Network (EVN). In the publication, the phase calibrator position of the EVN observations was already corrected to match that used by the VLBA. The 2020 position was referenced to the same phase calibrator, J2136+4301, and observed at 8.1 GHz. All star positions used in our analysis were referred to the ICRF3 position of this calibrator by subtracting 0.4542 mas in α and -0.5418 mas in δ to the positions (at 9 epochs) reported by Miller-Jones et al. (2013). As previously, these corrections come from Table E.3 of Lunz et al. (2023).

For $M_{new}$, the formal errors are reduced by 10% to 20% (for α and δ), 40% (for $\bar{\omega}$) and by 70% (for $\mu_{α*}$ and $\mu_{δ}$) compared to $M_{old}$. The additional observation uncertainty $E_{α*}$ is 184 µas, while in the δ direction no uncertainty was added, since the reduced $\chi^2$ was already below 1. For this star, the correlations between all parameters decrease when the new epoch in 2020 is included, especially those between α and δ, and between $\mu_{α*}$ and $\mu_{δ}$. Only the correlation between parallax and proper motion increases slightly. This indicates that further observations sensitive to the parallax determination are still needed to decorrelate those two parameters. This behavior contrasts with that for the other stars where adding the new epoch did not change the correlations significantly. The residuals at the X band and C band do not differ.

### 3.3.4. Brun 334

Brun 334 (VLBA 19, Parenago 1540) was observed by Kounkel et al. (2017) using the VLBA in 5 epochs between March 2014 and March 2016 at 5 GHz. The publication only gives the observation epochs to the precision of a full day. Therefore, the positions could be imprecise by up to 0.005 mas, considering the proper motion in Gaia EDR3 of 3.840 mas yr$^{-1}$. The star was observed in one epoch in January 2020 at 4.6 GHz relative to the same phase calibrator, J0529-0519. The shifts of -0.7522 mas in α and 0.1458 mas in δ from Lunz et al. (2023) were subtracted from the star positions from Kounkel et al. (2017) to reference them to ICRF3.

When the January 2020 position is included in the analysis, the formal errors of $\mu_{α*}$ and $\mu_{δ}$ improve by 60% and 40%, respectively, while those of the position and parallax degradation, especially for the δ coordinate, which has a formal error 50%

higher. Additionally, the residual scatter in each coordinate direction doubles and the additive variances (mostly in δ) need to be increased to obtain $\chi^2$ equal one. Although Brun 334 is a spectroscopic binary (Marschall & Mathieu 1988), the short-term orbital parameters, which may explain the degradation, were not considered in this work. Comparing our model estimates with those from Gaia EDR3 shows good agreement, except for $\mu_{δ}$ where the difference is at the $2\sigma$ level.

### 3.3.5. TYC 5346-538-1

Star TYC 5346-538-1 (VLBA 45) was observed in 5 epochs between March 2014 and March 2016 by Kounkel et al. (2017) at 5 GHz. The new observation in January 2020 was conducted at 4.6 GHz and the same phase calibrator as in the historical observations, J0542-0913, was used. The shifts of -0.9450 mas in α and -1.0675 mas in δ were subtracted to the star positions from Kounkel et al. (2017) to reference them to ICRF3. In addition, Kounkel et al. (2017) applied corrections of 0.256 mas in α and 0.771 mas in δ to the positions of the first epoch and of 0.204 mas in α and 0.659 mas in δ to the positions of the fourth epoch. It is believed that the reason is a pointing error due to calibration. These corrections were also applied in our analysis to obtain positions consistent with those shown in Fig. 4 in Kounkel et al. (2017).

The reduced $\chi^2$ in the α direction is already 0.07 without adding additional weights to the data in model $M_{old}$. The formal errors for $\mu_{α*}$ and $\mu_{δ}$ in $M_{new}$ decrease by about 40% compared to the formal errors in $M_{old}$, whereas those for $α*$, $δ$, and $\bar{\omega}$ increase by 50% to 70%. $M_{old}$ and $M_{new}$ differ by (0.299 ± 0.142) mas yr$^{-1}$ in $\mu_{δ}$, which is about 60% of its value, and by (-0.102 ± 0.049) mas yr$^{-1}$ in $\mu_{α*}$. Both parameters are closer to the Gaia EDR3 solution when $M_{new}$ is considered. $M_{old,orig}$ and $M_{new,orig}$ are solutions without applying the above corrections from Kounkel et al. (2017). The corresponding model estimates agree with those from solutions $M_{old}$, $M_{new}$, and Gaia EDR3 within the uncertainty bounds, except for $\mu_{δ}$ for which the difference between $M_{old}$ and $M_{old,orig}$ is at the $1.4\sigma$ level.

### 3.3.6. BH CVn

Star BH CVn (HR 5110) was observed at 15 epochs from May 1987 to May 1994 relative to J1317+3425 or relative to J1340+3754 in Lestrade et al. (1999). The 6 sessions relative to J1340+3754 (observed after 01 January 1992) were conducted at 8.4 GHz, while the frequency of the observations on J1340+3754 was 5.0 GHz for the four sessions before 01 January 1992 and 8.4 GHz for the sessions thereafter. The star was observed in January 2020 at 8.1 GHz relative to J1340+3754 and J1324+3622. A second reference source was chosen because J1340+3754 has significant structure and is relatively faint considering the sensitivity of the observations at this frequency band. J1324+3622 is closer to the target star than the second historical calibrator J1317+3425, though also relatively faint, and served as a backup calibrator. In Lestrade et al. (1999) the phase calibrator data were corrected for source structure only for J1340+3754 because of the significant structure of this source. Following the scheme described in Sect. 3.2, correction for the structure was applied to both calibrators (i.e., J1340+3754 and J1324+3622) in our 2020 data. The star positions were referenced to the ICRF3 position of the respective calibrators at all epochs. To this end, shifts of -0.0414 mas in α and 0.1481 mas in δ were subtracted from the star positions from Lestrade et al. (1999) referenced
to J1317+3425, while for those referenced to J1340+3754, the subtracted shifts were -0.8662 mas in $\alpha$ and 1.3768 mas in $\delta$.

As with the other stars in this study, several solutions with five astrometric parameters were determined based on three different selections of observations to show the reliability of the final solution. $M_{old}$ and $M_{new}$ are based on the observations relative to all calibrators, whereas $M_{J1340+3754}$ contains only positions relative to calibrator J1340+3754 and $M_{others}$ contains only positions referenced to the other two calibrators. Because $M_{J1340+3754}$ has the least residual scatter, J1340+3754 was used as the primary calibrator in solution $M_{w0.5,new}$, while J1317+3425 and J1324+3622 were considered as secondary calibrators. As noted above, the star positions relative to these two calibrators were thus down-weighted in the analysis. Comparing the various solutions, no significant differences can be found for $\mu_{\alpha*}$. On the other hand, the estimate for parameter $\mu_{\delta}$ varies when the 2020 position is added. The formal error of both parameters decreases by 80% when the 2020 position is added. Looking at specific solutions, the decrease is by more than 90% for solution $M_{J1340+3754}$, but only by 60% to 70% for solution $M_{others}$. Examination of the post-fit residuals, $RMS_{v,\alpha*}$ and $RMS_{v,\delta}$, indicates that even though J1340+3754 suffers from source structure, the individual absolute positions for BH CVn appear to be more consistent when using only this calibrator (and correcting for its structure), as in solution $M_{J1340+3754}$, than when referencing to several more compact calibrators (without consistent calibration of their structure), as in solution $M_{new}$, $M_{w0.5,new}$, and $M_{others}$.

The latter is the solution that determines the parallax worst in terms of formal error. Solution $M_{w0.5,new}$ is not significantly different from $M_{J1340+3754}$ and $M_{new}$, but $M_{J1340+3754}$ clearly has the lowest formal errors in all parameters. Thus, it was adopted for the analysis of the rotation parameters in Sect. 4.3 below, and the observations relative to the two other calibrators are thereby neglected. This is only possible due to the many positions available for BH CVn, allowing for a selection among these, which in general is not the case for the other stars in our sample, where solutions deteriorate if neglecting a subset of the positions.

### 3.3.7. Haro 1-6

Star Haro 1-6 (DoAr 21) was observed in 9 epochs from September 2005 to August 2006 relative to J1625-2527 at 8.42 GHz, and where 7 epochs were analyzed and published in Loinard et al. (2008). Another 8 epochs were conducted at the same frequency from July 2007 to September 2007, also relative to J1625-2527. In total, five of these epochs were discarded due to the influence of bad weather, thus leaving data for 12 epochs (Ortiz-León et al. 2017). Then Ortiz-León et al. (2017) observed the star in 2 epochs at 8 GHz and in 5 epochs at 5 GHz from August 2012 to October 2015 relative to J1627-2426². In the following, the reprocessed data from Ortiz-León et al. (2017), where all positions of Haro 1-6 were referenced to the position of a unique calibrator, J1627-2426, is used. In January 2020, the star was observed at 4.6 GHz relative to J1633-2557, since J1627-2426 is not in ICRF3 and probably too faint for the observation mode, and J1625-2527 has developed structure in the meantime.

To reference the star positions to the ICRF3 position of J1633-2557 at all epochs, the shifts of -0.345 mas in $\alpha$ and -1.2 mas in $\delta$ between J1627-2426 and J1625-2527, as determined by Ortiz-León et al. (2017) from observations between 2005 and 2007, were first subtracted from the star positions in Ortiz-León et al. (2017). Another shift of -0.03 mas in $\alpha$ and 0.06 mas in $\delta$ was further subtracted, corresponding to the offset between the positions of J1625-2527 from Ortiz-León et al. (2017) and from the rfc_2018b catalog used in the 2020 observation. Finally, a third shift of 0.7550 mas in $\alpha$ and -0.0207 mas in $\delta$ was subtracted, which represents the difference between the observed position of J1625-2527 in the 2020 experiment (when phase referenced to J1633-2557) and the rfc_2018b catalog position. This workaround is not ideal; however, it reduces the RMS$_{v,\delta}$ by 10%. Calibrator J1627-2426 would need to be observed with global astrometry and added to the next ICRF realization to resolve this issue.

To show the robustness of the final solution, several solutions based on three different subsets of the positions were obtained. $M_{all}$ includes all the above-mentioned positions as input data to the fitting process, while $M_{w07}$ excludes the position at the seventh epoch since that position showed the largest residuals (about 6 mas for $\alpha*$) and was considered as an outlier. $M_{>2012}$ includes only the latest more precise observations, similar to what was used for the final selection in Ortiz-León et al. (2017). For each of the three different selections, a solution with (new) or without (old) the January 2020 position was obtained. Down-weighting the 2020 position by 0.5 mas, because it is relative to a different calibrator (similar to what was tested for the observations of BH CVn), does not change the results significantly, which is to be expected as the positions for the different calibrators were previously homogenized.

Comparing the model estimates $M_{all}$ with $M_{w07}$, the additional noise required to achieve a $\chi^2$ of 1 and the scatter of residuals decrease significantly in $\alpha*$. The proper motion estimates do not change, but the formal errors of $\mu_{\alpha*}$ decrease by 60%. Furthermore, $\tau$ decreases by 8% and its formal error by 60%. While the parallax for the $M_{all}$ solution is consistent with that from Gaia EDR3, the estimate from the $M_{w07}$ solution instead deviates significantly due to the decreased formal error. Employing only the more precise positions in $M_{>2012}$ reduces the residual scatter by more than half. The estimate of $\tau$ in this solution matches the Gaia EDR3 result. On the other hand, $\mu_{\delta}$ is found to deviate significantly compared to the Gaia EDR3 value, even though the observation selection is closest in time to the Gaia EDR3 observations compared to the other solutions. Since the formal error in proper motion parameters is smallest for the $M_{w07,new}$ solution due to the long time interval between the first and last observation employed, this solution was adopted for the subsequent analysis in Sect. 4.3.

### 3.3.8. $\sigma^2$ CrB

Star $\sigma^2$ CrB was observed at 14 epochs from May 1987 to November 1994 relative to J1613+3412 in Lestrade et al. (1999). The observing frequency was 5.0 GHz for the first three epochs and 8.4 GHz thereafter. Of these data, only 12 epochs were found in the archive and thus used in our analysis. In January 2020 the star was observed at 8.1 GHz relative to the same calibrator with position in ICRF3. For the calculation of the new estimates for the model of stellar motion, the position of the phase calibrator uncorrected for source structure was employed, since the correction was also not made in Lestrade et al. (1999). The subtraction

<footnote n="2">² Positions $\alpha = 246:750025782$, $\delta = -24:444573598$ are given in the VLBA archive for radio source J1627-2427. It is assumed to be the radio source J1627-2426 ($\alpha = 246:750025121$, $\delta = -24:4445743194$) in catalog rfc_2018b, whose position only differs by 2.38 mas and 2.60 mas.</footnote>
of 0.1017 mas in $\alpha$ and 0.1965 mas in $\delta$ is required to reference the star positions from the historical observations to ICRF3.

Additional weights are only needed in case of $M_{new}$. Due to the large time difference between the historical observations and the 2020 position, small but significant linear acceleration terms (for parameterization see Loinard et al. 2007) are needed to properly model the star motion. This is evidenced by the large reduction of residual scatter between solutions $M_{new}$ and $M_{new,a}$ (see Table A.1). The derived acceleration terms are $(-0.046 \pm 0.004)$ mas yr$^{-2}$ in $\alpha*$ and $(-0.018 \pm 0.004)$ mas yr$^{-2}$ in $\delta$. The estimates for $M_{old}$ and $M_{new,a}$ differ by $(-1.137 \pm 0.060)$ mas yr$^{-1}$ in $\mu_{\alpha*}$ and $(-0.451 \pm 0.061)$ mas yr$^{-1}$ in $\mu_{\delta}$, while the formal errors for these parameters decrease by more than 20%. On the other hand, the formal error of $\bar{\omega}$ increases by about 15% and those of the coordinates double. This could be due to the reference epoch (2016.0) being away from the bulk of the observations. In all, the estimated proper motion and parallax parameters for $M_{new,a}$ are found to be closer to those from Gaia EDR3 compared to values in Lestrade et al. (1999) or $M_{new}$.

### 3.3.9. HD 199178

Star HD 199178 was observed at 8 epochs from September 1992 to September 1994 relative to J2102+4702 by Lestrade et al. (1999) at 8.4 GHz. However, the source was too weak for detection at two epochs, thus leaving 6 epochs for analysis. In January 2020, the star was observed at 8.1 GHz relative to the same calibrator with position in ICRF3. For the calculation of the new estimates for the model of stellar motion, the position of the phase calibrator uncorrected for source structure was employed, since the correction was also not made in Lestrade et al. (1999). The subtraction of 0.3760 mas in $\alpha$ and 1.3861 mas in $\delta$ is required to reference the star positions from the historical observations to ICRF3.

The scatter of the residuals does not change comparing $M_{old}$ and $M_{new}$. This shows that the observations fit each other well. Moreover, the formal errors decrease by about 97% for proper motion and by 30% for $\bar{\omega}$. The estimates from $M_{new}$ are closer to the Gaia EDR3 parameters than the ones from Lestrade et al. (1999); however, because of the smaller formal errors, the differences are significant for $M_{new}$, whereas they are not for the Lestrade et al. (1999) solution. No additional weights were added in either solution, and the $\chi^2$ was 0.5. This suggests that uncertainties for the individual positions were perhaps too large.

### 3.3.10. AR Lac

Star AR Lac was observed at 7 epochs from April 1989 to May 1994 relative to J2202+4216 in Lestrade et al. (1999), at 5.0 GHz in the first three epochs and at 8.4 GHz thereafter. However, the data from the first epoch was not available in the archive, thus leaving 6 epochs for our analysis. In January 2020, the star was observed at 8.1 GHz relative to the same calibrator with position in ICRF3. Additionally, the star was observed relative to J2153+4322 because J2202+4216 shows a lot of structure and future observations could benefit from a more compact calibrator, which can nowadays be used due to the increased sensitivity of the antenna network compared to historical observations. The position of the phase calibrator uncorrected for source structure was employed because the correction was also not made in Lestrade et al. (1999). The subtraction of 0.1730 mas in $\alpha$ and 0.0837 mas in $\delta$ is required to reference the star positions from the historical observations to ICRF3.

Adding the new position in 2020 to all available observations relative to J2202+4216, the formal errors of the proper motions decrease by 90% compared to solution $M_{old}$. Furthermore, leaving out the epoch in July 1990, which was affected by poor uv-coverage, was found to reduce the residual scatter and formal errors of all parameters by about 50%, compared to a solution including this epoch (see solution $M_{wo90,J2202+4216, new}$ in Table A.1). Adding the star position relative to J2153+4322 in 2020 in solution $M_{wo90,w0.5,new}$ with an additional uncertainty of 0.5 mas in quadrature, has almost no influence on the estimates, which is verified by the position with respect to J2202+4216 having residuals of zero, while the position with respect to J2153+4322 has residuals of about $\pm 0.28$ mas. If no additional uncertainty was added, the uncertainties of the model parameters would double. Solution $M_{wo90,w0.5,new}$ was selected for the analysis of the rotation parameters in Sect. 4.3. Compared to the solution from Lestrade et al. (1999), $\mu_{\alpha*}$ from $M_{wo90,w0.5,new}$ is closer to the value from Gaia EDR3. However, its formal error is considerably smaller, which is why the selected solution shows an offset relative to Gaia EDR3 that is larger than $3\sigma$, while it is less than $2\sigma$ in the solution of Lestrade et al. (1999). The same applies to $\mu_{\delta}$. In contrast, the parallax from solution $M_{wo90,w0.5,new}$ is significantly off that from Gaia EDR3, while the value from Lestrade et al. (1999) agrees with that from Gaia EDR3 within about $1\sigma$.

### 3.3.11. IM Peg

Star IM Peg was observed at 4 epochs from December 1991 to July 1994 relative to J2253+1608 in Lestrade et al. (1999) at 5.0 GHz in the first epoch and at 8.4 GHz thereafter. The subtraction of $-0.0318$ mas in $\alpha$ and $-0.0005$ mas in $\delta$ is required to reference the star positions to ICRF3. In addition, 35 positions between January 1997 and July 2005 from the Gravity Probe B experiment at 8.3 GHz could be used (Ratner et al. 2012). They are referenced to point C1 in J2253+1608 with coordinates $\alpha = 22h53m57.7479573$, $\delta = 16^\circ8'53.5'561281$ (Bartel et al. 2012), and the calibrator structure was corrected during the fringe-fit (Lebach et al. 2012). The subtraction of 0.2577 mas in $\alpha*$ and 0.4005 mas in $\delta$ is needed for these star positions to be referenced to ICRF3. To weigh this set of positions relative to the other data sets used in our study, the 0.06 mas WRMS scatter of residuals of the astrometric check source determined in Ratner et al. (2012) was used as position uncertainty in $\alpha*$ and $\delta$ for the star at each individual epoch. In January 2020, IM Peg was observed at 8.1 GHz relative to the same calibrator.

Adding the 2020 position, the formal errors of the proper motions in $\alpha*$ and $\delta$ decreased by 30%, while the formal error in $\bar{\omega}$ remained stable. The radio emission from star IM Peg has an orbital motion with a semi-major axis of 0.89 mas and an axis ratio of 0.30 (Ratner et al. 2012). The orbital motion has a period of 24.64877 days (Marsden et al. 2005). The reduction of the orbital pattern in the positions by subtracting the model values derived from the functional model for the linear orbital parameters in Table 3 of Ratner et al. (2012) resulted in a decrease in RMS$_{v,\alpha*}$ of about 0.15 mas and in RMS$_{v,\delta}$ of about 0.05 mas in the corresponding solution (labeled $M_{orb,new}$ in Table A.1). Smaller additive errors were needed for this solution and thus the formal errors of all parameters decreased. However, the astrometric parameters did not change significantly (they remain within $1\sigma$). Since the binary cannot be resolved by Gaia, the determination of the rotation parameters in Sect. 4.3 will use the estimates for the model without correction for the binary
orbit trajectory, $M_{new}$. The new model estimates are not significantly different from those presented in previous studies, but the formal error of the proper motion values is significantly smaller.

### 3.3.12. HD 22468

Star HD 22468 (HR 1099) was observed at 8 epochs from March 1991 to August 1994 relative to J0339-0146 in Lestrade et al. (1999), at 5.0 GHz in the first 3 sessions and at 8.4 GHz thereafter. The subtraction of -0.0200 mas in $\alpha$ and 0.2792 mas in $\delta$ is required for these star positions to be referenced to ICRF3. In addition, one position was obtained during experiment V515C (July 2018) relative to the same calibrator by the Long Baseline Array (LBA) supported by the ATNF (CSIRO), by Russian antennas operated by the Institute of Applied Astronomy of Russian Academy of Sciences (Shuygina et al. 2019), and by the HartRAO antenna in South Africa as described in Titov et al. (2020). The calibrator position was in ICRF3. Three additional positions in ICRF3 between May 2015 and July 2016 from absolute astrometry in S/X mode are published in Titov et al. (2020). Further details are given in Table A.2. In January 2020, the star was observed at 8.1 GHz relative to the same calibrator, J0339-0146, with the VLBA. For the calculation of the new estimates for the model of stellar motion, the position determined without correcting the phases of the calibrator for structure was employed, since the correction was not applied in the observations from the archive.

Several estimates of the astrometric model for HD 22468 were produced to test the impact of the new data compared to using only the Lestrade et al. (1999) data through solution $M_{old}$. The addition of the 2018 and 2020 positions (solution $M_{new}$) leads to a reduction in the formal errors of the proper motions of more than 90%, but also to a decrease of 40% of the formal error of $\omega$. The latter can be explained by the two positions being located at the local minimum and local maximum of the parallax pattern of the star.

As a final test, the three positions from absolute astrometry were included in the analysis, thus all 13 epochs were used. However, the first two such experiments (AOV003, AUA011) had to be discarded ($M_{abs, w09 & 10, new}$) because they show residuals several milliarcseconds large. Only the third experiment, AOV010, provides a good fit to the phase referenced data. The number of observations in the latter is 56 compared to 6 and 10 for the two rejected experiments (Titov et al. 2020). Adding this position to the phase referencing dataset does not change the results within the error bounds, but improves the formal errors by another 5% to 10% for all parameters. This solution was selected for comparison with *Gaia* and the determination of the rotation parameters in Sect. 4.3 below.

### 3.3.13. CoKu HP Tau G2

Star CoKu HP Tau G2 (HP Tau/G2) was observed at 8 epochs between September 2005 and December 2007 with the VLBA at 8.42 GHz by Torres et al. (2009) relative to the calibrator J0426+2327. In addition, Galli et al. (2018) detected the star in 5 epochs at 8.4 GHz and in 4 epochs at 5.0 GHz between December 2012 and October 2017 relative to J0438+2153. They corrected the 8 epochs from Torres et al. (2009) for the shift of the calibrator position and combined the time series. The homogenized data are used for astrometric fitting and the star position obtained relative to J0438+2153 at 4.6 GHz in January 2020, including correction for source structure, is added. The star positions from Galli et al. (2018) were corrected by subtracting

-0.4287 mas in $\alpha$ and -0.0492 mas in $\delta$, to refer them to the calibrator position in ICRF3, as also the case for the position from UL005.

The young star is the primary star in a gravitationally bound system with HP Tau G3 AB, which has a separation of about 10" (Harris et al. 2012; Nguyen et al. 2012). HP Tau G3 AB itself is a close binary (Richichi et al. 1994). Linear terms appear to be not sufficient to model the trajectory of CoKu HP Tau G2 given the above observations, as can be seen from the large residual scatter and the deviation of the estimated proper motion from the *Gaia* EDR3 values for the solution including all observations ($M_{new}$).

A nonlinear model is supported by the RUWE parameter (renormalized unit weight error, calculated from *Gaia* data; Lindegren et al. 2018) in *Gaia* EDR3, which is 3.85 for this star. A value larger than about 1.4 indicates that the star is not a single star as seen by *Gaia*. Galli et al. (2018) estimated orbital parameters in addition to position, proper motion, and parallax. In this work, orbital parameters are not considered. To compare the proper motion with *Gaia* EDR3, the VLBI time series were therefore trimmed to match the *Gaia* EDR3 observations between 25 July 2014 and 28 May 2017 (solution labeled $M_{Gaia}$ in Table A.1). In this solution, the residuals for X band and C band do not differ; therefore frequency dependent position offsets do not need to be considered. In order to best match the *Gaia* EDR3 model parameterization, $M_{Gaia}$ will be further used for the subsequent rotation parameter analysis in Sect. 4.4.

### 3.3.14. del Lib

Star del Lib (HD 132742) was observed at 3 epochs from July 2016 to July 2020. The positions in experiment UL005 (January 2020) were obtained relative to two different calibrators, J1456-0617 and J1510-0843, at 8.1 GHz. During experiment AOV010 (July 2016) it was observed relative to J1456-0617 with the Asia-Oceania VLBI network (AOV) and during experiment V583B (July 2020) it was observed relative to J1512-0905 by the LBA, both at about 8.4 GHz. The positions are listed in Table A.2. The star position derived without correcting the calibrator phases for structure was used for this fit because the other experiments did not use it either. To reference all observations to ICRF3, shifts of -0.4572 mas in $\alpha$ and -0.2092 mas in $\delta$ were subtracted from the position of J1512-0905, which was originally $\alpha = 15h12m50.5329, \delta = -9^\circ5'59.830$. All observations relative to J1456-0617 were already in ICRF3 and no correction was thus necessary.

As mentioned above, an additional uncertainty of 0.5 mas was applied to the star positions that are not relative to J1456-0617 (considered as the primary calibrator) to account for systematic errors. The resulting solution, labeled $M_{w0.5, new}$, is not significantly different from the solution without the additional uncertainties $M_{new}$. However, the VLBI model estimates do differ largely from those of *Gaia* EDR3, but for comparison, those of *Gaia* DR2 (also listed in Table A.1) are again different. Solution $M_{w0.5, new}$ was adopted for the subsequent analysis in Sect. 4.4.

### 3.3.15. HD 142184

Star HD 142184 (HR 5907, 1550-238) was not in the sample of sources observed by the VLBA in January 2020 because it is a Be-type star and therefore may exhibit radio emission from stellar winds leading to radio-optical offsets. It is one of the fastest rotating stars. Data from absolute astrometry performed with the LBA, VLBA, AOV, and other networks were collected.
During experiment V583B (July 2020) HD 142184 was further observed relative to J1553−2422 with the LBA at about 8.4 GHz. To reference the observation to ICRF3, shifts of −0.2098 mas in α and 0.0491 mas in δ were subtracted from the position of J1553−2422, which was originally α = 15h53m31.6278, δ = −24°22′6″036. All positions collected or derived from this work are listed in Table A.2.

The proper motion and parallax (solution M<sub>new</sub> in Table A.1) are in agreement with the *Gaia* EDR3 values within the error bounds. Because most of the data (i.e., the absolute positions) is not referenced to ICRF3 but to the aus2020b<fnref n="3" /> reference frame, the position information was not included in the analysis of the rotation parameters in Sect. 4.4 below, and only the proper motion and parallax were used<fnref n="4" />.

### 3.3.16. Additional single-epoch positions

For four additional stars (HD 167971, V 479 Sct, EI Eri, YY Men), observations at one or two epochs are available from experiments conducted at 8.4 GHz by the LBA, as listed in Table A.2. In experiments V583A (March 2020) and V583B (July 2020), star HD 167971 was observed, but relative to two different calibrators. Shifts of 42.8062 mas in α and −10.2367 mas in δ were subtracted for the star position relative to J1818−1108 and shifts of 1.0813 mas in α and 0.5267 mas in δ were subtracted for the star position relative to J1832−1035 to transfer them to ICRF3. Because the time interval between the first and last scans on HD 167971 in experiment V583A was only about 20–30 min, the uv-coverage was sparse. Thus, the side lobes were large and the uncertainty of the position was high. The same is true for stars V479 Sct and EI Eri, which were also observed in experiment V583A. For these stars, the shifts subtracted to reference positions to ICRF3 were −0.5414 mas in α and −0.4250 mas in δ (for V479 Sct) and −0.6007 mas in α and 0.4208 mas in δ (for EI Eri). Star YY Men was observed in absolute geodetic mode in two single-baseline experiments in 1990 and 1991, from which the position in 1990 could be determined. In addition, it was detected in experiment V583B relative to J0529−7245. In that case, shifts of 0.3641 mas in α and 0.5004 mas in δ were subtracted to reference the position to ICRF3.

### 3.3.17. Recapitulation

The RMS of the residual scatter shows that modern phase referencing observations can reach levels of 0.1 mas to 0.2 mas and below, whereas the inclusion of historical observations from around the 1990s (case of BH CVn, HD 199178, AR Lac, IM Peg, HD 22468) yields RMS levels of 0.2 mas to about 1.0 mas. At the present stage, accelerations are not considered in the rotation parameter analysis. However, six stars were found to have signs of nonlinear proper motion from their VLBI data verified by Student's t-test with a significance level of 5% (V410 Tau, Brun 334, σ² CrB, HD 22468, CoKu HP Tau G2, and HD 142184). For some stars, the model estimates better fit the values of the astrometric *Gaia* model when the same parameterization is applied (i.e., with no acceleration terms). This is usually the case when the VLBI mean epoch is close in time to the *Gaia* observations, such as for V410 Tau, Brun 334, and HD 142184. In the

case of σ² CrB, the proper motion estimates are closer to the *Gaia* estimates when linear accelerations are also estimated. This star is part of those objects for which the mean epoch of the VLBI observations is more distant from the *Gaia* observation time interval. Star CoKu HP Tau G2 has a more complicated trajectory. Restricting the VLBI time interval to the time interval of the *Gaia* EDR3 observations provides a better match to the *Gaia* data; however, there are no significant proper motion differences between models with and without acceleration terms. Star HD 22468 is a close binary with a period of about 2.8 days (Fekel 1983; García-Alvarez et al. 2003). It is located in another binary system ADS 2644A with an orbital period of 2101 yr and 6.2 orbital diameter (Lestrade et al. 1999). For this star, the acceleration terms were not found significant according to Student's t-test. The greater agreement of the parameter estimates with the *Gaia* EDR3 parameters and the reduction in E<sub>α*</sub> as well as RMS<sub>v,α*</sub> however suggests that a long-term acceleration in α should be tested again when more observations have been acquired in the future. The detailed investigation of the various star time series showed that for some stars more observations would help to prove the existence of long-term accelerations (V410 Tau, Brun 334, HD 22468), decorrelating the model parameters (SS Cyg, Haro 1-6), or investigating the reason for significant model differences between VLBI and *Gaia* (AR Lac).

Other reasons, such as orbital jitter (see Lestrade et al. 1999), may account for the differences between the estimated VLBI and *Gaia* models as well. This effect is specific to each source. For the present study, it is not considered, since long-term linear proper motions are of interest.

### 3.4. Absolute positions from VLBI and correction of parallax effect

In the calculations above, the absolute star position was lost by the application of source structure corrections during the fringe fit of the phase referencing calibrator data, as explained in Sect. 3.2. Thus, the mean positions α₀ and δ₀ from the five parameter fit at the mean epoch of the respective star time series have also no absolute reference. For this reason, the absolute single-epoch positions from Table 4 in Lunz et al. (2023) have been used instead for the frame tie. They are based on calibrator data that is not corrected for source structure, in agreement with the scheme used to derive the ICRF3 source positions. These single-epoch positions must be corrected for the parallax effect at the epoch of the observation, and the epoch *t* of the observation itself must be corrected for the Römer delay to match the barycentric time (Lindegren 2020a). For the stars considered in this study, the corrections employed the newly determined parallaxes listed in Table A.1. For the other stars, such as the additional stars in Sect. 3.3.16, the corrected *Gaia* parallax (Lindegren et al. 2021a) was used. All corrections are listed in Table 2. They are based on the barycentric coordinates of the Earth's center at the time of observation from the DE 421 ephemeris (Folkner et al. 2009) using the VLBI software VieVS@GFZ developed by GFZ in Potsdam, Germany.

The random position errors that come out from the phase referencing analysis are based on the beam shape and dynamic range of the image. However, these are only applicable to the relative calibrator-star positions. They do not represent a full error budget for an absolute star position from phase referencing at a single epoch. To this end, we consider a more realistic error budget by combining the random errors with additional uncertainties, using the same formulas as in Lunz et al. (2023). This

<footnote n="3">³ ftp://ivs.bkg.bund.de/pub/vlbi/ivsproducts/crf/aus2020b.crf.txt</footnote>

<footnote n="4">⁴ In future work, the position from V583B can be added to the rotation parameter analysis as an absolute position. However, for this work, it was not considered.</footnote>
<table><caption>**Table 2.** Römer delay, position corrections due to parallax effect, and radial velocity for the single-epoch positions used for the rotation parameter analysis.</caption><thead><tr><th>Star</th><th>Epoch t<br>(Julian year)</th><th>Roe (t)<br>(s)</th><th>Δα (t)<br>(mas)</th><th>Δδ (t)<br>(mas)</th><th>ν<sub>r</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>LS I +61 303</td><td>J2020.01520</td><td>251</td><td>+0.672</td><td>-0.043</td><td>-41.4</td></tr><tr><td>UX Ari</td><td>J2020.01542</td><td>316</td><td>+17.109</td><td>+1.681</td><td>50.7</td></tr><tr><td>HD 22468</td><td>J2020.01545</td><td>275</td><td>+25.132</td><td>+12.769</td><td>-15.3</td></tr><tr><td>HD 283447</td><td>J2020.01544</td><td>380</td><td>+6.020</td><td>+0.182</td><td>16.0</td></tr><tr><td>V410 Tau</td><td>J2020.01546</td><td>385</td><td>+5.450</td><td>+0.099</td><td>19.9</td></tr><tr><td>HD 283572</td><td>J2020.01547</td><td>389</td><td>+5.380</td><td>+0.081</td><td>14.2</td></tr><tr><td>T Tau</td><td>J2020.01548</td><td>384</td><td>+4.494</td><td>+0.890</td><td>19.2</td></tr><tr><td>Brun 334</td><td>J2020.01560</td><td>398</td><td>+0.930</td><td>+1.125</td><td>20.3</td></tr><tr><td>TYC 5346-538-1</td><td>J2020.01561</td><td>393</td><td>+0.803</td><td>+1.135</td><td>0.0</td></tr><tr><td>BH CVn</td><td>J2020.01370</td><td>71</td><td>-24.387</td><td>+7.892</td><td>6.4</td></tr><tr><td>σ<sup>2</sup> CrB</td><td>J2020.01416</td><td>-171</td><td>-33.110</td><td>+30.065</td><td>-14.7</td></tr><tr><td>Haro 1-6</td><td>J2020.01429</td><td>-395</td><td>-4.526</td><td>+0.401</td><td>-3.7</td></tr><tr><td>DoAr 51</td><td>J2020.01429</td><td>-402</td><td>-4.306</td><td>+0.322</td><td>0.0</td></tr><tr><td>HD 199178</td><td>J2020.01452</td><td>-160</td><td>+4.968</td><td>+7.388</td><td>-30.8</td></tr><tr><td>SS Cyg</td><td>J2020.01451</td><td>-128</td><td>+0.059</td><td>+0.068</td><td>-62.0</td></tr><tr><td>AR Lac</td><td>J2020.01460</td><td>-90</td><td>+22.514</td><td>+17.663</td><td>-33.8</td></tr><tr><td>IM Peg</td><td>J2020.01469</td><td>-189</td><td>+8.242</td><td>+5.305</td><td>-14.4</td></tr><tr><td>UV Psc</td><td>J2020.01506</td><td>38</td><td>+12.972</td><td>+5.242</td><td>6.5</td></tr><tr><td>HD 8357</td><td>J2020.01507</td><td>52</td><td>+20.280</td><td>+8.110</td><td>12.7</td></tr><tr><td>RZ Cas</td><td>J2020.01520</td><td>244</td><td>+36.601</td><td>-3.507</td><td>-39.4</td></tr><tr><td>B Per</td><td>J2020.01544</td><td>359</td><td>+8.367</td><td>-2.347</td><td>19.8</td></tr><tr><td>CoKu HP Tau G2</td><td>J2020.01549</td><td>403</td><td>+3.925</td><td>+0.441</td><td>16.6</td></tr><tr><td>SV Cam</td><td>J2020.01461</td><td>249</td><td>+9.617</td><td>-10.082</td><td>-13.8</td></tr><tr><td>54 Cam</td><td>J2020.01358</td><td>399</td><td>-4.082</td><td>-5.451</td><td>27.5</td></tr><tr><td>IL Hya</td><td>J2020.01348</td><td>268</td><td>-5.318</td><td>+6.125</td><td>-7.3</td></tr><tr><td>HU Vir</td><td>J2020.01371</td><td>75</td><td>-6.081</td><td>+2.755</td><td>-0.7</td></tr><tr><td>DK Dra</td><td>J2020.01362</td><td>211</td><td>-19.951</td><td>-0.544</td><td>-45.3</td></tr><tr><td>RS CVn</td><td>J2020.01370</td><td>106</td><td>-8.314</td><td>+2.315</td><td>-13.6</td></tr><tr><td>del Lib</td><td>J2020.01405</td><td>-243</td><td>-5.464</td><td>+2.082</td><td>-38.7</td></tr><tr><td>SZ Psc</td><td>J2020.01475</td><td>-211</td><td>+8.576</td><td>+4.264</td><td>12.0</td></tr><tr><td>HD 224085</td><td>J2020.01471</td><td>-36</td><td>+25.316</td><td>+12.056</td><td>-20.5</td></tr><tr><td>HD 142184</td><td>...</td><td>...</td><td>...</td><td>...</td><td>-9.20</td></tr><tr><td>HD 167971</td><td>J2020.21812</td><td>-36</td><td>-0.784</td><td>-0.022</td><td>14.2</td></tr><tr><td>HD 167971</td><td>J2020.53120</td><td>473</td><td>+0.226</td><td>-0.144</td><td>14.2</td></tr><tr><td>V479 Sct</td><td>J2020.21812</td><td>-53</td><td>-0.547</td><td>-0.023</td><td>0.0</td></tr><tr><td>EI Eri</td><td>J2020.21812</td><td>-231</td><td>+16.427</td><td>-1.238</td><td>17.6</td></tr><tr><td>YY Men</td><td>J2020.53120</td><td>87</td><td>-10.436</td><td>-3.707</td><td>-8.5</td></tr><tr><td>YY Men</td><td>J1990.94045</td><td>-70</td><td>+0.850</td><td>+4.580</td><td>-8.5</td></tr></tbody></table>

**Notes.** The Römer delay Roe needs to be subtracted from epoch *t*, which is the mean epoch between the first and last scans. In addition, the parallax effects, $\Delta\alpha$ and $\Delta\delta$, need to be added to the positions $\alpha(t)$ and $\delta(t)$ in Table A.2 and Table 4 in Lunz et al. (2023). The position of star HD 142184 was not used because it could not be transferred to ICRF3. The last column shows the radial velocity ($v_r$) taken from SIMBAD and used in the rotation parameter calculations.

comprises the celestial reference frame position uncertainty from the phase calibrator and an average for a combination of several types of errors (delay model errors, source structure errors, and the difference between positions obtained from phase and group-delays). The latter is about 0.28 mas in $\alpha*$ or 0.69 mas in $\delta$ from measurements to two different calibrators for a subset of the January 2020 VLBA positions in Lunz et al. (2023). Such more realistic uncertainties are given in Table 3 for the additional stars considered in this study (see Sect. 3.3.16). The positional uncertainty from absolute geodetic observations, such as that for YY Men in the last line of Table 3, remains unchanged.

## 4. Results

This section presents the results of the rotation parameter analysis for different scenarios in which the changes to the VLBI dataset described in the previous sections were tested. The analysis was done in the same way as in Lunz et al. (2023). In particular, correction of the Gaia EDR3 parallaxes was performed using the functions provided in Lindegren et al. (2021a). For details see Lunz et al. (2023) and Lindegren (2020a,b). The weighted mean WM, the weighted root mean square WRMS,

<table><caption>**Table 3.** Realistic error budget for the single-epoch positions of the additional stars.</caption><thead><tr><th>Star</th><th>Epoch t<br>(Julian year)</th><th>σ<sub>α*,absolute</sub><br>(mas)</th><th>σ<sub>δ,absolute</sub><br>(mas)</th></tr></thead><tbody><tr><td>HD 167971</td><td>J2020.21812</td><td>2.403</td><td>2.011</td></tr><tr><td>HD 167971</td><td>J2020.53120</td><td>4.896</td><td>5.138</td></tr><tr><td>V479 Sct</td><td>J2020.21812</td><td>0.320</td><td>0.763</td></tr><tr><td>EI Eri</td><td>J2020.21812</td><td>0.297</td><td>0.726</td></tr><tr><td>YY Men</td><td>J2020.53120</td><td>0.301</td><td>0.711</td></tr><tr><td>YY Men</td><td>J1990.94045</td><td>4.350</td><td>13.600</td></tr></tbody></table>

**Notes.** Realistic error budget $\sigma_{\alpha*,absolute}$ and $\sigma_{\delta,absolute}$ for the star positions $\alpha(t)$ and $\delta(t)$ at epoch *t* in Table A.2 when used as absolute positions in ICRF3 for comparison to the Gaia positions.

and the mean standard deviation ME of the rotation parameters using a representative range of iterations were used to assess the solutions and are listed in Table 4 for each scenario. These quantities were derived as described in Lunz et al. (2023). Similar to Lunz et al. (2023), the significance of the change in WM values between two scenarios is determined by a two-sample *t*-test with the null hypothesis that the mean values are equal (5% significance level). It is assumed that the rotation parameters are normally distributed with unknown and unequal variances (Behrens-Fisher problem). If the probability value is smaller than the significance level, the difference in the WM between the two scenarios is significant as the test cannot be accepted. Furthermore, as in Lindegren (2020a,b) and Lunz et al. (2023), an iteration solution was selected that best represents the rotation parameters of the respective scenario. It is referred to as the “baseline solution.” Table 5 summarizes the baseline solutions for the various scenarios

### 4.1. Rotation parameters including Galactocentric acceleration effect on positions and proper motions

We first reprocessed the results “55,EDR3” (where “55” represents the number of stars in the sample and “EDR3” represents the Gaia data release employed) from Lunz et al. (2023). The new solution includes the effect of Galactocentric acceleration as introduced in Sect. 2. The new rotation parameter results “55,EDR3,GA” (where “GA” indicates that the effect of Galactocentric acceleration was corrected) are presented in Fig. 1. Taking into account the effect of Galactocentric acceleration did not significantly change the WM, WRMS, and ME statistics of the rotation parameters (see Table 4). The baseline solution for scenario “55,EDR3,GA” was chosen at *k* = 13 number of rejected stars as for scenario “55,EDR3” because both orientation and spin parameters are stable for some iterations thereafter. The rotation parameters from this solution are very close to those of “55,EDR3”, both of which are reported in Table 5. Correlation coefficients for *k* = 13 and *T* = 2016.0 are

$$ \text{corr}[\epsilon_X(T), \epsilon_Y(T), \epsilon_Z(T), \omega_X, \omega_Y, \omega_Z] = \begin{bmatrix} +1.000 & +0.225 & +0.212 & +0.200 & +0.024 & -0.003 \\ \dots & +1.000 & +0.177 & +0.021 & +0.177 & +0.010 \\ \dots & \dots & +1.000 & -0.007 & +0.039 & +0.023 \\ \dots & \dots & \dots & +1.000 & +0.044 & +0.326 \\ \dots & \dots & \dots & \dots & +1.000 & +0.018 \\ \dots & \dots & \dots & \dots & \dots & +1.000 \end{bmatrix} \quad (5) $$
<table><caption>**Table 4.** WM, WRMS, and ME of the rotation parameters for various scenarios.</caption><thead><tr><th>Scenario</th><th></th><th>ε<sub>X</sub></th><th>ε<sub>Y</sub></th><th>ε<sub>Z</sub></th><th>ω<sub>X</sub></th><th>ω<sub>Y</sub></th><th>ω<sub>Z</sub></th></tr></thead><tbody><tr><td rowspan="3">55,EDR3</td><td>WM<sub>55,EDR3</sub></td><td>0.294</td><td>0.340</td><td>0.077</td><td>0.034</td><td>0.055</td><td>-0.020</td></tr><tr><td>WRMS<sub>55,EDR3</sub></td><td>0.049</td><td>0.074</td><td>0.084</td><td>0.016</td><td>0.014</td><td>0.012</td></tr><tr><td>ME<sub>55,EDR3</sub></td><td>0.079</td><td>0.099</td><td>0.060</td><td>0.011</td><td>0.013</td><td>0.012</td></tr><tr><td rowspan="3">55,EDR3,GA</td><td>WM<sub>55,EDR3,GA</sub></td><td>0.289</td><td>0.343</td><td>0.078</td><td>0.031</td><td>0.053</td><td>-0.017</td></tr><tr><td>WRMS<sub>55,EDR3,GA</sub></td><td>0.048</td><td>0.074</td><td>0.083</td><td>0.016</td><td>0.014</td><td>0.012</td></tr><tr><td>ME<sub>55,EDR3,GA</sub></td><td>0.078</td><td>0.099</td><td>0.060</td><td>0.011</td><td>0.013</td><td>0.012</td></tr><tr><td rowspan="3">50,EDR3,GA</td><td>WM<sub>50,EDR3,GA</sub></td><td>0.315</td><td>0.327</td><td>0.069</td><td>0.037</td><td>0.037</td><td>-0.004</td></tr><tr><td>WRMS<sub>50,EDR3,GA</sub></td><td>0.028</td><td>0.041</td><td>0.085</td><td>0.014</td><td>0.007</td><td>0.011</td></tr><tr><td>ME<sub>50,EDR3,GA</sub></td><td>0.081</td><td>0.102</td><td>0.061</td><td>0.012</td><td>0.014</td><td>0.012</td></tr><tr><td rowspan="3">55,EDR3,GA,NM</td><td>WM<sub>55,EDR3,GA,NM</sub></td><td>0.305</td><td>0.376</td><td>0.065</td><td>0.029</td><td>0.064</td><td>-0.014</td></tr><tr><td>WRMS<sub>55,EDR3,GA,NM</sub></td><td>0.062</td><td>0.081</td><td>0.082</td><td>0.010</td><td>0.013</td><td>0.014</td></tr><tr><td>ME<sub>55,EDR3,GA,NM</sub></td><td>0.092</td><td>0.115</td><td>0.068</td><td>0.010</td><td>0.012</td><td>0.010</td></tr><tr><td rowspan="3">49,EDR3,GA,NM</td><td>WM<sub>49,EDR3,GA,NM</sub></td><td>0.267</td><td>0.479</td><td>0.080</td><td>0.025</td><td>0.043</td><td>0.005</td></tr><tr><td>WRMS<sub>49,EDR3,GA,NM</sub></td><td>0.038</td><td>0.047</td><td>0.078</td><td>0.008</td><td>0.006</td><td>0.005</td></tr><tr><td>ME<sub>49,EDR3,GA,NM</sub></td><td>0.099</td><td>0.123</td><td>0.072</td><td>0.011</td><td>0.014</td><td>0.011</td></tr><tr><td rowspan="3">60,EDR3,GA,NM</td><td>WM<sub>60,EDR3,GA,NM</sub></td><td>0.347</td><td>0.360</td><td>0.060</td><td>0.031</td><td>0.065</td><td>-0.016</td></tr><tr><td>WRMS<sub>60,EDR3,GA,NM</sub></td><td>0.114</td><td>0.080</td><td>0.078</td><td>0.011</td><td>0.012</td><td>0.014</td></tr><tr><td>ME<sub>60,EDR3,GA,NM</sub></td><td>0.090</td><td>0.112</td><td>0.067</td><td>0.010</td><td>0.011</td><td>0.010</td></tr><tr><td rowspan="3">53,EDR3,GA,NM</td><td>WM<sub>53,EDR3,GA,NM</sub></td><td>0.282</td><td>0.458</td><td>0.100</td><td>0.025</td><td>0.041</td><td>0.006</td></tr><tr><td>WRMS<sub>53,EDR3,GA,NM</sub></td><td>0.046</td><td>0.054</td><td>0.100</td><td>0.009</td><td>0.009</td><td>0.008</td></tr><tr><td>ME<sub>53,EDR3,GA,NM</sub></td><td>0.098</td><td>0.122</td><td>0.071</td><td>0.011</td><td>0.014</td><td>0.011</td></tr></tbody></table>

**Notes.** The order of the scenarios is according to their appearance in the main text. For the calculation of the statistics iterations 12 to 47 were used for all scenarios including 55 stars, while iterations 11 to 42 were used for the scenario “50,EDR3,GA”, 11 to 41 for the scenario “49,EDR3,GA,NM”, 11 to 52 for the scenario “60,EDR3,GA,NM”, and 11 to 45 for the scenario “53,EDR3,GA,NM”. For the derivation of the ME values the last 10 iterations were also rejected, as standard deviations of the rotation parameters increase substantially if only few stars are available for the calculation. Units are milliarcseconds for $ε_X$, $ε_Y$, and $ε_Z$ and milliarcseconds per year for $ω_X$, $ω_Y$, and $ω_Z$.

These coefficients are identical to those in solution “55,EDR3” (see Lunz et al. 2023). In all, the baseline solution is very close to that of “55,EDR3”, hence meaning that adding Galactocentric acceleration has almost no impact on the results.

## 4.2. Impact of individual stars on rotation parameters

From Fig. 1, it can be seen that some stars that are not in the group of the first few stars that get rejected have a large impact on the spin parameters. Offsets of up to 0.03 mas yr⁻¹ are visible when just one star is rejected. This is the case for $σ^2$ CrB ($k = 12$), V410 Tau ($k = 19$), HD 22468 ($k = 23$) and LS I+61 303 ($k = 24$). All of these stars are in multiple star systems. As an example, V410 Tau has one of the lowest uncertainties in the VLBI proper motions in this scenario, with values of 0.017 mas yr⁻¹ in $μ_{α*}$ and 0.020 mas yr⁻¹ in $μ_δ$, and at the same time it is a star in a multiple star system. Excluding that star from the solution shows that the spin parameter in X becomes 0.025 mas yr⁻¹ higher after the rejection at $k = 19$, hence reflecting the large impact of this one star on the entire analysis. Inflating the uncertainty in the proper motion of V410 Tau to about 0.10 mas yr⁻¹ in each direction is required to achieve a similar effect as if the star were completely excluded from the analysis. This suggests that the proper motion uncertainties for this star were too optimistic or that acceleration parameters are significant and must be accounted for in the analysis.

In a similar way, an offset of about 0.15 mas appears in the orientation parameters along the Y and Z directions when

DoAr 51 gets rejected ($k = 13$). On the other hand, no offset in the spin parameters is visible. This star is an unresolved binary as seen by Gaia, while it shows two closeby components in the VLBI images (see Lunz et al. 2023). Since the offset appears only in the orientation parameters, it must be investigated which of the center of luminosity or barycenter works better as a counterpart to the Gaia position, which is equal to the photocenter. A solution where the barycenter was used instead of the center of luminosity changed the offset only by up to 2 µas per rotation axis, which is insignificant. Looking at the residuals of the January 2020 position for such a solution did not shed light on this question either, because the magnitude of the residuals is smaller in $α*$ for the center of luminosity and smaller in $δ$ for the barycenter. The other two close binaries in our sample, HD 283447 and UX Ari (see Lunz et al. 2023), are rejected within the first few iterations anyway – so they have no bearing on this analysis.

The statistics for a scenario excluding the five stars producing offsets in the rotation parameters are reported in Tables 4 and 5 under the label “50,EDR3,GA”. The WRMS of the orientation parameters in X and Y in that scenario is reduced by 40% compared to WRMS<sub>55,EDR3,GA</sub>. The WRMS in Z remained similar. For $ω_Y$, the decrease is by 50% and the difference in the parameter values was also found to be significant. For the other spin parameters, the decrease is by about 10%. The mean standard deviations slightly deteriorate (by less than 10%), which can be explained by the reduced number of objects in the dataset.
<table><caption>**Table 5.** Baseline solutions for various scenarios.</caption><thead><tr><th>Scenario</th><th>k</th><th>l</th><th>ε<sub>X</sub>(T)</th><th>ε<sub>Y</sub>(T)</th><th>ε<sub>Z</sub>(T)</th><th>ω<sub>X</sub></th><th>ω<sub>Y</sub></th><th>ω<sub>Z</sub></th><th>Q/n</th></tr><tr><td></td><td></td><td></td><td>σ<sub>ε<sub>X</sub></sub>(T)</td><td>σ<sub>ε<sub>Y</sub></sub>(T)</td><td>σ<sub>ε<sub>Z</sub></sub>(T)</td><td>σ<sub>ω<sub>X</sub></sub></td><td>σ<sub>ω<sub>Y</sub></sub></td><td>σ<sub>ω<sub>Z</sub></sub></td><td></td></tr></thead><tbody><tr><td>55,EDR3</td><td>13</td><td>42</td><td>+0.226<br>0.070</td><td>+0.327<br>0.091</td><td>+0.168<br>0.054</td><td>+0.022<br>0.010</td><td>+0.065<br>0.011</td><td>-0.016<br>0.010</td><td>5.58</td></tr><tr><td>55,EDR3,GA</td><td>13</td><td>42</td><td>+0.222<br>0.070</td><td>+0.327<br>0.091</td><td>+0.167<br>0.054</td><td>+0.020<br>0.010</td><td>+0.063<br>0.011</td><td>-0.013<br>0.010</td><td>5.58</td></tr><tr><td>50,EDR3,GA</td><td>10</td><td>40</td><td>+0.313<br>0.075</td><td>+0.314<br>0.097</td><td>+0.162<br>0.057</td><td>+0.042<br>0.012</td><td>+0.023<br>0.014</td><td>+0.009<br>0.012</td><td>6.88</td></tr><tr><td>55,EDR3,GA,NM</td><td>12</td><td>43</td><td>+0.285<br>0.081</td><td>+0.284<br>0.104</td><td>+0.169<br>0.061</td><td>+0.033<br>0.009</td><td>+0.073<br>0.010</td><td>-0.026<br>0.009</td><td>5.99</td></tr><tr><td>49,EDR3,GA,NM</td><td>13</td><td>36</td><td>+0.291<br>0.092</td><td>+0.499<br>0.118</td><td>+0.169<br>0.066</td><td>+0.026<br>0.011</td><td>+0.035<br>0.013</td><td>+0.010<br>0.010</td><td>4.30</td></tr><tr><td>60,EDR3,GA,NM</td><td>15</td><td>45</td><td>+0.322<br>0.080</td><td>+0.228<br>0.099</td><td>+0.163<br>0.061</td><td>+0.034<br>0.009</td><td>+0.072<br>0.010</td><td>-0.026<br>0.009</td><td>6.43</td></tr><tr><td>53,EDR3,GA,NM</td><td>12</td><td>41</td><td>+0.319<br>0.091</td><td>+0.423<br>0.111</td><td>+0.158<br>0.066</td><td>+0.030<br>0.011</td><td>+0.029<br>0.013</td><td>+0.010<br>0.010</td><td>8.27</td></tr></tbody></table>

**Notes.** The numbers *k* and *l* provide the number of rejected stars and the number of stars remaining in the sample for each baseline solution. The σ values indicate the formal errors of the parameter estimates from the baseline solution for each scenario. Final uncertainties are obtained by multiplying the σ values by $\sqrt{Q/n}$. Units are in milliarcseconds for ε<sub>X</sub>(T), ε<sub>Y</sub>(T), and ε<sub>Z</sub>(T) and milliarcseconds per year for ω<sub>X</sub>, ω<sub>Y</sub>, and ω<sub>Z</sub>.

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Results for scenario “55,EDR3,GA”, when using the VLBI data and Gaia EDR3 as in Fig. 3 in Lunz et al. (2023) but correcting for the effect of Galactocentric acceleration. The orientation and spin parameters for 51 different adjustment solutions are shown in the upper row. For each adjustment solution the star with the largest (Qᵢ/nᵢ) (which indicates the quality of the fit of the star data to a five parameter astrometric model, see Lindegren 2020a) was rejected in the following iteration. The respective (Qᵢ/nᵢ) along with the star name are shown in the lower left plot, whereas the lower right plot shows the quality of the overall fit Q/n, equivalent to χ² of the adjustment.</figcaption></figure>
### 4.3. Rotation parameters including the new models instead of one-epoch observations

In this section, we replace the VLBI data for 12 stars (HD 283572, V410 Tau, SS Cyg, Brun 334, TYC 5346-538-1, BH CVn, Haro 1-6, $\sigma^2$ CrB, HD 199178, AR Lac, IM Peg, and HD 22468) by the new models of stellar motion determined in Sect. 3.3 for these stars and replicate the rotation parameter analysis “55,EDR3,GA”. This way, the change of the rotation parameters due to the new model estimates can be directly assessed. In the analysis, absolute positions determined from the fringe fit of the calibrator to a point source model and corrected as described in Sect. 3.4, were used in place of the positions derived from the models of stellar motions to best connect to ICRF3. Furthermore, the new data was also corrected for the effect of Galactocentric acceleration as in described Sect. 4.1 and no stars were excluded. The resulting rotation parameter estimates for this scenario, denoted “55,EDR3,GA,NM” (where “NM” indicates that star data was replaced with the new models of stellar motion), are shown in Fig. A.2.

The new scenario still shows some offsets for iterations $10 \le k \le 35$ in both the orientation and spin parameters (as for “55,EDR3,GA”). This indicates that data from individual stars still has an impact on the derived rotation parameters. Small offsets in orientation occur at iterations 12, 22, and 24 when DoAr 51, LSI+61 303, and HD 283641 are rejected. Offsets in spin appear for the iterations when stars HD 22468 ($k=21$), LSI+61 303 ($k=22$), V410 Tau ($k=28$), and AR Lac ($k=32$) are rejected. Among these stars, only LSI+61 303 had no new model considered due to the complexity of its trajectory. It is worth noting that the introduction of the new stellar motion models reduced the offsets caused by $\sigma^2$ CrB. On the other hand, a new offset, not present before, was introduced for AR Lac. This suggests that more emphasis has to be put on the details of the VLBI-Gaia comparison for these objects.

When comparing “55,EDR3,GA,NM” to “55,EDR3,GA”, the largest deviations in WM occur in the Y direction (for both orientation and spin). The new estimates for models of stellar motion reduce the scatter in $\omega_X$ and increase it in $\epsilon_X$ and $\epsilon_Y$ (Figs. 1 and A.2), which is also reflected by the WRMS statistics (Table 4). On average, the ME decreased by about 15% for the spin, while for the orientation parameters it increased by about 15%. The “55,EDR3,GA” scenario has a lower ME for the orientation parameters because more VLBI positions are involved in the adjustment than in the “55,EDR3,GA,NM” scenario. At the same time, the uncertainties of the spin parameters are lower for the “55,EDR3,GA,NM” scenario because improved VLBI proper motion estimates were used.

The new baseline solution is chosen to be at $k=12$ rejected stars because both orientation and spin parameters are stable for some iterations thereafter. The corresponding rotation parameters are given in Table 5. The spin tends to be larger compared to the baseline solution of “55,EDR3,GA”. Correlation coefficients for $k=12$ and $T=2016.0$ are

$$ \text{corr}[\epsilon_X(T), \epsilon_Y(T), \epsilon_Z(T), \omega_X, \omega_Y, \omega_Z] $$

$$ = \begin{bmatrix} +1.000 & +0.199 & +0.193 & +0.001 & +0.004 & -0.050 \\ \dots & +1.000 & +0.161 & -0.002 & -0.030 & -0.018 \\ \dots & \dots & +1.000 & -0.052 & +0.007 & -0.093 \\ \dots & \dots & \dots & +1.000 & -0.034 & +0.334 \\ \dots & \dots & \dots & \dots & +1.000 & -0.098 \\ \dots & \dots & \dots & \dots & \dots & +1.000 \end{bmatrix}, \quad (6) $$

which shows that no strong correlations are present. Compared to the scenario with the old models, “55,EDR3,GA”, the correlation coefficients between the orientation and spin parameters in each axis have now vanished. Only a weak correlation between the X and Z spin components exists.

Another scenario “49,EDR3,GA,NM” was tested where the four stars HD 22468, LSI+61 303, V410 Tau, and AR Lac, which produce jumps in the iterative results of the spin parameters, and the two stars DoAr 51 and HD 283641, which produce jumps in the iterative results of the orientation parameters, were excluded from the beginning. The WRMS statistics for all parameters dropped by 5% to 65% compared to those for the scenario “55,EDR3,GA,NM”. At the same time, the WM values changed significantly for $\omega_Y$ and $\omega_Z$. The six rotation parameters of a new baseline solution at $k=13$ based on this scenario are provided in Table 5. As for the WM, some parameters changed significantly, namely $\epsilon_Y$, $\omega_Y$ and $\omega_Z$. In particular, the spin parameter along the Y direction was cut by half. Correlation coefficients for $k=13$ and $T=2016.0$ are

$$ \text{corr}[\epsilon_X(T), \epsilon_Y(T), \epsilon_Z(T), \omega_X, \omega_Y, \omega_Z] $$

$$ = \begin{bmatrix} +1.000 & +0.180 & +0.104 & -0.021 & -0.015 & -0.041 \\ \dots & +1.000 & +0.093 & -0.009 & -0.117 & +0.031 \\ \dots & \dots & +1.000 & -0.054 & +0.030 & -0.093 \\ \dots & \dots & \dots & +1.000 & -0.044 & +0.258 \\ \dots & \dots & \dots & \dots & +1.000 & -0.257 \\ \dots & \dots & \dots & \dots & \dots & +1.000 \end{bmatrix}, \quad (7) $$

which shows that the correlation coefficients between the orientation parameters decrease compared to scenario “55,EDR3,GA,NM”. Also, the correlation coefficient between $\omega_X$ and $\omega_Z$ decreases, while it increases (in the absolute sense) between $\omega_Y$ and $\omega_Z$.

### 4.4. Rotation parameters including the new models and additional observations

In another test, observations from five stars, HD 142184, EI Eri, HD 167971, V479 Sct, and YY Men, were added to the previous dataset. Furthermore, proper motion and parallax information was added to the two 2020 positions of CoKu HP Tau G2 and del Lib, respectively. The resulting plots for this scenario, labeled “60,EDR3,GA,NM”, are shown in Fig. A.3.

Star EI Eri is rejected at $k=7$, star HD 167971 at $k=8$, star V479 Sct at $k=10$, and star YY Men at $k=17$. For the first three stars and CoKu HP Tau G2 (rejected at $k=14$), small shifts occur in both orientation and spin when they are excluded. Star del Lib is rejected at $k=25$, and no shift can be identified. Star HD 142184 gets excluded as one of the last stars. Furthermore, the offsets in orientation for DoAr 51, LSI+61 303, and HD 283641 and in spin for HD 22468, LSI+61 303, V410 Tau, and AR Lac remain similar to those in scenario “55,EDR3,GA,NM”.

Comparing “60,EDR3,GA,NM” to “55,EDR3,GA,NM”, the MEs are slightly improved. The WMs and WRMSs do not significantly change, except for the WRMS of $\epsilon_X$, which almost doubles (increasing from 0.062 mas to 0.114 mas). This increase can be explained by CoKu HP Tau G2 now being included in the statistics.

The new baseline solution is chosen to be at $k=15$ rejected stars because both orientation and spin parameters remain stable for some iterations thereafter. The results of the rotation parameters for this scenario (see Table 5) are very close to those of the baseline solution for scenario “55,EDR3,GA,NM”.
Correlation coefficients for *k* = 15 and *T* = 2016.0 are not significantly different than for the baseline solution of scenario “55,EDR3,GA,NM”.

Another scenario “53,EDR3,GA,NM” was tested where the four stars HD 22468, LSI+61303, V410 Tau, and AR Lac, which produce jumps in the iterative results of the spin parameters, and the three stars DoAr 51, CoKu HP Tau G2, and HD 283641, which produce jumps in the iterative results of the orientation parameters, were excluded at the beginning. The iterative solutions are shown in Fig. A.4. The WRMS statistics dropped by 20% to 60% for all parameters except for $ε_Z$ compared to “60,EDR3,GA,NM”. Otherwise, the same conclusions can be drawn as when comparing “55,EDR3,GA,NM” and “49,EDR3,GA,NM”, so they are not repeated here. This shows that the new data has little impact on the results. The six rotation parameters of the corresponding baseline solution at *k* = 12 are listed in Table 5. The correlation coefficients are similar to those of “49,EDR3,GA,NM”.

## 5. Discussion

The different scenarios considered in the previous section characterize the changes in the analysis of the rotation parameters between ICRF3 and Gaia EDR3 due to the improved VLBI data. In all of them, the *Y* component of the rotation parameters is still less well determined than the *X* and *Z* components, as indicated by the formal errors. The $\chi^2$ of the iterative solutions in scenarios “50,EDR3,GA” and “55,EDR3,GA” drops to below unity when there are 15 stars left in the sample, and for “49,EDR3,GA,NM” and “55,EDR3,GA,NM” when there are 16 stars left. For “53,EDR3,GA,NM” and “60,EDR3,GA,NM”, the $\chi^2$ drops to below unity when 17 stars are still in the sample. This suggests that systematic errors were not reduced in the new solutions compared to scenario “55,EDR3” in Lunz et al. (2023), where the threshold was 16 stars. The addition of new data in scenario “60,EDR3,GA,NM” has little impact on the rotation parameter analysis as the additional noise that was required to be added to make the error budget realistic is currently too large.

All scenarios in which the stars that lead to offsets in the iterative rotation parameter results were excluded a priori (“50,EDR3,GA”, “49,EDR3,GA,NM”, and “53,EDR3,GA,NM”) show smaller $\omega_Y$ values compared to those from the scenarios in which all stars were included. To test the impact of the exclusion of other stars on the results, additional test scenarios were created in which the 11 most deviating stars of the “60,EDR3,GA,NM” scenario were excluded from the beginning. Then, a random selection of four other stars was excluded, so that in total 15 stars, as in the baseline solution of this scenario, are excluded. All possible combinations of four stars from the data set were tested. From each of the 211 876 individual scenarios, iterative solutions for the rotation parameters were obtained and the WM, WRMS, and ME statistics were calculated, discarding the last 10 iterations because otherwise, the formal errors become too large due to the small sample. The orientation offset WMs range from (0.117, 0.109, -0.030) mas to (0.550, 0.551, 0.222) mas, and the spin WMs range from (0.000, 0.021, -0.043) mas yr$^{-1}$ to (0.106, 0.097, 0.014) mas yr$^{-1}$. Their mean values ((0.339, 0.354, 0.069) mas for the orientation offset and (0.033, 0.095, -0.016) mas yr$^{-1}$ for the spin) are close to the WM values for scenario “60,EDR3,GA,NM” in Table 4. The $Q/n$ for the first iteration in these scenarios (without rejecting further stars) varies between 6.43 and 14.66, with a median of

13.46. The scenario where the sum of all $Q/n$ considered in the calculation of the WMs (i.e., including 35 stars) is minimum is the scenario where T Lep, DoAr 51, CoKu HP Tau G2, and SZ Psc were the four additional excluded stars. These are the same stars that were excluded in the baseline solution of “60,EDR3,GA,NM”, proving that the rejection process is sound. The four stars have RUWE parameters greater than 1.4, indicating that the Gaia data do not fit the standard model of stellar motion well.

The minimum magnitudes<fnref n="5" /> of the WRMSs are 0.100 mas for the orientation offset and 0.009 mas yr$^{-1}$ for the spin, and the minimum magnitudes of the MEs are 0.164 mas and 0.018 mas yr$^{-1}$, respectively. For the orientation offset, the minimum WRMS magnitude is reached when LSI+61303, HD 283641, DoAr51, and CoKu HP Tau G2 are excluded, while for the spin, it is reached when LSI+61303, HD 22468, CoKu HP Tau G2, and AR Lac are rejected. Except LSI+61303 for the orientation offset and CoKu HP Tau G2 for the spin they are the same stars as from the manual selection in “53,EDR3,GA,NM”. This overlap was expected because the stars that were rejected manually in the latter produce jumps and this also deteriorates the WRMSs. From the scatter of these results and the minimum MEs and WRMSs, it is concluded that the minimum uncertainty of WM is ~0.2 mas/√3 ≈ 0.12 mas for each of the three orientation offset components and ~0.02 mas yr$^{-1}$/√3 ≈ 0.01 mas yr$^{-1}$ for each of the three spin components for “60,EDR3,GA,NM”.

The solution adopted for this work is the baseline solution of scenario “60,EDR3,GA,NM”. This solution includes the most complete dataset and, as noted above, has the lowest sum of $Q/n$, thus indicating that it is superior to the other solutions. However, the value of $Q/n$ in this solution (6.43) remains larger than unity, meaning that the uncertainties of the input data are underestimated. To account for this, the formal errors of the rotation parameters have been multiplied by a factor of 2.5 (corresponding to $\sqrt{6.43}$), resulting in more realistic uncertainties. The rotation parameters then become (+0.322, +0.228, +0.163, +0.034, +0.072, -0.026)±(0.203, 0.251, 0.155, 0.023, 0.025, 0.023) mas or mas yr$^{-1}$. Comparing with the baseline solution of the best scenario in Lunz et al. (2023) ((+0.226, +0.327, +0.168, +0.022, +0.065, -0.016) ± (0.165, 0.215, 0.128, 0.024, 0.026, 0.024) mas or mas yr$^{-1}$, after scaling the uncertainties by $\sqrt{5.58}$), the values then do not differ significantly. In these solutions, the spin in the *Y* direction is the only parameter that is found significant at the $3\sigma$ level.

For the bright Gaia EDR3 reference frame, a correction was applied to its spin in one of the middle iterations of the processing, as explained in Sect. 1. The applied spin correction was (-0.0166, -0.0950, +0.0283) mas yr$^{-1}$ with an uncertainty of 0.024 mas yr$^{-1}$ per axis, and no orientation offset was considered. Reversing the applied correction by adding it back to the baseline solution from scenario “60,EDR3,GA,NM” results in an original, uncorrected spin in Gaia EDR3 of (0.016, -0.023, +0.002) ± (0.023, 0.025, 0.023) mas yr$^{-1}$. The spin in *Y* is then reduced to $1-\sigma$, which suggests that the uncorrected Gaia EDR3 bright frame is more consistent with the ICRF3 than the corrected frame. The results also show that the current level of accuracy in aligning the bright Gaia reference frame

<footnote n="5">⁵ Magnitudes in this context are defined as the square root of the quadratically added values for the three rotation axes. This value can be understood as a combined value for the total rotation and should provide a better comparison between the various solutions.</footnote>
spin to ICRF3 is the same as when employing the HIPPARCOS
positions.

The results for the residual spin for the scenarios “55,EDR3,GA”, “55,EDR3,GA,NM”, and “60,EDR3,GA,NM” (i.e., without manually excluding stars before the iterative adjustment process starts), align with the independent estimates of 80 µas yr⁻¹ for the total spin in the magnitude range G = 11–13 mag from Cantat-Gaudin & Brandt (2021), with the main component being the rotation about the Y axis with a value of 60 µas yr⁻¹ to 77 µas yr⁻¹. Their analysis is based on Gaia EDR3 data from open clusters and binaries with a bright and faint component, taking also into account a magnitude dependence of the optically bright frame. Furthermore, the results for the residual spin for the scenarios “50,EDR3,GA”, “49,EDR3,GA,NM”, and “53,EDR3,GA,NM” (with manually excluding stars that produce offsets in the rotation parameters before the iterative adjustment process starts), align with the independent spin estimates in the magnitude range G ≤ 10.5 mag of about 40 µas yr⁻¹ from the same authors. The largest component in these solutions, of 31 µas yr⁻¹ to 36 µas yr⁻¹, still is the rotation about the Y axis. The stars DoAr 51, CoKu HP Tau G2, and HD 283641, which were excluded a priori because they lead to offsets in the iterative orientation offset parameter results, are fainter than the threshold of G = 10.5 mag. This may indicate that the alignment results in this work are also sensitive to the magnitude. On the other hand, the RUWE parameters of these stars are also larger than 1.4, which means that the standard model of stellar motion does not fit their Gaia data well. Therefore, the effect seen could also well be a consequence of excluding the stars from the analysis. Fifteen stars in the sample have G ≥ 10.5 mag, and 28 stars have a RUWE ≥ 1.4, almost 50%. The stars producing the offsets in spin do not show any conspicuousness from the Gaia data but show acceleration or orbital motion from the VLBI data. In all, it is difficult to test the magnitude dependence of the alignment between Gaia and VLBI because the current VLBI dataset only includes 15 stars with G ≥ 10.5.

Orbital motion and accelerations should be considered in
future work in order to properly model the trajectories of some
of the stars. More accurate and precise VLBI data will help
to refine the model estimates, identify outliers (i.e., problem-
atic stars), and improve the overall reliability of the results.
The addition of only a few positions and model estimates of
a few new stars (Sect. 4.4) did not significantly improve the
rotation parameter analysis as systematic errors dominate over
the thermal errors. Furthermore, for some stars in our sam-
ple, the correlations between the astrometric parameters in the
model of stellar motion were large. This means these stars need
additional geometrically sensitive observations to decorrelate the
parameters.

The comparisons presented in this study show that the align-
ment of Gaia EDR3 with ICRF3 using VLBI observations of
radio stars will require more effort from the VLBI community
to reach the level of uncertainties expected for the final Gaia
data release (see Sect. 1). It should be emphasized that the align-
ment between VLBI and Gaia for individual objects is only as
good as shown by the comparison between the various VLBI
models and the astrometric Gaia model with five parameters, as
listed in Table A.1. The model values may differ due to, among
other things, orbital motions, differences between the barycen-
ter and center of luminosity of unresolved binary stars as seen
by Gaia, and radio-optical offsets (Brosche & Schuh 1999; Lunz
et al. 2023).

6. Conclusions

This study focused on further improving the determination of
residual orientation offsets and spin between the Gaia EDR3
bright (G ≤ 13 mag) frame and the ICRF3. VLBI observations
of radio stars consistently referenced to ICRF3 were used for
the alignment. It was shown that the influence of Galactocen-
tric acceleration on the rotation parameter analysis is negligible
with the currently available data.

In this work, new observations from January 2020 were
added to the historical time series whenever possible, allowing
large improvements in the stellar motion model estimates due
to the longer time span between the first and last observations.
Thanks to this additional data, it was possible to obtain improved
estimates for 12 stars and new estimates for three stars. Replacing
the proper motion and parallax information of these 12 stars in
the sample of 55 stars from Lunz et al. (2023) and using the Jan-
uary 2020 positions derived based on a point source model for
the calibrators reduced the mean standard deviation for the spin
parameters in the iterative rotation parameter analysis by 10% to
20% and the scatter of the spin parameter in the X direction by
35%. At the same time, the orientation offsets are less well deter-
mined, with an increase of the mean standard deviation by about
15% for these parameters, because a smaller number of VLBI
positions and not the more precise model positions are used in
the adjustment. Adding positions for four new stars and model
estimates for three stars reduced the mean formal errors in the
orientation offset parameters slightly, as expected.

We tested different scenarios in which selected stars were
excluded a priori. This test showed that excluding only a few
stars (namely 6 stars) besides the obvious outliers has a large
impact on the rotation parameter estimates.

A possible direction for further studies is to determine
whether the stars producing offsets in the iterative solu-
tion results should be excluded from the analysis, differently
weighted, or have their modeling extended (e.g., adding orbital
motion). New VLBI observations carried out in parallel with the
Gaia observations would clearly help to improve the analysis
and increase confidence in the estimated parameters. Addition-
ally, more stars in common between VLBI and Gaia would be
desirable to also study the possible magnitude dependence of
the rotation parameters. Finally, the systematic effects that lead
to inflating the error budget would need to be investigated as
well, and, if possible, reduced to give absolute positions a greater
influence on the determination of the spin.

The adopted scenario is the one that makes the sum of Q/n
over the representative range of the iterative spin solutions min-
imum. The selected baseline solution for this scenario shows
orientation offsets between Gaia EDR3 and ICRF3 on the order
of 0.2 mas to 0.3 mas and spins of about 0.03 mas yr⁻¹ in X and
Z and about 0.07 mas yr⁻¹ in Y. Among these, only the spin
component in Y is statistically significant (at the 3σ level).

Acknowledgements. We deeply thank the anonymous reviewer for the thorough feedback on our work. The authors acknowledge the use of the Very Long Base-line Array under the US Naval Observatory's time allocation. This work supports USNO's ongoing research into the celestial reference frame and geodesy. We thank also the Socorro correlator for reliably and quickly providing the correlated data. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. We acknowledge the use of two large radio telescopes: Tianma65 and Parkes. The Parkes radio telescope is part of the Australia Telescope National Facility which is funded by the Australian Government for operation as a
National Facility managed by CSIRO. The observations of radio stars were coordinated within the framework of the Asia-Oceania VLBI Group for Geodesy and Astrometry (AOV). The scheduling and data correlation were supported by three AOV member institutes which include Shanghai Astronomical Observatory of Chinese Academy of Sciences, University of Tasmania in Australia and Geospatial Information Authority of Japan. Experiments V515C and AUA020 were done with participation of the Scientific Equipment Sharing Center of the Quasar VLBI Network of the Institute of Applied Astronomy of the Russian Academy of Sciences (IAA RAS). We further thank Belén Tercero and Chris Phillips for their contribution. This project is supported by the DFG Grant Nos. SCHU 1103/7-2 and HE 5937/2-2. M.H. Xu was supported by the ERC StG Grant No. 101076060. This work has made use of the data from the European Space Agency (ESA) mission Gaia processed by the Gaia Data Processing and Analysis Consortium as well as from the mission HIPPARCOS. Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Calibrators were selected using the NRAO VLBA calibrator search tool (www vlba nrao edu/astro/calib/), the RFC calibrator search tool (http://astrogeo org/calib/search.html) and the Astrogeo VLBI FITS image database (http://astrogeo org/vlbi_images/). This research has made use of the VizieR catalog access tool and the SIMBAD database (Wenger et al. 2000), operated at CDS, Strasbourg, France. Calculations were made in MATLAB by The MathWorks, Inc.

## References

Arias, E. F., Charlot, P., Feissel, M., & Lestrade, J. F. 1995, A&A, 303, 604

Bartel, N., Bietenholz, M. F., Lebach, D. E., et al. 2012, ApJS, 201, 3

Brosche, P., & Schuh, H. 1999, ZFV – Zeitsch. Geodasie Geoinform. Landmanag., 124, 343

Cantat-Gaudin, T., & Brandt, T. D. 2021, A&A, 649, A124

Charlot, P., Jacobs, C. S., Gordon, D., et al. 2020, A&A, 644, A159

Condon, J. J. 1997, PASP, 109, 166

ESA 1997, ESA Special Publication, 1200, The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission

Fekel, F. C., J. 1983, ApJ, 268, 274

Folkner, W. M., Williams, J. G., & Boggs, D. H. 2009, IPN Progress Report, 178

Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2021a, A&A, 649, A1

Gaia Collaboration (Klioner, S. A., et al.) 2021b, A&A, 649, A9

Gaia Collaboration (Klioner, S. A., et al.) 2022, A&A, 667, A148

Galli, P. A. B., Loinard, L., Ortiz-Léon, G. N., et al. 2018, ApJ, 859, 33

García-Alvarez, D., Foing, B. H., Montes, D., et al. 2003, A&A, 397, 285

Harris, R. J., Andrews, S. M., Wilner, D. J., & Kraus, A. L. 2012, ApJ, 751, 115

Kounkel, M., Hartmann, L., Loinard, L., et al. 2017, ApJ, 834, 142

Kovalevsky, J., Lindegren, L., Perryman, M. A. C., et al. 1997, A&A, 323, 620

Lebach, D. E., Bartel, N., Bietenholz, M. F., et al. 2012, ApJS, 201, 4

Lestrade, J. F., Preston, R. A., Jones, D. L., et al. 1999, A&A, 344, 1014

Lindegren, L. 2020a, A&A, 633, A1

Lindegren, L. 2020b, A&A, 637, C5

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2

Lindegren, L., Bastian, U., Biermann, M., et al. 2021a, A&A, 649, A4

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021b, A&A, 649, A2

Loinard, L., Torres, R. M., Mioduszewski, A. J., et al. 2007, ApJ, 671, 546

Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29

Lunz, S., Anderson, J. M., Xu, M. H., et al. 2023, A&A, 676, A11

MacMillan, D. S., Fey, A., Gipson, J. M., et al. 2019, A&A, 630, A93

Marschall, L. A., & Mathieu, R. D. 1988, AJ, 96, 1956

Marsden, S. C., Berdyugina, S. V., Donati, J.-F., et al. 2005, ApJ, 634, L173

Miller-Jones, J. C. A., Sivakoff, G. R., Knigge, C., et al. 2013, Science, 340, 950

Nguyen, D. C., Brandeker, A., van Kerkwijk, M. H., & Jayawardhana, R. 2012, ApJ, 745, 119

Ortiz-León, G. N., Loinard, L., Kounkel, M. A., et al. 2017, ApJ, 834, 141

Ratner, M. I., Bartel, N., Bietenholz, M. F., et al. 2012, ApJS, 201, 5

Richichi, A., Leinert, C., Jameson, R., & Zinnecker, H. 1994, A&A, 287, 145

Secrest, N. J., Dudik, R. P., Dorland, B. N., et al. 2015, ApJS, 221, 12

Secrest, N. J., Dudik, R. P., Dorland, B. N., et al. 2016, VizieR Online Data Catalog: J/ApJS/221/12

Seidelmann, P. K. 1992, Explanatory supplement to the Astronomical almanac., rev. edn. (Mill Valley: University Science Books)

Shepherd, M. C. 1997, in Astronomical Society of the Pacific Conference Series, 125, Astronomical Data Analysis Software and Systems VI, eds. G. Hunt, & H. Payne, 77

Shuygina, N., Ivanov, D., Ipatov, A., et al. 2019, Geodesy Geodyn., 10, 150

Titov, O., Lambert, S. B., & Gontier, A. M. 2011, A&A, 529, A91

Titov, O., Shu, F., & Chen, W. 2020, in Astrometry, Earth Rotation, and Reference Systems in the GAIA era, ed. C. Bizouard, 173

Torres, R. M., Loinard, L., Mioduszewski, A. J., & Rodriguez, L. F. 2007, ApJ, 671, 1813

Torres, R. M., Loinard, L., Mioduszewski, A. J., & Rodríguez, L. F. 2009, ApJ, 698, 242

van Leeuwen, F. 2007, A&A, 474, 653

Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9
Appendix A: Tables and Figures

<table><caption>Table A.1: Newly derived astrometric models of stellar motion for 15 stars.</caption><thead><tr><th>model</th><th>E<sub>α*</sub><br>[μas]</th><th>E<sub>δ</sub><br>[mas]</th><th>RMS<sub>v,α*</sub><br>[mas]</th><th>RMS<sub>v,δ</sub><br>[Julian year]</th><th>Σn</th><th>t<sub>0</sub></th><th>α<sub>0</sub><br>[°]</th><th>σ<sub>α*0</sub><br>[mas]</th><th>δ<sub>0</sub><br>[°]</th><th>σ<sub>δ0</sub><br>[mas]</th><th>ω</th><th>σω</th><th>μ<sub>α*</sub><br>[mas yr<sup>-1</sup>]</th><th>σμ<sub>α*</sub><br>[mas yr<sup>-1</sup>]</th><th>μ<sub>δ</sub><br>[mas yr<sup>-1</sup>]</th><th>σμ<sub>δ</sub></th></tr></thead><tbody><tr><td colspan="16"><strong>HD 283572</strong></td></tr><tr><td>ref. (1,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2005.3600</td><td>65.495216803</td><td>0.032</td><td>28.3017697510.035</td><td>...</td><td>...</td><td>7.841</td><td>0.057</td><td>9.023</td><td>0.061</td><td>-26.445</td><td>0.077</td></tr><tr><td>ref. (1)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>7.722</td><td>0.057</td><td>8.853</td><td>0.096</td><td>-26.491</td><td>0.113</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>65.495246613</td><td>0.018</td><td>28.3016916180.009</td><td>...</td><td>...</td><td>7.873</td><td>0.019</td><td>8.837</td><td>0.027</td><td>-26.426</td><td>0.017</td></tr><tr><td>M<sub>old</sub></td><td>122</td><td>185</td><td>0.102</td><td>0.161</td><td>6</td><td>2005.36121</td><td>65.495216801</td><td>0.054</td><td>28.3017696800.080</td><td>...</td><td>...</td><td>7.736</td><td>0.068</td><td>8.865</td><td>0.114</td><td>-26.502</td><td>0.190</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>105</td><td>165</td><td>0.097</td><td>0.159</td><td>7</td><td>2007.45467</td><td>65.495222674</td><td>0.043</td><td>28.3017543390.068</td><td>...</td><td>...</td><td>7.745</td><td>0.059</td><td>8.892</td><td>0.008</td><td>-26.386</td><td>0.013</td></tr><tr><td colspan="16"><strong>V410 Tau</strong></td></tr><tr><td>ref. (1)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>7.751</td><td>0.027</td><td>8.703</td><td>0.017</td><td>-24.985</td><td>0.020</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>64.629661927</td><td>0.019</td><td>28.4543760790.009</td><td>...</td><td>...</td><td>7.730</td><td>0.021</td><td>8.846</td><td>0.025</td><td>-25.129</td><td>0.016</td></tr><tr><td>M<sub>old</sub></td><td>116</td><td>79</td><td>0.102</td><td>0.075</td><td>9</td><td>2015.75818</td><td>64.629661280</td><td>0.040</td><td>28.4543778750.029</td><td>...</td><td>...</td><td>7.749</td><td>0.043</td><td>8.705</td><td>0.028</td><td>-24.985</td><td>0.024</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>199</td><td>93</td><td>0.175</td><td>0.100</td><td>10</td><td>2016.18390</td><td>64.629662471</td><td>0.064</td><td>28.4543749150.033</td><td>...</td><td>...</td><td>7.707</td><td>0.067</td><td>8.786</td><td>0.033</td><td>-25.010</td><td>0.022</td></tr><tr><td colspan="16"><strong>SS Cyg</strong></td></tr><tr><td>ref. (3,2)</td><td>108</td><td>91</td><td>...</td><td>...</td><td>2011.5661</td><td>325.678846337</td><td>0.065</td><td>43.5861813920.070</td><td>...</td><td>...</td><td>8.800</td><td>0.120</td><td>112.420</td><td>0.070</td><td>33.380</td><td>0.070</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>325.679037286</td><td>0.022</td><td>43.5862224690.026</td><td>...</td><td>...</td><td>8.854</td><td>0.030</td><td>112.385</td><td>0.029</td><td>33.315</td><td>0.033</td></tr><tr><td>M<sub>old</sub></td><td>164</td><td>0</td><td>0.229</td><td>0.145</td><td>9</td><td>2011.22802</td><td>325.678831631</td><td>0.101</td><td>43.5861873930.050</td><td>...</td><td>...</td><td>8.841</td><td>0.094</td><td>112.4490</td><td>0.064</td><td>33.399</td><td>0.049</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>184</td><td>0</td><td>0.243</td><td>0.136</td><td>10</td><td>2012.10667</td><td>325.678869501</td><td>0.081</td><td>43.5861865350.043</td><td>...</td><td>...</td><td>8.896</td><td>0.054</td><td>112.415</td><td>0.019</td><td>33.353</td><td>0.014</td></tr><tr><td colspan="16"><strong>Brun 334</strong></td></tr><tr><td>ref. (4,2)<sup>(2)</sup></td><td>...</td><td>...</td><td>...</td><td>...</td><td>2015.1800</td><td>83.665666567</td><td>0.023</td><td>-5.4071122210.048</td><td>...</td><td>...</td><td>2.591</td><td>0.046</td><td>-4.010</td><td>0.080</td><td>-1.170</td><td>0.070</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>83.665665932</td><td>0.009</td><td>-5.4071124700.008</td><td>...</td><td>...</td><td>2.569</td><td>0.014</td><td>-3.736</td><td>0.012</td><td>-0.887</td><td>0.010</td></tr><tr><td>M<sub>old</sub></td><td>150</td><td>132</td><td>0.061</td><td>0.093</td><td>5</td><td>2015.17864</td><td>83.665666765</td><td>0.069</td><td>-5.4071122500.063</td><td>...</td><td>...</td><td>2.524</td><td>0.070</td><td>-3.943</td><td>0.098</td><td>-1.221</td><td>0.094</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>169</td><td>226</td><td>0.132</td><td>0.181</td><td>6</td><td>2015.98480</td><td>83.665665922</td><td>0.072</td><td>-5.4071124660.097</td><td>...</td><td>...</td><td>2.502</td><td>0.079</td><td>-3.766</td><td>0.037</td><td>-0.995</td><td>0.055</td></tr><tr><td colspan="16"><strong>TYC 5346-538-1</strong></td></tr><tr><td>ref. (4,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2015.2000</td><td>85.640319910</td><td>0.058</td><td>-8.120884301</td><td>0.140</td><td>...</td><td>2.348</td><td>0.069</td><td>0.680</td><td>0.090</td><td>-0.510</td><td>0.250</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>85.640320285</td><td>0.009</td><td>-8.120884072</td><td>0.010</td><td>...</td><td>2.363</td><td>0.013</td><td>0.580</td><td>0.012</td><td>-0.189</td><td>0.012</td></tr><tr><td>M<sub>old</sub></td><td>0</td><td>114</td><td>0.012</td><td>0.013</td><td>5</td><td>2015.19836</td><td>85.640320180</td><td>0.022</td><td>-8.1208840150.078</td><td>...</td><td>...</td><td>2.348</td><td>0.026</td><td>-0.686</td><td>0.043</td><td>-0.524</td><td>0.121</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>61</td><td>264</td><td>0.076</td><td>0.235</td><td>6</td><td>2016.00123</td><td>85.640320310</td><td>0.037</td><td>-8.120884063</td><td>0.126</td><td>...</td><td>2.328</td><td>0.039</td><td>0.584</td><td>0.024</td><td>-0.225</td><td>0.075</td></tr><tr><td>M<sub>old,orig</sub></td><td>139</td><td>398</td><td>0.111</td><td>0.302</td><td>5</td><td>2015.19836</td><td>85.640320157</td><td>0.070</td><td>-8.120884022</td><td>0.189</td><td>...</td><td>2.336</td><td>0.071</td><td>0.739</td><td>0.099</td><td>-0.092</td><td>0.274</td></tr><tr><td>M<sub>new,orig</sub></td><td>167</td><td>331</td><td>0.142</td><td>0.283</td><td>6</td><td>2016.00123</td><td>85.640320289</td><td>0.076</td><td>-8.120884058</td><td>0.151</td><td>...</td><td>2.336</td><td>0.083</td><td>0.603</td><td>0.041</td><td>-0.166</td><td>0.086</td></tr><tr><td colspan="16"><strong>BH CVn</strong></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1993.1088</td><td>203.698997828</td><td>0.373</td><td>37.1824333340.455</td><td>...</td><td>...</td><td>22.210</td><td>0.450</td><td>85.496</td><td>0.131</td><td>-9.220</td><td>0.160</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>203.699680762</td><td>0.044</td><td>37.1823719550.048</td><td>...</td><td>...</td><td>21.327</td><td>0.089</td><td>85.609</td><td>0.043</td><td>-9.549</td><td>0.051</td></tr><tr><td>M<sub>old</sub></td><td>1043</td><td>1373</td><td>1.656</td><td>1.649</td><td>14</td><td>1992.05054</td><td>203.698966582</td><td>0.580</td><td>37.1824358780.580</td><td>...</td><td>...</td><td>22.475</td><td>0.683</td><td>85.618</td><td>0.179</td><td>-9.351</td><td>0.228</td></tr><tr><td>M<sub>new</sub></td><td>869</td><td>1422</td><td>1.575</td><td>1.696</td><td>16</td><td>1995.54593</td><td>203.699070754</td><td>0.445</td><td>37.182426495</td><td>0.511</td><td>...</td><td>22.294</td><td>0.655</td><td>85.504</td><td>0.037</td><td>-9.673</td><td>0.044</td></tr><tr><td>M<sub>11340+3754</sub><sup>(1)</sup></td><td>0</td><td>0</td><td>0.687</td><td>0.468</td><td>7</td><td>1997.25288</td><td>203.699121926</td><td>0.282</td><td>37.1824217400.224</td><td>...</td><td>...</td><td>21.578</td><td>0.333</td><td>85.482</td><td>0.018</td><td>-9.696</td><td>0.012</td></tr><tr><td>M<sub>others</sub></td><td>1682053</td><td></td><td>1.282</td><td>2.015</td><td>9</td><td>1994.21831</td><td>203.699031029</td><td>0.614</td><td>37.1824300761.025</td><td>...</td><td>...</td><td>23.092</td><td>1.146</td><td>85.503</td><td>0.055</td><td>-9.618</td><td>0.085</td></tr><tr><td>M<sub>w0.5.new</sub></td><td>785</td><td>1348</td><td>1.575</td><td>1.698</td><td>16</td><td>1995.54593</td><td>203.699070762</td><td>0.443</td><td>37.1824264900.503</td><td>...</td><td>...</td><td>22.268</td><td>0.646</td><td>85.502</td><td>0.036</td><td>-9.675</td><td>0.044</td></tr><tr><td colspan="16"><strong>Haro 1-6</strong></td></tr><tr><td>ref. (5,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2007.9900</td><td>246.512566654</td><td>0.399</td><td>-24.3934464560.356</td><td>...</td><td>...</td><td>7.385</td><td>0.234</td><td>-19.630</td><td>0.190</td><td>-26.920</td><td>0.130</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>246.512518670</td><td>0.026</td><td>-24.393506245</td><td>0.014</td><td>...</td><td>7.344</td><td>0.031</td><td>-19.798</td><td>0.045</td><td>-27.146</td><td>0.033</td></tr><tr><td>M<sub>all,old</sub></td><td>1759</td><td>595</td><td>1.640</td><td>0.561</td><td>19</td><td>2009.45441</td><td>246.512557797</td><td>0.421</td><td>-24.393457135</td><td>0.139</td><td>...</td><td>7.592</td><td>0.435</td><td>-19.606</td><td>0.099</td><td>-27.047</td><td>0.037</td></tr><tr><td>M<sub>all,new</sub></td><td>1716</td><td>589</td><td>1.605</td><td>0.556</td><td>20</td><td>2009.98241</td><td>246.512554617</td><td>0.395</td><td>-24.393461095</td><td>0.134</td><td>...</td><td>7.521</td><td>0.417</td><td>-19.636</td><td>0.082</td><td>-27.034</td><td>0.031</td></tr><tr><td>M<sub>w0,old</sub></td><td>680</td><td>642</td><td>0.633</td><td>0.608</td><td>18</td><td>2009.61057</td><td>246.512556911</td><td>0.167</td><td>-24.393458297</td><td>0.153</td><td>...</td><td>6.961</td><td>0.187</td><td>-19.675</td><td>0.039</td><td>-27.049</td><td>0.041</td></tr><tr><td>M<sub>w0, new</sub><sup>(1)</sup></td><td>659</td><td>628</td><td>0.617</td><td>0.597</td><td>19</td><td>2010.15813</td><td>246.512553621</td><td>0.156</td><td>-24.393462404</td><td>0.146</td><td>...</td><td>6.947</td><td>0.177</td><td>-19.680</td><td>0.032</td><td>-27.037</td><td>0.034</td></tr><tr><td>M<sub>≥2012,old</sub></td><td>340</td><td>220</td><td>0.270</td><td>0.184</td><td>7</td><td>2014.15461</td><td>246.512529655</td><td>0.128</td><td>-24.3934924650</td><td>0.084</td><td>...</td><td>7.385</td><td>0.134</td><td>-19.631</td><td>0.108</td><td>-26.928</td><td>0.083</td></tr><tr><td>M<sub>≥2012, new</sub></td><td>340</td><td>190</td><td>0.286</td><td>0.178</td><td>8</td><td>2014.88707</td><td>246.512525240</td><td>0.120</td><td>-24.3934979500</td><td>0.072</td><td>...</td><td>7.381</td><td>0.136</td><td>-19.732</td><td>0.053</td><td>-26.949</td><td>0.036</td></tr><tr><td colspan="16"><strong>σ<sup>2</sup>CrB</strong></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1990.0014</td><td>243.671114605</td><td>0.104</td><td>33.8588538870</td><td>0.124</td><td>...</td><td>43.930</td><td>0.100</td><td>-267.048</td><td>0.037</td><td>-86.660</td><td>0.050</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>243.668787168</td><td>0.033</td><td>33.8582264310</td><td>0.041</td><td>...</td><td>44.057</td><td>0.046</td><td>-268.216</td><td>0.043</td><td>-87.282</td><td>0.058</td></tr><tr><td>M<sub>old</sub></td><td>0</td><td>0</td><td>0.341</td><td>0.280</td><td>12</td><td>1992.47265</td><td>243.670893774</td><td>0.110</td><td>33.8587943460</td><td>0.092</td><td>...</td><td>43.844</td><td>0.108</td><td>-267.078</td><td>0.047</td><td>-86.671</td><td>0.049</td></tr><tr><td>M<sub>new</sub></td><td>1596</td><td>344</td><td>1.527</td><td>0.543</td><td>13</td><td>1994.59123</td><td>243.670704024</td><td>0.472</td><td>33.8587431740</td><td>0.164</td><td>...</td><td>43.496</td><td>0.280</td><td>-267.725</td><td>0.049</td><td>-86.947</td><td>0.016</td></tr><tr><td>M<sub>new,a</sub><sup>(1)</sup></td><td>0</td><td>0</td><td>0.415</td><td>0.272</td><td>13</td><td>2016.00000</td><td>243.668787173</td><td>0.215</td><td>33.8582263710</td><td>0.206</td><td>...</td><td>43.809</td><td>0.124</td><td>-268.215</td><td>0.037</td><td>-87.122</td><td>0.036</td></tr><tr><td colspan="16"><strong>HD 199178</strong></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1993.7933</td><td>313.473486712</td><td>0.332</td><td>44.3864125650</td><td>397</td><td>...</td><td>8.590</td><td>0.330</td><td>26.595</td><td>0.407</td><td>-1.240</td><td>0.430</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>313.473714681</td><td>0.011</td><td>44.3864072430</td><td>0.013</td><td>...</td><td>8.891</td><td>0.015</td><td>26.451</td><td>0.015</td><td>-0.877</td><td>0.015</td></tr><tr><td>M<sub>old</sub></td><td>0</td><td>0</td><td>0.576</td><td>0.226</td><td>6</td><td>1993.79305</td><td>313.473486760</td><td>0.263</td><td>44.3864123510</td><td>200</td><td>...</td><td>8.637</td><td>0.285</td><td>26.555</td><td>0.389</td><td>-1.221</td><td>0.326</td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>0</td><td>0</td><td>0.548</td><td>0.277</td><td>7</td><td>1997.53897</td><td>313.473525218</td><td>0.210</td><td>44.3864114880</td><td>0.152</td><td>...</td><td>8.768</td><td>0.206</td><td>26.421</td><td>0.008</td><td>-0.819</td><td>0.011</td></tr><tr><td colspan="16"><strong>AR Lac</strong></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1992.4353</td><td>332.170232910</td><td>0.274</td><td>45.7421531880</td><td>361</td><td>...</td><td>23.970</td><td>0.370</td><td>-52.080</td><td>0.126</td><td>47.030</td><td>0.190</td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>332.169742736</td><td>0.018</td><td>45.742460766</td><td>0.017</td><td>...</td><td>23.525</td><td>0.023</td><td>-52.310</td><td>0.021</td><td>46.931</td><td>0.020</td></tr><tr><td>M<sub>old</sub></td><td>0</td><td>0</td><td>0.684</td><td>0.469</td><td>6</td><td>1992.70656</td><td>332.170227205</td><td>0.432</td><td>45.7421566810</td><td>339</td><td>...</td><td>24.125</td><td>0.409</td><td>-52.023</td><td>0.192</td><td>47.297</td><td>0.196</td></tr><tr><td>M<sub>w90.J2202+4216, new</sub></td><td>0</td><td>0</td><td>0.355</td><td>0.248</td><td>6</td><td>1997.62628</td><td>332.170124990</td><td>0.180</td><td>45.742221000</td><td>0.139</td><td>...</td><td>24.275</td><td>0.216</td><td>-52.233</td><td>0.010</td><td>47.032</td><td>0.011</td></tr><tr><td>M<sub>w90,w0.5.new</sub><sup>(1)</sup></td><td>0</td><td>0</td><td>0.346</td><td>0.249</td><td>7</td><td>2000.82461</td><td>332.170058496</td><td>0.164</td><td>45.742262784</td><td>0.121</td><td>...</td><td>24.275</td><td>0.238</td><td>-52.233</td><td>0.011</td><td>47.032</td><td>0.013</td></tr></tbody></table>
<table><caption>Table A.1: continued.</caption><thead><tr><th>model</th><th>E<sub>α*</sub><br>[µas]</th><th>E<sub>δ</sub><br>[mas]</th><th>RMS<sub>v,α*</sub><br>[mas]</th><th>RMS<sub>v,δ</sub><br>[mas]</th><th>Σn<br>[Julian year]</th><th>t<sub>0</sub><br>[°]</th><th>α<sub>0</sub><br>[°]</th><th>σ<sub>α*0</sub><br>[mas]</th><th>δ<sub>0</sub><br>[°]</th><th>σ<sub>δ0</sub><br>[mas]</th><th>σ<sub>ω</sub><br>[mas]</th><th>σ<sub>ω</sub><br>[mas]</th><th>μ<sub>α*</sub><br>[mas yr<sup>-1</sup>]</th><th>σ<sub>μα*</sub><br>[mas yr<sup>-1</sup>]</th><th>μ<sub>δ</sub><br>[mas yr<sup>-1</sup>]</th><th>σ<sub>μδ</sub><br>[mas yr<sup>-1</sup>]</th></tr></thead><tbody><tr><td><strong>IM Peg</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1992.9172</td><td>343.259484359</td><td>0.360</td><td>16.8412478610.392</td><td>10.280</td><td>0.620</td><td>-20.587</td><td>0.459</td><td>-27.530</td><td>0.400</td><td></td></tr><tr><td>ref.(7,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2005.0869</td><td>343.259410883</td><td>0.400</td><td>16.8411555690.390</td><td>10.370</td><td>0.074</td><td>-20.833</td><td>0.090</td><td>-27.267</td><td>0.095</td><td></td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>343.259344951</td><td>0.056</td><td>16.8410726250.049</td><td>10.166</td><td>0.065</td><td>-21.185</td><td>0.063</td><td>-27.444</td><td>0.060</td><td></td></tr><tr><td>M<sub>old</sub></td><td>572 635</td><td>0.558</td><td>0.640 39</td><td>2000.48153</td><td>343.259438677</td><td>0.096</td><td>16.841190373</td><td>0.106</td><td>10.424</td><td>0.117</td><td>-20.847</td><td>0.031</td><td>-27.306</td><td>0.034</td><td></td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>570 631</td><td>0.556</td><td>0.632 40</td><td>2000.96986</td><td>343.259435726</td><td>0.094</td><td>16.841186666</td><td>0.103</td><td>10.408</td><td>0.113</td><td>-20.826</td><td>0.021</td><td>-27.322</td><td>0.024</td><td></td></tr><tr><td>M<sub>orb.new</sub></td><td>382 545</td><td>0.401</td><td>0.584 40</td><td>2000.96986</td><td>343.259435710</td><td>0.064</td><td>16.8411866520.090</td><td>10.352</td><td>0.080</td><td>-20.805</td><td>0.015</td><td>-27.302</td><td>0.021</td><td></td></tr><tr><td><strong>HD 22468</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ref.(6,2)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1992.0000</td><td>54.197113376</td><td>0.406</td><td>0.588121131</td><td>0.401</td><td>33.880</td><td>0.470</td><td>-31.588</td><td>0.330</td><td>-161.690</td><td>0.310</td><td></td></tr><tr><td>ref. (8)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>-31.35</td><td>0.30</td><td>-160.9</td><td>0.3</td><td></td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>54.196900130</td><td>0.027</td><td>0.5870418110.026</td><td>33.978</td><td>0.035</td><td>-32.246</td><td>0.036</td><td>-162.073</td><td>0.032</td><td></td></tr><tr><td>M<sub>old</sub></td><td>0 662</td><td>0.630</td><td>0.967 8</td><td>1992.70183</td><td>54.197107163</td><td>0.312</td><td>0.5880898540.372</td><td>33.838</td><td>0.432</td><td>-31.574</td><td>0.303</td><td>-161.693</td><td>0.338</td><td></td></tr><tr><td>M<sub>new</sub></td><td>438 1125</td><td>0.781</td><td>1.131 10</td><td>1998.01762</td><td>54.197059930</td><td>0.282</td><td>0.587850689</td><td>0.411</td><td>34.307</td><td>0.280</td><td>-31.973</td><td>0.017</td><td>-161.966</td><td>0.033</td><td></td></tr><tr><td>M<sub>abs,wo9&amp;10.new</sub><sup>(1)</sup></td><td>428 1027</td><td>0.689</td><td>1.073 11</td><td>1999.70415</td><td>54.197044944</td><td>0.263</td><td>0.5877748110.366</td><td>34.220</td><td>0.256</td><td>-31.976</td><td>0.016</td><td>-161.965</td><td>0.028</td><td></td></tr><tr><td><strong>CoKu HP Tau G2</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ref. (1)</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>6.145</td><td>0.029</td><td>11.248</td><td>0.022</td><td>-15.686</td><td>0.013</td><td></td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>68.975704541</td><td>0.055</td><td>22.9036789230.032</td><td>5.979</td><td>0.059</td><td>13.459</td><td>0.074</td><td>-11.399</td><td>0.054</td><td></td></tr><tr><td>M<sub>new</sub></td><td>49865379</td><td>4.627</td><td>4.992 18</td><td>2011.81055</td><td>68.975691744</td><td>1.191</td><td>22.903697677</td><td>1.268</td><td>5.769</td><td>1.471</td><td>11.103</td><td>0.221</td><td>-16.085</td><td>0.257</td><td></td></tr><tr><td>M<sub>Gaia</sub><sup>(1)</sup></td><td>1333 1337</td><td>1.018</td><td>1.023 6</td><td>2015.98433</td><td>68.975704603</td><td>0.545</td><td>22.9036787750</td><td>0.547</td><td>6.353</td><td>0.583</td><td>11.830</td><td>0.586</td><td>-10.812</td><td>0.624</td><td></td></tr><tr><td><strong>del Lib</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gaia DR2</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2015.50000</td><td>225.242843762</td><td>0.269</td><td>-8.5189731950.238</td><td>8.328</td><td>0.334</td><td>-60.564</td><td>0.552</td><td>-3.651</td><td>0.457</td><td></td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>225.242834573</td><td>0.290</td><td>-8.5189741080.324</td><td>9.282</td><td>0.472</td><td>-63.051</td><td>0.385</td><td>-6.024</td><td>0.326</td><td></td></tr><tr><td>M<sub>new</sub></td><td>4722297</td><td>0.288</td><td>1.408 4</td><td>2019.28185</td><td>225.242775459</td><td>0.241</td><td>-8.518979228</td><td>1.146</td><td>5.918</td><td>0.313</td><td>-64.055</td><td>0.176</td><td>-8.398</td><td>0.726</td><td></td></tr><tr><td>M<sub>w0.5.new</sub><sup>(1)</sup></td><td>3242282</td><td>0.332</td><td>1.417 4</td><td>2019.28185</td><td>225.242775405</td><td>0.220</td><td>-8.518979190</td><td>1.144</td><td>6.754</td><td>0.345</td><td>-64.420</td><td>0.179</td><td>-8.255</td><td>0.719</td><td></td></tr><tr><td><strong>HD 142184</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gaia EDR3</td><td>...</td><td>...</td><td>...</td><td>...</td><td>2016.00000</td><td>238.482700570</td><td>0.069</td><td>-23.9782025420.050</td><td>6.990</td><td>0.074</td><td>-13.406</td><td>0.083</td><td>-24.111</td><td>0.064</td><td></td></tr><tr><td>M<sub>new</sub><sup>(1)</sup></td><td>1631 402</td><td>1.673</td><td>2.500 12</td><td>2019.18013</td><td>238.482688158</td><td>0.601</td><td>-23.9782239830.259</td><td>7.678</td><td>0.728</td><td>-13.663</td><td>0.536</td><td>-24.345</td><td>0.232</td><td></td></tr></tbody></table>

**Notes.** Estimates for astrometric models of stellar motions M with mean positions at mean epoch $t_0$. The number of individual positions employed for derivation of the respective model is given in column $\Sigma n$ for the fifteen stars discussed in this study. The parameters from literature references (ref.) and Gaia EDR3 are indicated on the first line for each star segment. In case two literature references are given in one line, the first one describes the original data, and the second one points to improved data as given in the respective reference (e.g., positions were shifted to ICRF3 or positions were determined as they were not provided in the first reference). In that case, the parameters reported on the corresponding line are those from the second reference. An ellipsis (...) indicates the omission of an entry. $^{(1)}$Chosen as the best solution for rotation parameter estimation in the following sections. $^{(2)}$Including orbital fit.

**References.** (1) Galli et al. (2018); (2) Lindegren (2020a); (3) Miller-Jones et al. (2013); (4) Kounkel et al. (2017); (5) Ortiz-León et al. (2017); (6) Lestrade et al. (1999); (7) Ratner et al. (2012); (8) Titov et al. (2020).
<table><caption>Table A.2: Star positions from absolute astrometry and alternate observations used in this study.</caption><thead><tr><th>Star</th><th>Calibrator</th><th>Epoch t<br>[Julian date]</th><th>α(t)<br>[h m s]</th><th>σ<sub>α,random</sub><br>[s]</th><th>δ(t)<br>[°′′]</th><th>σ<sub>δ,random</sub><br>[°′]</th><th>Calibrator α<br>[h m s]</th><th>Calibrator δ<br>[°′′]</th><th>Ref.</th><th>Project<br>code</th><th>Network</th></tr></thead><tbody><tr><td>HD 22468</td><td>J0339–0146</td><td>2458318.95833</td><td>3 36</td><td>47.25259</td><td>0.000006</td><td>0 35</td><td>12.9481</td><td>0.00008</td><td>3 39</td><td>30.93778633</td><td>-01 46 35.8041892</td><td>Titov et al. (2020)</td><td>V515C</td></tr><tr><td>HD 22468</td><td>absolute</td><td>2457160.00000</td><td>3 36</td><td>47.2573</td><td>0.00008</td><td>0 35</td><td>13.470</td><td>0.004</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV003</td></tr><tr><td>HD 22468</td><td>absolute</td><td>2457582.00000</td><td>3 36</td><td>47.2562</td><td>0.000034</td><td>0 35</td><td>13.308</td><td>0.002</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AUA011</td></tr><tr><td>HD 22468</td><td>absolute</td><td>2457597.00000</td><td>3 36</td><td>47.2569</td><td>0.000003</td><td>0 35</td><td>13.269</td><td>0.0003</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV010</td></tr><tr><td>E1Eri</td><td>J0408–0529</td><td>2458929.66667</td><td>4 9</td><td>40.9432071</td><td>0.0000007</td><td>-7 53</td><td>32.068041</td><td>0.000039</td><td>4 8</td><td>59.6499</td><td>-5 29</td><td>40.538</td><td>V583A</td></tr><tr><td>YY Men</td><td>J0529–7245</td><td>2459044.02083</td><td>4 58</td><td>17.9152711</td><td>0.00000271</td><td>-75 16</td><td>38.050921</td><td>0.000166</td><td>5 29</td><td>30.0422</td><td>-72 45</td><td>28.507</td><td>LBA</td></tr><tr><td>YY Men</td><td>absolute</td><td>2448236.00000</td><td>4 58</td><td>17.94792</td><td>0.00029</td><td>-75 16</td><td>38.0852</td><td>0.0136</td><td>...</td><td>...</td><td>...</td><td>...</td><td>LBA</td></tr><tr><td>del Lib</td><td>J1456–0617</td><td>2457596.50000</td><td>15 0</td><td>58.27745</td><td>0.000013</td><td>-8 31</td><td>8.3003</td><td>0.00012</td><td>14 56</td><td>41.39265522</td><td>-6 17</td><td>43.2039727</td><td>South-12</td></tr><tr><td>del Lib</td><td>J1512–0905</td><td>2459044.02083</td><td>15 0</td><td>58.2602504</td><td>0.0000031</td><td>-8 31</td><td>8.330971</td><td>0.000161</td><td>15 12</td><td>50.5329</td><td>-9 5</td><td>59.830</td><td>AOV010</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2457875.66760</td><td>15 53</td><td>55.8473903</td><td>0.0000217</td><td>-23 58</td><td>41.562171</td><td>0.000406</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV010</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458142.72890</td><td>15 53</td><td>55.8465556</td><td>0.0000466</td><td>-23 58</td><td>41.580231</td><td>0.000314</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV010</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458282.74890</td><td>15 53</td><td>55.8455544</td><td>0.0002387</td><td>-23 58</td><td>41.586771</td><td>0.008776</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV019</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458341.44780</td><td>15 53</td><td>55.8451101</td><td>0.0000348</td><td>-23 58</td><td>41.593428</td><td>0.001182</td><td>...</td><td>...</td><td>...</td><td>...</td><td>VLBA</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458473.21000</td><td>15 53</td><td>55.8456092</td><td>0.0000647</td><td>-23 58</td><td>41.597420</td><td>0.003519</td><td>...</td><td>...</td><td>...</td><td>...</td><td>VLBA</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458563.75110</td><td>15 53</td><td>55.8455989</td><td>0.0000405</td><td>-23 58</td><td>41.608241</td><td>0.000496</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV033</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458577.75080</td><td>15 53</td><td>55.8454197</td><td>0.0000353</td><td>-23 58</td><td>41.610263</td><td>0.000351</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV034</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458601.58380</td><td>15 53</td><td>55.8452023</td><td>0.0000334</td><td>-23 58</td><td>41.610419</td><td>0.001328</td><td>...</td><td>...</td><td>...</td><td>...</td><td>VLBA</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458862.46170</td><td>15 53</td><td>55.8449085</td><td>0.0000928</td><td>-23 58</td><td>41.624139</td><td>0.003391</td><td>...</td><td>...</td><td>...</td><td>...</td><td>VLBA</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458872.40590</td><td>15 53</td><td>55.8448337</td><td>0.0000232</td><td>-23 58</td><td>41.628742</td><td>0.000865</td><td>...</td><td>...</td><td>...</td><td>...</td><td>VLBA</td></tr><tr><td>HD 142184</td><td>absolute</td><td>2458968.72980</td><td>15 53</td><td>55.8441951</td><td>0.0003340</td><td>-23 58</td><td>41.629517</td><td>0.013609</td><td>...</td><td>...</td><td>...</td><td>...</td><td>AOV046</td></tr><tr><td>HD 142184</td><td>J1553–2422</td><td>2459044.02083</td><td>15 53</td><td>55.8433673</td><td>0.0000022</td><td>-23 58</td><td>41.638212</td><td>0.000133</td><td>15 53</td><td>31.6278</td><td>-24 22</td><td>6.036</td><td>V583B</td></tr><tr><td>HD 167971</td><td>J1818–1108</td><td>2458929.66667</td><td>18 18</td><td>5.8926205</td><td>0.0001457</td><td>-12 14</td><td>33.308782</td><td>0.001516</td><td>18 18</td><td>19.3166</td><td>-11 8</td><td>48.329</td><td>V583A</td></tr><tr><td>HD 167971</td><td>J1832–1035</td><td>2459044.02083</td><td>18 18</td><td>5.8798553</td><td>0.0000224</td><td>-12 14</td><td>33.360195</td><td>0.000321</td><td>18 32</td><td>20.8364</td><td>-10 35</td><td>11.200</td><td>LBA</td></tr><tr><td>V479 Set</td><td>J1838–1833</td><td>2458929.66667</td><td>18 26</td><td>15.0661506</td><td>0.0000003</td><td>-14 50</td><td>54.408599</td><td>0.000022</td><td>18 38</td><td>37.6916</td><td>-18 33</td><td>55.260</td><td>V583B</td></tr></tbody></table>

**Notes.** The star coordinates α(t) and δ(t) at epoch t and their respective uncertainties σ<sub>α,random</sub> and σ<sub>δ,random</sub>, information about the observing network as well as the position type (“absolute” in case of absolute geodetic VLBI and a calibrator name in case of phase referencing observations) are reported as collected from the respective references or derived from our work. The original calibrator coordinates “Calibrator α” and “Calibrator δ” taken from the references listed in column 10 are also given. (1) The network was enlarged with additional stations. Details can be found in the IVS master schedules available at https://ivscc.gsfc.nasa.gov/program/master.html.
<figure><img src="image_4.png" /><figcaption>(a) HD 283572.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>(b) V410 Tau.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>(c) SS Cyg.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>(d) Brun 334.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(e) TYC 5346-538-1.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(f) BH CVn.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(g) Haro 1-6.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(h) $\sigma^2$ CrB.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(i) HD 199178.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>(j) AR Lac.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>(k) IM Peg.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>(l) HD 22468.</figcaption></figure>

Fig. A.1: Sky motion for the 15 stars with astrometric models reported in Table A.1. The observed positions are marked with the red stars and their uncertainties are visualized by the black ellipses. The adjusted positions are labeled by the black dots and their uncertainties by the blue ellipses. The model is indicated by a black line.
<figure><img src="image_16.png" /><figcaption>(m) CoKu HP Tau G2.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>(n) del Lib.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>(o) HD 142184.</figcaption></figure>

Fig. A.1: continued.

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

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

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

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

Fig. A.2: Results for scenario “55,EDR3,GA,NM”, using the VLBI data and Gaia EDR3 as in Fig. 1, but replacing the VLBI data of stars HD 283572, V410 Tau, SS Cyg, Brun 334, TYC 5346-538-1, Haro 1-6, BH CVn, $\sigma^2$ CrB, HD 199178, AR Lac, IM Peg, and HD 22468 with newly determined models of stellar motion and newly corrected absolute positions as described in Sect. 4.3.
<figure><img src="image_23.png" /><figcaption>Fig. A.3: Results for scenario “60,EDR3,GA,NM”, using the VLBI data and Gaia EDR3 as in Fig. A.2, but adding the data of five new stars and astrometric information as described in Sect. 4.4.</figcaption></figure>
<figure><img src="image_24.png" /></figure>

Fig. A.4: Results for scenario “53,EDR3,GA,NM”, using the VLBI data and Gaia EDR3 as in Fig. A.3, but excluding seven stars (HD 22468, LSI+61 303, V410 Tau, AR Lac, DoAr 51, CoKu HP Tau G2, and HD 283641) from the beginning.