Planetary engulfment events have long been proposed as a lithium (Li) enrichment mechanism contributing to the population of Li-rich giants $(A(\mathrm{Li}) \geq 1.5$ dex). Using $GALAH$ survey data and MESA stellar models, we calculate the strength and duration of the Li enrichment signature produced in the convective envelope of a host star that has engulfed a hot Jupiter (HJ) companion. We consider solar-metallicity stars in the mass range of $1.0-2.0~\mathrm{M_{\odot}}$ and the Li supplied by a HJ of $1.0~\mathrm{M_{J}}$. We explore engulfment events that occur near the main sequence turn-off (MSTO) and out to orbital separations of $R_{\star}{\sim}~0.1~\mathrm{AU}= 22~\mathrm{R_{\odot}}$. We map our results onto the Hertzsprung-Russell (H-R) Diagram, revealing the parameter space where planetary engulfment events produce significant Li enrichment signatures. We also map the associated survival times of these signatures, which range across 9 orders of magnitude. Our calculations indicate that if the HJ engulfment event occurs near the MSTO or on the subgiant branch, Li enrichment can be measured at a $5\sigma$ confidence level and with strengths that exceed meteoritic abundance measurements. Moreover, for stars of $1.4~\mathrm{M_{\odot}}$, these signatures are predicted to survive for up to 1 Gyr. We determine that Li enrichment beyond the subgiant branch must be produced by other mechanisms, such as the Cameron-Fowler process or accretion of material from an AGB companion.

We employ hydrodynamical simulations to investigate the underlying mechanism responsible for the low levels of gas and dust in globular clusters. Our models examine the competing effects of energy and mass supply from the various components of the evolved stellar population for globular clusters 47 Tucanae, M15, NGC 6440, and NGC 6752. Ignoring all other gas evacuation processes, we find that the energy output from the stars that have recently turned off the main sequence are capable of effectively clearing the evolved stellar ejecta and producing intracluster gas densities consistent with current observational constraints. This result distinguishes a viable gas and dust evacuation mechanism that is ubiquitous in globular clusters. In addition, we extend our analysis to probe the efficiency of pulsar wind feedback in globular clusters. The detection of intracluster ionized gas in 47~Tucanae allows us to place particularly strict limits on pulsar wind thermalization efficiency, which must be extremely low in the cluster's core in order to be in accordance with the observed density constraints.

We present a catalog of 1,143 periodic variables, compiled from our image-subtracted photometric analysis of the K2 Campaign-0 super stamp. This super stamp is centered on the open clusters M35 and NGC 2158. Approximately 46% of our periodic variables were previously unreported. Of the catalog variables, we find that 331 are members of M35 and 56 are members of NGC 2158 (${P_{\mathrm{m}}~>~0.5}$). Our catalog contains two new transiting exoplanet candidates, both of which orbit field stars. The smaller planet candidate has a radius of $0.35~\pm~0.04 \ \rm R_{J}$ and orbits a K dwarf ($K_{\mathrm{p}}=15.4$ mag) with a transit depth of 2.9 millimag. The larger planet candidate has a radius of $0.72~\pm~0.02 \ \rm R_{J}$ and orbits a late G type star ($K_{\mathrm{p}}=15.7$ mag) with a transit depth of 2.2 millimag. The larger planet candidate maybe an unresolved binary or a false alarm. Our catalog includes 44 eclipsing binaries, including ten new detections. Of the eclipsing binaries, one is an M35 member and five are NGC 2158 members. Our catalog contains a total of 1,097 non-transiting variable stars, including a field $\delta$ Cepheid exhibiting double mode pulsations, 561 rotational variables, and 251 pulsating variables (primarily $\gamma$ Doradus and $\delta$ Scuti types). The periods of our catalog sources range between 43 minutes to 24 days. The known ages of our reported cluster variables will facilitate investigations of a variety of stellar evolutionary processes.

We have discovered a young M star of mass $0.16\,M_\odot$ and radius $0.63\,R_\odot$, likely in the Upper Sco Association, that exhibits only a single $80\%$ deep occultation of 1-day duration. The star has frequent flares and a low-amplitude rotational modulation, but is otherwise quiet over 160 days of cumulative observation during K2 Campaigns C2 and C15. We discuss how such a deep eclipse is not possible by one star crossing another in any binary or higher-order stellar system in which no mass transfer has occurred. The two possible explanations we are left with are (1) orbiting dust or small particles (e.g., a disk bound to a smaller orbiting body, or unbound dust that emanates from such a body); or (2) a transient accretion event of dusty material near the corotation radius of the star. In either case, the time between such occultation events must be longer than $\sim$80 days. We model a possible orbiting occulter both as a uniform elliptically shaped surface (e.g., an inclined circular disk) and as a `dust sheet' with a gradient of optical depth behind its leading edge. The required masses in such dust features are then $\gtrsim 3 \times 10^{19}$ g and $\gtrsim 10^{19}$ g, for the two cases, respectively.

Planets accompany most sun-like stars. The orbits of many are sufficiently close that they will be engulfed when their host stars ascend the giant branch. This Letter compares the power generated by orbital decay of an engulfed planet to the intrinsic stellar luminosity. Orbital decay power is generated by drag on the engulfed companion by the surrounding envelope. As stars ascend the giant branch their envelope density drops and so does the power injected through orbital decay, scaling approximately as $L_{\rm decay} \propto R_*^{-9/2}$. Their luminosity, however, increases along the giant branch. These opposed scalings indicate a crossing, where $L_{\rm decay}= L_*$. We consider the engulfment of planets along isochrones in the Hertzsprung-Russell (H-R) diagram. We find that the conditions for such a crossing occur around $L_*\approx 10^2$~$L_\odot$ (or $a\approx 0.1$~au) for Jovian planetary companions. The consumption of closer-in giant planets, such as hot Jupiters, leads to $L_{\rm decay}\gg L_*$, while more distant planets such as warm Jupiters, $a\approx 0.5$~au, lead to minor perturbations of their host stars with $L_{\rm decay} \ll L_*$. Our results map out the parameter space along the giant branch in the H-R Diagram where interaction with planetary companions leads to significant energetic disturbance of host stars.

The reduction of the \emph{K2}'s Campaign 9 (K2C9) microlensing data is challenging mostly because of the very crowded field and the unstable pointing of the spacecraft. In this work, we present the first method that can extract microlensing signals from this K2C9 data product. The raw light curves and the astrometric solutions are first derived, using the techniques from Soares-Furtado et al. and Huang et al. for \emph{K2} dense field photometry. We then minimize and remove the systematic effect by performing simultaneous modeling with the microlensing signal. We also derive precise $(K_p-I)$ vs. $(V-I)$ color-color relations that can predict the microlensing source flux in the \emph{Kepler} bandpass. By implementing the color-color relation in the light curve modeling, we show that the microlensing parameters can be better constrained. In the end, we use two example microlensing events, OGLE-2016-BLG-0980 and OGLE-2016-BLG-0940, to test our method.

Observations were made of the open clusters M35 and NGC 2158 during the initial K2 campaign (C0). Reducing these data to high-precision photometric time-series is challenging due to the wide point spread function (PSF) and the blending of stellar light in such dense regions. We developed an image-subtraction-based K2 reduction pipeline that is applicable to both crowded and sparse stellar fields. We applied our pipeline to the data-rich C0 K2 super-stamp, containing the two open clusters, as well as to the neighboring postage stamps. In this paper, we present our image subtraction reduction pipeline and demonstrate that this technique achieves ultra-high photometric precision for sources in the C0 super-stamp. We extract the raw light curves of 3960 stars taken from the UCAC4 and EPIC catalogs and de-trend them for systematic effects. We compare our photometric results with the prior reductions published in the literature. For detrended, TFA-corrected sources in the 12--12.25 $\rm K_{p}$ magnitude range, we achieve a best 6.5 hour window running rms of 35 ppm falling to 100 ppm for fainter stars in the 14--14.25 $ \rm K_{p}$ magnitude range. For stars with $\rm K_{p}> 14$, our detrended and 6.5 hour binned light curves achieve the highest photometric precision. Moreover, all our TFA-corrected sources have higher precision on all time scales investigated. This work represents the first published image subtraction analysis of a K2 super-stamp. This method will be particularly useful for analyzing the Galactic bulge observations carried out during K2 campaign 9. The raw light curves and the final results of our detrending processes are publicly available at \url{http://k2.hatsurveys.org/archive/}.

We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHE; E > 100 GeV) by VERITAS with a significance of $8.9\sigma$ and showed clear variability on time scales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of $3.6 \pm 0.4_{\mathrm{stat}} \pm 0.3_{\mathrm{syst}}$ with an integral flux above 200 GeV of $(8.0 \pm 0.9_{\mathrm{stat}} \pm 3.2_{\mathrm{syst}}) \times 10^{-12}\, \mathrm{cm}^{-2} \mathrm{s}^{-1}$. No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneous with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift-UVOT), X-ray (Swift-XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be well reproduced with model parameters typical for VHE-detected BL Lacs.