1. The Sample
We performed high spatial resolution observations six Be stars during Cycle 29 using HST/FGS (HST proposal #16868 PI: S. Caballero-Nieves). We identified each of the six stars (Table 1) spectroscopically as Be stars via BeSS (Be Star Spectral catalog: Neiner 2018) showing emission lines in the visible or infrared. This study focuses on wide companions of Be stars that are possible triple systems where the inner binary has undergone some transfer of angular momentum in order to produce the decretion disk on the Be star.
2. Literature Search
None of the stars in our sample have been previously identified as having companions. Of note, HD 50737, HD 52244, and HD 206135 are identified as having a proper motion anomaly (Kervella et al. Reference Kervella, Arenou, Mignard and Thévenin2019). Furthermore, Cruzalèbes et al. (Reference Cruzalèbes, Petrov, Robbe-Dubois, Varga, Burtscher, Allouche, Berio, Hofmann, Hron, Jaffe, Lagarde, Lopez, Matter, Meilland, Meisenheimer, Millour and Schertl2019) observed HD 50737 and HD 206135 interferometricaly using the VLTI and detected no companion, Bodensteiner, J. et al. (2020) via literature review determined HD 35345 is presumably single, Peter, D. et al. (2012) used a luck imaging technique to observe HD 239618 and did not detect a companion.
2.1. Gaia RUWE
Castro-Ginard et al. (Reference Castro-Ginard, Penoyre, Casey, Brown, Belokurov, Cantat-Gaudin, Drimmel, Fouesneau, Khanna, Kurbatov, Price-Whelan, Rix and Smart2024) suggest a Gaia DR3 RUWE value of
$\gt 1.4$
to indicate a positional anomaly due to a binary companion. The decretion disk around Be stars can add complications to the RUWE value due to their intrinsic variability. Fitton et al. (Reference Fitton, Tofflemire and Kraus2022) used a NEOWISE photometry criterion of W1-W3
$\gt1$
, to identify stars with a disk and adopts a value of RUWE
$\gt$
2.5 as the threshold for potential binarity for such stars. Five out of six stars in our sample meet the NEOWISE color criterion for having a disk and have a RUWE
$\lt2.5$
suggesting no positional anomaly due to binarity. The star that does not meet the NEOWISE color criterion, HD 206135, has a RUWE
$\lt1.4$
. Therefore, none of our targets have a RUWE value that suggests the presence of a companion.
3. Multiplicity Results
HST/FGS is a a white light shearing interferometer predominantly used for the guiding of HST. This style of interferometer scans and produces two orthogonal S-curves seen in the center plots of Fig. 1 (Solid Line). Of the six stars observed one showed a non-single S-curve. In order to fit the parameters of HD 52244, we used two main methods. The first is a
$\chi^2$
fitting algorithm from the Space Telescope Science Institute (STScI) (Nelan and et al. Reference Nelan2011) where two copies of a single star “calibrator” S-curve and scales and shifts/scales them to fit the observed S-curve. The other algorithm we used employed the use of a cross-correlation function, where the “calibrator” is cross correlated with the observed S-curve. These two methods agree for the y-axis. For the x-axis the separation is small and the CCF fails to identify a companion so we used a second derivative test to check for the presence of a companion but due to the nature of this method we do not receive numerical results (separation and differential magnitude) for the companion. The results of these algorithms are presented in table 1.
Cross-correlation Function (CCF) and second derivative test for HD 52244 in x-axis (left) and y-axis (right). The top row is the mean residuals (grey) from the CCF routine and the observed S-curve with the detected component removed (black). Vertical lines show the location of a detected component. The middle row is the fit S-curve (dashed line) and the observed S-curve (solid line) from the CCF. The bottom row shows the results from the second derivative test with the black line being observed minus residuals, the grey region the standard deviation of the calibrators, and the red line being the second derivative.

Figure 1 Long description
Panel A: The left column contains three line graphs. The top graph shows the mean residuals from the cross-correlation function (CCF) routine in grey and the observed S-curve with the detected component removed in black. The vertical dashed line indicates the location of a detected component. The middle graph displays the fit S-curve as a dashed line and the observed S-curve as a solid line from the CCF. The bottom graph presents the results from the second derivative test, with the black line representing observed minus residuals, the grey region indicating the standard deviation of the calibrators, and the red line showing the second derivative. Panel B: The right column contains three line graphs similar to those in Panel A but for the y-axis. The top graph shows the mean residuals from the CCF routine in grey and the observed S-curve with the detected component removed in black. The vertical dashed line marks the location of a detected component. The middle graph displays the fit S-curve as a dashed line and the observed S-curve as a solid line from the CCF. The bottom graph presents the results from the second derivative test, with the black line representing observed minus residuals, the grey region indicating the standard deviation of the calibrators, and the red line showing the second derivative.
Plx. from Gaia Collaboration (Reference Collaboration2023) (A) Levenhagen and Leister (Reference Levenhagen and Leister2006) (B) Guetter (Reference Guetter1968) (C) Morgan et al. (Reference Morgan, Whitford and Code1953) (D) Boulon and Fehrenbach (Reference Boulon and Fehrenbach1959) (E) Racine (Reference Racine1968)

Table 1 Long description
A table listing astronomical data for various stars. The table has 6 rows and 9 columns. Column headers are: HD Number, RA (J2000), Dec (J2000), Spectral Class, V (mag), B-V (mag), Plx. mas, RUWE, W1-W3 (mag). Row 1: HD Number, 33232; RA (J2000), 05 10 48.2; Dec (J2000), +41 00 10.4; Spectral Class, B2Vne(B); V (mag), 8.16; B-V (mag), 0.07; Plx. mas, 0.76; RUWE, 0.86; W1-W3 (mag), 2.00. Row 2: HD Number, 35345; RA (J2000), 05 25 44.8; Dec (J2000), +35 38 49.9; Spectral Class, B1Vpe(A); V (mag), 8.43; B-V (mag), 0.11; Plx. mas, 0.88; RUWE, 0.96; W1-W3 (mag), 1.60. Row 3: HD Number, 50737; RA (J2000), 06 53 52.2; Dec (J2000), -13 11 09.1; Spectral Class, B2Vnne(A); V (mag), 8.68; B-V (mag), -0.04; Plx. mas, 0.79; RUWE, 1.02; W1-W3 (mag), 1.51. Row 4: HD Number, 52244; RA (J2000), 06 59 46.5; Dec (J2000), -16 12 02.7; Spectral Class, B2IV/Vnde(A); V (mag), 9.68; B-V (mag), -0.06; Plx. mas, 0.58; RUWE, 1.84; W1-W3 (mag), 1.15. Row 5: HD Number, 239618; RA (J2000), 21 14 45.5; Dec (J2000), +59 45 39.7; Spectral Class, B2Ve(C); V (mag), 8.65; B-V (mag), 0.40; Plx. mas, 1.01; RUWE, 1.85; W1-W3 (mag), 1.20. Row 6: HD Number, 206135; RA (J2000), 21 36 57.1; Dec (J2000), +68 11 07.3; Spectral Class, B3V(E); V (mag), 8.27; B-V (mag), 0.12; Plx. mas, 1.10; RUWE, 0.92; W1-W3 (mag), 0.40.
Results of the STScI and Cross-correlation algorithms to check for binary companions Notes:
$^a$
This
$\Delta \text{m}$
is fit in the y-axis and then held for the x-axis.
$^b$
These results use the
$\Delta x$
from the STScI code since it was not resolved with the CCF or second derivative test.

Table 2 Long description
The table presents a comparison of results from the STScI and Cross-correlation algorithms for checking binary companions. It has two rows and eight columns. The columns are labeled as follows: Delta x [mas], Delta y [mas], rho [mas], rho [AU], theta, Delta m_x, Delta m_y, and Delta m. The first row corresponds to the STScI algorithm and the second row to the Cross-correlation algorithm. Row 1: Delta x, -13.8 ± 0.2; Delta y, -39.8 ± 0.2; rho [mas], 42.2; rho [AU], 73.3; theta, 144.5°; Delta m_x, 1.91^a; Delta m_y, 1.91; Delta m, 1.91^a. Row 2: Delta x, -41.0 ± 0.2; Delta y, not applicable; rho [mas], 43.3; rho [AU], 75.1; theta, 143.9°^b; Delta m_x, 1.90; Delta m_y, not applicable; Delta m, not applicable.
4. Summary
Of the six Be stars we observed, HD 52244, had a detection of an unknown companion with a separation of 74 AU (42 mas) at a position angle of
$144^\circ\pm0.3$
E of N and a differential magnitude in the Hubble F583W filter of 1.91. This observation used the F583W filter, that is much wider than the standard Johnson V filter. We assume a comparable differentiation magnitude in the V-Band and that the Be star is on the main sequence. Using
$\Delta m_V =1.91$
and a
$M_V=-1.5$
from Pecaut and Mamajek (Reference Pecaut and Mamajek2013), we arrive at the companion being a late B7-8 V star with
$M_V=0.41$
. With a separation of 74 AU and a total mass of the system of 9.17
$M_\odot$
. A circular orbit for this system would have a period of about 210 years, and follow-up observations with FGS could detect such motion.







a
Δm
b
Δx