Dynamical Masses and Radiative Transfer Modeling of HD 698: a Be Binary in Evolutionary Transition

TL;DR

HD 698 is resolved as a Be primary (M_Be ≈ 7.48 M_⊙) accompanied by a bloated, hydrogen-poor stripped companion (M_comp ≈ 1.23 M_⊙) in a near-circular orbit with P ≈ 55.93 d. By combining high-resolution spectroscopy, long-baseline interferometry, and three-component radiative-transfer modeling, the study derives a dynamical distance of d ≈ 888 pc and constrains the companion to T_eff ≈ 10.0 kK and R ≈ 13.1 R_⊙, with the disk characterized by ρ_0 ≈ 5×10^−12 g cm^−3 and n ≈ 3.0. The SED-fitting supports a post-mass-transfer transitional phase for HD 698, while Balmer/metal-line diagnostics indicate a hydrogen-poor, CNO-processed atmosphere; discrepancies point to non-solar abundances and complex disk irradiation effects. These results place HD 698 among the growing class of Be+bloated binaries, offering rare empirical constraints on the immediate aftermath of binary mass transfer and informing the evolutionary pathways toward Be+sdO/B systems like φ Persei. The work also demonstrates a framework for jointly constraining Be binaries with stripped companions through spectro-interferometric and radiative-transfer techniques.

Abstract

We present a detailed analysis of the early post-mass-transfer binary HD 698 (V742 Cas) combining high-resolution optical spectroscopy, long-baseline interferometry, and radiative-transfer modeling. Counter-phased radial-velocity curves yield a circular orbit with P=55.927+/-0.001 d and component masses M_Be=7.48+/-0.07 M_sun and M_comp=1.23+/-0.02 M_sun. The Be primary is traced by broad H alpha wings, while narrow metallic absorption lines arise from a slowly rotating companion. The interferometric separation implies a dynamical distance of 888+/-5 pc. The spectral energy distribution is reproduced with E(B-V)=0.321+/-0.016 and a viscous decretion disk of base density rho_0~5x10^-12 g cm^-3 at r=R_eq, declining radially as rho(r)~r^-n with n=3.0. The companion is luminous and inflated, with T_eff=10.0(+0.2,-0.1) kK, R_comp=13.1+/-0.2 R_sun, and log(L/L_sun)=3.19, contributing significantly to the flux (L_comp/L_Be~0.3). Spectral line mismatches further suggest a hydrogen-poor, CNO-processed atmosphere, consistent with a stripped-envelope star. HD 698 thus adds to the emerging class of Be+bloated OB binaries, capturing a brief post-mass-transfer phase when the donor remains spectroscopically detectable prior to the subdwarf stage.

Figures (16)

• Figure 1: Top panels: 78 spectra (grey) of the 4450--4545 Å (upper plot), 4815--4890 Å (middle plot), and 6520--6690 Å regions (lower plot), overlaid with an averaged spectrum (blue). Bottom panels: phase-resolved spectra of the same regions. Blue dashed line demonstrates motion of the companion, red dashed line---motion of the Be star.
• Figure 2: Phase-folded RV curves: companion's RVs (circles), Be star's RVs (triangles), with fits overlaid (solid blue line for the companion and solid red line for the Be star). The average error of the companion's RV measurements is 1.67 km s$^{-1}$, of the Be star's---2.07 km s$^{-1}$.
• Figure 3: Surface density map for a Be binary model with a 50-day orbital period, $q = 0.16$, and disk viscosity parameter $\alpha = 1.0$. The white lines trace the Roche equipotential contours, with the two large black circles representing the Be star (in x, y=0.0) and the companion (x=25.6 Req, y=0.0 Req). The small black dots trace the two spiral arms formed in the Be disk due to the presence of the companion. More details about the simulation can be found in 2025AA...698A.309R.
• Figure 4: Top panel: The total modeled SED (purple solid line) is the sum of the Be star and its disk (yellow solid line) and companion's (green dashed line) contributions, both generated with HDUST. Observational data is shown in black. Yellow dashed line represents flux contribution from the disk. Bottom panel: residuals distribution across wavelengths.
• Figure 5: Distribution of SED $\chi^2$ values as a function of base density for HDUST models at two distances (703 pc and 888 pc) and two density slope parameters (n=3.0 and n=3.5). Each panel shows the same data colored by different free parameters: companion's effective temperature (left panel), companion's radius (central panel), and interstellar extinction $E(B-V)$ (right panel). Diamonds represent models at 703 pc, circles represent models at 888 pc. Larger symbols correspond to n=3.0, smaller symbols to n=3.5. Gray lines connect models with identical distance and slope parameters to highlight density-dependent trends. Lower $\chi^2$ values indicate better fits to the observed SED.