We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star V410 Tau based on data collected mostly with SPIRou, the near-infrared (NIR) spectropolarimeter recently installed at the Canada-France-Hawaii Telescope, as part of the SPIRou Legacy Survey large programme, and with TESS between October and December 2019. Using Zeeman-Doppler Imaging (ZDI), we obtained the first maps of photospheric brightness and large-scale magnetic field at the surface of this young star derived from NIR spectropolarimetric data. For the first time, ZDI is also simultaneously applied to high-resolution spectropolarimetric data and very-high-precision photometry. V410 Tau hosts both dark and bright surface features and magnetic regions similar to those previously imaged with ZDI from optical data, except for the absence of a prominent dark polar spot. The brightness distribution is significantly less contrasted than its optical equivalent, as expected from the difference in wavelength. The large-scale magnetic field (~410 G), found to be mainly poloidal, features a dipole of ~390 G, again compatible with previous studies at optical wavelengths. NIR data yield a surface differential rotation slightly weaker than that estimated in the optical at previous epochs. Finally, we measured the radial velocity of the star and filtered out the stellar activity jitter using both ZDI and Gaussian Process Regression down to a precision of ~0.15 and 0.08 $\mathrm{km\,s^{-1}}$ RMS, respectively, confirming the previously published upper limit on the mass of a potential close-in massive planet around V410 Tau.

Intermediate Mass Black Hole (IMBH) mergers with masses $10^4 - 10^6$ $M_{\odot}$ are expected to produce gravitational waves (GWs) detectable by the Laser Interferometer Space Antenna (LISA) with high signal to noise ratios out to redshift 20. IMBH mergers are expected to take place within dwarf galaxies, however, the dynamics, timescales, and effect on their hosts are largely unexplored. In a previous study, we examined how IMBHs would pair and merge within nucleated dwarf galaxies. IMBHs in nucleated hosts evolve very efficiently, forming a binary system and coalescing within a few hundred million years. Although the fraction of dwarf galaxies ($10^7$ M$_{\odot} \leq$ $M_{\star} \leq 10^{10}$ M$_{\odot}$) hosting nuclear star clusters is between 60-100\%, this fraction drops to 20-70\% for lower-mass dwarfs ($M_{\star}\approx 10^7$ M$_{\odot}$), with the largest drop in low-density environments. Here, we extend our previous study by performing direct $N-$body simulations to explore the dynamics and evolution of IMBHs within non-nucleated dwarf galaxies, under the assumption that IMBHs exist within these dwarfs. To our surprise, none of IMBHs in our simulation suite merge within a Hubble time, despite many attaining high eccentricities $e \sim 0.7-0.95$. We conclude that extremely low stellar density environments in the centers of non-nucleated dwarfs do not provide an ample supply of stars to interact with IMBHs binary resulting in its stalling, in spite of triaxiality and high eccentricity, common means to drive a binary to coalescence. Our findings underline the importance of considering all detailed host properties to predict IMBH merger rates for LISA.

In this paper we present the "Small Bodies: Near and Far" Infrared Database, an easy-to-use tool intended to facilitate the modeling of thermal emission of small Solar System bodies. Our database collects thermal emission measurements of small Solar Systems targets that are otherwise available in scattered sources and gives a complete description of the data, with all information necessary to perform direct scientific analyses and without the need to access additional, external resources. This public database contains representative data of asteroid observations of large surveys (e.g. AKARI, IRAS and WISE) as well as a collection of small body observations of infrared space telescopes (e.g. the Herschel Space Observatory) and provides a web interface to access this data (this https URL). We also provide an example for the direct application of the database and show how it can be used to estimate the thermal inertia of specific populations, e.g. asteroids within a given size range. We show how different scalings of thermal inertia with heliocentric distance (i.e. temperature) may affect our interpretation of the data and discuss why the widely-used radiative conductivity exponent ($\alpha$=-3/4) might not be adequate in general, as hinted by previous studies.

The Rosette complex is a well studied region of the galactic plane which presents the apparent characteristics of a triggered star forming region. This is however still debated as no strong evidence corroborates this statement. We focus on characterizing the young stellar population in the Rosette complex to improve our understanding of the processes that regulate the star formation in this region. We propose an original method that relies on the joint analysis of the star color and density in the near-infrared. It leads to mapping the molecular cloud spatial distribution and detecting the embedded clusters with their characterization in terms of member number and age estimation. We have identified 13 clusters, 2 of which are new discoveries, and we estimate that the total number of young stellar objects in the Rosette ranges between 4000 and 8000 members. We find that the age distribution of the young clusters is not consistent with a general triggered scenario for the star formation in this molecular cloud. This study proves that the Rosette complex evolution is not governed by the influence of its OB star population. It suggests that the simple morphological appearance of an active region is not sufficient to conclude much about the triggering role in the star formation process. Our method of constraining the cluster properties using UKIDSS and WISE data has proven efficient, and studies of other regions of the galactic plane would definitely benefit from this approach.

The composition of the solar corona differs from that of the photosphere, with the plasma thought to fractionate in the solar chromosphere according to the First Ionisation Potential (FIP) of the different elements. This produces a FIP bias, wherein elements with a low FIP are preferentially enhanced in the corona compared to their photospheric abundance, but direct observations of this process remain elusive. Here we use a series of spectroscopic observations of Active Region AR 12759 as it transited the solar disc over a period of 6 days from 2-7 April 2020 taken using the Hinode Extreme ultraviolet Imaging Spectrometer (EIS) and Interface Region Imaging Spectrograph (IRIS) instruments to look for signatures of plasma fractionation in the solar chromosphere. Using the Si X/S X and Ca XIV/Ar XIV diagnostics, we find distinct differences between the FIP bias of the leading and following polarities of the active region. The widths of the IRIS Si IV lines exhibited clear differences between the leading and following polarity regions, indicating increased unresolved wave activity in the following polarity region compared to the leading polarity region, with the chromospheric velocities derived using the Mg II lines exhibiting comparable, albeit much weaker, behaviour. These results are consistent with plasma fractionation via resonant/non-resonant waves at different locations in the solar chromosphere following the ponderomotive force model, and indicate that IRIS could be used to further study this fundamental physical process.

Context. The dust and gas temperature in proto-planetary disks play critical roles in determining their chemical evolution and influencing planet formation processes. Aims. We attempted an accurate measurement of the dust and CO temperature profile in the edge-on disk of the Flying Saucer. Methods. We used the unique properties of the Flying Saucer, its edge-on geometry and its fortunate position in front of CO clouds with different brightness temperatures to provide independent constraints on the dust temperature. We compared it with the dust temperature derived using the radiative transfer code DiskFit and the CO gas temperature. Results. We find clear evidence for a substantial gas temperature vertical gradient, with a cold (10 K) disk mid-plane and a warmer CO layer where T(r) is 27 K at 100 au, dropping with exponent 0.3. Direct evidence for CO depletion in the mid-plane, below about 1 scale height, is also found. At this height, the gas temperature is 15-20 K, consistent with the expected CO freeze out temperature. The dust disk appears optically thin at 345 GHz, and exhibits moderate settling.

The K2 mission of the Kepler Space Telescope allowed the observations of light curves of small solar system bodies throughout the whole Solar system. In this paper we present the results of a collection of K2 transneptunian object observations, between Campaigns C03 (November 2014 -- February 2015) to C19 (August -- September, 2018), which includes 66 targets. Due to the faintness of our targets the detectability rate of a light curve period is $\sim$56%, notably lower than in the case of other small body populations, like Hildas or Jovian trojans. We managed to obtain light curve periods with an acceptable confidence for 37 targets; the majority of these cases are new identifications. We were able to give light curve amplitude upper limits for the other 29 targets. Several of the newly detected light curve periods are longer than $\sim$24 h, in many cases close to $\sim$80 h, i.e., these targets are slow rotators. This relative abundance of slowly rotating objects is similar to that observed among Hildas, Jovian trojans and Centaurs in the K2 mission, and also among main belt asteroids measured with the TESS Space Telescope. Transneptunian objects show notably higher light curve amplitudes at large (D $\gtrsim$ 300 km) sizes than that found among large main belt asteroids, in contrast to the general expectation that due to their lower compressive strength they reach hydrostatic equlibrium at smaller sizes than their inner solar system counterparts.

Long-lived radioactive nuclides, such as $^{40}$K, $^{232}$Th, $^{235}$U and $^{238}$U, contribute to persistent heat production in the mantle of terrestrial-type planets. As refractory elements, the concentrations of Th and U in a terrestrial exoplanet are implicitly reflected in the photospheric abundances in the stellar host. However, a robust determination of these stellar abundances is difficult in practice owing to the general paucity and weakness of the relevant spectral features. We draw attention to the refractory, $r-$process element europium, which may be used as a convenient and practical proxy for the population analysis of radiogenic heating in exoplanetary systems. As a case study, we present a determination of Eu abundances in the photospheres of $\alpha$ Cen A and B. We find that europium is depleted with respect to iron by $\sim$ 0.1 dex and to silicon by $\sim$ 0.15 dex compared to solar in both binary components. To first order, the measured Eu abundances can be converted to the abundances of $^{232}$Th, $^{235}$U and $^{238}$U with observational constraints while the abundance of $^{40}$K is approximated independently with a Galactic chemical evolution model. We find that the radiogenic heat budget in an $\alpha$-Cen-Earth is $73.4^{+8.3}_{-6.9}$ TW upon its formation and $8.8^{+1.7}_{-1.3}$ TW at the present day, respectively $23\pm5$ % and $54\pm5$ % lower than that in the Hadean and modern Earth. As a consequence, mantle convection in an $\alpha$-Cen-Earth is expected to be overall weaker than that of the Earth (assuming other conditions are the same) and thus such a planet would be less geologically active, suppressing its long-term potential to recycle its crust and volatiles. With Eu abundances being available for a large sample of Sun-like stars, the proposed approach can extend our ability to make predictions about the nature of other rocky worlds.

The optical observations of Ic-4 supernova (SN) 2016coi/ASASSN-16fp, from $\sim 2$ to $\sim450$ days after explosion, are presented along with analysis of its physical properties. The SN shows the broad lines associated with SNe Ic-3/4 but with a key difference. The early spectra display a strong absorption feature at $\sim 5400$ \AA\ which is not seen in other SNe~Ic-3/4 at this epoch. This feature has been attributed to He I in the literature. Spectral modelling of the SN in the early photospheric phase suggests the presence of residual He in a C/O dominated shell. However, the behaviour of the He I lines are unusual when compared with He-rich SNe, showing relatively low velocities and weakening rather than strengthening over time. The SN is found to rise to peak $\sim 16$ d after core-collapse reaching a bolometric luminosity of Lp $\sim 3\times10^{42}$ \ergs. Spectral models, including the nebular epoch, show that the SN ejected $2.5-4$ \msun\ of material, with $\sim 1.5$ \msun\ below 5000 \kms, and with a kinetic energy of $(4.5-7)\times10^{51}$ erg. The explosion synthesised $\sim 0.14$ \msun\ of 56Ni. There are significant uncertainties in E(B-V)host and the distance however, which will affect Lp and MNi. SN 2016coi exploded in a host similar to the Large Magellanic Cloud (LMC) and away from star-forming regions. The properties of the SN and the host-galaxy suggest that the progenitor had $M_\mathrm{ZAMS}$ of $23-28$ \msun\ and was stripped almost entirely down to its C/O core at explosion.

Aims. We introduce a novel way to identify new compact hierarchical triple stars by exploiting the huge potential of Gaia DR3 and also its future data releases. We aim to increase the current number of compact hierarchical triples significantly. Methods. We utilize several eclipsing binary catalogs from different sky surveys totaling more than 1 million targets for which we search for Gaia DR3 Non-single Star orbital solutions with periods substantially longer than the eclipsing periods of the binaries. Those solutions in most cases should belong to outer orbits of tertiary stars in those systems. We also try to validate some of our best-suited candidates using TESS eclipse timing variations. Results. We find 403 objects with suitable Gaia orbital solutions of which 27 are already known triple systems. This makes 376 newly identified hierarchical triple system candidates in our sample. We analyze the cumulative probability distribution of the outer orbit eccentricities and find that it is very similar to the ones found by earlier studies based on the observations of the Kepler and OGLE missions. We found measurable non-linear eclipse timing variations or third-body eclipses in the TESS data for 192 objects which we also consider to be confirmed candidates. Out of these, we construct analytical light-travel time effect models for the eclipse timing variations of 22 objects with well-sampled TESS observations. We compare the outer orbital parameters from our solutions with the ones from the Gaia solutions and find that the most reliable orbital parameter is the orbital period, while the values of the other parameters should be used with caution.

Blue large-amplitude pulsators (BLAPs) are a recently discovered group of hot pulsating stars whose evolutionary status remains uncertain. Their supposed progenitors are either $\simeq 0.3M_{\odot}$ shell H-burning stars or $\simeq 1.0M_{\odot}$ core He-burning stars, both relying on mass loss or a merger event in a (rarely observed) close interacting binary system. With the goal to understand the stellar masses of BLAPs, we therefore carried out a linear non-adiabatic analysis of a grid of models computed using mesa-rsp, with appropriate input stellar parameters $ZXMLT_{\rm eff}$ and convection parameter sets. We discuss the impact of stellar mass, metallicity, helium abundance and convection parameters on the theoretical instability regions of BLAPs. We also derive new theoretical period relations; our theoretical period relations using low stellar masses seem to be in better agreement with the observed period relations. Although only two BLAPS have been observed to be multi-periodic oscillator so far, we analyse theoretical $P_{1O}/P_F$ ratios and compare these values with other classical pulsators. Furthermore, we provide the first asteroseismic mass estimate for the triple-mode pulsator, OGLE-BLAP-030 which seems to be well-constrained in the range of $0.62-0.64 M_{\odot}$ with a high metallicity of $Z=0.07$, albeit with a few sources of uncertainty involved. This would place the BLAP star intermediate to the two proposed mass scenarios so far.

We present a new technique to generate the light curves of RRab stars in different photometric bands ($I$ and $V$ bands) using Artificial Neural Networks (ANN). A pre-computed grid of models was used to train the ANN, and the architecture was tuned using the $I$ band light curves. The best-performing network was adopted to make the final interpolators in the $I$ and $V$ bands. The trained interpolators were used to predict the light curve of RRab stars in the Magellanic Clouds, and the distances to the LMC and SMC were determined based on the reddening independent Wesenheit index. The estimated distances are in good agreement with the literature. The comparison of the predicted and observed amplitudes, and Fourier amplitude ratios showed good agreement, but the Fourier phase parameters displayed a few discrepancies. To showcase the utility of the interpolators, the light curve of the RRab star EZ Cnc was generated and compared with the observed light curve from the Kepler mission. The reported distance to EZ Cnc was found to be in excellent agreement with the updated parallax measurement from Gaia EDR3. Our ANN interpolator provides a fast and efficient technique to generate a smooth grid of model light curves for a wide range of physical parameters, which is computationally expensive and time-consuming using stellar pulsation codes.

Few spectra of directly-imaged exoplanets have been obtained in the mid-infrared (> 3 $\mu$m). This region is particularly rich in molecular spectral signatures, whose measurements can help recover atmospheric parameters and provide a better understanding of giant planet formation and atmospheric dynamics. In the past years, exoplanet interferometry with the VLTI/GRAVITY instrument has provided medium-resolution spectra of a dozen substellar companions in the near infrared. The 100-meter interferometric baselines allow for the stellar and planetary signals to be efficiently disentangled at close angular separations (< 0.3''). We aim to extend this technique to the mid-infrared using MATISSE, the VLTI's mid-infrared spectro-interferometer. We take advantage of the fringe tracking and off-axis pointing capabilities recently brought by the GRA4MAT upgrade. Using this new mode, we observed the giant planet $\beta$ Pictoris b in L and M bands (2.75-5 $\mu$m) at a spectral resolution of 500. We developed a method to correct chromatic dispersion and non-common paths effects in the fringe phase and modelled the planet astrometry and stellar contamination. We obtained a high-signal-to-noise spectrum of $\beta$ Pictoris b, showing the planet continuum in L (for the first time) and M bands, which contains broad absorption features of H$_2$O and CO. In conjunction with a new GRAVITY spectrum, we modelled it with the ForMoSA nested sampling tool and the Exo-REM grid of atmospheric models, and found a solar carbon-to-oxygen ratio in the planet atmosphere. This study opens the way to the characterization of fainter and closer-in planets with MATISSE, which could complement the JWST at angular separations too close for it to obtain exoplanet spectra. Starting in 2025, the new adaptive optics system brought by the GRAVITY+ upgrade will further extend the detection limits of MATISSE.

We have discovered a triply eclipsing triple-star system, TIC 290061484, with the shortest known outer period, Pout, of only 24.5 days. This "eclipses" the previous record set by lambda Tauri at 33.02 days, which held for 68 yr. The inner binary, with an orbital period of Pin = 1.8 days, produces primary and secondary eclipses and exhibits prominent eclipse timing variations with the same periodicity as the outer orbit. The tertiary star eclipses, and is eclipsed by, the inner binary with pronounced asymmetric profiles. The inclinations of both orbits evolve on observable timescales such that the third-body eclipses exhibit dramatic depth variations in TESS data. A photodynamical model provides a complete solution for all orbital and physical parameters of the triple system, showing that the three stars have masses of 6.85, 6.11, and 7.90 MSun, radii near those corresponding to the main sequence, and Teff in the range of 21,000-23,700 K. Remarkably, the model shows that the triple is in fact a subsystem of a hierarchical 2+1+1 quadruple with a distant fourth star. The outermost star has a period of ~3200 days and a mass comparable to the stars in the inner triple. In ~20 Myr, all three components of the triple subsystem will merge, undergo a Type II supernova explosion, and leave a single remnant neutron star. At the time of writing, TIC 290061484 is the most compact triple system and one of the tighter known compact triples (i.e., Pout/Pin = 13.7).

Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$\sigma$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared.

Variable accretion in young stellar objects reveals itself photometrically and spectroscopically over a continuum of timescales and amplitudes. Most dramatic are the large outbursts (e.g., FU Ori, V1647 Ori, and EX Lup type events), but more frequent are the less coherent, smaller burst-like variations in accretion rate. Improving our understanding of time-variable accretion directly addresses the fundamental question of how stars gain their masses. We review variability phenomena, as characterized from observations across the wavelength spectrum, and how those observations probe underlying physical conditions. The diversity of observed lightcurves and spectra at optical and infrared wavelengths defies a simple classification of outbursts and bursts into well-defined categories. Mid-infrared and submillimeter wavelengths are sensitive to lower-temperature phenomena and more embedded, younger sources, and it is currently unclear if observed flux variations probe similar or distinct physics relative to the shorter wavelengths. We highlight unresolved issues and emphasize the value of spectroscopy, multiwavelength studies, and ultimately patience in using variable accretion to understand stellar mass assembly.

Context. Cepheids are pulsating stars that play a crucial role in several astrophysical contexts. Among the different types, the Classical Cepheids are fundamental tools for the calibration of the extragalactic distance ladder. They are also powerful stellar population tracers in the context of Galactic studies. The Gaia Third Data Release (DR3) publishes improved data on Cepheids collected during the initial 34 months of operations. Aims. We present the Gaia DR3 catalogue of Cepheids of all types, obtained through the analysis carried out with the Specific Object Study (SOS) Cep&RRL pipeline. Methods. We discuss the procedures adopted to clean the Cepheid sample from spurious objects, to validate the results, and to re-classify sources with a wrong outcome from the SOS Cep&RRL pipeline. Results. The Gaia DR3 includes multi-band time-series photometry and characterisation by the SOS Cep&RRL pipeline for a sample of 15,006 Cepheids of all types. The sample includes 4,663, 4,616, 321 and 185 pulsators, distributed in the LMC, SMC, M31 and M33, respectively, as well as 5 221 objects in the remaining All Sky sub-region which includes stars in the MW field/clusters and in a number of small satellites of our Galaxy. Among this sample, 327 objects were known as variable stars in the literature but with a different classification, while, to the best of our knowledge, 474 stars have not been reported before to be variable stars and therefore they likely are new Cepheids discovered by Gaia.

How multiple close-in super-Earths form around stars with masses lower than that of the Sun is still an open issue. Several recent modeling studies have focused on planet formation around M-dwarf stars, but so far no studies have focused specifically on K dwarfs, which are of particular interest in the search for extraterrestrial life. We aim to reproduce the currently known population of close-in super-Earths observed around K-dwarf stars and their system characteristics. We performed 48 high-resolution N-body simulations of planet formation via planetesimal accretion using the existing GENGA software running on GPUs. In the simulations we varied the initial disk mass and the solid and gas surface density profiles. Each simulation began with 12000 bodies with radii of between 200 and 2000 km around two different stars, with masses of 0.6 and 0.8 $M_{\odot}$. Most simulations ran for 20 Myr, with several simulations extended to 40 or 100 Myr. The mass distributions for the planets with masses between 2 and 12 $M_\oplus$ show a strong preference for planets with masses M_p&lt;6 $M_\oplus$ and a lesser preference for planets with larger masses, whereas the mass distribution for the observed sample increases almost linearly. However, we managed to reproduce the main characteristics and architectures of the known planetary systems and produce mostly long-term angular-momentum-deficit-stable, nonresonant systems, but we require an initial disk mass of 15 $M_\oplus$ or higher and a gas surface density value at 1 AU of 1500 g cm$^{-2}$ or higher. Our simulations also produce many low-mass planets with M&lt;2 $M_\oplus$, which are not yet found in the observed population, probably due to the observational biases. The final systems contain only a small number of planets, which could possibly accrete substantial amounts of gas, and these formed after the gas had mostly dissipated.

We outline the impact of a small two-band UV-photometry satellite mission on the field of stellar physics, magnetospheres of stars, binaries, stellar clusters, interstellar matter, and exoplanets. On specific examples of different types of stars and stellar systems, we discuss particular requirements for such satellite missions in terms of specific mission parameters such as bandpass, precision, cadence, and mission duration. We show that such a mission may provide crucial data not only for hot stars that emit most of their light in UV, but also for cool stars, where UV traces their activity. This is important, for instance, for exoplanetary studies, because the level of stellar activity influences habitability. While the main asset of the two-band UV mission rests in time-domain astronomy, an example of open clusters proves that such a mission would be important also for the study of stellar populations. Properties of the interstellar dust are best explored when combining optical and IR information with observations in UV. It is well known that dust absorbs UV radiation efficiently. Consequently, we outline how such a UV mission can be used to detect eclipses of sufficiently hot stars by various dusty objects and study disks, rings, clouds, disintegrating exoplanets or exoasteroids. Furthermore, UV radiation can be used to study the cooling of neutron stars providing information about the extreme states of matter in the interiors of neutron stars and used for mapping heated spots on their surfaces.

Tidal disruption of stars in dense nuclear star clusters containing supermassive central black holes (SMBH) is modeled by high-accuracy direct N-body simulation. Stars getting too close to the SMBH are tidally disrupted and a tidal disruption event (TDE) happens. TDEs probe properties of SMBH, their accretion disks, and the surrounding nuclear stellar cluster. In this paper we compare rates of full tidal disruption events (FTDE) with partial tidal disruption events (PTDE). Since a PTDE does not destroy the star, a leftover object emerges; we use the term 'leftover star' for it; two novel effects occur in the simulation: (1) variation of the leftover star's mass and radius, (2) variation of the leftover star's orbital energy. After switching on these two effects in our simulation, the number of FTDEs is reduced by roughly 28%, and the reduction is mostly due to the ejection of the leftover stars from PTDEs coming originally from relatively large distance. The number of PTDEs is about 75% higher than the simple estimation given by Stone et al. (2020), and the enhancement is mainly due to the multiple PTDEs produced by the leftover stars residing in the diffusive regime. We compute the peak mass fallback rate for the PTDEs and FTDEs recorded in the simulation, and find 58% of the PTDEs have peak mass fallback rate exceeding the Eddington limit, and the number of super-Eddington PTDEs is 2.3 times the number of super-Eddington FTDEs.