Companion Star Research Articles

Planetary nebulae seen with TESS: New and revisited short-period binary central star candidates from Cycles 1 to 4

Context. High-precision and high-cadence photometric surveys such as Kepler or TESS are making huge progress not only in the detection of new extrasolar planets but also in the study of a great number of variable stars. This is the case for central stars of planetary nebulae (PNe), which have similarly benefited from the capabilities of these missions, increasing the number of known binary central stars and helping us to constrain the relationship between binarity and the complex morphologies of their host PNe. Aims. In this paper, we analyse the TESS light curves of a large sample of central stars of PNe with the aim of detecting signs of variability that may hint at the presence of short-period binary nuclei. This will have important implications in understanding PN formation evolution as well as the common envelope phase. Methods. We analysed 62 central stars of true, likely, or possible PNe and modelled the detected variability through an MCMC approach accounting for three effects: reflection, ellipsoidal modulations due to tidal forces, and the so-called Doppler beaming. Among the 62 central stars, only 38 are amenable for this study. The remaining 24 show large contamination from nearby sources preventing an optimal analysis. Also, eight targets are already known binary central stars, which we revisit here with the new high precision of the TESS data. Results. In addition to recovering the eight already known binaries in our sample, we find that 18 further central stars show clear signs of periodic variability in the TESS data, probably resulting from different physical effects compatible with the binary scenario. We propose them as new candidate binary central stars. We also discuss the origin of the detected variability in each particular case by using the TESS_localize algorithm. Finally, 12 targets show no or only weak evidence of variability at the sensitivity of TESS. Conclusions. Our study demonstrates the power of space-based photometric surveys in searching for close binary companions of central stars of PNe. Although our detections can only be catalogued as candidate binaries, we find a large percentage of possible stellar pairs associated with PNe, supporting the hypothesis that binarity plays a key role in shaping these celestial structures.

Short-term variability of the transitional pulsar candidate CXOU J110926.4−650224 from X-rays to infrared

CXOU J110926.4−650224 is a candidate transitional millisecond pulsar (tMSP) with X-ray and radio emission properties reminiscent of those observed in confirmed tMSPs in their X-ray ‘subluminous’ disc state. We present the results of observing campaigns that, for the first time, characterise the optical and near-infrared variability of this source and establish a connection with the mode-switching phenomenon observed in X-rays. The optical emission exhibited flickering activity, frequent dipping episodes where it appeared redder, and a multi-peaked flare where it was bluer. The variability pattern was strongly correlated with that of the X-ray emission. Each dip matched an X-ray low-mode episode, indicating that a significant portion of the optical emission originates from nearly the same region as the X-ray emission. The near-infrared emission also displayed remarkable variability, including a dip of 20 min in length during which it nearly vanished. Time-resolved optical spectroscopic observations reveal significant changes in the properties of emission lines from the disc and help infer the spectral type of the companion star to be between K0 and K5. We compare the properties of CXOU J110926.4−650224 with those of other tMSPs in the X-ray subluminous disc state and discuss our findings within the context of a recently proposed scenario that explains the phenomenology exhibited by the prototypical tMSP PSR J1023+0038.

VELOcities of CEpheids (VELOCE)

Classical Cepheids provide valuable insights into the evolution of stellar multiplicity among intermediate-mass stars. Here, we present a systematic investigation of single-lined spectroscopic binaries (SB1s) based on high-precision velocities measured by the VELOcities of CEpheids (VELOCE) project. We detected 76 (29%) SB1 systems among the 258 Milky Way Cepheids in the first VELOCE data release, 32 (43%) of which were not previously known to be SB1 systems. We determined 30 precise and three tentative orbital solutions, 18 (53%) of which are reported for the first time. This large set of Cepheid orbits provides a detailed view of the eccentricity e and orbital period Porb distribution among evolved intermediate-mass stars, ranging from e ∈ [0.0, 0.8] and Porb ∈ [240, 9000] d. The orbital motion on timescales exceeding the 11 yr VELOCE baseline was investigated using a template-fitting technique applied to literature data. Particularly interesting objects include (a) R Cru, the Cepheid with the shortest orbital period in the Milky Way (∼238 d); (b) ASAS J103158−5814.7, a short-period overtone Cepheid exhibiting time-dependent pulsation amplitudes as well as orbital motion; and (c) 17 triple systems with outer visual companions, among other interesting objects. Most VELOCE Cepheids (21/23) that exhibit evidence of a companion based on a Gaia proper motion anomaly are also spectroscopic binaries, whereas the remaining do not exhibit significant (> 3σ) orbital radial velocity variations. Gaia quality flags, notably the renormalized unit weight error (RUWE), do not allow Cepheid binaries to be identified reliably although statistically the average RUWE of SB1 Cepheids is slightly higher than that of non-SB1 Cepheids. A comparison with Gaia photometric amplitudes in G-, Bp, and Rp also does not allow one to identify spectroscopic binaries among the full VELOCE sample, indicating that the photometric amplitudes in this wavelength range are not sufficiently informative of companion stars.

A Precessing Stellar Disk Model for Superorbital Modulations of the Gamma-Ray Binary LS I+61° 303

Gamma-ray binary LS I+61° 303 consists of a neutron star orbiting around a Be star with a period of P orb ≃ 26.5 days. Apart from orbital modulations, the binary shows long-term flux variations with a superorbital period of Psup≃4.6yr as seen in nearly all wavelengths. The origin of this superorbital modulation is still not well understood. Under the pulsar wind–stellar outflow interaction scenario, we propose that the superorbital modulations of LS I+61° 303 could be caused by the precession of the Be disk. Assuming X-rays arise from synchrotron radiation of the intrabinary shock, we develop an analytical model to calculate expected flux modulations over the orbital and superorbital phases. The asymmetric two-peak profiles in orbital light curves and sinusoidal-like long-term modulations are reproduced under the precessing disk scenario. The observed orbital phase drifting of the X-ray peak and our fitting of long-term X-ray data indicate that the neutron star is likely orbiting around the star with a small eccentricity and periastron phase around Φp ∼ 0.6. We compare the Corbet diagrams of LS I+61° 303 with other Be/X-ray binaries, and the linear correlation in the Psup–Porb diagram suggests that the precession of the Be disk in LS I+61° 303 is induced by the tidal torque of its neutron star companion.

X-Ray Flashes on Helium Novae

A helium nova occurs on a white dwarf (WD) accreting hydrogen-deficient matter from a helium star companion. When the mass of a helium envelope on the WD reaches a critical value, unstable helium burning ignites to trigger a nova outburst. A bright soft X-ray phase appears in an early outbursting phase of a helium nova before it optically rises toward maximum. Such an X-ray bright phase is called the X-ray flash. We present theoretical light curves of X-ray flashes for 1.0, 1.2, and 1.35 M ☉ helium novae with mass accretion rates of (1.6–7.5) × 10−7 M ☉ yr−1. Long durations of the X-ray flashes (100 days–10 yr) and high X-ray luminosities (∼1038 erg s−1) indicate that X-ray flashes are detectable as a new type of X-ray transient or persistent X-ray sources. An X-ray flash is a precursor of optical brightening, so that the detection of X-ray flashes on helium novae enables us to plan arranged observation for optical premaximum phases that have been one of the frontiers of the study of novae. We found a candidate object of helium-burning X-ray flash from the literature on extragalactic X-ray surveys. This X-ray transient source is consistent with our X-ray flash model of a 1.35 M ⊙ WD.

Spectroscopic Observations of the GALEX Nearby Young Star Survey Sample. I. Nearby Moving Group Candidates

The GALEX Nearby Young Star Search (GALNYSS) yielded the identification of more than 2000 late-type stars that, based on their ultraviolet and infrared colors and pre-Gaia proper motions, are potentially of age<200 Myr and lie within ~120 pc of Earth. We present the results of a campaign of medium- and high-resolution optical spectroscopy of 471 GALNYSS stars aimed at confirming their youth and their potential membership in nearby young stellar moving groups. We present radial velocity (RV), Li absorption, and H-alpha emission measurements for these spectroscopically observed GALNYSS stars, and assess their Li absorption and optical emission-line properties and infrared excesses. Our RV measurements are combined with literature and Gaia DR3 RV measurements and Gaia DR3 astrometry and photometry to obtain the spatiokinematics and color-magnitude positions of GALNYSS stars. We use these results to assess membership in the TW Hya, Tuc-Hor, Carina, Columba, and Argus Associations and the β Pic and AB Dor moving groups. We have identified 132 stars as candidate members of these seven groups; roughly half of these candidates are newly identified on the basis of data presented here. At least one-third of the 132 candidates are spectroscopic and/or photometric binaries and/or have comoving (visual) binary companions in Gaia DR3. The contingent of young, low-mass stars in the solar vicinity we identify here should provide excellent subjects for future direct imaging exoplanet surveys and studies of the early evolution of low-mass stars and their planetary progeny.

BSN: The First Light Curve Analysis of the Total Eclipse Binary System EL Tuc

We conducted the first light curve study of the binary star EL Tuc within the Binary Systems of South and North project’s framework. The photometric observations were made using standard multiband BVR c I c filters at an observatory in Argentina. We presented a new ephemeris for EL Tuc and a linear fit to the O – C diagram, utilizing our extracted times of minima and additional literature. We employed the PHysics Of Eclipsing BinariEs Python code and the Markov chain Monte Carlo approach for the system’s light curve analysis. The target system’s light curve solution required a cold starspot on the hotter component. We conclude that EL Tuc is a total contact binary system with a low mass ratio of q = 0.172 ± 0.002, an orbital inclination of i = 83.°74 ± 0.°40, and a fillout factor of f = 53.7% ± 1.6%. We used the P–a relationship and the Gaia Data Release 3 parallax method to determine the absolute parameters of EL Tuc to compare the precision of our results. This system was classified as W-type based on the mass and effective temperature of the companion stars. The positions of the systems were depicted on the M–L, M–R, T–M, and q–L ratio diagrams. The relationship between the spectroscopic and photometric mass ratios of binaries was discussed.

Merger Precursor: Year-long Transients Preceding Mergers of Low-mass Stripped Stars with Compact Objects

Binary mass transfer can occur at high rates due to rapid expansion of the donor’s envelope. In the case where mass transfer is unstable, the binary can rapidly shrink its orbit and lead to a merger. In this work we consider the appearance of the system preceding merger, specifically for the case of a low-mass ( ≈2.5- 3M⊙) helium star with a neutron star (NS) companion. Modeling the mass transfer history as well as the wind launched by super-Eddington accretion onto the NS, we find that such systems can power slowly rising transients with timescales as long as years, and luminosities of ∼1040- 1041 erg s −1 from optical to UV. The final explosion following the merger (or core-collapse of the helium star in some cases) leads to an interaction-powered transient with properties resembling Type Ibn supernovae (SNe), possibly with a bright early peak powered by shock cooling emission for merger-powered explosions. We apply our model to the Type Ibn SN 2023fyq, that displayed a long-term precursor activity from years before the terminal explosion.

Constraining Wind-driven Accretion onto Gaia BH3 with Chandra

Gaia BH3 is the most massive known stellar-origin black hole in the Milky Way, with a mass M • ≈ 33 M ⊙. Detected from Gaia’s astrometry, this black hole is in the mass range of those observed via gravitational waves, whose nature is still highly debated. Hosted in a binary system with a companion giant star that is too far away for Roche-lobe mass transfer, this black hole could nonetheless accrete at low levels due to wind-driven mass loss from its companion star, thus accreting in advection-dominated accretion flow, or ADAF, mode. Using stellar wind models, we constrain its Eddington ratio in the range 10−9 < f Edd < 10−7, corresponding to radiative efficiencies 5 × 10−5 < ϵ < 10−3, compatible with radiatively inefficient accretion modes. Chandra ACIS-S observed this object and obtained the most sensitive upper bound of its [2–10] keV flux: F X < 3.25 × 10−15 erg s−1 cm−2 at 90% confidence level corresponding to L [2–10] < 2.10 × 1029 erg s−1. Using ADAF emission models, we constrained its accretion rate to f Edd < 4.91 × 10−7 at the apastron, in agreement with our theoretical estimate. At the periastron, we expect fluxes ∼50 times larger. Because of the inferred low rates, accretion did not significantly contribute to black hole growth over the system’s lifetime. Detecting the electromagnetic emission from Gaia BH3 will be fundamental to informing stellar wind and accretion disk models.

The Evolution of Massive Binary Stars

Massive stars play a major role in the evolution of their host galaxies and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and interact with their companion stars during their lifetimes. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also responsible for the production of interesting and exotic interaction products. ▪Extensive observation campaigns with well-understood detection sensitivities have enabled the conversion of observed properties into intrinsic characteristics, facilitating a direct comparison to theory.▪Studies of large samples of massive stars in our Galaxy and the Magellanic Clouds have unveiled new types of interaction products, providing critical constraints on the mass transfer phase and the formation of compact objects.▪The direct detection of gravitational waves has revolutionized the study of stellar mass compact objects, providing a new window to study massive star evolution. Their formation processes are, however, still unclear. The known sample of compact object mergers will increase by orders of magnitude in the coming decade, which is vastly outgrowing the number of stellar-mass compact objects detected through electromagnetic radiation.

The true nature of HE 0057-5959, the most metal-poor, Li-rich star

The Li-rich stars are a class of rare objects with a surface lithium abundance, A(Li), that exceeds that of other stars in the same evolutionary stage. The origin of these stars is still debated, and valuable routes to look at include the Cameron-Fowler mechanism, a mass-transfer process in a binary system, or the engulfment of rocky planets or brown dwarfs. Metal-poor ([Fe/H]&lt;−1 dex) stars are only a small fraction of the entire population of Li-rich stars. We observed the metal-poor ([Fe/H]=−3.95±0.11 dex) giant star HE 0057–5959 with MIKE at the Magellan Telescope, deriving A(Li)NLTE=+2.09±0.07 dex. Such an Li abundance is significantly higher, by about 1 dex, than that of other stars at the same evolutionary stage. A previous analysis of the same target suggested that its high A(Li) reflects an ongoing first-dredge-up process. We revised the nature of HE 0057-5959 by comparing its stellar parameters and A(Li) with appropriate stellar evolution models describing Li depletion due to the deepening of the convective envelope. This comparison rules out that HE 0057-5959 is caught during its first dredge-up, the latter having already ended according to the parameters of this star. Its A(Li), remarkably higher than the typical lithium plateau drawn by similar giant stars, demonstrates that HE 0057-5959 joins the class of the rare metal-poor, Li-rich stars. HE 0057-5959 is the most metal-poor, Li-rich star discovered so far. We considered different scenarios to explain this star also comparing it with the other metal-poor, Li-rich stars. No internal mixing able to activate the Cameron-Fowler mechanism is known for metal-poor stars at this evolutionary stage. The engulfment of planets is also disfavoured because such metal-poor stars should not host planets. Finally, HE 0057-5959 is one of the most Na-rich among the Li-rich stars, and we found that a strong excess of Na abundance is common to all three Li-rich stars with [Fe/H]&lt;–3 dex. This finding could support the scenario of mass transfer from a massive companion star (able to simultaneously produce large amounts of both elements) in a binary system, even if we found no evidence of radial velocity variations.

An Alternative Model for the Orbital Decay of M82 X-2: The Anomalous Magnetic Braking of a Bp Star

Recently, the first pulsating ultraluminous X-ray source M82 X-2 was reported to be experiencing a rapid orbital decay at a rate of Ṗ=−(5.69±0.24)×10−8ss−1 based on 7 yr NuSTAR data. To account for the observed orbital-period derivative, it requires a mass transfer rate of ∼200Ṁedd ( Ṁedd is the Eddington accretion rate) from the donor star to the accreting neutron star. However, other potential models cannot be completely excluded. In this work, we propose an anomalous magnetic braking (AMB) model to interpret the detected orbital decay of M82 X-2. If the donor star is an Ap/Bp star with an anomalously strong magnetic field, the magnetic coupling between strong surface magnetic field and irradiation-driven wind from the surface of the donor star could cause an efficient angular-momentum loss, driving a rapid orbital decay observed in M82 X-2. The AMB mechanism of an Ap/Bp star with a mass of 5.0–15.0 M ⊙ and a surface magnetic field of 3000–4500 G could produce the observed Ṗ of M82 X-2. We also discuss the possibility of other alternative models including the companion star expansion and a surrounding circumbinary disk.

Astrometry and Precise Radial Velocities Yield a Complete Orbital Solution for the Nearby Eccentric Brown Dwarf LHS 1610 b

The LHS 1610 system consists of a nearby (d = 9.7 pc) M5 dwarf hosting a candidate brown dwarf companion in a 10.6 days, eccentric (e ∼ 0.37) orbit. We confirm this brown dwarf designation and estimate its mass ( 49.5−3.5+4.3 M Jup) and inclination (114.5° −10.0+7.4 ) by combining discovery radial velocities (RVs) from the Tillinghast Reflector Echelle Spectrograph and new RVs from the Habitable-zone Planet Finder with the available Gaia astrometric two-body solution. We highlight a discrepancy between the measurement of the eccentricity from the Gaia two-body solution (e = 0.52 ± 0.03) and the RV-only solution (e = 0.3702 ± 0.0003). We discuss possible reasons for this discrepancy, which can be further probed when the Gaia astrometric time series become available as part of Gaia Data Release 4. As a nearby mid-M star hosting a massive short-period companion with a well-characterized orbit, LHS 1610 b is a promising target to look for evidence of sub-Alfvénic interactions and/or auroral emission at optical and radio wavelengths. LHS 1610 has a flare rate (0.28 ± 0.07 flares per day) on the higher end for its rotation period (84 ± 8 days), similar to other mid-M dwarf systems such as Proxima Cen and YZ Ceti that have recent radio detections compatible with star–planet interactions. While available Transiting Exoplanet Survey Satellite photometry is insufficient to determine an orbital phase dependence of the flares, our complete orbital characterization of this system makes it attractive to probe star–companion interactions with additional photometric and radio observations.

The Completeness of Accreting Neutron Star Binary Candidates from the Chinese Space Station Telescope

Abstract A neutron star (NS) has many extreme physical conditions, and one may obtain some important information about an NS via accreting neutron star binary (ANSB) systems. The upcoming Chinese Space Station Telescope (CSST) provides an opportunity to search for a large sample of ANSB candidates. Our goal is to check the completeness of the potential ANSB samples from CSST data. In this paper, we generate some ANSBs and normal binaries under the CSST photometric system by binary evolution and binary population synthesis method and use a machine learning method to train a classification model. Although the Precision (94.56%) of our machine learning model is as high as before study, the Recall is only about 63.29%. The Precision/Recall is mainly determined by the mass transfer rate between the NSs and their companions. In addition, we also find that the completeness of ANSB samples from CSST photometric data by the machine learning method also depends on the companion mass and the age of the system. ANSB candidates with a low initial mass companion star (0.1 M ⊙ to 1 M ⊙) have a relatively high Precision (94.94%) and high Recall (86.32%), whereas ANSB candidates with a higher initial mass companion star (1.1 M ⊙ to 3 M ⊙) have similar Precision (93.88%) and quite low Recall (42.67%). Our results indicate that although the machine learning method may obtain a relatively pure sample of ANSBs, a completeness correction is necessary for one to obtain a complete sample.

Radiative plasma simulations of black hole accretion flow coronae in the hard and soft states

Stellar-mass black holes in x-ray binary systems are powered by mass transfer from a companion star. The accreted gas forms an accretion disk around the black hole and emits x-ray radiation in two distinct modes: hard and soft state. The origin of the states is unknown. We perform radiative plasma simulations of the electron-positron-photon corona around the inner accretion flow. Our simulations extend previous efforts by self-consistently including all the prevalent quantum electrodynamic processes. We demonstrate that when the plasma is turbulent, it naturally generates the observed hard-state emission. In addition, we show that when soft x-ray photons irradiate the system—mimicking radiation from an accretion disk—the turbulent plasma transitions into a new equilibrium state that generates the observed soft-state emission. Our findings demonstrate that turbulent motions of magnetized plasma can power black-hole accretion flow coronae and that quantum electrodynamic processes control the underlying state of the plasma.

A Pilot Search for Gravitational Self-Lensing Binaries with the Zwicky Transient Facility

Binary systems containing a compact object may exhibit periodic brightening episodes due to gravitational lensing as the compact object transits the companion star. Such “self-lensing” signatures have been detected before for white dwarf binaries. We attempt to use these signatures to identify detached stellar-mass neutron star and black hole binaries using data from the Zwicky Transient Facility (ZTF). We present a systematic search for self-lensing signals in Galactic binaries from a subset of high-cadence ZTF data taken in 2018. We identify 12 plausible candidates from the search, although because each candidate is observed to only brighten once, other origins such as stellar flares are more likely. We discuss prospects for more comprehensive future searches of the ZTF data.

An Alternative Channel to Black Hole Low-mass X-Ray Binaries: Dynamical Friction of Dark Matter?

Both the anomalous magnetic braking of Ap/Bp stars and the surrounding circumbinary disk models can account for the formation of black hole (BH) low-mass X-ray binaries (LMXBs), while the simulated effective temperatures of the donor stars are significantly higher than the observed values. Therefore, the formation of BH LMXBs is still not completely understood. In this work, we diagnose whether the dynamical friction between dark matter and the companion stars can drive BH binaries to evolve toward the observed BH LMXBs and alleviate the effective temperature problem. Assuming that there exists a density spike of dark matter around BH, the dynamical friction can produce an efficient angular momentum loss, driving BH binaries with an intermediate-mass companion star to evolve into BH LMXBs for a spike index higher than γ = 1.58. Our detailed stellar evolution models show that the calculated effective temperatures can match the observed value of most BH LMXBs for a spike index range of γ = 1.7–2.1. However, the simulated mass-transfer rates when γ = 2.0 and 2.1 are too high to be consistent with the observed properties showing that BH LMXBs appear as soft X-ray transients. Therefore, the dynamical friction of dark matter can only alleviate the effective temperature problem of those BH LMXBs with a relatively short orbital period.

To Be or Not To Be: The Role of Rotation in Modeling Galactic Be X-Ray Binaries

Be X-ray binaries (Be-XRBs) are one of the largest subclasses of high-mass X-ray binaries, comprised of a rapidly rotating Be star and neutron star companion in an eccentric orbit, intermittently accreting material from a decretion disk around the donor. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of supernovae, common envelope, and mass transfer (MT). Using the state-of-the-art BPS code, POSYDON, which relies on precomputed grids of detailed binary stellar evolution models, we investigate the Galactic Be-XRB population. POSYDON incorporates stellar rotation self-consistently during MT phases, enabling detailed examination of the rotational distribution of Be stars in multiple phases of evolution. Our fiducial BPS and Be-XRB model aligns well with the orbital properties of Galactic Be-XRBs, emphasizing the role of rotational constraints. Our modeling reveals a rapidly rotating population (ω/ω crit ≳ 0.3) of Be-XRB-like systems with a strong peak at intermediate rotation rates (ω/ω crit ≃ 0.6) in close alignment with observations. All Be-XRBs undergo an MT phase before the first compact object forms, with over half experiencing a second MT phase from a stripped helium companion (Case BB). Computing rotationally limited MT efficiencies and applying them to our population, we derive a physically motivated MT efficiency distribution, finding that most Be-XRBs have undergone highly nonconservative MT ( β¯rot≃0.15 ). Our study underscores the importance of detailed angular momentum modeling during MT in interpreting Be-XRB populations, emphasizing this population as a key probe for the stability and efficiency of MT in interacting binaries.

Whispering in the dark: Faint X-ray emission from black holes with OB star companions

Recently, astrometric and spectroscopic surveys of OB stars revealed a few stellar-mass black holes (BHs) with orbital periods of as low as 10 days. Contrary to wind-fed BH high-mass X-ray binaries (HMXBs), no X-ray counterpart was detected, probably because of the absence of a radiatively efficient accretion disc around the BH. Nevertheless, dissipative processes in the hot, dilute, and strongly magnetised plasma around the BH (so-called BH corona) can still lead to non-thermal X-ray emission (e.g. synchrotron). We determine the X-ray luminosity distribution from BH+OB star binaries up to orbital periods of a few thousand days. odot $ and orbital periods of 1-3000\,d. We computed the X-ray luminosity for a broad range of radiative efficiencies that depend on the mass accretion rate and flow geometry. odot $, which aligns well with our predictions and may be interesting sources for follow-up observations. The predicted luminosities of the OB companions to these X-ray-emitting BHs are 10$^ odot These findings advocate for prolonged X-ray observations of the stellar-mass black hole candidates identified in the vicinity of OB stars. Such long exposures could reveal the underlying population of X-ray-faint BHs and provide constraints for the evolution from single to double degenerate binaries and identify the progenitors of gravitational wave mergers.

Bumpy Superluminous Supernovae Powered by a Magnetar–Star Binary Engine

Wolf–Rayet stars in close binary systems can be tidally spun up by their companions, potentially leaving behind fast-spinning, highly magnetized neutron stars, known as “magnetars,” after core collapse. These newborn magnetars can transfer rotational energy into heating and accelerating the ejecta, producing hydrogen-poor superluminous supernovae (SLSNe). In this Letter, we propose that the magnetar wind of the newborn magnetar could significantly evaporate its companion star, typically a main-sequence or helium star, if the binary system is not disrupted by the abrupt mass loss and supernova (SN) kick. The subsequent heating and acceleration of the evaporated star material along with the SN ejecta by the magnetar wind can produce a postpeak bump in the SLSN lightcurve. Our model can reproduce the primary peaks and postpeak bumps of four example observed multiband SLSN lightcurves, revealing that the mass of the evaporated material could be ∼0.4–0.6 M ⊙ if the material is hydrogen-rich. We propose that the magnetar could induce strongly enhanced evaporation from its companion star near the pericenter if the orbit of the post-SN binary is highly eccentric, ultimately generating multiple postpeak bumps in the SLSN lightcurves. This “magnetar–star binary engine” model may offer a possible explanation for the evolution of polarization, along with the origin and velocity broadening of late-time hydrogen or helium broad spectral features observed in some bumpy SLSNe. The diversity in the lightcurves and spectra of SLSNe may be attributed to the wide variety of companion stars and post-SN binary systems.