In Kawahara et al. (2018) and Masuda et al. (2019), we reported the discovery of four self-lensing binaries consisting of F/G-type stars and (most likely) white dwarfs whose masses range from 0.2 to 0.6 solar masses. Here we present their updated system parameters based on new radial velocity data from the Tillinghast Reflector Echelle Spectrograph at the Fred Lawrence Whipple Observatory, and the Gaia parallaxes and spectroscopic parameters of the primary stars. We also briefly discuss the astrophysical implications of these findings.

We investigate the roles of magnetically driven disk wind (MDW) and thermally driven photoevaporative wind (PEW) in the long-time evolution of protoplanetary disks. We start simulations from the early phase in which the disk mass is $0.118\,{\mathrm{M}_{\odot}}$ around a $1\,{\mathrm{M}_{\odot}}$ star and track the evolution until the disk is completely dispersed. We incorporate the mass loss by PEW and the mass loss and magnetic braking (wind torque) by MDW, in addition to the viscous accretion, viscous heating, and stellar irradiation. We find that MDW and PEW respectively have different roles: magnetically driven wind ejects materials from an inner disk in the early phase, whereas photoevaporation has a dominant role in the late phase in the outer ($\gtrsim1\,$au) disk. The disk lifetime, which depends on the combination of MDW, PEW, and viscous accretion, shows a large variation of $\sim1$-$20\,$Myr; the gas is dispersed mainly by the MDW and the PEW in the cases with a low viscosity and the lifetime is sensitive to the mass-loss rate and torque of the MDW, whereas the lifetime is insensitive to these parameters when the viscosity is high. Even in disks with very weak turbulence, the cooperation of MDW and PEW enables the disk dispersal within a few Myr.

We present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substantial amount of gas is also associated with 49 Ceti. The [C I] emission map shows two peaks inside the dust ring, and its overall extent is comparable to that of the dust continuum emission and the CO emission. We find that the [C I]/CO(J=3-2) intensity ratio significantly varies along the major axis. The ratio takes the minimum value of 1.8 around the dust peak position, and increases inwards and outwards. The enhanced ratio around the central star (~3) likely originates from the stellar UV radiation, while that in the outer disk (~10) from the interstellar UV radiation. Such complex distributions of the [C I] and CO(J=3-2) emission will be a key to understand the origin of the gas in 49 Ceti, and will also provide a stringent constraint on physical and chemical models of gaseous debris disks.

We report the discovery of the fifth self-lensing binary in which a low-mass white dwarf (WD) gravitationally magnifies its 15th magnitude G-star companion, KIC 8145411, during eclipses. The system was identified from a pair of such self-lensing events in the Kepler photometry, and was followed up with the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5m telescope at the Fred Lawrence Whipple Observatory and the High-Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. A joint analysis of the TRES radial velocities, the HDS spectrum, and the Kepler photometry of the primary star determines the WD mass $0.20\pm0.01\,M_\odot$, orbital semi-major axis $1.28\pm0.03\,\mathrm{au}$, and orbital eccentricity $0.14\pm0.02$. Because such extremely low-mass WDs cannot be formed in isolation within the age of the Galaxy, their formation is believed to involve binary interactions that truncated evolution of the WD progenitor. However, the observed orbit of the KIC 8145411 system is at least ten times wider than required for this scenario to work. The presence of this system in the Kepler sample, along with its similarities to field blue straggler binaries presumably containing WDs, may suggest that some 10% of post-AGB binaries with Sun-like primaries contain such anomalous WDs.

We want to investigate how planet formation is imprinted on stellar surface composition using up-to-date stellar evolution models. We simulate the evolution of pre-main-sequence stars as a function of the efficiency of heat injection during accretion, the deuterium mass fraction, and the stellar mass. For simplicity, we assume that planet formation leads to the late accretion of zero-metallicity gas, diluting the surface stellar composition as a function of the mass of the stellar outer convective zone. We adopt $150\,{\mathrm{M}_\oplus}(M_\star/\mathrm{M}_\odot)(Z/\mathrm{Z}_\odot)$ as an uncertain but plausible estimate of the mass of heavy elements that is not accreted by stars with giant planets, including our Sun. By combining our stellar evolution models to these estimates, we evaluate the consequences of planet formation on stellar surface composition. We show that after the first $\sim0.1$ Myr, the evolution of the convective zone follows classical evolutionary tracks within a factor of two in age. We find that planet formation should lead to a scatter in stellar surface composition that is larger for high-mass stars than for low-mass stars. We predict a spread in [Fe/H] of approximately $0.02$ dex for stars with $T_\mathrm{eff}\sim 5500\,$K, marginally compatible with differences in metallicities observed in some binary stars with planets. Stars with $T_\mathrm{eff}\geq 7000\,$K may show much larger [Fe/H] deficits, by 0.6 dex or more, compatible with the existence of refractory-poor $\lambda$ Boo stars. We also find that planet formation may explain the lack of refractory elements seen in the Sun as compared to solar twins, but only if the ice-to-rock ratio in the solar-system planets is less than $\approx0.4$ and planet formation began less than $\approx1.3$ Myr after the beginning of the formation of the Sun. (abbreviated)

The recent discovery of hot dust grains in the vicinity of main-sequence stars has become a hot issue among the scientific community of debris disks. Hot grains must have been enormously accumulated near their sublimation zones, but it is a mystery how such a high concentration of hot grains is sustained. The most difficult conundrum is that the size of hot dust grains is estimated to lie in the submicrometer range, while submicrometer-sized grains are instantly swept away from near-stellar environments by stellar radiation pressure. One and only mechanism proposed for prolonging the residence time of hot grains in the near-stellar environments is trapping of charged nanoparticles by stellar magnetic fields. We revisit the model of magnetic grain trapping around main-sequence stars of various spectral classes by taking into account sublimation and electric charging of the grains. The model of magnetic grain trapping predicts that hot dust grains are present in the vicinity of main-sequence stars with high rotation velocities and intermediate magnetic-field strengths. On the contrary, we find that the detection of hot dust grains has no correlation with the rotation velocities of central stars nor the magnetic field strengths of the stars. Our numerical evaluation of electric grain charging indicates that the surface potential of submicrometer-sized grains in the vicinity of main-sequence stars is typically 4-5 V, which is one order of magnitude smaller than the value assumed by the model of magnetic grain trapping. On the basis of our numerical simulation on sublimation of dust grains in the vicinity of a star, it turns out that their lives end due to sublimation in a timescale much shorter than the period of one revolution at the gyroradius. It is, therefore, infeasible to dynamically extend the dwell time of hot grains inside the sublimation zone by magnetic trapping, ...

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.)

Formation, evolution, and dispersal processes of protoplanetary disks are investigated and the disk lifetime is estimated. Gravitational collapse of a pre-stellar core forms both a central star and a protoplanetary disk. The central star grows by accretion from the disk, and irradiation by the central star heats up the disk and generates thermal wind, which results in the disk dispersal. We calculate the evolution of protoplanetary disks from their parent pre-stellar cores to dispersal of the disks. We find that the disk lifetimes of typical pre-stellar cores are around 2--4 million years (Myr). A pre-stellar core with high angular momentum forms a larger disk whose lifetime is long, while a disk around a X-ray luminous star has a short lifetime. Integrating the disk lifetimes under various mass and angular velocity of prestellar cores and X-ray luminosities of young stellar objects, we obtain disk fraction at a given stellar age and mean lifetime of the disks. Our model indicates that the mean lifetime of protoplanetary disks is 3.7 Myr, which is consistent with the observational estimate from young stellar clusters. We also find that the dispersion of X-ray luminosity is needed to reproduce the observed disk fraction.

The aim of this work is to search for Li-rich giants in a sample of clusters where planets have been searched, thus we can study the planet engulfment scenario to explain Li replenishment using a proper comparison sample of stars without detected giant planets. We derived Li abundances for a sample of 67 red giant stars in 12 different open clusters using standard spectral synthesis techniques and high resolution spectra (from HARPS and UVES). We also determined masses, ages and radius from PARSEC stellar isochrones to constrain the evolutionary stage of these stars. We find three stars in different clusters with clearly enhanced Li abundances compared to other stars within the cluster. Interestingly, the only two stars with a detected substellar companion in our sample belong to that group. One of the planet hosts, NGC2423No3, might lie close to the luminosity bump on the HR diagram, a phase where Li production by the Cameron-Fowler process is supported by extra-mixing to bring fresh Li up to the surface. On the other hand, NGC4349No127 is a more massive and more evolved giant that does not seem to be in the evolutionary phase where other Li-rich stars are found. We discuss the possibility that the Li enhancement of this star is triggered by the engulfment of a planet, considering that close-in planets hardly survive the RGB-tip and the early-AGB phases.

Recent radial-velocity surveys for GK clump giants have revealed that planets also exist around ~1.5-3 Msun stars. However, no planets have been found inside 0.6 AU around clump giants, in contrast to solar-type main-sequence stars, many of which harbor short-period planets such as hot Jupiters. In this study we examine the possibility that planets were engulfed by host stars evolving on the red-giant branch (RGB). We integrate the orbital evolution of planets in the RGB and helium burning (HeB) phases of host stars, including the effects of stellar tide and stellar mass loss. Then we derive the critical semimajor axis (or the survival limit) inside which planets are eventually engulfed by their host stars after tidal decay of their orbits. Especially, we investigate the impact of stellar mass and other stellar parameters on the survival limit in more detail than previous studies. In addition, we make detailed comparison with measured semimajor axes of planets detected so far, which no previous study did. We find that the critical semimajor axis is quite sensitive to stellar mass in the range between 1.7 and 2.1 Msun, which suggests a need for careful comparison between theoretical and observational limits of existence of planets. Our comparison demonstrates that all those planets are beyond the survival limit, which is consistent with the planet-engulfment hypothesis. However, on the high-mass side (> 2.1 Msun), the detected planets are orbiting significantly far from the survival limit, which suggests that engulfment by host stars may not be the main reason for the observed lack of short-period giant planets. To confirm our conclusion, the detection of more planets around clump giants, especially with masses > 2.5 Msun, is required.

Around low- and intermediate-mass (1.5-3 M_sun) red giants, no planets have been found inside 0.6 AU. Such a paucity is not seen in the case of 1 M_sun main sequence stars. In this study, we examine the possibility that short-period planets were engulfed by their host star evolving off the main sequence. To do so, we have simulated the orbital evolution of planets, including the effects of stellar tide and mass loss, to determine the critical semimajor axis, a_crit, beyond which planets survive the RGB expansion of their host star. We have found that a_crit changes drastically around 2 M_sun: In the lower-mass range, a_crit is more than 1 AU, while a_crit is as small as about 0.2 AU in the higher-mass range. Comparison with measured semimajor axes of known planets suggests that there is a lack of planets that only planet engulfment never accounts for in the higher-mass range. Whether the lack is real affects our understanding of planet formation. Therefore, increasing the number of planet samples around evolved intermediate-mass stars is quite meaningful to confirm robustness of the lack of planets.