The enigmatic dance of the HD 189733A system: A quest for accretion^★

INAF – Osservatorio Astronomico di Palermo, 90123 Palermo PA, Italy
e-mail: salvatore.colombo@inaf.it

Received: 16 October 2023
Accepted: 6 January 2024

Abstract

Context. Several studies suggest that the emission properties of a star can be significantly affected by its interaction with a nearby planet through magnetic fields or interaction between their respective winds. However, the actual observability of these effects remains a subject of debate. An illustrative example is the HD 189733A system: certain characteristics of its emissions have been interpreted as indicative of ongoing interactions between the star and its associated planet. Other studies attribute these characteristics to the coronal activity of the star.

Aims. In this study we aimed to investigate whether the observed stellar X-ray flare events, which appear to be in phase with the planetary period in the HD 189733A system, could be attributed to the accretion of the planetary wind onto the stellar surface or if they resulted from an interaction between the planetary and stellar winds.

Methods. We developed a 3D magnetohydrodynamic model with the PLUTO code that describes the system HD 189733A , including the central host star and its hot Jupiter along with their respective winds. The effects of gravity and the magnetic fields of both the star and the planet are taken into account.

Results. Our analysis reveals that, in the cases examined in this study, the accretion scenario is only viable when the stellar magnetic field strength is at 5 G and the planetary magnetic field strength is at 1 G. In this scenario, the Rayleigh-Taylor instabilities lead to the formation of an accretion column that connects the star to the planet. Once formed the accretion column remains stable for the duration of the simulation. The accretion column produces an accretion rate of the order of 10¹² g s⁻¹ and shows an average density of about 10⁷ cm⁻³. In the other case explored, the accretion column does not form because the Rayleigh-Taylor instability is suppressed by the stronger magnetic field intensities assumed for both the star and the planet. We synthesized the emission resulting from the shocked planetary wind and found that the total X-ray emission ranges from 5 × 10²³ to 10²⁴ erg s⁻¹.

Conclusions. In the case of accretion, the emission originating from the hotspot cannot be distinguished from the coronal activity. Also, the interaction between the planetary and stellar winds cannot be responsible for the X-ray emission, as the total emission produced is about four orders of magnitude lower than the average X-ray luminosity of the star.