Stellar Evolution Code Research Articles

The first two-dimensional stellar structure and evolution models of rotating stars

Contact. Rotation is a key ingredient in the theory of stellar structure and evolution. Until now, stellar evolution codes operate in a one-dimensional framework for which the validity domain in regards to the rotation rate is not well understood. Aims. In this Letter, we present the first results of self-consistent stellar models in two spatial dimensions, which compute the time evolution of a star and its rotation rate along the main sequence (MS). We also present a comparison to observations. Methods. We make use of an extended version of the ESTER code, which solves the stellar structure of a rotating star in two dimensions with time evolution, including chemical evolution, and an implementation of rotational mixing. We computed evolution tracks for a 12 M⊙ model, once for an initial rotation rate equal to 15% of the critical frequency, and once for 50%. Results. We first show that our model initially rotating at 15% of the critical frequency is able to reproduce all the observations of the β Cephei star HD 192575, which was recently studied with asteroseismology. Beyond the classical surface parameters, such as effective temperature or luminosity, our model also reproduces the core mass along with the rotation rate of the core and envelope at the estimated age of the star. This particular model also shows that the meridional circulation has a negligible influence on the transport of chemical elements such as nitrogen, for which the abundance may be increased at the stellar surface. Furthermore, it shows that in the late MS, nuclear evolution is faster than the relaxation time needed to reach a steady state of stellar angular momentum distribution. Conclusions. We demonstrate that we have successfully taken a new step towards two-dimensional evolutionary modelling of rotating stars. This opens new perspectives on the understanding of the dynamics of fast rotating stars and on the way rotation impacts stellar evolution.

Angular momentum transport by magnetic fields in main-sequence stars with Gamma Doradus pulsators

Context. Asteroseismic studies show that cores of post-main-sequence stars rotate more slowly than theoretically predicted by stellar models with purely hydrodynamical transport processes. Recent studies of main-sequence stars, particularly Gamma Doradus (γ Dor) stars, have revealed the internal rotation rates for hundreds of stars, offering a counterpart on the main sequence for studies of angular momentum transport. Aims. We investigate whether such a disagreement between observed and predicted internal rotation rates is present in main-sequence stars by studying angular momentum transport in γ Dor stars. Furthermore, we test whether models of rotating stars with internal magnetic fields can reproduce their rotational properties. Methods. We computed rotating models with the Geneva stellar evolution code taking into account meridional circulation and shear instability. We also computed models with internal magnetic fields using a general formalism for transport by the Tayler-Spruit dynamo. We then compared these models to observational constraints for γ Dor stars that we compiled from the literature, thus combining the core rotation rates, projected rotational velocities from spectroscopy, and constraints on their fundamental parameters. Results. We show that combining the different observational constraints available for γ Dor stars enable us to clearly distinguish the different scenarios for internal angular momentum transport. Stellar models with purely hydrodynamical processes are in disagreement with the data, whereas models with internal magnetic fields can reproduce both core and surface constraints simultaneously. Conclusions. Similarly to results obtained for subgiant and red giant stars, angular momentum transport in radiative regions of γ Dor stars is highly efficient, in good agreement with predictions of models with internal magnetic fields.

The Imprint of Convection on Type I X-Ray Bursts: Pauses in Photospheric Radius Expansion Lightcurves

Motivated by the recent observation by NICER of a type I X-ray burst from SAX J1808.4–3658 with a distinct “pause” feature during its rise, we show that bursts which ignite in a helium layer underneath a hydrogen-rich shell naturally give rise to such pauses, as long as enough energy is produced to eject the outer layers of the envelope by super-Eddington winds. The length of the pause is determined by the extent of the convection generated after ignition, while the rate of change of luminosity following the pause is set by the hydrogen gradient left behind by convection. Using the MESA stellar evolution code, we simulate the accumulation, nuclear burning, and convective mixing prior to and throughout the ignition of the burst, followed by the hydrodynamic wind. We show that the results are sensitive to the treatment of convection adopted within the code. In particular, the efficiency of mixing at the H/He interface plays a key role in determining the shape of the lightcurve. The data from SAX J1808.4–3658 favor strong mixing scenarios. Multidimensional simulations will be needed to properly model the interaction between convection and nuclear burning during these bursts, which will then enable a new way to use X-ray burst lightcurves to study neutron star surfaces.

Convective dynamos of black widow companions

ABSTRACT Black widows and redbacks are binary millisecond pulsars with close low-mass companions that are irradiated and gradually ablated by the pulsar’s high-energy luminosity Lirr. These binaries evolve primarily through magnetic braking, which extracts orbital angular momentum and pushes the companion to overflow its Roche lobe. Here, we use the stellar evolution code mesa to examine how the irradiation modifies the companion’s structure. Strong Lirr inhibits convection to the extent that otherwise fully convective stars become almost fully radiative. By computing the convective velocities and assuming a dynamo mechanism, we find that the thin convective envelopes of such strongly irradiated companions ($L_{\rm irr}\gtrsim 3\, {\rm L}_\odot$) generate much weaker magnetic fields than previously thought – halting binary evolution. With our improved magnetic braking model, we explain most observed black widow and redback companions as remnants of main-sequence stars. We also apply our model (with Lirr) to evolved companions that overflow their Roche lobe close to the end of their main-sequence phase. The evolutionary tracks of such companions bifurcate, explaining the shortest period systems (which are potential gravitational wave sources) as well as the longest period ones (which are the progenitors of common pulsar–white dwarf binaries). The variety of black widow structures and evolutionary trajectories may be utilized to calibrate the dependence of magnetic braking on the size of the convective layer and on the existence of a radiative–convective boundary, with implications for single stars as well as other binaries, such as cataclysmic variables and AM Canum Venaticorum stars.

Understanding Binary Systems—a Comparison between COSMIC and MESA

We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our Universe.

Overview and Validation of the Asteroseismic Modeling Portal v2.0

The launch of NASA’s Kepler space telescope in 2009 revolutionized the quality and quantity of observational data available for asteroseismic analysis. While Kepler was able to detect solar-like oscillations in hundreds of main-sequence and subgiant stars, the Transiting Exoplanet Survey Satellite is now making similar observations for thousands of the brightest stars in the sky. The Asteroseismic Modeling Portal (AMP) is an automated and objective stellar model-fitting pipeline for asteroseismic data, which was originally developed to use models from the Aarhus Stellar Evolution Code. We briefly summarize an updated version of the AMP pipeline that uses Modules for Experiments in Stellar Astrophysics, and we present initial modeling results for the Sun and several solar analogs to validate the precision and accuracy of the inferred stellar properties.

Impact of central mixing scheme and nuclear reaction network on the extent of convective cores

Convective cores are the hydrogen reservoirs of main sequence stars that are more massive than around 1.2 solar masses. The characteristics of the cores have a strong impact on the evolution and structure of the star. However, such results rely on stellar evolution codes, in which simplistic assumptions are often made on the physics in the core. Indeed, mixing is commonly considered to be instantaneous and the most basic nuclear networks assume beryllium at its equilibrium abundance. Those assumptions lead to significant differences in the central composition of the elements for which the timescale to reach nuclear equilibrium is lower than the convective timescale. In this work, we show that those discrepancies impact the nuclear energy production and, therefore, the size of convective cores in models computed with overshoot. We find that cores computed with instantaneous mixing are up to 30% bigger than those computed with diffusive mixing. Similar differences are found when using basic nuclear networks. Additionally, we observed an extension of the duration of the main sequence due to those core size differences. We then investigated the impact of those structural differences on the seismic modeling of solar-like oscillators. Modeling two stars observed by Kepler, we find that the overshoot parameter of the best models computed with a basic nuclear network is significantly lower, compared to models computed with a full nuclear network. This work is a necessary step in improving the modeling of convective cores, which is key to determining accurate ages in the framework of future space missions such as Plato.

Modelling stellar evolution in mass-transferring binaries and gravitational-wave progenitors with metisse

ABSTRACT Massive binaries are vital sources of various transient processes, including gravitational-wave mergers. However, large uncertainties in the evolution of massive stars, both physical and numerical, present a major challenge to the understanding of their binary evolution. In this paper, we upgrade our interpolation-based stellar evolution code metisse to include the effects of mass changes, such as binary mass transfer or wind-driven mass loss, not already included within the input stellar tracks. metisse’s implementation of mass loss (applied to tracks without mass loss) shows excellent agreement with the sse fitting formulae and with detailed mesa tracks, except in cases where the mass transfer is too rapid for the star to maintain equilibrium. We use this updated version of metisse within the binary population synthesis code bse to demonstrate the impact of varying stellar evolution parameters, particularly core overshooting, on the evolution of a massive (25 and 15 M⊙) binary system with an orbital period of 1800 d. Depending on the input tracks, we find that the binary system can form a binary black hole or a black hole–neutron star system, with primary (secondary) remnant masses ranging between 4.47 (1.36) and 12.30 (10.89) M⊙, and orbital periods ranging from 6 d to the binary becoming unbound. Extending this analysis to a population of isolated binaries uniformly distributed in mass and orbital period, we show that the input stellar models play an important role in determining which regions of the binary parameter space can produce compact binary mergers, paving the way for predictions for current and future gravitational-wave observatories.

Long-duration Gamma-Ray Burst Progenitors and Magnetar Formation

Millisecond magnetars produced in the center of dying massive stars are one prominent model to power gamma-ray bursts (GRBs). However, their detailed nature remains a mystery. To explore the effects of the initial mass, rotation rate, wind mass loss, and metallicity on the GRB progenitors and the newborn magnetar properties, we evolve 227 of 10–30 M ⊙ single star models from the pre-main sequence to core collapse by using the stellar evolution code MESA. The presupernova properties, the compactness parameter, and the magnetar characteristics of models with different initial parameters are presented. The compactness parameter remains a nonmonotonic function of the initial mass and initial rotation rate when the effects of varying metallicity and the “Dutch” wind scale factor are taken into account. We find that the initial rotation rate and mass play the dominant roles in whether a star can evolve into a GRB progenitor. The minimum rotation rate necessary to generate a magnetar gradually reduces as the initial mass increases. The greater the initial metallicity and “Dutch” wind scale factor, the larger the minimum rotation rate required to produce a magnetar. In other words, massive stars with low metallicity are more likely to harbor magnetars. Furthermore, we present the estimated period, magnetic field strength, and masses of magnetars in all cases. The typical rotational energy of these millisecond magnetars is sufficient to power long-duration GRBs.

Modeling fast radio burst heating a main-sequence companion star in a close binary using MESA

The exact origin of fast radio bursts (FRBs) remains a mystery. The repeating fast radio burst source, FRB20200120E, was discovered in a globular cluster containing old stellar populations. Yang (2021) suggested that this FRB might be in close binaries with low-mass main-sequence (MS) stars. They analytically investigated the observational consequences caused by the heating of FRB radio radiation onto the low-mass MS companion star in a close binary, suggesting that the radio radiation emitted by FRB could make the MS companion star more luminous and detectable in future multi-wavelength follow-up observations for a Galactic FRB. We revisited the study of Yang (2021) by numerically modeling the detailed process of FRB heating onto an MS companion with 1D stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). Our results are consistent with the trends derived from the analytical model of Yang (2021), except that the typical re-emission luminosities of our main sequence (MS) models, caused by the heating from FRBs, are generally dimmer by about two orders of magnitude compared to his findings, and our models have a longer re-emission timescale. This may indicate that the searches of the optical transients caused by the radio radiation heating companion star are more likely to be successful within a distance of 0.3 Mpc.

Mass-stream trajectories with non-synchronously rotating donors

ABSTRACT Mass-transfer interactions in binary stars can lead to accretion disc formation, mass-loss from the system, and spin-up of the accretor. To determine the trajectory of the mass-transfer stream, and whether it directly impacts the accretor, or forms an accretion disc, requires numerical simulations. The mass-transfer stream is approximately ballistic, and analytic approximations based on such trajectories are used in many binary population synthesis codes as well as in detailed stellar evolution codes. We use binary population synthesis to explore the conditions under which mass transfer takes place. We then solve the reduced three-body equations to compute the trajectory of a particle in the stream for systems with varying system mass ratio, donor synchronicity, and initial stream velocity. Our results show that, on average, both more mass and more time are spent during mass transfer from a sub-synchronous donor than from a synchronous donor. Moreover, we find that at low initial stream velocity, the asynchronous rotation of the donor leads to self-accretion over a large range of mass ratios, especially for supersynchronous donors. The stream (self-)intersects in a narrow region of parameter space where it transitions between accreting on to the donor or the accretor. Increasing the initial stream velocity leads to larger areas of the parameter space where the stream accretes on to the accretor, but also more (self-)intersection. The radii of closest approach generally increase, but the range of specific angular momenta that these trajectories carry at the radius of closest approach gets broader. Our results are made publicly available.

Observational predictions for Thorne–Żytkow objects

ABSTRACT Thorne–Żytkow objects (TŻO) are potential end products of the merger of a neutron star with a non-degenerate star. In this work, we have computed the first grid of evolutionary models of TŻOs with the MESA stellar evolution code. With these models, we predict several observational properties of TŻOs, including their surface temperatures and luminosities, pulsation periods, and nucleosynthetic products. We expand the range of possible TŻO solutions to cover $3.45 \lesssim \rm {\log \left(T_{eff}/K\right)}\lesssim 3.65$ and $4.85 \lesssim \rm {\log \left(L/L_{\odot }\right)}\lesssim 5.5$. Due to the much higher densities our TŻOs reach compared to previous models, if TŻOs form we expect them to be stable over a larger mass range than previously predicted, without exhibiting a gap in their mass distribution. Using the GYRE stellar pulsation code we show that TŻOs should have fundamental pulsation periods of 1000–2000 d, and period ratios of ≈0.2–0.3. Models computed with a large 399 isotope fully coupled nuclear network show a nucleosynthetic signal that is different to previously predicted. We propose a new nucleosynthetic signal to determine a star’s status as a TŻO: the isotopologues $\mathrm{^{44}Ti} \rm {O}_2$ and $\mathrm{^{44}Ti} \rm {O}$, which will have a shift in their spectral features as compared to stable titanium-containing molecules. We find that in the local Universe (∼SMC metallicities and above) TŻOs show little heavy metal enrichment, potentially explaining the difficulty in finding TŻOs to-date.

Feasibility of structure inversions for gravity-mode pulsators

Context. Gravity-mode asteroseismology has significantly improved our understanding of mixing in intermediate mass stars. However, theoretical pulsation periods of stellar models remain in tension with observations, and it is often unclear how the models of these stars should be further improved. Inversions provide a path forward by directly probing the internal structure of these stars from their pulsation periods, quantifying which parts of the model are in need of improvement. This method has been used with success in the case of solar-like pulsators, but has not yet been applied to main-sequence gravity-mode pulsators. Aims. Our aim is to determine whether structure inversions for gravity-mode pulsators are feasible. We focus on the case of slowly rotating slowly pulsating B-type (SPB) stars. Methods. We computed and analyzed dipole mode kernels for three variables pairs: (ρ, c), (N2, c), and (N2, ρ). We assessed the potential of these kernels by predicting the oscillation frequencies of a model after perturbing its structure. We then tested two inversion methods, regularized least squares (RLS) and subtractive optimally localized averages (SOLA), using a model grid computed with the MESA stellar evolution code and the GYRE pulsation code. Results. We find that changing the stellar structure affects the oscillation frequencies in a nonlinear way. The oscillation modes for which this nonlinear dependency is the strongest are in resonance with the near-core peak in the buoyancy frequency. The near-core region of the star can be probed with SOLA, while RLS requires fine tuning to obtain accurate results. Both RLS and SOLA are strongly affected by the nonlinear dependencies on the structure differences, as these methods are based on a first-order approximation. These inversion methods need to be modified for meaningful applications of inversions to SPB stars. Conclusions. Our results show that inversions of gravity-mode pulsators are possible in principle, but that the typical inversion methods developed for solar-like oscillators are not applicable. Future work should focus on developing nonlinear inversion methods.

Neutrino Spectrum and Energy Loss Rates Due to Weak Processes on Hot 56Fe in Pre-Supernova Environment

Applying TQRPA calculations of Gamow–Teller strength functions in hot nuclei, we compute the (anti)neutrino spectra and energy loss rates arising from weak processes on hot 56Fe under pre-supernova conditions. We use a realistic pre-supernova model calculated by the stellar evolution code MESA. Taking into account both charged and neutral current processes, we demonstrate that weak reactions with hot nuclei can produce high-energy (anti)neutrinos. We also show that, for hot nuclei, the energy loss via (anti)neutrino emission is significantly larger than that for nuclei in their ground state. It is found that the neutral current de-excitation via the νν¯-pair emission is presumably a dominant source of antineutrinos. In accordance with other studies, we confirm that the so-called single-state approximation for neutrino spectra might fail under certain pre-supernova conditions.

A theory of mass transfer in binary stars

ABSTRACT Calculation of the mass transfer (MT) rate $\dot{M}_\text{d}$ of a Roche lobe overflowing star is a fundamental task in binary star evolution theory. Most of the existing MT prescriptions are based on a common set of assumptions that combine optically thick and optically thin regimes with different flow geometries. In this work, we develop a new model of MT based on the assumption that the Roche potential sets up a nozzle converging on the inner Lagrangian point and that the gas flows mostly along the axis connecting both stars. We derive a set of 1D hydrodynamic equations governing the gas flow with $\dot{M}_\text{d}$ determined as the eigenvalue of the system. The inner boundary condition directly relates our model to the structure of the donor obtained from 1D stellar evolution codes. We obtain an algebraic solution for the polytropic equation-of-state (EOS). This gives $\dot{M}_\text{d}$ within a factor of 0.9 to 1.0 of existing optically thick prescriptions and reduces to the existing optically thin prescription for isothermal gas. For a realistic EOS, we find that $\dot{M}_\text{d}$ differs by up to a factor of 4 from existing models. We illustrate the effects of our new MT model on a $30\, M_\odot$ low-metallicity star undergoing intensive thermal time-scale MT and find that it is more likely to become unstable to L2 overflow and common-envelope evolution than expected according to MT prescriptions. Our model provides a framework to include additional physics such as radiation or magnetic fields.

Carbon–oxygen ultra-massive white dwarfs in general relativity

ABSTRACT We employ the La Plata stellar evolution code, lpcode, to compute the first set of constant rest-mass carbon–oxygen ultra-massive white dwarf evolutionary sequences for masses higher than 1.29 M⊙ that fully take into account the effects of general relativity on their structural and evolutionary properties. In addition, we employ the lp-pul pulsation code to compute adiabatic g-mode Newtonian pulsations on our fully relativistic equilibrium white dwarf models. We find that carbon–oxygen white dwarfs more massive than 1.382 M⊙ become gravitationally unstable with respect to general relativity effects, being this limit higher than the 1.369 M⊙ we found for oxygen–neon white dwarfs. As the stellar mass approaches the limiting mass value, the stellar radius becomes substantially smaller compared with the Newtonian models. Also, the thermo-mechanical and evolutionary properties of the most massive white dwarfs are strongly affected by general relativity effects. We also provide magnitudes for our cooling sequences in different passbands. Finally, we explore for the first time the pulsational properties of relativistic ultra-massive white dwarfs and find that the period spacings and oscillation kinetic energies are strongly affected in the case of most massive white dwarfs. We conclude that the general relativity effects should be taken into account for an accurate assessment of the structural, evolutionary, and pulsational properties of white dwarfs with masses above ∼1.30 M⊙.

Asteroseismic age constraints on the open cluster NGC 2477 using oscillating stars identified with TESS FFI

Context. Asteroseismology is one of the few methods to derive ages of individual stars due to the high precision of their observations. Isochrone fitting is a powerful alternative method for deriving ages by studying clusters of stars. Pulsating stars in clusters should therefore allow for detailed studies of the stellar models. Aims. Our objectives are to exploit the NASA TESS data along with ESA Gaia data to search for and detect oscillations in cluster member stars. We analyse the intermediate-age open cluster NGC 2477, known to suffer from differential extinction, to explore if asteroseismology and cluster characteristics can help us understand the metallicity and extinction, as well as result in better age determinations than isochrone-fitting alone. Methods. We combined a multitude of recent observations from Gaia, high-resolution spectroscopy, and extinction maps to analyse the cluster and then search for and detect variability in the member stars using TESS full frame images (FFIs) data. To interpret all of these data, we used stellar structure, evolution and oscillation codes. Results. We conducted an in-depth analysis of the extinction and metallicity of NGC 2477, using the most recent spectroscopic, photometric, and extinction observations for the cluster. Analysis of dust and extinction maps confirmed that the differential extinction in the direction of the cluster is not due to the background. The cluster’s metallicity from high-resolution spectroscopy varies from 0.06 to 0.16 dex. We performed an isochrone fitting to the cluster using publically available isochrones (BASTI, MIST, and PARSEC), which provides a cluster age of between 0.6 to 1.1 Ga. Then using TESS FFI, we analysed the time dimension of the members of this cluster. We created optimised pixel light curves using the tessipack package which allows us to consider possible contamination by nearby stars. Using these light curves, we identified many interesting levels of variability of stars in this cluster, including binaries and oscillating stars. For the asteroseismic analysis, we selected a few uncontaminated A–F type oscillating stars and used the MESA and GYRE codes to interpret the frequency signals. By comparing the theoretical and the observed spectra, we identified frequency separations, Δν, for four stars. Then using the identified Δν and imposing that the best matched theoretical models have the same age, metallicity, and background extinction, we constrained the cluster’s age to 1.0 ± 0.1 Ga. Conclusions. We conclude that using the TESS FFI data, we can identify oscillating stars in clusters and constrain the age of the cluster using asteroseismology.

Lazarus Stars: Numerical investigations of stellar evolution with star-lifting as a life extension strategy

Abstract The aging and gradual brightening of the Sun will challenge Earth’s habitability in the next few billion years. If life exists elsewhere in the Universe, the aging of its host star similarly poses an existential threat. One solution, which we dub a Lazarus star, is for an advanced civilization to remove (or star-lift) mass from their host star at a rate that offsets the increase in luminosity, keeping the flux on the habitable planet(s) constant and extending the lifetime of their star. While this idea has existed since 1985 when it was first proposed by Criswell, numerical investigations of star-lifting have been lacking. Here, we use the stellar evolution code MESA to find mass vs. age and $\dot{M}$vs. age relations which would hold the flux on surrounding planets constant. We explore initial mass ranging from 0.2 M⊙ to 1.2 M⊙. For most stars with a mass initially below about 0.4M⊙, we find that star-lifting increases their main-sequence lifetimes up to 500 Gyr until they approach the hydrogen burning limit and star-lifting is no longer possible. For more massive stars, star-lifting increase main-sequence lifetimes by 1 Gyr to 100 Gyr, though they still enter the red-giant phase. For example, a Sun-like star has a main-sequence lifetime which can be increased by up to 3 Gyr. This requires a mass-loss rate of about 0.05 MCeres per year. We compare star-lifting to other survival strategies and briefly discuss methods for detecting these engineered stars.

Internal gravity waves in massive stars

Stars that are over 1.6 solar masses are generally known to possess convective cores and radiative envelopes, which allows for the propagation of outwardly travelling internal gravity waves (IGWs). Here, we study the generation and propagation of IGWs in such stars using two-dimensional, fully non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics. Compared to previous similar works, this study utilises radius-dependent thermal diffusivity profiles for five different stellar masses at the middle of the main sequence: 3 M⊙, 5 M⊙, 7 M⊙, 10 M⊙, and 13 M⊙. From the simulations, we find that the surface perturbations are larger for higher masses, but no noticeable trends are observed for the frequency slopes with different stellar masses. The slopes are also similar to the results from previous works. We compared our simulation results with stellar photometric data from a recent survey and we found that for frequency intervals above 8 μHz, there is a good agreement between the temperature frequency slopes from the simulations and the surface brightness variations of these observed stars. This indicates that the brightness variations are caused by core-generated IGWs.

Extreme evaporation of planets in hot thermally unstable protoplanetary discs: the case of FU Ori

ABSTRACT Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 yr ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures of ≳ 30 000 K. Such temperatures are reached in the inner 0.1 AU during thermal instability bursts. In our 1D time-dependent code the disc and an embedded planet interact through gravity, heat, and mass exchange. We use disc viscosity constrained by simulations and observations of dwarf novae instabilities, and we constrain planet properties with a stellar evolution code. We show that dusty gas giants born in the outer self-gravitating disc reach the innermost disc in a ∼O(104) yr with radius of ∼10RJ. We show that their EE rates are $\gtrsim 10^{-5} {\rm {\rm M}_{\odot }}$ yr−1; if this exceeds the background disc accretion activity then the system enters a planet-sourced mode. Like a stellar secondary in mass-transferring binaries, the planet becomes the dominant source of matter for the star, albeit for ∼O(100) yr. We find that a ∼6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of $\sim 10^{-6} {\rm {\rm M}_{\odot }}$ yr−1 appears to explain all that we currently know about FU Ori accretion outburst. More massive planets and/or planets in older less massive discs do not experience EE process. Future FUOR modelling may constrain planet internal structure and evolution of the earliest discs.