skip to main content


Title: Sub-stellar engulfment by a main-sequence star: Where is the lithium?
Context. Since the discovery of exoplanetary systems, questions have been raised as to the sub-stellar companions that can survive encounters with their host star, and how this interaction may affect the internal structure and evolution of the hosting star, and particularly its surface chemical composition. Aims. We study whether the engulfment of a brown dwarf (BD) by a solar-like main-sequence (MS) star can significantly alter the structure of the star and the Li content on its surface. Methods. We performed 3D smoothed particle hydrodynamics simulations of the engulfment of a BD with masses 0.01 and 0.019 M ⊙ , on an MS star of 1 M ⊙ and solar composition, in three different scenarios: a head-on collision, a grazing collision with an impact parameter η  = 0.5  R ⊙ , and a merger. We studied the dynamics of the interaction in detail, and the relevance of the type of interaction and the mass of the BD on the final fate of the sub-stellar object and the host star in terms of mass loss of the system, angular momentum transfer, and changes in the Li abundance on the surface of the host star. Results. In all the studied scenarios, most of the BD mass is diluted in the denser region of the MS star. Only in the merger scenario a significant fraction (∼40%) of the BD material would remain in the outer layers. We find a clear increase in the surface rotational velocity of the host star after the interaction, ranging between 25 km s −1 (grazing collision) to 50 km s −1 (merger). We also find a significant mass loss from the system (in the range 10 −4  − 10 −3   M ⊙ ) due to the engulfment, which in the case of the merger may form a circumstellar disk-like structure. Assuming that neither the depth of the convective envelope of the host star nor its mass content are modified during the interaction, a small change in the surface Li abundance in the head-on and grazing collisions is found. However, in the merger we find large Li enhancements, by factors of 20 − 30, depending on the BD mass. Some of these features could be detected observationally in the host star, provided they remained for a long enough time. Conclusions. In our 3D simulations, a sizable fraction of the BD survives long enough to be mixed with the inner core of the MS star. This is at odds with previous suggestions based on 1D simulations. In some cases the final surface rotational velocity is very high, coupled with enough mass loss that may form a circumstellar disk. Merger scenarios tend to dilute considerably more BD material on the surface of the MS star, which could be detected as a Li-enhancement. The dynamic of the simulated scenarios suggests the development of asymmetries in the structure of the host star that can only be tackled with 3D codes, including the long-term evolution of the system.  more » « less
Award ID(s):
1927130
NSF-PAR ID:
10464895
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
670
ISSN:
0004-6361
Page Range / eLocation ID:
A155
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Aims: We present a detailed long-term study of the single M6 III giant RZ Ari to obtain direct and simultaneous measurements of the magnetic field, activity indicators, and radial velocity in order to infer the origin of its activity. We study its magnetic activity in the context of stellar evolution, and for this purpose, we also refined its evolutionary status and Li abundance. In general, for the M giants, little is known about the properties of the magnetic activity and its causes. RZ Ari possess the strongest surface magnetic field of the known Zeeman-detected M giants and is bright enough to allow a deep study of its surface magnetic structure. The results are expected to shed light on the activity mechanism in these stars.

    Methods: We used the spectropolarimeter Narval at the Télescope Bernard Lyot (Observatoire du Pic du Midi, France) to obtain a series of Stokes I and V profiles for RZ Ari. Using the least-squares deconvolution technique, we were able to detect the Zeeman signature of the magnetic field. We measured its longitudinal component by means of the averaged Stokes I and V profiles. In addition, we also applied Zeeman-Doppler imaging (ZDI) to search for the rotation period of the star, and we constructed a tentative magnetic map. It is the first magnetic map for a star that evolved at the tip of red giant branch (RGB) or even on the asymptotic giant branch (AGB). The spectra also allowed us to monitor chromospheric emission lines, which are well-known indicators of stellar magnetic activity. From the observations obtained between September 2010 and August 2019, we studied the variability of the magnetic field of RZ Ari. We also redetermined the initial mass and evolutionary status of this star based on current stellar evolutionary tracks and on the angular diameter measured from CHARA interferometry. Results: Our results point to an initial mass of 1.5Mso that this giant is more likely an early-AGB star, but a lotaction at the tip of the RGB is not completely excluded. With a v sin i of 6.0 ±0.5 km s−1, the upper limit for the rotation period is found to be 909 days. On the basis of our dataset and AAVSO photometric data, we determined periods longer than 1100 days for the magnetic field and photometric variability, and 704 days for the spectral line activity indicators. The rotation period determined on the basis of the Stokes V profiles variability is 530 days. A similar period of 544 days is also found for the photometric data. When we take this rotation period and the convective turnover time into account, an effective action of an α-ω type dynamo seems to be unlikely, but other types of dynamo could be operating there. The star appears to lie outside the two magnetic strips on the giant branches, where the α-ω-type dynamo is expected to operate effectively, and it also has a much higher lithium content than the evolutionary model predicts. These facts suggest that a planet engulfment could speed up its rotation and trigger dynamo-driven magnetic activity. On the other hand, the period of more than 1100 days cannot be explained by rotational modulation and could be explained by the lifetime of large convective structures. The absence of linear polarization at the time the magnetic field was detected, however, suggests that a local dynamo probably does not contribute significantly to the magnetic field, at least for that time interval. 
    more » « less
  2. Aims. We present a detailed visible and near-infrared spectro-interferometric analysis of the Be-shell star o Aquarii from quasi-contemporaneous CHARA/VEGA and VLTI/AMBER observations. Methods. We analyzed spectro-interferometric data in the H α (VEGA) and Br γ (AMBER) lines using models of increasing complexity: simple geometric models, kinematic models, and radiative transfer models computed with the 3D non-LTE code HDUST. Results. We measured the stellar radius of o Aquarii in the visible with a precision of 8%: 4.0 ± 0.3 R ⊙ . We constrained the circumstellar disk geometry and kinematics using a kinematic model and a MCMC fitting procedure. The emitting disk sizes in the H α and Br γ lines were found to be similar, at ~10–12 stellar diameters, which is uncommon since most results for Be stars indicate a larger extension in H α than in Br γ . We found that the inclination angle i derived from H α is significantly lower (~15°) than the one derived from Br γ : i ~ 61.2° and 75.9°, respectively. While the two lines originate from a similar region of the disk, the disk kinematics were found to be near to the Keplerian rotation (i.e., β = −0.5) in Br γ ( β ~ −0.43), but not in H α ( β ~ −0.30). After analyzing all our data using a grid of HDUST models (BeAtlas), we found a common physical description for the circumstellar disk in both lines: a base disk surface density Σ 0 = 0.12 g cm −2 and a radial density law exponent m = 3.0. The same kind of discrepancy, as with the kinematic model, is found in the determination of i using the BeAtlas grid. The stellar rotational rate was found to be very close (~96%) to the critical value. Despite being derived purely from the fit to interferometric data, our best-fit HDUST model provides a very reasonable match to non-interferometric observables of o Aquarii: the observed spectral energy distribution, H α and Br γ line profiles, and polarimetric quantities. Finally, our analysis of multi-epoch H α profiles and imaging polarimetry indicates that the disk structure has been (globally) stable for at least 20 yr. Conclusions. Looking at the visible continuum and Br γ emission line only, o Aquarii fits in the global scheme of Be stars and their circumstellar disk: a (nearly) Keplerian rotating disk well described by the viscous decretion disk (VDD) model. However, the data in the H α line shows a substantially different picture that cannot fully be understood using the current generation of physical models of Be star disks. The Be star o Aquarii presents a stable disk (close to the steady-state), but, as in previous analyses, the measured m is lower than the standard value in the VDD model for the steady-state regime ( m = 3.5). This suggests that some assumptions of this model should be reconsidered. Also, such long-term disk stability could be understood in terms of the high rotational rate that we measured for this star, the rate being a main source for the mass injection in the disk. Our results on the stellar rotation and disk stability are consistent with results in the literature showing that late-type Be stars are more likely to be fast rotators and have stable disks. 
    more » « less
  3. ABSTRACT

    The detection of the binary black hole merger GW190521, with primary black hole mass $85^{+21}_{-14} {\rm M}_{\odot }$, proved the existence of black holes in the theoretically predicted pair-instability gap ($\sim 60-120 \, {\rm M}_{\odot }$) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon–oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code StarSmasher, finding that a head-on collision can remove up to 12 per cent of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product.

     
    more » « less
  4. The detection of the binary black hole merger GW190521, with primary mass 85+21−14 M⊙ , proved the existence of black holes in the theoretically predicted pair-instability gap ( ∼60−120M⊙ ) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon-oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code StarSmasher, finding that a head-on collision can remove up to 12% of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the carbon-oxygen core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product. 
    more » « less
  5. Abstract

    Progenitor models for the “luminous” subclass of Fast Blue Optical Transients (LFBOTs; prototype: AT2018cow) are challenged to simultaneously explain all of their observed properties: fast optical rise times of days or less; peak luminosities ≳1044erg s−1; low yields ≲0.1Mof56Ni; aspherical ejecta with a wide velocity range (≲3000 km s−1to ≳0.1–0.5cwith increasing polar latitude); presence of hydrogen-depleted-but-not-free dense circumstellar material (CSM) on radial scales from ∼1014cm to ∼3 × 1016cm; embedded variable source of non-thermal X-ray/γ-rays, suggestive of a compact object. We show that all of these properties are consistent with the tidal disruption and hyper-accretion of a Wolf-Rayet (WR) star by a black hole or neutron star binary companion. In contrast with related previous models, the merger occurs with a long delay (≳100 yr) following the common envelope (CE) event responsible for birthing the binary, as a result of gradual angular momentum loss to a relic circumbinary disk. Disk-wind outflows from the merger-generated accretion flow generate the56Ni-poor aspherical ejecta with the requisite velocity range. The optical light curve is powered primarily by reprocessing X-rays from the inner accretion flow/jet, though CSM shock interaction also contributes. Primary CSM sources include WR mass loss from the earliest stages of the merger (≲1014cm) and the relic CE disk and its photoevaporation-driven wind (≳1016cm). Longer delayed mergers may instead give rise to supernovae Type Ibn/Icn (depending on the WR evolutionary state), connecting these transient classes with LFBOTs.

     
    more » « less