Search

We present new stellar evolutionary sequences of very metal-rich stars evolved with the Monash Stellar Structure code and with MESA. The Monash models include masses of $1-8M_{\odot}$ with metallicities $Z=0.04$ to $Z=0.1$ and are evolved from the main sequence to the thermally-pulsing asymptotic giant branch (AGB). These are the first $Z=0.1$ AGB models in the literature. The MESA models include intermediate-mass models with $Z=0.06$ to $Z=0.09$ evolved to the onset of the thermally-pulsing phase. Third dredge-up only occurs in intermediate-mass models $Z \le 0.08$. Hot bottom burning (HBB) shows a weaker dependence on metallicity, with the minimum mass increasing from 4.5$M_{\odot}$ for $Z=0.014$ to $\approx 5.5 M_{\odot}$ for Z = 0.04, $6M_{\odot}$ for $ 0.05 \le Z \le 0.07$ and above 6.5$ M_{\odot}$ for $Z\ge 0.08$. The behaviour of the $Z=0.1$ models is unusual; most do not experience He-shell instabilities owing to rapid mass-loss on the early part of the AGB. Turning off mass-loss produces He-shell instabilities, however thermal pulses are weak and result in no third dredge-up. The minimum mass for carbon ignition is reduced from 8$M_{\odot}$ for $Z=0.04$ to 7$M_{\odot}$ for $Z=0.1$, which implies a reduction in the minimum mass for core-collapse supernovae. MESA models of similarly high metallicity ($Z=0.06 - 0.09$) show the same lowering of the minimum mass for carbon ignition: carbon burning is detected in a $6 M_{\odot}$ model at the highest metallicity ($Z=0.09$) and in all $7 M_{\odot}$ models with $Z \ge 0.06$. This demonstrates robustness of the lowered carbon burning threshold across codes., 18 pages, 10 figures, accepted for publication in MNRAS

We present a supervised machine learning methodology to classify stellar populations in the Local Group dwarf-irregular galaxy NGC6822. Near-IR colours (J-H, H-K, and J-K), K-band magnitudes and far-IR surface brightness (at 70 and 160 micron) measured from Spitzer and Herschel images are the features used to train a Probabilistic Random Forest (PRF) classifier. Point-sources are classified into eight target classes: young stellar objects (YSOs), oxygen- and carbon-rich asymptotic giant branch stars, red giant branch and red super-giant stars, active galactic nuclei, massive main-sequence stars and Galactic foreground stars. The PRF identifies sources with an accuracy of ~90 percent across all target classes rising to ~96 percent for YSOs. We confirm the nature of 125 out of 277 literature YSO candidates with sufficient feature information, and identify 199 new YSOs and candidates. Whilst these are mostly located in known star forming regions, we have also identified new star formation sites. These YSOs have mass estimates between ~15-50 Msun, representing the most massive YSO population in NGC 6822. Another 82 out of 277 literature candidates are definitively classified as non-YSOs by the PRF analysis. We characterise the star formation environment by comparing the spatial distribution of YSOs to those of gas and dust using archival images. We also explore the potential of using (unsupervised) t-distributed stochastic neighbour embedding maps for the identification of the same stellar population classified by the PRF., Comment: 20 pages, 16 figures. Accepted for publication in MNRAS

We study the evolution of six exoplanetary systems with the stellar evolutionary code MESA and conclude that they will likely spin-up the envelope of their parent stars on the red giant branch (RGB) or later on the asymptotic giant branch (AGB) to the degree that the mass loss process might become non-spherical. We choose six observed exoplanetary systems where the semi-major axis is ~1-2AU, and use the binary mode of MESA to follow the evolution of the systems. In four systems the star engulfs the planet on the RGB, and in two systems on the AGB, and the systems enter a common envelope evolution (CEE). In two systems where the exoplanet masses are Mp~10MJ, where MJ is Jupiter mass, the planet spins-up the envelope to about 10% of the break-up velocity. Such envelopes are likely to have significant non-spherical mass loss geometry. In the other four systems where Mp~MJ the planet spins-up the envelope to values of ~1-2% of break-up velocity. Magnetic activity in the envelope that influences dust formation might lead to a small departure from spherical mass loss even in these cases. In the two cases of CEE on the AGB the planet deposits energy to the envelope that amounts to >10% of the envelope binding energy. We expect this to cause a non-spherical mass loss that will shape an elliptical planetary nebula in each case., Accepted for publication in MNRAS

We present the abundance analyses of 7 Carbon enhanced metal-poor (CEMP) stars to understand the origin of carbon in them. We used high-resolution optical spectra to derive abundances of various elements. We also used low-resolution Near-Infrared (NIR) spectra to derive the abundance of O and 12C/13C from the CO molecular band and compared their values with those derived from high-resolution optical spectra. We identified a good agreement between the values. Thus, in cool CEMP stars, the NIR observations complement the high-resolution optical observations to derive the oxygen abundance and the 12C/13C ratio. This enables us to probe fainter cool CEMP stars using NIR spectroscopy. C, N, O abundances of all the program stars in this study show abundances that are consistent with binary mass transfer from a low-mass low-metallicity Asymptotic Giant Branch (AGB) companion which is further supported by the presence of enhancement in neutron-capture elements and detection of radial velocity variation. One of the stars show abundance patterns similar to a CEMP-s star whereas the abundance pattern of the rest of the stars satisfy the criteria required to classify them as CEMP-r/s stars. The sub-classification of some of the stars studied here is revisited. The abundance of neutron capture elements in these CEMP-r/s stars resembles to that of i-process models where proton ingestion episodes in the companion low-mass low-metallicity AGB stars produce the necessary neutron density required for the onset of i-process., 18 pages, 14 figures, 5 tables, accepted for publication in MNRAS

We use a combined binary evolution code including dynamical effects to study nova eruptions in a symbiotic system. Following the evolution, over $\sim10^5$ years, of multiple consecutive nova eruptions on the surface of a $1.25M_\odot$ white dwarf (WD) accretor, we present a comparison between simulations of two types of systems. The first is the common, well known, cataclysmic variable (CV) system in which a main sequence donor star transfers mass to its WD companion via Roche-lobe overflow. The second is a detached, widely separated, symbiotic system in which an asymptotic giant branch donor star transfers mass to its WD companion via strong winds. For the latter we use the Bondi-Hoyle-Lyttleton prescription along with orbital dynamics to calculate the accretion rate. We use the combined stellar evolution code to follow the nova eruptions of both simulations including changes in mass, accretion rate and orbital features. We find that while the average accretion rate for the CV remains fairly constant, the symbiotic system experiences distinct epochs of high and low accretion rates. The examination of epochs for which the accretion rates of both simulations are similar, shows that the evolutionary behaviors are identical. We obtain that for a given WD mass, the rate that mass is accreted ultimately determines the development, and that the stellar class of the donor is of no significance to the development of novae. We discuss several observed systems and find that our results are consistent with estimated parameters of novae in widely separated symbiotic systems., Comment: Accepted for publication in MNRAS

Only a small number of exoplanets has been identified in stellar cluster environments. We initiated a high angular resolution direct imaging search using the Hubble Space Telescope (HST) and its NICMOS instrument for self-luminous giant planets in orbit around seven white dwarfs in the 625 Myr old nearby (45 pc) Hyades cluster. The observations were obtained with NIC1 in the F110W and F160W filters, and encompass two HST roll angles to facilitate angular differential imaging. The difference images were searched for companion candidates, and radially averaged contrast curves were computed. Though we achieve the lowest mass detection limits yet for angular separations >0.5 arcsec, no planetary mass companion to any of the seven white dwarfs, whose initial main sequence masses were >2.8 Msun, was found. Comparison with evolutionary models yields detection limits of 5 to 7 Jupiter masses according to one model, and between 9 and 12 Mjup according to another model, at physical separations corresponding to initial semimajor axis of >5 to 8 A.U. (i.e., before the mass loss events associated with the red and asymptotic giant branch phase of the host star). The study provides further evidence that initially dense cluster environments, which included O- and B-type stars, might not be highly conducive to the formation of massive circumstellar disks, and their transformation into giant planets (with m>6 Mjup and a>6 A.U.). This is in agreement with radial velocity surveys for exoplanets around G- and K-type giants, which did not find any planets around stars more massive than about 3 Msun., 6 pages, 3 figures; https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/staa3422/5956540

Stellar winds of cool carbon stars enrich the interstellar medium with significant amounts of carbon and dust. We present a study of the influence of two-fluid flow on winds where we add descriptions of frequency-dependent radiative transfer. Our radiation hydrodynamic models in addition include stellar pulsations, grain growth and ablation, gas-to-dust drift using one mean grain size, dust extinction based on both the small particle limit and Mie scattering, and an accurate numerical scheme. We calculate models at high spatial resolution using 1024 gridpoints and solar metallicities at 319 frequencies, and we discern effects of drift by comparing drift models to non-drift models. Our results show differences of up to 1000 per cent in comparison to extant results. Mass-loss rates and wind velocities of drift models are typically, but not always, lower than in non-drift models. Differences are larger when Mie scattering is used instead of the small particle limit. Amongst other properties, the mass-loss rates of the gas and dust, dust-to-gas density ratio, and wind velocity show an exponential dependence on the dust-to-gas speed ratio. Yields of dust in the least massive winds increase by a factor four when drift is used. We find drift velocities in the range 10-67 km/s, which is drastically higher than in our earlier works that use grey radiative transfer. It is necessary to include an estimate of drift velocities to reproduce high yields of dust and low wind velocities., 32 pages, 19 figures, accepted for publication in MNRAS

We present extensive observations of SN 2018zd covering the first $\sim450$\,d after the explosion. This SN shows a possible shock-breakout signal $\sim3.6$\,hr after the explosion in the unfiltered light curve, and prominent flash-ionisation spectral features within the first week. The unusual photospheric temperature rise (rapidly from $\sim 12,000$\,K to above 18,000\,K) within the earliest few days suggests that the ejecta were continuously heated. Both the significant temperature rise and the flash spectral features can be explained with the interaction of the SN ejecta with the massive stellar wind ($0.18^{+0.05}_{-0.10}\, \rm M_{\odot}$), which accounts for the luminous peak ($L_{\rm max} = [1.36\pm 0.63] \times 10^{43}\, \rm erg\,s^{-1}$) of SN 2018zd. The luminous peak and low expansion velocity ($v \approx 3300$ km s$^{-1}$) make SN 2018zd to be like a member of the LLEV (luminous SNe II with low expansion velocities) events originated due to circumstellar interaction. The relatively fast post-peak decline allows a classification of SN 2018zd as a transition event morphologically linking SNe~IIP and SNe~IIL. In the radioactive-decay phase, SN 2018zd experienced a significant flux drop and behaved more like a low-luminosity SN~IIP both spectroscopically and photometrically. This contrast indicates that circumstellar interaction plays a vital role in modifying the observed light curves of SNe~II. Comparing nebular-phase spectra with model predictions suggests that SN 2018zd arose from a star of $\sim 12\,\rm M_{\odot}$. Given the relatively small amount of $^{56}$Ni ($0.013 - 0.035 \rm M_{\odot}$), the massive stellar wind, and the faint X-ray radiation, the progenitor of SN 2018zd could be a massive asymptotic giant branch star which collapsed owing to electron capture., Accepted for publication in MNRAS, 20 pages, 11 figures

We investigate the impact of bursts in star formation on the predictions of one-zone chemical evolution models, adopting oxygen (O), iron (Fe), and strontium (Sr), as representative $\alpha$, iron-peak, and s-process elements, respectively. To this end, we develop the Versatile Integrator for Chemical Evolution (VICE), a python package. Starbursts driven by a temporary boost of gas accretion rate create loops in [O/Fe]-[Fe/H] evolutionary tracks and a peak in the stellar [O/Fe] distribution at intermediate values. Bursts driven by a temporary boost of star formation efficiency have a similar effect, and they also produce a population of $\alpha$-deficient stars during the depressed star formation phase that follows the burst. This $\alpha$-deficient population is more prominent if the outflow rate is tied to a time-averaged star formation rate (SFR) instead of the instantaneous SFR. Theoretical models of Sr production predict a strong metallicity dependence of supernova and asymptotic giant branch (AGB) star yields, though comparison to data suggests an additional source that is nearly metallicity-independent. Evolution of [Sr/Fe] and [Sr/O] during a starburst is complex because of the yield metallicity dependence and the multiple timescales in play. Moderate amplitude (10-20\%) sinusoidal oscillations in SFR produce loops in [O/Fe]-[Fe/H] tracks and multiple peaks in [O/Fe] distributions, which could be one source of intrinsic scatter in observed sequences. We investigate models that have a factor of ~2 enhancement of SFR at t = 12 Gyr, as suggest by some recent Milky Way observations. A late episode of enhanced star formation could help explain the existence of young stars with moderate $\alpha$-enhancements and the surprisingly young median age found for solar metallicity stars in the solar neighborhood, while also raising the possibility that this starburst has not fully decayed., Comment: 17 pages; 11 figures; submitted to MNRAS; comments welcome

We study the dust evolution in galaxies by implementing a detailed dust prescription in the SAGE semi-analytical model for galaxy formation. The new model, called Dusty SAGE, follows the condensation of dust in the ejecta of type II supernovae and asymptotic giant branch (AGB) stars, grain growth in the dense molecular clouds, destruction by supernovae shocks, and the removal of dust from the ISM by star formation, reheating, inflows and outflows. Our model successfully reproduces the observed dust mass function at redshift z = 0 and the observed scaling relations for dust across a wide range of redshifts. We find that the dust mass content in the present Universe is mainly produced via grain growth in the interstellar medium (ISM). By contrast, in the early Universe, the primary production mechanism for dust is the condensation in stellar ejecta. The shift of the significant production channel for dust characterises the scaling relations of dust-to-gas (DTG) and dust-to-metal (DTM) ratios. In galaxies where the grain growth dominates, we find positive correlations for DTG and DTM ratios with both metallicity and stellar mass. On the other hand, in galaxies where dust is produced primarily via condensation, we find negative or no correlation for DTM and DTG ratios with either metallicity or stellar mass. In agreement with observation showing that the circumgalactic medium (CGM) contains more dust than the ISM, our model also shows the same trend for z < 4. Our semi-analytic model is publicly available at https: //github.com/dptriani/dusty-sage., 19 pages, 14 figures, accepted for publication in MNRAS

We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the $Gaia$ DR2 catalogue, the Sloan Digital Sky Survey (SDSS) and the $Galaxy~Evolution~Explorer$ (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20000K to 9000K we find that 22% of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5$\sigma$. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs - which include the full overshoot region - lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60% of white dwarfs must have a hydrogen mass larger than $M_{\rm H}/M_{\rm WD} = 10^{-10}$, another 25% have masses in the range $M_{\rm H}/M_{\rm WD} = 10^{-14}-10^{-10}$, and 15% have less hydrogen than $M_{\rm H}/M_{\rm WD} = 10^{-14}$. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch (AGB) and accretion of planetesimals onto white dwarfs., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

Barium stars are one of the important probes to understand the origin and evolution of slow neutron-capture process elements in the Galaxy. These are extrinsic stars, where the observed s-process element abundances are believed to have an origin in the now invisible companions that produced these elements at their asymptotic giant branch (AGB) phase of evolution. We have attempted to understand the s-process nucleosynthesis, as well as the physical properties of the companion stars through a detailed comparison of observed elemental abundances of 10 barium stars with the predictions from AGB nucleosynthesis models, FRUITY. For these stars, we have presented estimates of abundances of several elements, C, N, O, Na, Al, α-elements, Fe-peak elements, and neutron-capture elements Rb, Sr, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, and Eu. The abundance estimates are based on high resolution spectral analysis. Observations of Rb in four of these stars have allowed us to put a limit to the mass of the companion AGB stars. Our analysis clearly shows that the former companions responsible for the surface abundance peculiarities of these stars are low-mass AGB stars. Kinematic analysis has shown the stars to be members of Galactic disc population.

We performed detailed spectroscopic analyses of a young C-rich planetary nebula (PN) Jonckheere900 (J900) in order to characterise the properties of the central star and nebula. Of the derived 17 elemental abundances, we present the first determination of eight elemental abundances. We present the first detection of the [F IV] 4059.9 A, [F V] 13.4 um, and [Rb IV] 5759.6 A lines in J900. J900 exhibits a large enhancement of F and neutron-capture elements Se, Kr, Rb, and Xe. We investigated the physical conditions of the H2 zone using the newly detected mid-IR H2 lines while also using the the previously measured near-IR H2 lines, which indicate warm (~670 K) and hot (~3200 K) temperature regions. We built the spectral energy distribution (SED) model to be consistent with all the observed quantities. We found that about 67 % of all dust and gas components (4.5x10^-4 Msun and 0.83 Msun, respectively) exists beyond the ionisation front, indicating critical importance of photodissociation regions in understanding stellar mass loss. The best-fitting SED model indicates that the progenitor evolved from an initially ~2.0 Msun star which had been in the course of the He-burning shell phase. Indeed, the derived elemental abundance pattern is consistent with that predicted by a asymptotic giant branch star nucleosynthesis model for a 2.0 Msun star with Z = 0.003 and partial mixing zone mass of 6.0x10^-3 Msun. Our study demonstrates how accurately determined abundances of C/F/Ne/neutron-capture elements and gas/dust masses help us understand the origin and the internal evolution of the PN progenitors., Comment: 22 pages, 11 figures, accepted for publication by MNRAS. The title is corrected for accordance with the paper published by MNRAS

Many possible scenarios for the formation of multiple stellar populations (MSP) in globular clusters (GCs) have been discussed so far, including the involvement of asymptotic giant branch stars, fast rotating main sequence stars, very massive main sequence stars and mass-transferring massive binaries based on stellar evolution modelling. But self-consistent, dynamical simulations of very young GCs are usually not considered. In this work, we perform direct $N$-body modelling such systems with total masses up to $3.2\times10^5$ M$_\odot$, taking into account the observationally constrained primordial binary properties, and discuss the stellar-mergers driven both by binary stellar evolution and dynamical evolution of GCs. The occurrence of stellar mergers is enhanced significantly in binary-rich clusters such that stars forming from the gas polluted by mergers-driven ejection/winds would appears as MSPs. We thus emphasize that stellar mergers can be an important process that connects MSP formation with star cluster dynamics, and that multiple MSP formation channels can naturally work together. The scenario studied here, also in view of a possible top-heavy IMF, may be particularly relevant for explaining the high mass fraction of MSPs (the mass budget problem) and the absence of MSPs in young and low-mass star clusters., Comment: 16 pages, 11 figures, accepted for MNRAS

Using ALMA observations of $^{12}$CO(2-1), $^{28}$SiO(5-4) and $^{32}$SO$_2$(16$_{6,10}$-17$_{5,13}$) emissions of the circumstellar envelope of AGB star EP Aqr, we describe the morpho-kinematics governing the nascent wind. Main results are: 1) Two narrow polar structures, referred to as jets, launched from less than 25 au away from the star, build up between $\sim$ 20 au and $\sim$ 100 au to a velocity of $\sim$ 20 \kms. They fade away at larger distances and are barely visible in CO data. 2) SO$_2$, SiO and CO emissions explore radial ranges reaching respectively $\sim$30 au, 250 au and 1000 au from the star, preventing the jets to be detected in SO$_2$ data. 3) Close to the star photosphere, rotation (undetected in SiO and CO data) and isotropic radial expansion combine with probable turbulence to produce a broad SO$_2$ line profile ($\sim$ 7.5 \kms\ FWHM). 4) A same axis serves as axis of rotation close to the star, as jet axis and as axi-symmetry axis at large distances. 5) A radial wind builds up at distances up to $\sim$ 300 au from the star, with larger velocity near polar than equatorial latitudes. 6) A sharp depletion of SiO and CO emissions, starting near the star, rapidly broadens to cover the whole blue-western quadrant, introducing important asymmetry in the CO and particularly SiO observations. 7) The $^{12}$C/$^{13}$C abundance ratio is measured as 9$\pm$2. 8) Plausible interpretations are discussed, in particular assuming the presence of a companion., 18 pages, 16 figures, MNRAS accepted

The material lost through stellar winds of Asymptotic Giant Branch (AGB) stars is one of the main contributors to the chemical enrichment of galaxies. The general hypothesis of the mass loss mechanism of AGB winds is a combination of stellar pulsations and radiative pressure on dust grains, yet current models still suffer from limitations. Among others, they assume chemical equilibrium of the gas, which may not be justified due to rapid local dynamical changes in the wind. This is important as it is the chemical composition that regulates the thermal structure of the wind, the creation of dust grains in the wind, and ultimately the mass loss by the wind. Using a self-consistent hydrochemical model, we investigated how non-equilibrium chemistry affects the dynamics of the wind. This paper compares a hydrodynamical and a hydrochemical dust-free wind, with focus on the chemical heating and cooling processes. No sustainable wind arises in a purely hydrodynamical model with physically reasonable pulsations. Moreover, temperatures are too high for dust formation to happen, rendering radiative pressure on grains impossible. A hydrochemical wind is even harder to initiate due to efficient chemical cooling. However, temperatures are sufficiently low in dense regions for dust formation to take place. These regions occur close to the star, which is needed for radiation pressure on dust to sufficiently aid in creating a wind. Extending this model self-consistently with dust formation and evolution, and including radiation pressure, will help to understand the mass loss by AGB winds., Accepted for publication in MNRAS. 30 pages (incl. Appendix), 19 figures

ALMA observations of CO(1-0) and CO(2-1) emissions of the circumstellar envelope of EP Aqr, an oxygen-rich AGB star, are reported. A thorough analysis of their properties is presented using an original method based on the separation of the data-cube into a low velocity component associated with an equatorial outflow and a faster component associated with a bipolar outflow. A number of important and new results are obtained concerning the distribution in space of the effective emissivity, the temperature, the density and the flux of matter. A mass loss rate of (1.6$\pm$0.4)10$^{-7}$ solar masses per year is measured. The main parameters defining the morphology and kinematics of the envelope are evaluated and uncertainties inherent to de-projection are critically discussed. Detailed properties of the equatorial region of the envelope are presented including a measurement of the line width and a precise description of the observed inhomogeneity of both morphology and kinematics. In particular, in addition to the presence of a previously observed spiral enhancement of the morphology at very small Doppler velocities, a similarly significant but uncorrelated circular enhancement of the expansion velocity is revealed, both close to the limit of sensitivity. The results of the analysis place significant constraints on the parameters of models proposing descriptions of the mass loss mechanism, but cannot choose among them with confidence., Comment: 26 pages, 31 figures, accepted for publication in MNRAS

We present a near-infrared imaging polarimetric study of the pre-planetary nebula: Frosty Leo. The observations were carried out in J, H and K' bands using the new polarimeter POLICAN mounted on the 2.1m telescope of the Guillermo Haro Astrophysical Observatory, Sonora, Mexico. The most prominent result observed in the polarization maps is a large and well defined dusty envelope (35\arcsec\ diameter in H-band). The polarization position angles in the envelope are particularly well ordered and nearly parallel to the equator of the nebula (seen in J and H bands). The nebula presents a known bipolar outflow and the envelope completely wraps around it. Within the bipolar lobes, we find high polarization levels ranging from $60\%$ (J band) to $90\%$ (K' band) and the polarization angles trace a centrosymmetric pattern. We found the remnants of superwind shells at the edges of the bipolar lobes and the duration of this phase is around 600 yrs. The origin of polarization features in the nebula is most likely due to a combination of single and multiple scattering. Our results clearly demonstrate new structures that provide new hints on the evolution of Frosty Leo from its previous asymptotic giant branch phase., Comment: 5 figures, 9 pages, accepted in MNRAS

We report the detection of variable stars within a 11.5' x 11.5' region near the Galactic centre (GC) that includes the Arches and Quintuplet clusters, as revealed by the VISTA Variables in the Via Lactea (VVV) survey. There are 353 sources that show Ks-band variability, of which the large majority (81%) correspond to red giant stars, mostly in the asymptotic giant branch (AGB) phase. We analyze a population of 52 red giants with long-term trends that cannot be classified into the typical pulsating star categories. Distances and extinctions are calculated for 9 Mira variables, and we discuss the impact of the chosen extinction law on the derived distances. We also report the presence of 48 new identified young stellar object (YSO) candidates in the region., Comment: Resubmitted to MNRAS after a positive referee report, 21 pages, 23 figures and 3 tables

We present a new grid of non-adiabatic, linear pulsation models of Long-Period Variables (LPVs), including periods and growth rates for radial modes from the fundamental to the fourth overtone. The models span a wide range in mass, luminosity, metallicity, C/O ratio and helium abundance, effectively covering the whole thermally-pulsing asymptotic giant branch (TP-AGB) evolution, and representing a significant update with respect to previous works. The main improvement is the inclusion of detailed atomic and molecular opacities, consistent with the models chemical mixture, that makes the present set of models the first to systematically account for variability in C-stars. We examine periods and growth rates in the models, and find that, while the fundamental mode is affected by the structure of the envelope, overtones are less sensitive to the interior and largely determined by the global properties. In the models, the frequency of the overtone with the largest degree of excitation is found to scale with the acoustic cut-off frequency at the stellar surface, a behaviour similar to that observed for the frequency of maximum oscillation power for solar-like oscillations in less evolved red giants. This allows us to provide a simple analytic prescription to predict the most-likely dominant mode as a function of stellar parameters. Best-fit relations for periods are also provided. By applying results of pulsation models to evolutionary tracks, we present a general picture of the evolution of long-period variability during the TP-AGB, that we find consistent with observations. Models are made public through a dedicated web interface., 23 pages, 24 figures, accepted for publication in MNRAS --- web interface service under construction

In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium., Comment: 21 pages, 13 figures (plus 6 page, 7 figure appendix). Accepted to MNRAS

We present a study of the high-resolution optical spectrum for the hot post-asymptotic giant branch (post-AGB) star, Hen 3-1013, identified as the optical counterpart of the infrared sourceIRAS 14331-6435. For the first time the detailed identifications of the observed absorption and emission features in the wavelength range 3700-9000 \AA\ is carried out. Absorption lines of HI, HeI, CI, NI, OI, NeI, CII, NII, OII, SiII, SII, ArII, FeII, MnII, CrII, TiII, CoII, NiII, SIII, FeIII and SIV were detected. From the absorption lines, we derived heliocentric radial velocities of $V_r=-29.6\pm0.4$ km/s. We have identified emission permitted lines of OI, NI, FeII, MgII, SiII and AlII. The forbidden lines of [NI], [FeII], [CrII] and [NiII] have been identified also. Analysis of [NiII] lines in the gaseous shell gives an estimate for the electron density $N_e\sim10^7$ cm$^{-3}$ and the expansion velocity of the nebula $V_{exp}=12$ km/s. The mean radial velocity as measured from emission features of the envelope is $V_r=-36.0\pm0.4$ km/s. The Balmer lines H$\alpha$, H$\beta$ and H$\gamma$ show P Cyg behaviour which indicate ongoing post-AGB mass-loss. Based on ASAS and ASAS-SN data, we have detected rapid photometric variability in Hen 3-1013 with an amplitude up to 0.2 mag in the V band. The star's low-resolution spectrum underwent no significant changes from 1994 to 2012. Based on archival data, we have traced the photometric history of the star over more than 100 years. No significant changes in the star brightness have been found., Comment: 19 pages, 6 figures, 8 tables, accepted for publication in MNRAS

We present the results of an analysis aimed at testing the accuracy and precision of the PARSEC v1.2S library of stellar evolution models, in a Bayesian framework, to infer stellar parameters. We mainly employ the online DEBCat catalogue by Southworth, a compilation of detached eclipsing binary systems with published measurements of masses and radii to $\sim$ 2 per cent precision. We select a sample of 318 binary components, with masses between 0.10 and 14.5 Msun, at distances between 1.3 pc and ~ 8 kpc for Galactic objects and ~ 44-68 kpc for extragalactic ones. The Bayesian analysis applied takes as input effective temperature, radius, and [Fe/H], and their uncertainties, returning theoretical predictions for other stellar parameters. From the comparison with dynamical masses, we conclude that the inferred masses are precisely derived for stars on the main-sequence and in the core-helium-burning phase, with uncertainties of 4 per cent and 7 per cent, respectively, on average. Masses for subgiants and red giants are predicted within 14 per cent, and those for early asymptotic giant branch stars within 24 per cent. These results are helpful to further improve the models, in particular for advanced evolutionary stages for which our understanding is limited. We obtain distances and ages for the binary systems and compare them, whenever possible, with precise literature estimates, finding excellent agreement. We discuss evolutionary effects and challenges for inferring stellar ages from evolutionary models. We also provide useful polynomial fittings to theoretical zero-age main-sequence relationships., Comment: 25 pages, 26 figures

We calculate the outward energy transport time by convection and photon diffusion in an inflated common envelope and find this time to be shorter than the envelope expansion time. We conclude therefore that most of the hydrogen recombination energy ends in radiation rather than in kinetic energy of the outflowing envelope. We use the stellar evolution code MESA and inject energy inside the envelope of an asymptotic giant branch star to mimic energy deposition by a spiraling-in stellar companion. During 1.7 years the envelope expands by a factor of more than 2. Along the entire evolution the convection can carry the energy very efficiently outwards, to the radius where radiative transfer becomes more efficient. The total energy transport time stays within several months, shorter than the dynamical time of the envelope. Had we included rapid mass loss, as is expected in the common envelope evolution, the energy transport time would have been even shorter. It seems that calculations that assume that most of the recombination energy ends in the outflowing gas might be inaccurate., Comment: To be published in MNRAS

By following the evolution of several observed exoplanetary systems we show that by lowering the mass loss rate of single solar-like stars during their two giant branches, these stars will swallow their planets at the tip of their asymptotic giant branch (AGB) phase. This will most likely lead the stars to form elliptical planetary nebulae (PNe). Under the traditional mass loss rate these stars will hardly form observable PNe. Stars with a lower mass loss rate as we propose, about 15 per cent of the traditional mass loss rate of single stars, leave the AGB with much higher luminosities than what traditional evolution produces. Hence, the assumed lower mass loss rate might also account for the presence of bright PNe in old stellar populations. We present the evolution of four exoplanetary systems that represent stellar masses in the range of 0.9-1.3 Mo. The justification for this low mass loss rate is our assumption that the stellar samples that were used to derive the traditional average single-star mass loss rate were contaminated by stars that suffer binary interaction., accepted to MNRAS

We argue that outward transport of energy by convection and photon diffusion in a common envelope evolution (CEE) of giant stars substantially reduces the fraction of the recombination energy of hydrogen and helium that is available for envelope removal. We base our estimate on the properties of an unperturbed asymptotic giant branch (AGB) spherical model, and on some simple arguments. Since during the CEE the envelope expands and energy removal by photon diffusion becomes more efficient, our arguments underestimate the escape of recombination energy. We hence strengthen earlier claims that recombination energy does not contribute much to common envelope removal. A large fraction of the energy that jets deposit to the envelope, on the other hand, might be in the form of kinetic energy of the expanding and buoyantly rising hot bubbles. These rapidly rising bubbles remove mass from the envelope. We demonstrate this process by conducting a three-dimensional hydrodynamical simulation where we deposit hot gas in the location of a secondary star that orbits inside the envelope of a giant star. Despite the fact that we do not include the large amount of gravitational energy that is released by the in-spiraling secondary star, the hot bubbles alone remove mass at a rate of about 0.1 Mo/yr, which is much above the regular mass loss rate., accepted to MNRAS

We report an intrinsically precise relation of the barium-to-iron enrichment as a function of age for a local, volume-complete (N = 30) sample of ancient Population II (tau >= 12 Gyr) and intermediate-disc stars (tau similar or equal to 10 Gyr), which suggests a common, r-process-dominated nucleosynthesis site for both elements in the early stages of the Milky Way. Deviants from this empirical relation are to a large extent identified as formerly known or new blue straggler stars. We report in particular the striking case of the Population II star HD159062, whose barium overabundance is difficult to explain without wind accretion of s-process material from a former asymptotic giant branch (AGB) primary that very likely survived as a white dwarf companion. The weak but significant barium enhancement that we measure for HR3578 and 104 Tau also suggests that both may be accompanied by faint degenerate companions. If confirmed through precision astrometry or direct imaging observations, this would mean a very efficient method to uncover ancient stellar remnant companions around solar-type stars.

We present a detailed analysis of IC 4776, a planetary nebula displaying a morphology believed to be typical of central star binarity. The nebula is shown to comprise a compact hourglass-shaped central region and a pair of precessing jet-like structures. Time-resolved spectroscopy of its central star reveals periodic radial velocity variability consistent with a binary system. While the data are insufficient to accurately determine the parameters of the binary, the most likely solutions indicate that the secondary is probably a low-mass main sequence star. An empirical analysis of the chemical abundances in IC 4776 indicates that the common-envelope phase may have cut short the AGB evolution of the progenitor. Abundances calculated from recombination lines are found to be discrepant by a factor of approximately two relative to those calculated using collisionally excited lines, suggesting a possible correlation between low abundance discrepancy factors and intermediate-period post-common-envelope central stars and/or Wolf-Rayet central stars. The detection of a radial velocity variability associated with binarity in the central star of IC 4776 may be indicative of a significant population of (intermediate-period) post-common-envelope binary central stars which would be undetected by classic photometric monitoring techniques., Comment: Accepted for publication in MNRAS

Recent observations have revealed that at least several old globular clusters (GCs) in the Galaxy have discrete distributions of stars along the Mg-Al anti-correlation. In order to discuss this recent observation, we construct a new one-zone GC formation model in which the maximum stellar mass (m_max) in the initial mass function (IMF) of stars in a forming GC depends on the star formation rate (SFR), as deduced from independent observations. We investigate the star formation histories of forming GCs. The principal results are as follows. About 30 Myr after the formation of the first generation (1G) of stars within a particular GC, new stars can be formed from ejecta from asymptotic giant branch (AGB) stars of 1G. However, the formation of this second generation (2G) of stars can last only for [10-20] Myr, because the most massive SNe of 2G expel all of the remaining gas. The third generation (3G) of stars are then formed from AGB ejecta ~ 30$ Myr after the truncation of 2G star formation. This cycle of star formation followed by its truncation by SNe can continue until all AGB ejecta is removed from the GC by some physical process. Thus, it is inevitable that GCs have discrete multiple stellar populations in the [Mg/Fe]-[Al/Fe] diagram. Our model predicts that low-mass GCs are unlikely to have discrete multiple stellar populations, and young massive clusters may not have massive OB stars owing to low m_max ([20-30] M_sun) during the secondary star formation., 15 pages, 11 figures, accepted for publication in MNRAS

Gravitationally focused wind accretion in binary systems consisting of an evolved star with a gaseous envelope and a compact accreting companion is a possible mechanism to explain mass transfer in symbiotic binaries. We study the mass accretion around the secondary caused by the strong wind from the primary late-type component using global three-dimensional hydrodynamic numerical simulations during quiescence and outburst stages. In particular, the dependence of the mass accretion rate on the mass-loss rate, wind parameters and phases of wind outburst development is considered. For a typical wind from an asymptotic giant branch star with a mass-loss rate of 1e-6 Msun/year and wind speeds of 20-50 km/s, the mass transfer through a focused wind results in efficient infall on to the secondary. Accretion rates onto the secondary of 5-20 per cent of the mass-loss from the primary are obtained during quiescence and outburst periods where the wind velocity and mass-loss rates are varied, about 20-50 per cent larger than in the standard Bondi-Hoyle-Lyttleton approximation. This mechanism could be an important method for explaining observed accretion luminosities and periodic modulations in the accretion rates for a broad range of interacting binary systems., Comment: 11 pages, 8 figures, to be published in MNRAS

Dust formation and resulting mass loss around Asymptotic Giant Branch (AGB) stars with initial metallicity in the range of $0 \leq Z_{\rm ini} \leq 10^{-4}$ and initial mass $2\leq M_{\rm ini}/M_{\odot} \leq 5$ are explored by the hydrodynamical calculations of dust-driven wind (DDW) along the AGB evolutionary tracks. We employ the MESA code to simulate the evolution of stars, assuming an empirical mass-loss rate in the post-main sequence phase, and considering the three types of low-temperature opacities (scaled-solar, CO-enhanced, and CNO-enhanced opacities) to elucidate the effect on the stellar evolution and the DDW. We find that the treatment of low-temperature opacity strongly affects the dust formation and resulting DDW; in the carbon-rich AGB phase, the maximum $\dot{M}$ of $M_{\rm ini} \geq$ 3 $M_{\odot}$ star with the CO-enhanced opacity is at least one order of magnitude smaller than that with the CNO-enhanced opacity. A wide range of stellar parameters being covered, a necessary condition for driving efficient DDW with $\dot{M} \ge 10^{-6}$ $M_{\odot}$ yr$^{-1}$ is expressed as the effective temperature $T_{\rm eff} \lesssim 3850$ K and $\log(\delta_{\rm C}L/\kappa_{\rm R} M) \gtrsim 10.43\log T_{\rm eff}-32.33 $ with the carbon excess $\delta_{\rm C}$ defined as $\epsilon_{\rm C} - \epsilon_{\rm O}$ and the Rosseland mean opacity $\kappa_{\rm R}$ in units of cm$^2$g$^{-1}$ in the surface layer, and the stellar mass (luminosity) $M$ $(L)$ in solar units. The derived fitting formulae of gas and dust mass-loss rates in terms of input stellar parameters could be useful for investigating the dust yield from AGB stars in the early Universe being consistent with the stellar evolution calculations., Comment: 26 pages, 7 figures, 4 tables, accepted for publication in MNRAS

We present a joint analysis of the four most prominent sodium-sensitive features (NaD, NaI8190, NaI1.14, and NaI2.21), in the optical and Near-Infrared spectral range, of two nearby, massive (sigma~300km/s), early-type galaxies (named XSG1 and XSG2). Our analysis relies on deep VLT/X-Shooter long-slit spectra, along with newly developed stellar population models, allowing for [Na/Fe] variations, up to 1.2dex, over a wide range of age, total metallicity, and IMF slope. The new models show that the response of the Na-dependent spectral indices to [Na/Fe] is stronger when the IMF is bottom heavier. For the first time, we are able to match all four Na features in the central regions of massive early-type galaxies, finding an overabundance of [Na/Fe], in the range 0.5-0.7dex, and a bottom-heavy IMF. Therefore, individual abundance variations cannot be fully responsible for the trends of gravity-sensitive indices, strengthening the case towards a non-universal IMF. Given current limitations of theoretical atmosphere models, our [Na/Fe] estimates should be taken as upper limits. For XSG1, where line strengths are measured out to 0.8Re, the radial trend of [Na/Fe] is similar to [Mg/Fe] and [C/Fe], being constant out to 0.5Re, and decreasing by 0.2-0.3dex at 0.8Re, without any clear correlation with local metallicity. Such a result seems to be in contrast with the predicted increase of Na nucleosynthetic yields from AGB stars and TypeII SNe. For XSG1, the Na-inferred IMF radial profile is consistent, within the errors, with that derived from TiO features and the Wing-Ford band, presented in a recent paper., Comment: 22 pages, 8 figure, accepted for publication in MNRAS. The new Na-enhanced models will be available soon at http://miles.iac.es/

By means of grid-based, 3D hydrodynamical simulations we study the formation of second generation (SG) stars in a young globular cluster (GC) of mass 10^7 Msun, the possible progenitor of an old GC with a present mass ~(1-5) * 10^6 Msun. The cluster accretes external gas as its first generation (FG) asymptotic giant branch (AGB) stars release their ejecta and SG stars form. We consider two models characterised by different densities of the external gas. In both cases, we find that a very compact SG subsystem with central density >10^5M sun/pc^3 forms in the innermost regions of the cluster. The low-density model forms a population of extreme SG stars with high helium enhancement, followed by the formation of another SG group out of a mix of pristine gas and AGB ejecta and characterised by a modest helium enhancement. On the other hand, the high-density model forms in prevalence SG stars with modest helium enhancement. Our simulations illustrate the dynamical processes governing the formation of SG populations in GCs and shed light on the structural properties emerging at the end of this phase. The newly born SG groups have different concentrations, with more extreme SG stars more centrally concentrated than those with less extreme chemical abundances. The very high density of the SG subsystems implies that SG massive stars, if formed, might suffer frequent close encounters, collisions and gas stripping, thus possibly contributing further gas to the SG formation., MNRAS, accepted for publication. 17 pages, 7 figures

H2O (22 GHz) and SiO masers (43, 86, 129 GHz) in the bipolar proto-planetary nebula OH231.8+4.2 were simultaneously monitored using the 21-m antennas of the Korean VLBI Network in 2009-2015. Both species exhibit periodic flux variations that correlate with the central star's optical light curve, with a phase delay of up to 0.15 for the maser flux variations with respect to the optical light curve. The flux densities of SiO v = 2, J = 1-0 and H2O masers decrease with time, implying that they may disappear in 10-20 years. However, there seems to have been a transient episode of intense H2O maser emission around 2010. We also found a systematic behaviour in the velocity profiles of these masers. The velocities of the H2O maser components appear to be remarkably constant, suggesting ballistic motion for the bipolar outflow in this nebula. On the other hand, those of the SiO maser clumps show a systematic radial acceleration of the individual clumps, converging to the outflow velocity of the H2O maser clumps. Measuring the full widths at zero power of the detected lines, we estimated the expansion velocities of the compact bipolar outflow traced by H2O maser and SiO thermal line, and discussed the possibility of the expanding SiO maser region in the equatorial direction. All of our analyses support that the central host star of OH231.8 is close to the tip of the AGB phase, and that the mass-loss rate recently started to decrease because of incipient post-AGB evolution., Comment: 21 pages, 16 figures

The origin of the progenitors of type Ia supernovae (SNe Ia) is still uncertain. The core-degenerate (CD) scenario has been proposed as an alternative way for the production of SNe Ia. In this scenario, SNe Ia are formed at the final stage of common-envelope evolution from a merger of a carbon-oxygen white dwarf (CO WD) with the CO core of an asymptotic giant branch companion. However, the birthrates of SNe Ia from this scenario are still not well determined. In this work, we performed a detailed investigation on the CD scenario based on a binary population synthesis approach. The SN Ia delay times from this scenario are basically in the range of 90Myr-2500Myr, mainly contributing to the observed SNe Ia with short and intermediate delay times although this scenario can also produce some old SNe Ia. Meanwhile, our work indicates that the Galactic birthrates of SNe Ia from this scenario are no more than 20% of total SNe Ia due to more careful treatment of mass transfer. Although the SN Ia birthrates in the present work are lower than those in Ilkov & Soker, the CD scenario cannot be ruled out as a viable mechanism for the formation of SNe Ia. Especially, SNe Ia with circumstellar material from this scenario contribute to 0.7-10% of total SNe Ia, which means that the CD scenario can reproduce the observed birthrates of SNe Ia like PTF 11kx. We also found that SNe Ia happen systemically earlier for a high value of metallicity and their birthrates increase with metallicity., 13 pages, 8 figures, 1 table, accepted for publication in MNRAS

We have used the tidal equations of Zahn to determine the maximum orbital distance at which companions are brought into Roche lobe contact with their giant primary, when the primary expands during the giant phases. This is a key step when determining the rates of interaction between giants and their companions. Our stellar structure calculations are presented as maximum radii reached during the red and asymptotic giant branch (RGB and AGB, respectively) stages of evolution for masses between 0.8 and 4.0 Mo (Z=0.001 - 0.04) and compared with other models to gauge the uncertainty on radii deriving from details of these calculations. We find overall tidal capture distances that are typically 1-4 times the maximum radial extent of the giant star, where companions are in the mass range from 1 Jupiter mass to a mass slightly smaller than the mass of the primary. We find that only companions at initial orbital separations between ~320 and ~630 Ro will be typically captured into a Roche lobe-filling interaction or a common envelope on the AGB. Comparing these limits with the period distribution for binaries that will make PN, we deduce that in the standard scenario where all ~1-8 Mo stars make a PN, at most 2.5 per cent of all PN should have a post-common envelope central star binary, at odds with the observational lower limit of 15-20 per cent. The observed over-abundance of post-interaction central stars of PN cannot be easily explained considering the uncertainties. We examine a range of explanations for this discrepancy., Comment: 19 pages, 16 figures, accepted by Monthly Notices of the Royal Astronomical Society

Recent observations of the L2 Puppis system suggest that this Mira-like variable may be in the early stages of forming a bipolar planetary nebula (PN). As one of nearest and brightest AGB stars, thought be a binary, L2 Puppis serves as a benchmark object for studying the late-stages of stellar evolution. We perform global, three-dimensional, adaptive-mesh-refinement hydrodynamic simulations of the L2 Puppis system with AstroBEAR. We use the radiative transfer code RADMC-3D to construct the broad-band spectral-energy-distribution (SED) and synthetic observational images from our simulations. Given the reported binary parameters, we are able to reproduce the current observational data if a short pulse of dense material is released from the AGB star with a velocity sufficient to escape the primary but not the binary. Such a situation could result from a thermal pulse, be induced by a periastron passage of the secondary, or could be launched if the primary ingests a planet., 6 pages, 12 figures, submitted to MNRAS

We used models of thermally-pulsing asymptotic giant branch (AGB) stars, that also describe the dust-formation process in the wind, to interpret the combination of near- and mid-infrared photometric data of the dwarf galaxy IC 1613. This is the first time that this approach is extended to an environment different from the Milky Way and the Magellanic Clouds (MCs). Our analysis, based on synthetic population techniques, shows a nice agreement between the observations and the expected distribution of stars in the colour-magnitude diagrams obtained with JHK and Spitzer bands. This allows a characterization of the individual stars in the AGB sample in terms of mass, chemical composition, and formation epoch of the progenitors. We identify the stars exhibiting the largest degree of obscuration as carbon stars evolving through the final AGB phases, descending from 1-1.25Msun objects of metallicity Z=0.001 and from 1.5-2.5Msun stars with Z=0.002. Oxygen-rich stars constitute the majority of the sample (65%), mainly low mass stars (, 14 pages, 6 figures, accepted for publication on MNRAS

We analyse the planetary nebulae (PNe) population of the Small Magellanic Cloud (SMC), based on evolutionary models of stars with metallicities in the range $10^{-3} \leq Z \leq 4\times 10^{-3}$ and mass $0.9 M\odot < M < 8M\odot$, evolved through the asymptotic giant branch (AGB) phase. The models used account for dust formation in the circumstellar envelope. To characterise the PNe sample of the SMC, we compare the observed abundances of the various species with the final chemical composition of the AGB models: this study allows us to identify the progenitors of the PNe observed, in terms of mass and chemical composition. According to our interpretation, most of the PNe descend from low-mass ($M < 2 M\odot$) stars, which become carbon rich, after experiencing repeated third dredge-up episodes, during the AGB phase. A fraction of the PNe showing the signature of advanced CNO processing are interpreted as the progeny of massive AGB stars, with mass above $\sim 6 M\odot$, undergoing strong hot bottom burning. The differences with the chemical composition of the PNe population of the Large Magellanic Cloud (LMC) is explained on the basis of the diverse star formation history and age-metallicity relation of the two galaxies. The implications of the present study for some still highly debated points regarding the AGB evolution are also commented., Accepted for publication in MNRAS, 11 pages, 4 figures

We present the detection of a variable star with an $18.05$ minute period in far-ultraviolet (FUV) images of the globular cluster NGC 1851 taken with the Hubble Space Telescope (HST). A candidate optical counterpart lies on the red horizontal branch or the asymptotic giant branch star of the cluster, but it is statistically possible that this is a chance superposition. This interpretation is supported by optical spectroscopt obtained with HST/STIS: the spectrum contains none of the strong emission lines that would be expected if the object was a symbiotic star (i.e. a compact accretor fed by a giant donor). We therefore consider two other possibilities for the nature of FUV variable: (i) an intermediate polar (i.e. a compact binary containing an accreting magnetic white dwarf), or (ii) an AM CVn star (i.e. an interacting double-degenerate system). In the intermediate polar scenario, the object is expected to be an X-ray source. However, no X-rays are detected at its location in $\simeq 65$~ksec of {\em Chandra} imaging, which limits the X-ray luminosity to $L_X \leqslant 10^{32}$~erg~s$^{-1}$. We therefore favour the AM CVn interpretation, but a FUV spectrum is needed to distinguish conclusively between the two possibilities. If the object is an AM CVn binary, it would be the first such system known in any globular cluster., 10 pages, 10 figures, Accepted to MNRAS 2016 May 16

The possibility that mass ejected during Asymptotic Giant Branch (AGB) stellar evolution phases falls back towards the star has been suggested in applications ranging from the formation of accretion disks to the powering of late-thermal pulses. In this paper, we seek to explicate the properties of fallback flow trajectories from mass-loss events. We focus on a transient phase of mass ejection with sub-escape speeds, followed by a phase of a typical AGB wind. We solve the problem using both hydrodynamic simulations and a simplified one-dimensional analytic model that matches the simulations. For a given set of initial wind characteristics, we find a critical shell velocity that distinguishes between "shell fallback" and "shell escape". We discuss the relevance of our results for both single and binary AGB stars. In particular, we discuss how our results help to frame further studies of fallback as a mechanism for forming the substantial population of observed post-AGB stars with dusty disks., 7 pages, 13 figures, accepted by MNRAS on Jan. 4th 2016

We evolve stellar models to study the common envelope (CE) interaction of an early asymptotic giant branch star of initial mass $5\,\rm M_{\odot}$ with a companion star of mass ranging from $0.1$ to $2\,\rm M_{\odot}$. We model the CE as a fast stripping phase in which the primary experiences rapid mass loss and loses about 80 per cent of its mass. The post-CE remnant is then allowed to thermally readjust during a Roche-lobe overflow (RLOF) phase and the final binary system and its orbital period are investigated. We find that the post-CE RLOF phase is long enough to allow nuclear burning to proceed in the helium shell. By the end of this phase, the donor is stripped of both its hydrogen and helium and ends up as carbon-oxygen white dwarf of mass about $0.8\,\rm M_{\odot}$. We study the sensitivity of our results to initial conditions of different companion masses and orbital separations at which the stripping phase begins. We find that the companion mass affects the final binary separation and that helium-shell burning causes the star to refill its Roche lobe leading to post-CE RLOF. Our results show that double mass transfer in such a binary interaction is able to strip the helium and hydrogen layers from the donor star without the need for any special conditions or fine tuning of the binary parameters., Comment: Accepted to MNRAS 10 pages, 8 figures

We use cosmological, hydrodynamical simulations from the EAGLE and OWLS projects to assess the significance of recycled stellar ejecta as fuel for star formation. The fractional contributions of stellar mass loss to the cosmic star formation rate (SFR) and stellar mass densities increase with time, reaching $35 \%$ and $19 \%$, respectively, at $z=0$. The importance of recycling increases steeply with galaxy stellar mass for $M_{\ast} < 10^{10.5}$ M$_{\odot}$, and decreases mildly at higher mass. This trend arises from the mass dependence of feedback associated with star formation and AGN, which preferentially suppresses star formation fuelled by recycling. Recycling is more important for satellites than centrals and its contribution decreases with galactocentric radius. The relative contribution of AGB stars increases with time and towards galaxy centers. This is a consequence of the more gradual release of AGB ejecta compared to that of massive stars, and the preferential removal of the latter by star formation-driven outflows and by lock up in stellar remnants. Recycling-fuelled star formation exhibits a tight, positive correlation with galaxy metallicity, with a secondary dependence on the relative abundance of alpha elements (which are predominantly synthesized in massive stars), that is insensitive to the subgrid models for feedback. Hence, our conclusions are directly relevant for the origin of the mass-metallicity relation and metallicity gradients. Applying the relation between recycling and metallicity to the observed mass-metallicity relation yields our best estimate of the mass-dependent contribution of recycling. For centrals with a mass similar to that of the Milky Way, we infer the contributions of recycled stellar ejecta to the SFR and stellar mass to be $35 \%$ and $20 \%$, respectively., 20 pages, 12 figures (not including appendices), accepted for publication in MNRAS

The study of extragalactic integrated light can yield partial information on stellar population ages, abundances, and the initial mass function (IMF). The power-law slope of the IMF has been studied in recent investigations with gravity-sensitive spectral indicators that hopefully measure the ratio between KM dwarfs and giants. We explore two additional effects that might mimic the effects of the IMF slope in integrated light, the low mass cutoff (LMCO) and a variable contribution of light from the asymptotic giant branch (AGB). We show that the spectral effects of these three (IMF slope, LMCO, AGB strength) are subtle compared to age-abundance effects. We illustrate parameter degeneracies and covariances and conclude that the three effects can be disentangled, but only in the regime of very accurate observations, with enhanced effectiveness if high-precision photometry is combined with spectroscopy.

We present high spectral resolution mid-IR observations of SiS towards the C-rich AGB star IRC+10216 carried out with the Texas Echelon-cross-Echelle Spectrograph mounted on the NASA Infrared Telescope Facility. We have identified 204 ro-vibrational lines of 28Si32S, 26 of 29Si32S, 20 of 28Si34S, and 15 of 30Si32S in the frequency range 720–790 cm−1. These lines belong to bands v = 1–0, 2–1, 3–2, 4–3, and 5–4, and involve rotational levels with Jlow ≲ 90. About 30 per cent of these lines are unblended or weakly blended and can be partially or entirely fitted with a code developed to model the mid-IR emission of a spherically symmetric circumstellar envelope composed of expanding gas and dust. The observed lines trace the envelope at distances to the star ≲35R⋆ (≃0.7 arcsec). The fits are compatible with an expansion velocity of 1 + 2.5(r/R⋆ − 1) km s−1 between 1 and 5R⋆, 11 km s−1 between 5 and 20R⋆, and 14.5 km s−1 outwards. The derived abundance profile of 28Si32S with respect to H2 is 4.9 × 10−6 between the stellar photosphere and 5R⋆, decreasing linearly down to 1.6 × 10−6 at 20R⋆ and to 1.3 × 10−6 at 50R⋆. 28Si32S seems to be rotationally under local thermodynamic equilibrium (LTE) in the region of the envelope probed with our observations and vibrationally out of LTE in most of it. There is a red-shifted emission excess in the 28Si32S lines of band v = 1−0 that cannot be found in the lines of bands v = 2−1, 3–2, 4–3, and 5–4. This excess could be explained by an enhancement of the vibrational temperature around 20R⋆ behind the star. The derived isotopic ratios 28Si/29Si, and 32S/34S are 17 and 14, compatible with previous estimates., JPF was supported during part of this study by the UNAM through a postdoctoral fellowship. We thank the Spanish MINECO/MICINN for funding support through grants AYA2009-07304, AYA2012-32032, the ASTROMOL Consolider project CSD2009-00038 and the European Research Council (ERC Grant 610256: NANOCOSMOS). MJR is supported by grant AST 03-07497. Development of TEXES was supported by grants from the NSF and USRA. MJR, JHL, and others want to thank IRTF, which is operated by the University of Hawaii under Cooperative Agreement NCC 5-538 with the National Aeronautics and Space Administration, Office of Space Science, Planetary Astronomy Program.

We present a detailed study based on infrared photometry of all Galactic RV Tauri stars from the General Catalogue of Variable Stars (GCVS). RV Tauri stars are the brightest among the population II Cepheids. They are thought to evolve away from the asymptotic giant branch (AGB) towards the white dwarf domain. IRAS detected several RV Tauri stars because of their large IR excesses and it was found that they occupy a specific region in the [12] - [25], [25] - [60] IRAS two-colour diagram. We used the all sky survey of WISE to extend these studies and compare the infrared properties of all RV Tauri stars in the GCVS with a selected sample of post-AGB objects with the goal to place the RV Tauri pulsators in the context of post-AGB evolution. Moreover, we correlated the IR properties of both the RV Tauri stars and the comparison sample with other observables like binarity and the presence of a photospheric chemical anomaly called depletion. We find that Galactic RV Tauri stars display a range of infrared properties and we differentiate between disc sources, objects with no IR excess and objects for which the spectral energy distribution (SED) is uncertain. We obtain a clear correlation between disc sources and binarity. RV Tauri stars with a variable mean magnitude are exclusively found among the disc sources. We also find evidence for disc evolution among the binaries. Furthermore our studies show that the presence of a disc seems to be a necessary but not sufficient condition for the depletion process to become efficient., Comment: 14 pages, 13 figures, 6 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

We performed a detailed analysis of elemental abundances, dust features, and polycyclic aromatic hydrocarbons (PAHs) in the C-rich planetary nebula (PN) Wray16-423 in the Sagittarius dwarf spheroidal galaxy, based on a unique dataset taken from the Subaru/HDS, MPG/ESO FEROS, HST/WFPC2, and Spitzer/IRS. We performed the first measurements of Kr, Fe, and recombination O abundance in this PN. The extremely small [Fe/H] implies that most Fe atoms are in the solid phase, considering into account the abundance of [Ar/H]. The Spitzer/IRS spectrum displays broad 16-24 um and 30 um features, as well as PAH bands at 6-9 um and 10-14 um. The unidentified broad 16-24 um feature may not be related to iron sulfide (FeS), amorphous silicate, or PAHs. Using the spectral energy distribution model, we derived the luminosity and effective temperature of the central star, and the gas and dust masses. The observed elemental abundances and derived gas mass are in good agreement with asymptotic giant branch nucleosynthesis models for an initial mass of 1.90 Msun and a metallicity of Z=0.004. We infer that respectively about 80 %, 50 %, and 90 % of the Mg, S, and Fe atoms are in the solid phase. We also assessed the maximum possible magnesium sulfide (MgS) and iron-rich sulfide (Fe50S) masses and tested whether these species can produce the band flux of the observed 30 um feature. Depending on what fraction of the sulfur is in sulfide molecules such as CS, we conclude that MgS and Fe50S could be possible carriers of the 30 um feature in this PN., Comment: 25 pages, 17 figures, 16 tables, Accepted for publication by MNRAS

There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $\omega$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $\omega$ Centauri., Comment: 21 pages, 16 figures, accepted for publication by MNRAS. Stellar yields included as online data tables

NGC~300 X-1 is a Wolf Rayet + black hole binary that exhibits periodic decreases in X-ray flux. We present two new observations of NGC~300 X-1 from the Chandra X-ray Observatory (totaling $\sim$130 ks) along with ACS imaging data from the Hubble Space Telescope. We observe significant short-term variability in the X-ray emission that is inconsistent with an occultation by the donor star, but is consistent with structure in the outer accretion disk or the wind of the donor star. We simultaneously fit a partially-covered disk blackbody and Comptonized corona model to the eclipse egress and non-eclipsing portions of the X-ray spectrum. We find that the only model parameters that varied between the eclipse egress and non-eclipsing portions of the spectra were the partial covering fraction ($\sim$86% during eclipse egress and $\sim$44% during non-eclipse) and absorbing column ($\sim$12.3$\times10^{22}$ cm$^{-2}$ during eclipse egress, compared to $\sim$1.4$\times10^{22}$ cm$^{-2}$ during non-eclipse). The X-ray spectra are consistent with the movement of the X-ray source through the dense stellar winds of the companion star. From our new HST imaging, we find the WR star within the X-ray error circle, along with additional optical sources including an AGB star and an early-type main sequence star. Finally, we use our egress measurement to rephase previous radial velocity measurements reported in the literature, and find evidence that the velocities are strongly affected by the ionization of the wind by the compact object. Thus, we argue the inferred mass of the black hole may not be reliable.