We present the study of the south-west group in the Fornax cluster centred on the brightest group galaxy (BGG) Fornax A, observed as part of the Fornax Deep Survey (FDS). This includes the analysis of the bright group members (mB < 16 mag) and the intra-group light (IGL). The main objective of this work is to investigate the assembly history of the Fornax A group and to compare its physical quantities as a function of the environment to that of the Fornax cluster core. For all galaxies, we extract the azimuthally averaged surface brightness profiles in three optical bands (g, r, i) by modelling the galaxy's isophotes. We derive their colour profiles and structural parameters in all respective bands. The long integration time and large covered area of the FDS allow us to also estimate the amount of IGL. The majority of galaxies in the Fornax A group are late-type galaxies (LTGs), spanning a range of stellar mass of $8 < log (M_* M_{\odot}) < 10.5$. Six out of nine LTGs show a Type III (up-bending) break in their light profiles, which is either suggestive of strangulation halting star-formation in their outskirts or their HI-richness causing enhanced star-formation in their outer-discs. The estimated luminosity of the IGL is $6 \pm 2 \times 10^{10} L_{\odot}$ in g-band, which corresponds to about 16% of the total light in the group. The Fornax A group appears to be in an early stage of assembly with respect to the cluster core. The environment of the Fornax A group is not as dense as that of the cluster core, with all galaxies except the BGG showing similar morphology, comparable colours and stellar masses, and Type III disc-breaks, without any clear trend of these properties with group-centric distances. The main contribution to the IGL is from the minor merging in the outskirts of the BGG NGC1316 and, probably, the disrupted dwarf galaxies close to the group centre.

In early-type stars a fossil magnetic field may be generated during the star formation process or be the result of a stellar merger event. Surface magnetic fields are thought to be erased by (sub)surface convection layers, which typically leave behind weak disordered fields. However, if the fossil field is strong enough it can prevent the onset of (sub)surface convection and so be preserved onto the main sequence. We calculate the critical field strength at which this occurs, and find that it corresponds well with the lower limit amplitude of observed fields in strongly magnetised Ap/Bp stars ($\approx$ 300 G). The critical field strength is predicted to increase slightly during the main sequence evolution, which could also explain the observed decline in the fraction of magnetic stars. This supports the conclusion that the bimodal distribution of observed magnetic fields in early-type stars reflects two different field origin stories: strongly magnetic fields are fossils fields inherited from star formation or a merger event, and weak fields are the product of on-going dynamo action.

This paper continues the series of the Fornax Deep Survey (FDS). Following the previous studies dedicated to extended Fornax cluster members, we present the catalogs of compact stellar systems in the Fornax cluster as well as extended background sources and point-like sources. We derive ugri photometry of ~1.7 million sources over the $\sim$21 sq. degree area of FDS centered on NGC1399. For a wider area, of $\sim$27 sq. degs extending in the direction of NGC1316, we provide gri data for ~3.1 million sources. To improve the morphological characterization of sources we generate multi-band image stacks by coadding the best seeing gri-band single exposures with a cut at FWHM<=0.9 arcsec. We use the multi-band stacks as detection frames. The identification of compact sources is obtained from a combination of photometric and morphometric selection criteria taking as reference the properties of sources with well-defined classification from the literature. We present a preliminary analysis of globular cluster (GC) distributions in the Fornax area. The study confirms and extends further previous results. We observe the inter-galactic population of GCs, a population of mainly blue GCs centered on NGC1399, extends over $\sim$0.9Mpc, with an ellipticity $\sim$0.65. Several sub-structures extend over $\sim$0.5Mpc along various directions. Two of these structures do not cross any bright galaxy; one of them appears to be connected to NGC1404, a bright galaxy close to the cluster core and particularly poor of GCs. Using the gri catalogs we analyze the GC distribution over the extended FDS area, and do not find any obvious GC sub-structure bridging the two brightest cluster galaxies, NGC1316 and NGC1399. Although NGC1316 is twice brighter of NGC1399 in optical bands we estimate a factor of 3-4 richer GC population around NGC1399 compared to NGC1316, out to galactocentric distances of 40 arcmin

This work is based on deep multi-band (g, r, i) data from the Fornax Deep Survey with VST. We analyse the surface brightness profiles of the 19 bright ETGs inside the virial radius of the Fornax cluster. The main aim of this work is to identify signatures of accretion onto galaxies by studying the presence of outer stellar halos, and understand their nature and occurrence. Our analysis also provides a new and accurate estimate of the intra-cluster light inside the virial radius of Fornax. We performed multi-component fits to the azimuthally averaged surface brightness profiles available for all sample galaxies. This allows to quantify the relative weight of all components in the galaxy structure that contribute to the total light. In addition, we derived the average g-i colours in each component identified by the fit, as well as the azimuthally averaged g-i colour profiles, to correlate them with the stellar mass of each galaxy and the location inside the cluster. We find that in the most massive and reddest ETGs the fraction of light in, probably accreted, halos is much larger than in the other galaxies. Less-massive galaxies have an accreted mass fraction lower than 30%, bluer colours and reside in the low-density regions of the cluster. Inside the virial radius of the cluster, the total luminosity of the intra-cluster light, compared with the total luminosity of all cluster members, is about 34%. Inside the Fornax cluster there is a clear correlation between the amount of accreted material in the stellar halos of galaxies and the density of the environment in which those galaxies reside. By comparing this quantity with theoretical predictions and previous observational estimates, there is a clear indication that the driving factor for the accretion process is the total stellar mass of the galaxy, in agreement with the hierarchical accretion scenario.

Using deep g,r,i imaging from the VEGAS survey, we have searched for ultra diffuse galaxies (UDGs) in the IC 1459 group. Assuming they are group members, we identify 9 galaxies with physical sizes and surface brightnesses that match the UDG criteria within our measurement uncertainties. They have mean colours of g--i = 0.6 and stellar masses of $\sim$10$^8$ M$_{\odot}$. Several galaxies appear to have associated systems of compact objects, e.g. globular clusters. Two UDGs contain a central bright nucleus, with a third UDG revealing a remarkable double nucleus. This appears to be the first reported detection of a double nucleus in a UDG - its origin is currently unclear.

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 report on the search for the optical counterpart of the gravitational event GW170814, which was carried out with the VLT Survey Telescope (VST) by the GRAvitational Wave Inaf TeAm (GRAWITA). Observations started 17.5 hours after the LIGO/Virgo alert and we covered an area of 99 deg$^2$ that encloses $\sim 77\%$ and $\sim 59\%$ of the initial and refined localization probability regions, respectively. A total of six epochs were secured over nearly two months. The survey reached an average limiting magnitude of 22 AB mag in the $r-$band. After assuming the model described in Perna et al. 2019 that derives as possible optical counterpart of a BBH event a transient source declining in about one day, we have computed a survey efficiency of about $5\%$. This paper describes the VST observational strategy and the results obtained by our analysis pipelines developed to search for optical transients in multi-epoch images. We report the catalogue of the candidates with possible identifications based on light-curve fitting. We have identified two dozens of SNe, nine AGNs, one QSO. Nineteen transients characterized by a single detection were not classified. We have restricted our analysis only to the candidates that fall into the refined localization map. None out of 39 left candidates could be positively associated with GW170814. This result implies that the possible emission of optical radiation from a BBH merger had to be fainter than r $\sim$ 22 ($L_{optical}$ $\sim$ $1.4 \times 10^{42}$ erg/s) on a time interval ranging from a few hours up to two months after the GW event.

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

In the last years, Astroinformatics has become a well defined paradigm for many fields of Astronomy. In this work we demonstrate the potential of a multidisciplinary approach to identify globular clusters (GCs) in the Fornax cluster of galaxies taking advantage of multi-band photometry produced by the VLT Survey Telescope using automatic self-adaptive methodologies. The data analyzed in this work consist of deep, multi-band, partially overlapping images centered on the core of the Fornax cluster. In this work we use a Neural-Gas model, a pure clustering machine learning methodology, to approach the GC detection, while a novel feature selection method ($\Phi$LAB) is exploited to perform the parameter space analysis and optimization. We demonstrate that the use of an Astroinformatics based methodology is able to provide GC samples that are comparable, in terms of purity and completeness with those obtained using single band HST data (Brescia et al. 2012) and two approaches based respectively on a morpho-photometric (Cantiello et al. 2018b) and a PCA analysis (D'Abrusco et al. 2015) using the same data discussed in this work.

The formation of Low mass X-ray binaries (LMXB) is favored within dense stellar systems such as Globular Clusters (GCs). The connection between LMXB and Globular Clusters has been extensively studied in the literature, but these studies have always been restricted to the innermost regions of galaxies. We present a study of LMXB in GCs within the central 1.5 deg^2 of the Fornax cluster with the aim of confirming the existence of a population of LMXB in intra-cluster GCs and understand if their properties are related to the host GCs, to the environment or/and to different formation channels.

We present the study of a magnitude limited sample (mB < 16.6 mag) of 13 late type galaxies (LTGs), observed inside the virial radius, Rvir 0.7 Mpc, of the Fornax cluster within the Fornax Deep Survey (FDS). The main objective is to use surface brightness (SB) profiles and g-i colour maps to obtain information on the internal structure of these galaxies and find signatures of the mechanisms that drive their evolution in high-density environment, which is inside the virial radius of the cluster. By modelling galaxy isophotes, we extract the azimuthally averaged surface brightness profiles in four optical bands. We also derive g-i colour profiles, and structural parameters like total magnitude and effective radius. For 10 of the galaxies in this sample, we observe a clear discontinuity in their SB profiles, derive their break radius (Br), and classify their disc-breaks into Type-II (down-bending) or Type-III (up-bending). We find that Type-II galaxies have bluer average (g-i) colour in their outer discs while Type-III galaxies are redder. Br increases with stellar mass and molecular gas mass while decreases with molecular gas-fractions. The inner and outer scale-lengths increase monotonically with absolute magnitude, as found in other works. In Fornax, galaxies with morphological type 5< T< 9 (~60 % of the sample) are located beyond the high-density, ETG-dominated regions, but no correlation found between T and the disc-break type. The main results of this work suggest that the disc-breaks of LTGs inside the virial radius of the Fornax cluster seem to have arisen through a variety of mechanisms, which is evident in their outer-disc colours and the absence of molecular gas beyond their break radius in some cases. This can result in a variety of stellar populations inside and outside the break radii.

Prior to the detection of black holes (BHs) via the gravitational waves (GWs) they generate at merger, the presence of BHs was inferred in X-ray binaries, mostly via dynamical measurements, with masses in the range between $\sim 5-20~M_\odot$. The LIGO discovery of the first BHs via GWs was surprising in that the two BHs that merged had masses of $35.6^{+4.8}_{-3.0}$ and $30.6^{+3.0}_{-4.4}\,M_\odot$, which are both above the range inferred from X-ray binaries. With 20 BH detections from the O1/O2 runs, the distribution of masses remains generally higher than the X-ray inferred one, while the effective spins are generally lower, suggesting that, at least in part, the GW-detected population might be of dynamical origin rather than produced by the common evolution of field binaries. Here we perform high-resolution N-body simulations of a cluster of isolated BHs with a range of initial mass spectra and upper mass cut-offs, and study the resulting binary mass spectrum resulting from the dynamical interactions. Our clusters have properties similar to those of the massive remnants in an OB association $\sim 10 \, \mathrm{Myr}$ after formation. We perform a likelihood analysis for each of our dynamically-formed binary population against the data from the O1 and O2 LIGO/Virgo runs. We find that an initial mass spectrum $M_{\rm BH}\propto M^{-2.35}$ with an upper mass cutoff $M_{\rm max}\sim 50M_\odot$ is favored by the data, together with a slight preference for a merger rate that increases with redshift.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

Many ultra diffuse galaxies (UDGs) have now been identified in clusters of galaxies. However, the number of nearby UDGs suitable for detailed follow-up remain rare. Our aim is to begin to identify UDGs in the environments of nearby bright early-type galaxies from the VEGAS survey. Here we use a deep g band image of the NGC 5846 group, taken as part of the VEGAS survey, to search for UDGs. We found one object with properties of a UDG if it associated with the NGC 5846 group, which seems likely. The galaxy, we name NGC 5846$\_$UDG1, has an absolute magnitude of M$_g$ = -14.2, corresponding to a stellar mass of $\sim$10$^8$ M$_{\odot}$. It also reveals a system of compact sources which are likely globular clusters. Based on the number of globular clusters detected we estimate a halo mass that is greater than 8$\times$10$^{10}$ M$_{\odot}$ for UDG1.

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS.

The Fornax Deep Survey Dwarf galaxy Catalog (FDSDC) includes 564 dwarf galaxies in the Fornax cluster and the in-falling Fornax A subgroup. We use the FDSDC galaxies for statistical analysis of the structural and stellar population differences in the range of galactic environments within the Fornax cluster. We present the standard scaling relations for the dwarfs and analyze trends as a function of cluster-centric radius. We find a different behavior for the bright dwarfs (-18.5 mag < M$_r$ < -16 mag) as compared to the fainter ones (M$_r$ > -16 mag): While considering galaxies in the same magnitude-bins, we find that, while for fainter dwarfs the g'-r' color is redder for lower surface brightness objects (as expected from fading stellar populations), for brighter dwarfs the color is redder for the higher surface brightness and higher S\'ersic n objects. The trend of the bright dwarfs might be explained by those galaxies being affected by harassment and by slower quenching of star formation in their inner parts. As the fraction of early-type dwarfs with respect to late-types increases toward the central parts of the cluster, the color-surface brightness trends are also manifested in the cluster-centric trends, confirming that it is indeed the environment that changes the galaxies. We also estimate the strengths of the ram-pressure stripping, tidal disruption, and harassment in the Fornax cluster, and find that our observations are consistent with the theoretically expected ranges of galaxy properties where each of those mechanisms dominate. We furthermore find that the luminosity function, color-magnitude relation, and axis-ratio distribution of the dwarfs in the center of the Fornax cluster are similar to those in the center of the Virgo cluster.

Weak magnetic fields have recently been detected in a number of A-type stars, including Vega and Sirius. At the same time, space photometry observations of A- and late B-type stars from Kepler and TESS have highlighted the existence of rotational modulation of surface features akin to stellar spots. Here we explore the possibility that surface magnetic spots might be caused by the presence of small envelope convective layers at or just below the stellar surface, caused by recombination of H and He. Using 1D stellar evolution calculations and assuming an equipartition dynamo, we make simple estimates of field strength at the photosphere. For most models the largest effects are caused by a convective layer driven by second helium ionization. While it is difficult to predict the geometry of the magnetic field, we conclude that the majority of intermediate-mass stars should have dynamo-generated magnetic fields of order a few gauss at the surface. These magnetic fields can appear at the surface as bright spots, and cause photometric variability via rotational modulation, which could also be wide-spread in A-stars. The amplitude of surface magnetic fields and their associated photometric variability is expected to decrease with increasing stellar mass and surface temperature, so that magnetic spots and their observational effects should be much harder to detect in late B-type stars.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

Light curves and periodograms of 160 B stars observed by the TESS space mission and 29 main-sequence B stars from Kepler and K2 were used to classify the variability type. There are 114 main-sequence B stars in the TESS sample, of which 45 are classified as possible rotational variables. This confirms previous findings that a large fraction (about 40 percent) of A and B stars may exhibit rotational modulation. Gaia DR2 parallaxes were used to estimate luminosities, from which the radii and equatorial rotational velocities can be deduced. It is shown that observed values of the projected rotational velocities are lower than the estimated equatorial velocities for nearly all the stars, as they should be if rotation is the cause of the light variation. We conclude that a large fraction of main-sequence B stars appear to contain surface features which cannot likely be attributed to abundance patches.

The halo masses $M_{halo}$ of low surface brightness (LSB) galaxies are critical measurements for understanding their formation processes. One promising method to estimate a galaxy's $M_{halo}$ is to exploit the empirical scaling relation between $M_{halo}$ and the number of associated globular clusters ($N_{\mathrm{GC}}$). We use a Bayesian mixture model approach to measure $N_{\mathrm{GC}}$ for 175 LSB ($23\leq\left\langle \mu_{e,r} \right\rangle [\mathrm{mag\ arcsec}^{-2}]\leq 28$) galaxies in the Fornax cluster using the Fornax Deep Survey (FDS) data; this is the largest sample of low mass galaxies so-far analysed for this kind of study. The proximity of the Fornax cluster means that we can measure galaxies with much smaller physical sizes ($0.3\leq r_{e,r}\ [\mathrm{kpc}]\leq 9.5$) compared to previous studies of the GC systems of LSB galaxies, probing stellar masses down to $M_{*}\sim10^{5}\mathrm{M_{\odot}}$. The sample also includes \nudg\ ultra-diffuse galaxies (UDGs), with projected $r$-band half-light radii greater than 1.5 kpc. Our results are consistent with an extrapolation of the $M_{*}-M_{halo}$ relation predicted from abundance matching. In particular, our UDG measurements are consistent with dwarf sized halos, having typical masses between $10^{10}$ and $10^{11}\mathrm{M_{\odot}}$. Overall, our UDG sample is statistically indistinguishable from smaller LSB galaxies in the same magnitude range. We do not find any candidates likely to be as rich as some of those found in the Coma cluster. We suggest that environment might play a role in producing GC-rich LSB galaxies.

Uncertainties in stellar structure and evolution theory are largest for stars undergoing core convection on the main sequence. A powerful way to calibrate the free parameters used in the theory of stellar interiors is asteroseismology, which provides direct measurements of angular momentum and element transport. We report the detection and classification of new variable O and B stars using high-precision short-cadence (2-min) photometric observations assembled by the Transiting Exoplanet Survey Satellite (TESS). In our sample of 154 O and B stars, we detect a high percentage (90%) of variability. Among these we find 23 multiperiodic pulsators, 6 eclipsing binaries, 21 rotational variables, and 25 stars with stochastic low-frequency variability. Several additional variables overlap between these categories. Our study of O and B stars not only demonstrates the high data quality achieved by TESS for optimal studies of the variability of the most massive stars in the Universe, but also represents the first step towards the selection and composition of a large sample of O and B pulsators with high potential for joint asteroseismic and spectroscopic modeling of their interior structure with unprecedented precision.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

This paper is based on the multi-band (ugri) Fornax Deep Survey (FDS) with the VLT Survey Telescope (VST). We study bright early-type galaxies (m_B<15 mag) in the 9 square degrees around the core of the Fornax cluster, which covers the virial radius. The main goal of the present work is to provide the analysis of the light distribution for all galaxies out to unprecedented limits (in radius and surface brightness) and to release the main products resulting from this analysis in all FDS bands. From the isophote fit, we derive the azimuthally averaged surface brightness profiles, the position angle and ellipticity profiles as a function of the semi-major axis. In each band, we derive the total magnitudes, effective radii, integrated colors and stellar mass-to-light ratios. The long integration times, the arcsec-level angular resolution of OmegaCam@VST and the large covered area of FDS allow us to map the light and color distributions out to large galactocentric distances (up to about 10-15 R_e) and surface brightness levels beyond mu_r = 27 mag/arcsec^2 (mu_B > 28 mag/arcsec^2). Therefore, the new FDS data allow us to explore in great detail the morphology and structure of cluster galaxies out to the region of the stellar halo. The observations suggest that the Fornax cluster is not completely relaxed inside the virial radius. The bulk of the gravitational interactions between galaxies happens in the W-NW core region of the cluster, where most of the bright early-type galaxies are located and where the intra-cluster baryons (diffuse light and GCs) are found. We suggest that the W-NW sub-clump of galaxies results from an infalling group onto the cluster, which has modified the structure of the galaxy outskirts (making asymmetric stellar halos) and has produced the intra-cluster baryons (ICL and GCs), concentrated in this region of the cluster.

The Fornax Deep Survey (FDS), an imaging survey in the u', g', r', and i'-bands, has a supreme resolution and image depth compared to the previous spatially complete Fornax Cluster Catalog (FCC). Our new data allows us to study the galaxies down to r'-band magnitude m$_{r'}\approx$21 mag (M$_{r'}\approx$-10.5 mag). These data provide an important legacy dataset to study the Fornax cluster. We aim to present the Fornax Deep Survey (FDS) dwarf galaxy catalog, focusing on explaining the data reduction and calibrations, assessing the quality of the data, and describing the methods used for defining the cluster memberships for the catalog objects. As a first step we used the SExtractor fine-tuned for dwarf galaxy detection, to find galaxies from the FDS data, covering a 26 deg$^2$ area of the main cluster, and the area around the Fornax A substructure. We made 2D-decompositions of the identified galaxies using GALFIT. We used color-magnitude, luminosity-radius and luminosity-concentration relations to separate the cluster galaxies from the background galaxies. We then divided the cluster galaxies into early- and late-type galaxies according to their morphology and gave first order morphological classifications. Our final catalog includes 14,095 galaxies. We classify 590 galaxies as being likely Fornax cluster galaxies, of which 564 are dwarfs (M$_{r'}$ > -18.5 mag) consisting our Fornax dwarf catalog. Of the cluster dwarfs we classify 470 as early-types, and 94 as late-type galaxies. Our final catalog reaches its 50% completeness limit at magnitude M$_{r'}$ = -10.5 mag and surface brightness $\bar{\mu}_{e,r'}$ = 26 mag arcsec-2, which is approximately three magnitudes deeper than the FCC. Based on previous works and comparison with a spectroscopically confirmed subsample, we estimate that our final Fornax dwarf galaxy catalog has < 10% contamination from the background objects.

(Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

The conclusion that the diffuse galaxy NGC1052-DF2 is deficient in dark matter depends critically on the distance, as noted in the original publication and subsequent works. The results of our preliminary assessment of the SBF distance have already been referenced on the web. Here we present the details of a more careful analysis, supply some context, and comment on the remaining uncertainty.

Most of the galaxies in the Universe at present day are in groups, which are key to understanding the galaxy evolution. In this work we present a new deep mosaic of 1.2 x 1.0 square degrees of the group of galaxies centered on NGC 5018, acquired at the ESO VLT Survey Telescope. We use u, g, r images to analyse the structure of the group members and to estimate the intra-group light. Taking advantage of the deep and multiband photometry and of the large field of view of the VST telescope, we studied the structure of the galaxy members and the faint features into the intra-group space and we give an estimate of the intragroup diffuse light in the NGC 5018 group of galaxies. We found that ~ 41% of the total g-band luminosity of the group is in the form of intragroup light (IGL). The IGL has a (g - r) color consistent with those of other galaxies in the group, indicating that the stripping leading to the formation of IGL is ongoing. From the study of this group we can infer that there are at least two different interactions involving the group members: one between NGC 5018 and NGC 5022, which generates the tails and ring-like structures detected in the light, and another between NGC 5022 and MCG-03-34-013 that have produced the HI tail. A minor merging event also happened in the formation history of NGC 5018 that have perturbed the inner structure of this galaxy.

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multi-wavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma-rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubts on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multi-wavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30 degrees away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.

It has been demonstrated that the time variability of a star's brightness at different frequencies can be used to infer its surface gravity, radius, mass, and age. With large samples of light curves now available from Kepler and K2, and upcoming surveys like TESS, we wish to quantify the overall information content of this data and identify where the information resides. As a first look into this question we ask which stellar parameters we can predict from the brightness variations in red-giant stars data and to what precision, using a data-driven model. We demonstrate that the long-cadence (30-minute) Kepler light curves for 2000 red-giant stars can be used to predict their $T_{\rm eff}$ and $\log g$. Our inference makes use of a data-driven model of a part of the autocorrelation function (ACF) of the light curve, where we posit a polynomial relationship between stellar parameters and the ACF pixel values. We find that this model, trained using 1000 stars, can be used to recover the temperature $T_{\rm eff}$ to $<$100 K, the surface gravity to $<$ 0.1 dex, and the asteroseismic power-spectrum parameters $\rm \nu_{max}$ and $\rm \Delta{\nu}$ to $<11$ $\mu$Hz and $<0.9$ $\mu$Hz ($\lesssim$ 15\%). We recover $T_{\rm eff}$ from range of time-lags 0.045 $<$ $T_{\rm lag}$ $<$ 370 days and the $\log g$, $\rm \nu_{max}$ and $\rm \Delta{\nu}$ from the range 0.045 $<$ $T_{\rm lag}$ $<$ 35 days. We do not discover any information about stellar metallicity. The information content of the data about each parameter is empirically quantified using this method, enabling comparisons to theoretical expectations about convective granulation.

Ultra-long gamma ray bursts (ULGRBs) are a distinct class of GRBs characterized by durations of several thousands of seconds, about two orders of magnitude longer than those of standard long GRBs (LGRBs). The driving engine of these events has not been uncovered yet, and ideas range from magnetars, to tidal disruption events, to extended massive stars, such as blue super giants (BSG). BSGs, a possible endpoint of stellar evolution, are attractive for the relatively long free-fall times of their envelopes, allowing accretion to power a long-lasting central engine. At the same time, their large radial extension poses a challenge to the emergence of a jet. Here we perform an end-to-end simulation aimed at assessing the viability of BSGs as ULGRB progenitors. The evolution to core collapse of a BSG star model is calculated with the MESA code. We then compute the accretion rate for the fraction of envelope material with enough angular momentum to circularize and form an accretion disk, and input the corresponding power into a jet which we evolve through the star envelope with the FLASH code. Our simulation shows that the jet can emerge, and the resulting light curves resemble those observed in ULGRBs, with durations $T_{90}$ ranging from $\approx$4000~s to $\approx10^4$ s depending on the viewing angle.

We present the largest and most spatially extended planetary nebulae (PNe) catalog ever obtained for the Fornax cluster. We measured velocities of 1452 PNe out to 200 kpc in the cluster core using a counter-dispersed slitless spectroscopic technique with data from FORS2 on the VLT. With such extended spatial coverage, we can study separately the stellar halos of some of the cluster main galaxies and the intracluster light. In this second paper of the Fornax Cluster VLT Spectroscopic Survey (FVSS), we identify and classify the emission-line sources, describe the method to select PNe and calculate their coordinates and velocities from the dispersed slitless images. From the PN 2D velocity map we identify stellar streams that are possibly tracing the gravitational interaction of NGC1399 with NGC1404 and NGC1387. We also present the velocity dispersion profile out to $\sim 200$ kpc radii, which shows signatures of a superposition of the bright central galaxy and the cluster potential, with the latter clearly dominating the regions outside R $\sim 1000$" ($\sim 100$ kpc).

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

The joint detection of gravitational waves and electromagnetic radiation from the binary neutron star (BNS) merger GW170817 has provided unprecedented insight into a wide range of physical processes: heavy element synthesis via the $r$-process; the production of relativistic ejecta; the equation of state of neutron stars and the nature of the merger remnant; the binary coalescence timescale; and a measurement of the Hubble constant via the "standard siren" technique. In detail, all of these results depend on the distance to the host galaxy of the merger event, NGC4993. In this paper we measure the surface brightness fluctuation (SBF) distance to NGC4993 in the F110W and F160W passbands of the Wide Field Camera 3 Infrared Channel on the Hubble Space Telescope (HST). For the preferred F110W passband we derive a distance modulus of $m{-}M=33.05\pm0.08\pm0.10$ mag, or a linear distance $d=40.7\pm1.4\pm1.9$ Mpc (random and systematic errors, respectively); a virtually identical result is obtained from the F160W data. This is the most precise distance to NGC4993 available to date. Combining our distance measurement with the corrected recession velocity of NGC4993 implies a Hubble constant $H_0=71.9\pm 7.1~km~s^{-1}~Mpc^{-1}$. A comparison of our result to the GW-inferred value of $H_0$ indicates a binary orbital inclination of $i\,{\gtrsim}\,137~\deg$. The SBF technique can be applied to early-type host galaxies of BNS mergers to ${\sim\,}100$ Mpc with HST and possibly as far as ${\sim\,}300$ Mpc with the James Webb Space Telescope, thereby helping to break the inherent distance-inclination degeneracy of the GW signals at distances where many future BNS mergers are likely to be detected.

We describe a program to measure surface brightness fluctuation (SBF) distances to galaxies observed in the Next Generation Virgo Cluster Survey (NGVS), a photometric imaging survey covering $104~deg^2$ of the Virgo cluster in the ${u}^*,g,i,z$ bandpasses with the Canada-France Hawaii Telescope. We describe the selection of the sample galaxies, the procedures for measuring the apparent $i$-band SBF magnitude $\bar{i}$, and the calibration of the absolute $\bar{M}_i$ as a function of observed stellar population properties. The multi-band NGVS data set provides multiple options for calibrating the SBF distances, and we explore various calibrations involving individual color indices as well as combinations of two different colors. Within the color range of the present sample, the two-color calibrations do not significantly improve the scatter with respect to wide-baseline, single-color calibrations involving $u^{*}$. We adopt the ${u}^*{-}z$ calibration as reference for the present galaxy sample, with an observed scatter of 0.11 mag. For a few cases that lack good ${u}^*$ photometry, we use an alternative relation based on a combination of $g{-}i$ and $g{-}z$ colors, with only a slightly larger observed scatter of 0.12 mag. The agreement of our measurements with the best existing distance estimates provides confidence that our measurements are accurate. We present a preliminary catalog of distances for 89 galaxies brighter than $B_T\approx13.0$ mag within the survey footprint, including members of the background M and W Clouds at roughly twice the distance of the main body of the Virgo cluster. The extension of the present work to fainter and bluer galaxies is in progress.

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.

In the core of the Fornax cluster, on the West side of NGC1399, we have detected a previously unknown region of intra-cluster light (ICL). It is made up by several faint ($\mu_r \simeq 28 - 29$~mag/arcsec$^2$) {\it patches} of diffuse light. The bulk of the ICL is located in between the three bright galaxies in the core, NGC1387, NGC1379 and NGC1381, at $10\leq R \leq40$~arcmin ($\sim 58 - 230$~kpc) from the central galaxy NGC~1399. We show that the ICL is the counterpart in the diffuse light of the known over-density in the population of blue globular clusters (GCs). The total g-band luminosity of the ICL is $L_g\simeq 8.3 \times 10^{9}$ $L_{\odot}$, which is $\sim5\%$ of the total luminosity of NGC1399. This is consistent with the fraction of the blue GCs in the same region of the cluster. The ICL has $g-r \sim 0.7$~mag, which is similar to the colors in the halo of the bright galaxies in the cluster core. The new findings were compared with theoretical predictions for the ICL formation and they support a scenario in which the intra-cluster population detected in the core of the Fornax cluster is build up by the tidal stripping of material (stars and GCs) from galaxy outskirts in a close passage with the cD. Moreover, the diffuse form of the ICL and its location close to the core of the cluster is expected in a dynamically evolved cluster as Fornax.

We analyze the globular cluster (GC) systems in two very different galaxies, N3115 and N1399. With the papers of this series, we aim at highlighting common and different properties in the GC systems in galaxies covering a wide range of parameter space. We compare the GCs in N3115 and N1399 as derived from the analysis of u, g, and i band images taken with the VST telescope as part of the VEGAS and FDS surveys. We selected GC candidates using as reference the morpho-photometric and color properties of confirmed GCs. The surface density of GCs in N3115 reveal a morphology similar to the profile of field stars; the same is true when blue and red GCs are taken separately. The GC maps for N1399 are richer in structure and confirm the existence of an intracluster GC component. We confirm the presence of a spatial offset in the N1399 GC centroid and find that the centroid of the GCs for N3115 coincides well with the galaxy center. Both GC systems show unambiguous color bimodality in g-i and u-i. The azimuthal average of the radial density profiles in both galaxies reveals a larger spatial extent for the total GCs population with respect to the galaxy surface brightness profile. As for the specific frequency of GCs, Sn, we find it is a factor of two higher in N1399 than for N3115. By inspecting the radial behavior of the specific frequency, Sn(<r), for the total, blue, and red GCs, we find notable similarities between the trends for red GCs in the two targets. In spite of extremely different host environments, the red GCs in both cases appear closely linked to the light distribution of field stars. Blue GCs extend to larger galactocentric scales than red GCs, marking a significant difference between the two galaxies: the blue/red GCs and field stellar components of N3115 appear well thermalized with each other and the blue GCs in N1399 appear to fade into an unrelaxed intracluster GC population.

Aims. We analyze the properties of the sources in the NGC253 to define an up to date catalog of GC candidates in the galaxy. Methods. Our analysis is based on the science verification data of two ESO survey telescopes, VST and VISTA. Using ugri photometry from VST and JKs from VISTA, GC candidates were selected using the morpho-photometric and color properties of spectroscopically confirmed GCs available in the literature. The strength of the results was verified against available archival HST/ACS data from the GHOSTS survey. Results. The adopted GC selection leads to the definition of a sample of ~350 GC candidates. At visual inspection, we find that 82 objects match all the requirements for selecting GC candidates and 155 are flagged as uncertain GC candidate; 110 are unlikely GCs, most likely background galaxies. Furthermore, our analysis shows that four of the previously spectroscopically confirmed GCs, i.e., ~20% of the total spectroscopic sample, are more likely either background galaxies or high-velocity Milky Way stars. The radial density profile of the selected best candidates shows the typically observed r1/4-law radial profile. The analysis of the color distributions reveals only marginal evidence of the presence of color bimodality, which is normally observed in galaxies of similar luminosity. The GC luminosity function does not show the typical symmetry, mainly because of the lack of bright GCs. Part of the bright GCs missing might be at very large galactocentric distances or along the line of sight of the galaxy dusty disk. Conclusions. Using ugriJKs photometry we purged the list of GCs with spectroscopic membership and photometric GC candidates in NGC 253. Our results show that the use of either spectroscopic or photometric data only does not generally ensure a contaminant-free sample and a combination of both spectroscopy and photometry is preferred.

Studies of low surface brightness (LSB) galaxies in nearby clusters have revealed a sub-population of extremely diffuse galaxies with central surface brightness $\mu_{0,g'}$ > 24 mag arcsec$^{-2}$ and effective radius between 1.5 kpc < R$_{e}$ < 10 kpc. The origin of these Ultra Diffuse Galaxies (UDGs) is still unclear, although several theories have been suggested. We exploit the deep g', r' and i'-band images of the Fornax Deep Survey (FDS), in order to identify LSB galaxies in the center of the Fornax cluster. We identified visually all extended structures having r'-band central surface brightness of $\mu_{0,r'}$ > 23 mag arcsec$^{-2}$. We classified the objects based on their appearance and performed 2D S\'ersic model fitting with GALFIT. We analyzed their distribution and orientations in the cluster, and studied their colors and compared the LSB galaxies in Fornax with those in other environments. Our sample consists of 205 galaxies of which 196 are LSB dwarfs (with R$_e$ < 1.5kpc) and nine are UDGs (R$_e$ > 1.5 kpc). We show that the UDGs have g'-r' colors similar to those of LSB dwarfs. The largest UDGs in our sample appear different from the other LSB galaxies, in that they are significantly more elongated and extended, whereas the smaller UDGs differ from the LSB dwarfs only by their effective radii. We do not find clear differences between the structural parameters of the UDGs in our sample and those of UDGs in other galaxy environments. We find that the dwarf LSB galaxies in our sample are less concentrated in the cluster center than the galaxies with higher surface brightness, and that their number density drops in the core of the cluster. Our findings are consistent with the small UDGs forming the tail of a continuous distribution of less extended LSB galaxies. However, the elongated and distorted shapes of the large UDGs could imply that they are tidally disturbed galaxies.

We present the results of numerical simulations of the prompt emission of short-duration gamma-ray bursts. We consider emission from the relativistic jet, the mildly relativistic cocoon, and the non-relativistic shocked ambient material. We find that the cocoon material is confined between off-axis angles 15<theta<45 degrees and gives origin to X-ray transients with a duration of a few to ~10 seconds, delayed by a few seconds from the time of the merger. We also discuss the distance at which such transients can be detected, finding that it depends sensitively on the assumptions that are made about the radiation spectrum. Purely thermal cocoon transients are detectable only out to a few Mpc, Comptonized transients can instead be detected by the FERMI GBM out to several tens of Mpc.

The arXiv is the most popular preprint repository in the world. Since its inception in 1991, the arXiv has allowed researchers to freely share publication-ready articles prior to formal peer review. The growth and the popularity of the arXiv emerged as a result of new technologies that made document creation and dissemination easy, and cultural practices where collaboration and data sharing were dominant. The arXiv represents a unique place in the history of research communication and the Web itself, however it has arguably changed very little since its creation. Here we look at the strengths and weaknesses of arXiv in an effort to identify what possible improvements can be made based on new technologies not previously available. Based on this, we argue that a modern arXiv might in fact not look at all like the arXiv of today.

Observations of diffuse starlight in the outskirts of galaxies are thought to be a fundamental source of constraints on the cosmological context of galaxy assembly in the $\Lambda$CDM model. Such observations are not trivial because of the extreme faintness of such regions. In this work, we investigate the photometric properties of six massive early type galaxies (ETGs) in the VEGAS sample (NGC 1399, NGC 3923, NGC 4365, NGC 4472, NGC 5044, and NGC 5846) out to extremely low surface brightness levels, with the goal of characterizing the global structure of their light profiles for comparison to state-of-the-art galaxy formation models. We carry out deep and detailed photometric mapping of our ETG sample taking advantage of deep imaging with VST/OmegaCAM in the g and i bands. By fitting the light profiles, and comparing the results to simulations of elliptical galaxy assembly, we identify signatures of a transition between "relaxed" and "unrelaxed" accreted components and can constrain the balance between in situ and accreted stars. The very good agreement of our results with predictions from theoretical simulations demonstrates that the full VEGAS sample of $\sim 100$ ETGs will allow us to use the distribution of diffuse light as a robust statistical probe of the hierarchical assembly of massive galaxies.

As part of the Fornax Deep Survey with the ESO VLT Survey Telescope, we present new $g$ and $r$ bands mosaics of the SW group of the Fornax cluster. It covers an area of $3 \times 2$ square degrees around the central galaxy NGC1316. The deep photometry, the high spatial resolution of OmegaCam and the large covered area allow us to study the galaxy structure, to trace stellar halo formation and look at the galaxy environment. We map the surface brightness profile out to 33arcmin ($\sim 200$kpc $\sim15R_e$) from the galaxy centre, down to $\mu_g \sim 31$ mag arcsec$^{-2}$ and $\mu_r \sim 29$ mag arcsec$^{-2}$. This allow us to estimate the scales of the main components dominating the light distribution, which are the central spheroid, inside 5.5 arcmin ($\sim33$ kpc), and the outer stellar envelope. Data analysis suggests that we are catching in act the second phase of the mass assembly in this galaxy, since the accretion of smaller satellites is going on in both components. The outer envelope of NGC1316 still hosts the remnants of the accreted satellite galaxies that are forming the stellar halo. We discuss the possible formation scenarios for NGC1316, by comparing the observed properties (morphology, colors, gas content, kinematics and dynamics) with predictions from cosmological simulations of galaxy formation. We find that {\it i)} the central spheroid could result from at least one merging event, it could be a pre-existing early-type disk galaxy with a lower mass companion, and {\it ii)} the stellar envelope comes from the gradual accretion of small satellites.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.

Asteroseismology probes the interiors of stars by studying oscillation modes at a star's surface. Although pulsation spectra are well understood for solar-like oscillators, a substantial fraction of red giant stars observed by Kepler exhibit abnormally low-amplitude dipole oscillation modes. Fuller et al. (2015) suggests this effect is produced by strong core magnetic fields that scatter dipole internal gravity waves (IGWs) into higher multipole IGWs or magnetic waves. In this paper, we study the interaction of IGWs with a magnetic field to test this mechanism. We consider two background stellar structures: one with a uniform magnetic field, and another with a magnetic field that varies both horizontally and vertically. We derive analytic solutions to the wave propagation problem and validate them with numerical simulations. In both cases, we find perfect conversion from IGWs into magnetic waves when the IGWs propagate into a region exceeding a critical magnetic field strength. Downward propagating IGWs cannot reflect into upward propagating IGWs because their vertical wavenumber never approaches zero. Instead, they are converted into upward propagating slow (Alfvenic) waves, and we show they will likely dissipate as they propagate back into weakly magnetized regions. Therefore, strong internal magnetic fields can produce dipole mode suppression in red giants, and gravity modes will likely be totally absent from the pulsation spectra of sufficiently magnetized stars.

We use three dimensional radiation magneto-hydrodynamic simulations to study the effects of magnetic fields on the energy transport and structure of radiation pressure dominated main sequence massive star envelopes at the region of the iron opacity peak. We focus on the regime where the local thermal timescale is shorter than the dynamical timescale, corresponding to inefficient convective energy transport. We begin with initially weak magnetic fields relative to the thermal pressure, from 100-1000G in differing geometries. The unstable density inversion amplifies the magnetic field, increasing the magnetic energy density to values close to equipartition with the turbulent kinetic energy density. By providing pressure support, the magnetic field's presence significantly increases the density fluctuations in the turbulent envelope, thereby enhancing the radiative energy transport by allowing photons to diffuse out through low density regions. Magnetic buoyancy brings small scale magnetic fields to the photosphere and increases the vertical energy transport with the energy advection velocity proportional to the Alfv\'en velocity, although in all cases we study photon diffusion still dominates the energy transport. The increased radiative and advective energy transport causes the stellar envelope to shrink by several scale heights. We also find larger turbulent velocity fluctuations compared to the purely hydrodynamic case, reaching $\approx$ 100 km/s at the stellar photosphere. The photosphere also shows vertical oscillations with similar averaged velocities and periods of a few hours. The increased turbulent velocity and oscillations will have strong impacts on the line broadening and periodic signals in massive stars.

The Next Generation Virgo Cluster Survey is a deep (with a $2\sigma$ detection limit $\mu_g$ = 29~mag~arcsec$^{-2}$ in the $g-$band) optical panchromatic survey targeting the Virgo cluster from its core to virial radius, for a total areal coverage of 104 square degrees. As such, the survey is well suited for the study of galaxies' outskirts, haloes and low surface brightness features that arise from dynamical interactions within the cluster environment. We report the discovery of extremely faint ($\mu_g$ $>$ 25 mag arcsec$^{-2}$) shells in three Virgo cluster early-type dwarf galaxies, VCC~1361, VCC~1447 and VCC~1668. Among them, VCC~1447 has an absolute magnitude M$_{g}$ = -11.71 mag and is {\it the least massive galaxy with a shell system discovered to date}. We present a detailed study of these low surface brightness features. We detect between three and four shells in each of our galaxies. Within the uncertainties, we find no evidence of a color difference between the galaxy main body and shell features. The observed arcs of the shells are located upto several effective radii of the galaxies. We further explore the origin of these low surface brightness features with the help of idealized numerical simulations. We find that a near equal mass merger is best able to reproduce the main properties of the shells, including their quite symmetric appearance and their alignment along the major axis of the galaxy. The simulations provide support for a formation scenario in which a recent merger, between two near-equal mass, gas-free dwarf galaxies forms the observed shell systems.

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.

Deep observations of the dwarf elliptical (dE) galaxy NGC 1396 (M$_V = -16.60$, Mass $\sim 4\times10^8$ M$_\odot$), located in the Fornax cluster, have been performed with the VLT/ MUSE spectrograph in the wavelength region from $4750-9350$ \AA{}. In this paper we present a stellar population analysis studying chemical abundances, the star formation history (SFH) and the stellar initial mass function (IMF) as a function of galacto-centric distance. Different, independent ways to analyse the stellar populations result in a luminosity-weighted age of $\sim$ 6 Gyr and a metallicity [Fe/H]$\sim$ $-0.4$, similar to other dEs of similar mass. We find unusually overabundant values of [Ca/Fe] $\sim +0.1$, and under-abundant Sodium, with [Na/Fe] values around $-0.1$, while [Mg/Fe] is overabundant at all radii, increasing from $\sim+0.1$ in the centre to $\sim +0.2$ dex. We notice a significant metallicity and age gradient within this dwarf galaxy. To constrain the stellar IMF of NGC 1396, we find that the IMF of NGC 1396 is consistent with either a Kroupa-like or a top-heavy distribution, while a bottom-heavy IMF is firmly ruled out. An analysis of the abundance ratios, and a comparison with galaxies in the Local Group, shows that the chemical enrichment history of NGC 1396 is similar to the Galactic disc, with an extended star formation history. This would be the case if the galaxy originated from a LMC-sized dwarf galaxy progenitor, which would lose its gas while falling into the Fornax cluster.

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

We report the discovery of a complex extended density enhancement in the Globular Clusters (GCs) in the central $\sim 0.5(^{\circ})^2$ ($\sim 0.06$ Mpc$^2$) of the Fornax cluster, corresponding to $\sim 50\%$ of the area within 1 core radius. This overdensity connects the GC system of NGC1399 to most of those of neighboring galaxies within $\sim 0.6^{\circ}$ ($\sim 210$ kpc) along the W-E direction. The asymmetric density structure suggests that the galaxies in the core of the Fornax cluster experienced a lively history of interactions that have left a clear imprint on the spatial distribution of GCs. The extended central dominant structure is more prominent in the distribution of blue GCs, while red GCs show density enhancements that are more centrally concentrated on the host galaxies. We propose that the relatively small-scale density structures in the red GCs are caused by galaxy-galaxy interactions, while the extensive spatial distribution of blue GCs is due to stripping of GCs from the halos of core massive galaxies by the Fornax gravitational potential. Our investigations is based on density maps of candidate GCs extracted from the multi-band VLT Survey Telescope (VST) survey of Fornax (FDS), identified in a three-dimensional color space and further selected based on their $g$-band magnitude and morphology.

The asteroseismology of red giant stars has continued to yield surprises since the onset of high-precision photometry from space-based observations. An exciting new theoretical result shows that the previously observed suppression of dipole oscillation modes in red giants can be used to detect strong magnetic fields in the stellar cores. A fundamental facet of the theory is that nearly all the mode energy leaking into the core is trapped by the magnetic greenhouse effect. This results in clear predictions for how the mode visibility changes as a star evolves up the red giant branch, and how that depends on stellar mass, spherical degree, and mode lifetime. Here, we investigate the validity of these predictions with a focus on the visibility of different spherical degrees. We find that mode suppression weakens for higher degree modes with an average reduction in the quadrupole mode visibility of up to 49% for the least evolved stars in our sample, and no detectable suppression of octupole modes, in agreement with the theoretical predictions. We furthermore find evidence for the influence of increasing mode lifetimes on the measured visibilities along the red giant branch, in agreement with previous independent observations. These results support the theory that strong internal magnetic fields are responsible for the observed suppression of non-radial modes in red giants. We also find preliminary evidence that stars with suppressed dipole modes on average have slightly lower metallicity than normal stars.

Recent asteroseismic analyses have revealed the presence of strong (B $\gtrsim 10^5$ G) magnetic fields in the cores of many red giant stars. Here, we examine the implications of these results for the evolution of stellar magnetic fields, and we make predictions for future observations. Those stars with suppressed dipole modes indicative of strong core fields should exhibit moderate but detectable quadrupole mode suppression. The long magnetic diffusion times within stellar cores ensure that dynamo-generated fields are confined to mass coordinates within the main sequence convective core, and the observed sharp increase in dipole mode suppression rates above $1.5 \, M_\odot$ may be explained by the larger convective core masses and faster rotation of these more massive stars. In clump stars, core fields of $\sim10^5 \, {\rm G}$ can suppress dipole modes, whose visibility should be equal to or less than the visibility of suppressed modes in ascending red giants. High dipole mode suppression rates in low-mass ($M \lesssim 2 \, M_\odot$) clump stars would indicate that magnetic fields generated during the main sequence can withstand subsequent convective phases and survive into the compact remnant phase. Finally, we discuss implications for observed magnetic fields in white dwarfs and neutron stars, as well as the effects of magnetic fields in various types of pulsating stars.

[Abrigded] We have started a new deep, multi-imaging survey of the Fornax cluster, dubbed Fornax Deep Survey (FDS), at the VLT Survey Telescope. In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC1399 in the core of the cluster. We found a very extended and diffuse envelope surrounding the luminous galaxy NGC1399: we map the surface brightness out to 33 arcmin (~ 192 kpc) from the galaxy center and down to about 31 mag/arcsec^2 in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region between NGC1399 and NGC1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC1387 on its east side is stripped away. By fitting the light profile, we found that it exists a physical break radius in the total light distribution at R=10 arcmin (~58 kpc) that sets the transition region between the bright central galaxy and the outer exponential stellar halo. We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive BCGs, we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with a stellar mass in the range 10^8 - 10^11 M_sun. This might suggest that the halo of NGC1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, differently from the Virgo cluster, the extended stellar halo around NGC1399 is characterised by a more diffuse and well-mixed component, including the ICL.

We present multi-wavelength observations of SN2014C during the first 500 days. These observations represent the first solid detection of a young extragalactic stripped-envelope SN out to high-energy X-rays. SN2014C was the explosion of an H-stripped progenitor star with ordinary explosion parameters. However, over the time scale of ~1yr, SN2014C experienced a complete metamorphosis and evolved from an ordinary H-poor supernova of type Ib into a strongly interacting, H-rich supernova of type IIn. Signatures of the SN shock interacting with a dense medium are observed across the spectrum. Coordinated observations with Swift, Chandra and NuSTAR have captured the evolution in detail and revealed the presence of a massive shell of ~1 Msun of hydrogen-rich material at ~6d16 cm from the explosion site. We estimate that the shell was ejected by the progenitor star in the decades to centuries before core collapse. This result poses significant challenges to current theories of massive star evolution, as it requires a physical mechanism responsible for the ejection of the deepest hydrogen layer of H-poor SN progenitors synchronized with the onset of stellar collapse. Theoretical investigations point at binary interactions and/or instabilities during the last stages of nuclear burning in massive stars as potential triggers of the time-dependent mass loss. We constrain these scenarios utilizing the sample of 183 SNe Ib/c with public radio observations. Our analysis identifies SN2014C-like signatures in ~10% of SNe with constraining radio data. This fraction is somewhat larger but reasonably consistent with the expectation from the theory of recent envelope ejection due to binary evolution IF the ejected material can survive in the close environment for 1000-10000 yrs. Alternatively, nuclear burning instabilities extending all the way to the core C-burning phase might also play a critical role.

Magnetic fields play a role in almost all stages of stellar evolution. Most low-mass stars, including the Sun, show surface fields that are generated by dynamo processes in their convective envelopes. Intermediate-mass stars do not have deep convective envelopes, although 10% exhibit strong surface fields that are presumed to be residuals from the stellar formation process. These stars do have convective cores that might produce internal magnetic fields, and these might even survive into later stages of stellar evolution, but information has been limited by our inability to measure the fields below the stellar surface. Here we use asteroseismology to study the occurrence of strong magnetic fields in the cores of low- and intermediate-mass stars. We have measured the strength of dipolar oscillation modes, which can be suppressed by a strong magnetic field in the core, in over 3600 red giant stars observed by Kepler. About 20% of our sample show mode suppression but this fraction is a strong function of mass. Strong core fields only occur in red giants above 1.1 solar masses (1.1Msun), and the occurrence rate is at least 60% for intermediate-mass stars (1.6--2.0Msun), indicating that powerful dynamos were very common in the convective cores of these stars.

Internal stellar magnetic fields are inaccessible to direct observations and little is known about their amplitude, geometry and evolution. We demonstrate that strong magnetic fields in the cores of red giant stars can be identified with asteroseismology. The fields can manifest themselves via depressed dipole stellar oscillation modes, which arises from a magnetic greenhouse effect that scatters and traps oscillation mode energy within the core of the star. The Kepler satellite has observed a few dozen red giants with depressed dipole modes which we interpret as stars with strongly magnetized cores. We find field strengths larger than $\sim\! 10^5 \,{\rm G}$ may produce the observed depression, and in one case we infer a minimum core field strength of $\approx \! \! 10^7 \,{\rm G}$.

The spin of a number of black holes (BHs) in binary systems has been measured. In the case of BHs found in low-mass X-ray binaries (LMXBs) the observed values are in agreement with some theoretical predictions based on binary stellar evolution. However, using the same evolutionary models, the calculated spins of BHs in high-mass X-ray binaries (HMXBs) fall short compared to the observations. A possible solution to this conundrum is the accretion of high-specific-angular-momentum material after the formation of the BH, although this requires accretion above the Eddington limit. Another suggestion is that the observed high values of the BHs spin could be the result of an asymmetry during Core Collapse (CC). The only available energy to spin up the compact object during CC is its binding energy. A way to convert it to rotational kinetic energy is by using a Standing Accretion Shock Instability (SASI), which can develop during CC and push angular momentum into the central compact object through a spiral mode ($m = 1$). Here we study the CC-SASI scenario and discuss, in the case of LMXBs and HMXBs, the limits for the spin of a stellar-mass BHs. Our results predict a strong dichotomy in the maximum spin of low-mass compact objects and massive BHs found in HMXBs. The maximum spin value ($|a_\star|$) for a compact object near the mass boundary between BHs and NSs is found to be somewhere between 0.27 and 0.38, depending on whether secular or dynamical instabilities limit the efficiency of the spin up process. For more massive BHs, such as those found in HMXBs, the natal spin is substantially smaller and for $M_{\rm BH}\!>\!8~M_\odot$ spin is limited to values $|a_\star|\!\lesssim\!0.05$. Therefore we conclude that the observed high spins of BHs in HMXBs cannot be the result of a CC-SASI spin up.

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $\tau_0$, to $\tau_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $\tau_0 \gg \tau_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $\tau_0 \ll \tau_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 30 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 5 \times 10^3 \, {\rm s}$ in typical core-collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 500 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.

We present the VST Early-type GAlaxy Survey (VEGAS), which is designed to obtain deep multiband photometry in g, r, i, of about one hundred nearby galaxies down to 27.3, 26.8, and 26 mag/arcsec^2 respectively, using the ESO facility VST/OmegaCAM. The goals of the survey are 1) to map the light distribution up to ten effective radii, r_e, 2) to trace color gradients and surface brightness fluctuation gradients out to a few r_e for stellar population characterization, and 3) to obtain a full census of the satellite systems (globular clusters and dwarf galaxies) out to 20% of the galaxy virial radius. The external regions of galaxies retain signatures of the formation and evolution mechanisms that shaped them, and the study of nearby objects enables a detailed analysis of their morphology and interaction features. To clarify the complex variety of formation mechanisms of early-type galaxies (ETGs), wide and deep photometry is the primary observational step, which at the moment has been pursued with only a few dedicated programs. The VEGAS survey has been designated to provide these data for a volume-limited sample with exceptional image quality. In this commissioning photometric paper we illustrate the capabilities of the survey using g- and i-band VST/OmegaCAM images of the nearby galaxy NGC 4472 and of smaller ETGs in the surrounding field. Our surface brightness profiles reach rather faint levels and agree excellently well with previous literature. Genuine new results concern the detection of an intracluster light tail in NGC 4472 and of various substructures at increasing scales. We have also produced extended (g-i) color profiles. The VST/OmegaCAM data that we acquire in the context of the VEGAS survey provide a detailed view of substructures in the optical emission from extended galaxies, which can be as faint as a hundred times below the sky level.

We present new calibrations of the near-infrared surface brightness fluctuation (SBF) distance method for the F110W (J) and F160W (H) bandpasses of the Wide Field Camera 3 Infrared Channel (WFC3/IR) on the Hubble Space Telescope. The calibrations are based on data for 16 early-type galaxies in the Virgo and Fornax clusters observed with WFC3/IR and are provided as functions of both the optical (g-z) and near-infrared (J-H) colors. The scatter about the linear calibration relations for the luminous red galaxies in the sample is approximately 0.10 mag, corresponding to a statistical error of 5% in distance. Our results imply that the distance to any suitably bright elliptical galaxy can be measured with this precision out to about 80 Mpc in a single-orbit observation with WFC3/IR, making this a remarkably powerful instrument for extragalactic distances. The calibration sample also includes much bluer and lower-luminosity galaxies than previously used for IR SBF studies, revealing interesting population differences that cause the calibration scatter to increase for dwarf galaxies.

We present a study of globular clusters (GCs) and other small stellar systems (SSSs) in the field of NGC3115, observed as part of the VEGAS imaging survey, carried out with the VST telescope. We use deep g and i data of NGC3115, a well-studied lenticular galaxy with excellent scientific literature. This is fundamental to test the methodologies, verify the results, and probe the capabilities of the VEGAS-SSS. Leveraging the large field of view of the VST allow us to accurately study of the distribution and properties of SSSs as a function of galactocentric distance Rgc, well beyond ~20 galaxy effective radii, in a way not often possible. Our analysis of colors, magnitudes and sizes of SSS candidates confirms the results from existing studies, some of which carried out with 8-10m class telescopes, and further extends them to unreached Rgc distances, with similar accuracy. We find a color bimodality for the GC population and a r1/4 profile for the surface density of GCs as for the galaxy light profile. The radial color gradient of blue and red GCs previously found, for instance by the SLUGGS survey, is further extended out to the largest Rgc inspected, ~65 kpc. The surface density profiles of blue and red GCs taken separately are well approximated by a r1/4 density profile, with the fraction of blue GCs being slightly larger at larger Rgc. We do not find hints of a trend for the red GC and for the GC turnover magnitude to vary with radius, but we observe a ~0.2 mag difference in the turnover magnitude of the blue and red GCs subpopulations. Inspecting SSS sizes and colors we obtained a list of UCDs and GC candidates suitable for future spectroscopic follow-up. In conclusion, our study shows the reliability of the methodologies developed to study SSSs in the field of bright early-type galaxies, and the great potential of the VEGAS survey to produce original results on SSSs science.

We present the first deep photometry of a good candidate for a forming polar ring galaxy at redshift z~0.05. This object, named FCSS J033710.0-354727, is a background galaxy in the VST deep field of the Fornax cluster. The deep exposures combined with the high angular resolution of the OmegaCAM at VST allow us to carry out the first detailed photometric analysis for this system in the g and i bands to derive the galaxy structure and colors. Results show that the central object resembles a disk galaxy, surrounded by a ring-like structure 2 times more extended than the central disk. The warped geometry and the presence of bright knots observed along the polar direction, as well as the several debris detected on the NW side with colors comparable to those of the galaxy, suggest that the polar structure is still forming. We argue that the wide polar ring/disk is the result of the ongoing disruption of a companion galaxy in the potential of the central object, which is 2-3 times more massive than the accreting galaxy.

Recent asteroseismic advances have allowed for direct measurements of the internal rotation rates of many sub-giant and red giant stars. Unlike the nearly rigidly rotating Sun, these evolved stars contain radiative cores that spin faster than their overlying convective envelopes, but slower than they would in the absence of internal angular momentum transport. We investigate the role of internal gravity waves in angular momentum transport in evolving low mass stars. In agreement with previous results, we find that convectively excited gravity waves can prevent the development of strong differential rotation in the radiative cores of Sun-like stars. As stars evolve into sub-giants, however, low frequency gravity waves become strongly attenuated and cannot propagate below the hydrogen burning shell, allowing the spin of the core to decouple from the convective envelope. This decoupling occurs at the base of the sub-giant branch when stars have surface temperatures of roughly 5500 K. However, gravity waves can still spin down the upper radiative region, implying that the observed differential rotation is likely confined to the deep core near the hydrogen burning shell. The torque on the upper radiative region may also prevent the core from accreting high-angular momentum material and slow the rate of core spin-up. The observed spin-down of cores on the red giant branch cannot be totally attributed to gravity waves, but the waves may enhance shear within the radiative region and thus increase the efficacy of viscous/magnetic torques.

Hercules is the prototype of the ultra-faint dwarf (UFD) galaxies. To date, there are still no firm constraints on its total luminosity, due to the difficulty of disentangling Hercules bona-fide stars from the severe Galactic field contamination. In order to better constrain Hercules properties we aim at removing foreground and background contaminants in the galaxy field using the proper motions of the Milky Way stars and the colour-colour diagram. We have obtained images of Hercules in the rSloan, BBessel and Uspec bands with the Large Binocular Telescope (LBT) and LBC-BIN mode capabilities. The rSloan new data-set combined with data from the LBT archive span a time baseline of about 5 yr, allowing us to measure for the first time proper motions of stars in the Hercules direction. The Uspec data along with existing LBT photometry allowed us to use colour-colour diagram to further remove the field contamination. Thanks to a highly-accurate procedure to derive the rSloan-filter geometric distortion solution for the LBC-red, we were able to measure stellar relative proper motions to a precision of better than 5 mas yr^-1 down to rSloan=22 mag and disentangle a significant fraction (\>90\%) of Milky Way contaminants. We ended up with a sample of 528 sources distributed over a large portion of the galaxy body (0.12 deg^2). Of these sources, 171 turned out to be background galaxies and additional foreground stars, from the analysis of the Uspec - BBessel vs. BBessel - rSloan colour-colour diagram. This leaves us with a sample of 357 likely members of the Hercules UFD. We compared the cleaned colour-magnitude diagram (CMD) with evolutionary models and synthetic CMDs, confirming the presence in Hercules of an old population (t=12\pm 2 Gyr), with a wide spread in metallicity (-3.3\<[Fe/H]\<-1.8).

STREGA (STRucture and Evolution of the GAlaxy) is a Guaranteed Time survey being performed at the VST (the ESO VLT Survey Telescope) to map about 150 square degrees in the Galactic halo, in order to constrain the mechanisms of galactic formation and evolution. The survey is built as a five-year project, organized in two parts: a core program to explore the surrounding regions of selected stellar systems and a second complementary part to map the southern portion of the Fornax orbit and extend the observations of the core program. The adopted stellar tracers are mainly variable stars (RR~Lyraes and Long Period Variables) and Main Sequence Turn-off stars for which observations in the g,r,i bands are obtained. We present an overview of the survey and some preliminary results for three observing runs that have been completed. For the region centered on $\omega$~Cen (37 deg^2), covering about three tidal radii, we also discuss the detected stellar density radial profile and angular distribution, leading to the identification of extratidal cluster stars. We also conclude that the cluster tidal radius is about 1.2 deg, in agreement with values in the literature based on the Wilson model.

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

STEP (the SMC in Time: Evolution of a Prototype interacting late-type dwarf galaxy) is a Guaranteed Time Observation survey being performed at the VST (the ESO VLT Survey Telescope). STEP will image an area of 74 deg$^2$ covering the main body of the Small Magellanic Cloud (32 deg$^2$), the Bridge that connects it to the Large Magellanic Cloud (30 deg$^2$) and a small part of the Magellanic Stream (2 deg$^2$). Our $g,r,i,H_{\alpha}$ photometry is able to resolve individual stars down to magnitudes well below the main-sequence turnoff of the oldest populations. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products of the survey. We also present preliminary results for the first two fields for which data acquisition is completed, including some detailed analysis of the two stellar clusters IC\,1624 and NGC\,419.

We combined new near-infrared VLT/HAWK-I data of the globular clusters (GCs) in the isolated edge-on S0 galaxy NGC 3115 with optical and spectroscopic ones taken from the literature, with the aim of analyzing the multiband GC color distributions. A recent study from the SLUGGS survey has shown that the GCs in this galaxy follow a bimodal distribution of Ca II triplet indices. Thus, NGC 3115 presents a critical example of a GC system with multiple, distinct, metallicity subpopulations, and this may argue against the "projection" scenario, which posits that the ubiquitous color bimodality mainly results from nonlinearities in the color-metallicity relations. Using optical, NIR, and spectroscopic data, we found strong and consistent evidence of index bimodality, which independently confirms the metallicity bimodality in NGC 3115 GCs. At the same time, we also found evidence for some color-color nonlinearity. Taken in the broader context of previous studies, the multicolor consistency of the GC bimodality in NGC 3115 suggests that in cases where GC systems exhibit clear differences between their optical and optical-NIR color distributions (as in some giant ellipticals), the apparent inconsistencies most likely result from nonlinearities in the color-metallicity relations.

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

The distance of NGC1316, the brightest galaxy in Fornax, is an interesting test for the cosmological distance scale. First, because Fornax is the 2nd largest cluster of galaxies at <~25 Mpc after Virgo and, in contrast to Virgo, has a small line-of-sight depth; and second, because NGC1316 is the galaxy with the largest number of detected SNeIa, giving the opportunity to test the consistency of SNeIa distances internally and against other indicators. We measure SBF mags in NGC1316 from ground and space-based imaging data, providing a homogeneous set of measurements over a wide wavelength interval. The SBF, coupled with empirical and theoretical calibrations, are used to estimate the distance to the galaxy. We present the first B-band SBF measurements of NGC1316 and use them together with the optical and near-IR SBF data to analyze the properties of field stars. Our distance modulus m-M=31.59 +-0.05(stat) +-0.14(sys), when placed in a consistent Cepheid distance scale, agrees with the results from other indicators. However, our result is ~17% larger than the most recent estimate based on SNeIa. Possible explanations for this disagreement are the uncertainties on internal extinction, or calibration issues. Concerning the stellar population analysis, we confirm earlier results from other indicators: the field stars in NGC1316 are dominated by a solar metallicity, intermediate age component. A substantial mismatch exists between B-band SBF models and data, a behavior that can be accounted for by an enhanced percentage of hot horizontal branch stars. Our study of the SBF distance to NGC1316, and the comparison with distances from other indicators, raises some concern about the homogeneity between the calibrations of different indicators. If not properly placed in the same reference scale, significant differences can occur, with dramatic impact on the cosmological distance ladder.

We predict the flux and surface velocity perturbations produced by convectively excited gravity modes (g-modes) in main sequence stars. Core convection in massive stars can excite g-modes to sufficient amplitudes to be detectable with high precision photometry by Kepler and CoRoT, if the thickness of the convective overshoot region is < 30 per cent of a pressure scale height. The g-modes manifest as excess photometric variability, with amplitudes of ~ 10 micromagnitudes at frequencies ~ 10 microHz (0.8 1/d) near the solar metallicity zero-age main sequence. The flux variations are largest for stars with M > 5 solar masses, but are potentially detectable down to M ~ 2 - 3 solar masses. During the main sequence evolution, radiative damping decreases such that ever lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 100 micromagnitudes at the terminal-age main sequence. Using the same convective excitation model, we confirm previous predictions that solar g-modes produce surface velocity perturbations of < 0.3 mm/s. This implies that stochastically excited g-modes are more easily detectable in the photometry of massive main sequence stars than in the Sun.

We present HST/STIS time-series spectroscopy of the central star of the "Cat's Eye" planetary nebula NGC 6543. Intensive monitoring of the UV lines over a 5.8 hour period reveals well defined details of large-scale structure in the fast wind, which are exploited to provide new constraints on the rotation rate of the central star. We derive characteristics of the line profile variability that support a physical origin due to co-rotating interaction regions (CIRs) that are rooted at the stellar surface. The recurrence time of the observed spectral signatures of the CIRs is used to estimate the rotation period of the central star and, adopting a radius between 0.3 and 0.6 \Rsun constrains the rotational velocity to the range 54 \leq v_{rot} \leq 108 \kms. The implications of these results for single star evolution are discussed based on models calculated here for low-mass stars. Our models predict a sub-surface convective layer in NGC 6543 which we argue to be causally connected to the occurrence of structure in the fast wind.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

Weak magnetic fields have recently been detected in Vega and Sirius. Here, we explore the possibility that these fields are the remnants of some field inherited or created during or shortly after star formation and, unlike true fossil fields, are still evolving as we observe them. The timescale of this evolution is given in terms of the Alfven timescale and the rotation frequency by tau_evol ~ tau_A^2 Omega, which would be comparable to the age of the star. It is shown that it is likely that all intermediate- and high-mass stars contain fields of at least the order of the strength found so far in Vega and Sirius. Faster rotators are expected to have stronger magnetic fields. Stars may experience an increase in field strength during their early main-sequence, but for most of their lives field strength will decrease slowly. The length scale of the magnetic structure on the surface may be small in very young stars but should quickly increase to at least very approximately a fifth of the stellar radius.

The interpretation that bimodal colour distributions of globular clusters (GCs) reflect bimodal metallicity distributions has been challenged. Non-linearities in the colour to metallicity conversions caused by the horizontal branch (HB) stars may be responsible for transforming a unimodal metallicity distribution into a bimodal (optical) colour distribution. We study optical/near-infrared (NIR) colour distributions of the GC systems in 14 E/S0 galaxies. We test whether the bimodal feature, generally present in optical colour distributions, remains in the optical/NIR ones. The latter colour combination is a better metallicity proxy than the former. We use KMM and GMM tests to quantify the probability that different colour distributions are better described by a bimodal, as opposed to a unimodal distribution. We find that double-peaked colour distributions are more commonly seen in optical than in optical/NIR colours. For some of the galaxies where the optical (g-z) distribution is clearly bimodal, the (g-K) and (z-K) distributions are better described by a unimodal distribution. The two most cluster-rich galaxies in our sample, NGC4486 and NGC4649, show some interesting differences. The (g-K) distribution of NGC4649 is better described by a bimodal distribution, while this is true for the (g-K) distribution of NGC4486 GCs only if restricted to a brighter sub-sample with small K-band errors (< 0.05 mag). Formally, the K-band photometric errors cannot be responsible for blurring bimodal metallicity distributions to unimodal (g-K) colour distributions. However, simulations including the extra scatter in the colour-colour diagrams (not fully accounted for in the photometric errors) show that such scatter may contribute to the disappearance of bimodality in (g-K) for the full NGC4486 sample. For the less cluster-rich galaxies results are inconclusive due to poorer statistics. [Abridged]

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively.

The Cepheid mass discrepancy, the difference between mass estimates from stellar pulsation and stellar evolution models, is a long standing challenge for the understanding of stellar astrophysics. We discuss the current state of the mass discrepancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass loss is a significant contributor to the mass discrepancy, but it is not clear if the metallicity dependence of Cepheid mass loss is consistent with the measured metallicity dependence of the mass discrepancy.

Hot luminous stars show a variety of phenomena in their photospheres and winds which still lack clear physical explanation. Among these phenomena are photospheric turbulence, line profile variability (LPV), non-thermal emission, non-radial pulsations, discrete absorption components (DACs) and wind clumping. Cantiello et al. (2009) argued that a convection zone close to the stellar surface could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity associated with iron-group elements and is referred to as the "iron convection zone" (FeCZ). Assuming dynamo action producing magnetic fields at equipartition in the FeCZ, we investigate the occurrence of subsurface magnetism in OB stars. Then we study the surface emergence of these magnetic fields and discuss possible observational signatures of magnetic spots. Simple estimates are made using the subsurface properties of massive stars, as calculated in 1D stellar evolution models. We find that magnetic fields of sufficient amplitude to affect the wind could emerge at the surface via magnetic buoyancy. While at this stage it is difficult to predict the geometry of these features, we show that magnetic spots of size comparable to the local pressure scale height can manifest themselves as hot, bright spots. Localized magnetic fields could be widespread in those early type stars that have subsurface convection. This type of surface magnetism could be responsible for photometric variability and play a role in X-ray emission and wind clumping.

We study the two galaxies NGC4621 and NGC4374 in the Virgo cluster to derive their distances and stellar population properties. The targets have hosted three type Ia Supernova events allowing to investigate the correlations between the SNeIa and their host stellar systems. Using deep BVR data, obtained with FORS2 at the VLT, we analyse the Surface Brightness Fluctuations (SBF) properties of the targets. We adopt our measurements and existing calibrations to estimate the distance of NGC4621 and NGC4374. For stellar population analysis, we measured SBF amplitudes in different galaxy regions. We present a detailed comparison between data and models to constrain the characteristics of the dominant stellar components at i) various galactic radii, and ii) in the regions where SNeIa were recorded. Our V and R SBF measures provide distances in agreement with literature estimates. The median of our and literature SBF-based distances agrees with the one from non-SBF methods. Comparing data with models we find that stellar populations properties do not change significantly along galactic radius, with a dominant population having old age and solar chemical composition. The galaxies appear similar in all properties analysed, except for B-band SBF. Since the SBF magnitudes in this band are sensitive to the properties of a hot stellar component, we speculate that such behaviour is a consequence of different diffuse hot components in the galaxies. We find that the presence of a percentage of hot-HB stars in old and metal rich stellar populations could be at the origin of the observed differences. We find a good uniformity in the V and R SBF and integrated colours in the regions where the three SNeIa exploded. On the other hand, the B-band SBF signal shows intriguing differences.

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass discrepancy. Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in the stellar models. Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the Cepheid mass discrepancy.

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

We explore the role of mass loss and convective core overshoot in the evolution of Classical Cepheids. Stellar evolution models are computed with a recipe for pulsation-driven mass loss and it is found that mass loss alone is unable to account for the long-standing Cepheid mass discrepancy. However, the combination of mass loss and moderate convective core overshooting does provide a solution, bringing the amount of convective core overshooting in Cepheids closer to that found for other stars.

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

Rotation can have severe consequences for the evolution of massive stars. It is now considered as one of the main parameters, alongside mass and metallicity that determine the final fate of single stars. In massive, fast rotating stars mixing processes induced by rotation may be so efficient that helium produced in the center is mixed throughout the envelope. Such stars evolve almost chemically homogeneously. At low metallicity they remain blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into progenitors of long gamma-ray bursts. In binaries this type of evolution may occur because of (I) tides in very close binaries, as a result of (II) spin up by mass transfer, as result of (III) a merger of the two stars and (IV) when one of the components in the binary was born with a very high initial rotation rate. As these stars stay compact, the evolutionary channels are very different from what classical binary evolutionary models predict. In this contribution we discuss examples of nearly chemically homogeneous evolution in very close tidally-locked binaries. Even in such very close massive binaries, the stars may remain compact and avoid mass transfer, while Roche lobe overflow and a merger would be inevitable in the classical picture. This type of evolution may provide an alternative path to form tight Wolf-Rayet binaries and massive black hole binaries.

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

We measure surface brightness fluctuation (SBF) magnitudes in the F814W filter and g-I colors for nine bright early-type Fornax cluster galaxies imaged with the Hubble Space Telescope Advanced Camera for Surveys (ACS). The goal is to achieve the first systematic SBF calibration for the ACS/F814W bandpass. Because of its much higher throughput, F814W is more efficient for SBF studies of distant galaxies than the ACS/F850LP bandpass that has been used to study nearby systems. Over the color range spanned by the sample galaxies, 1.06<g-I<1.32 (AB mag), the dependence of SBF magnitude mbar_I on g-I is linear to a good approximation, with slope $\sim2$. When the F850LP SBF distance measurements from the ACS Fornax Cluster Survey are used to derive absolute Mbar_I magnitudes, the dependence on g-I becomes extremely tight, with a slope of $1.8\pm0.2$ and scatter of 0.03 mag. The small observed scatter indicates both that the estimated random errors are correct, and that the intrinsic deviations from the SBF-color relation are strongly correlated between the F814W and F850LP bandpasses, as expected. The agreement with predictions from stellar population models is good, both in slope and zero point, indicating that our mean Fornax distance of 20 Mpc is accurate. The models predict curvature in the relation beyond the color limits of our sample; thus, the linear calibration should not be extrapolated naively. In the Appendices, we reconsider the Tonry ground-based and Jensen NICMOS SBF distance catalogues; we provide a correction formula to ameliorate the small apparent bias in the former, and the offset needed to make the latter consistent with other SBF studies. We also tabulate two new SBF distances to galaxies observed in the ACS Virgo Cluster Survey.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

Thermohaline mixing has recently been proposed to occur in low-mass red giants, with large consequence for the chemical yields of low-mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We use a stellar evolution code which includes rotational mixing, internal magnetic fields and thermohaline mixing. We confirm that during the red giant stage, thermohaline mixing has the potential to decrease the abundance of ^3He which is produced earlier on the main sequence. In our models we find that this process is working on the RGB only in stars with initial mass M \simle 1.5Msun. Moreover we report that thermohaline mixing is present also during core helium burning and beyond, and has the potential to change the surface abundances of AGB stars. While we find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants. To explain the surface abundances observed at the bump in the luminosity function, the speed of the mixing process needs to be more than two orders of magnitude higher than in our models. However it is not clear if thermohaline mixing is the only physical process responsible for these surface abundance anomalies. Therefore, at this stage, it is not possible to calibrate the efficiency of thermohaline mixing against the observations.

Pulsations driven by partial ionization of hydrogen in the envelope are often considered important for driving winds from red supergiants (RSGs). In particular, it has been suggested by some authors that the pulsation growth rate in a RSG can be high enough to trigger an unusually strong wind (or a super-wind), when the luminosity to mass ratio becomes sufficiently large. Using both hydrostatic and hydrodynamic stellar evolution models with initial masses ranging from 15 to 40 \Msun, we investigate 1) how the pulsation growth rate depends on the global parameters of supergiant stars, and 2) what would be the consequences of a pulsation-driven super-wind, if it occurred, for the late stages of massive star evolution. We suggest that such a super-wind history would be marked by a runaway increase, followed by a sudden decrease, of the winds mass loss rate. The impact on the late evolution of massive stars would be substantial, with stars losing a huge fraction of their H-envelope even with a significantly lower initial mass than previously predicted. This might explain the observed lack of Type II-P supernova progenitors having initial mass higher than about 17 \Msun. We also discuss possible implications for a subset of Type IIn supernovae.

We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance in massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence stars.

Two recent empirical developments in the study of extragalactic globular cluster (GC) populations are the color-magnitude relation of the blue GCs (the "blue tilt") and the nonlinearity of the dependence of optical GC colors on metallicity. The color-magnitude relation, interpreted as a mass-metallicity relation, is thought to be a consequence of self-enrichment. Nonlinear color-metallicity relations have been shown to produce bimodal color distributions from unimodal metallicity distributions. We simulate GC populations including both a mass-metallicity scaling relation and nonlinear color-metallicity relations motivated by theory and observations. Depending on the assumed range of metallicities and the width of the GC luminosity function (GCLF), we find that the simulated populations can have bimodal color distributions with a "blue tilt" similar to observations, even though the metallicity distribution appears unimodal. The models that produce these features have the relatively high mean GC metallicities and nearly equal blue and red peaks characteristic of giant elliptical galaxies. The blue tilt is less apparent in the models with metallicities typical of dwarf ellipticals; the narrower GCLF in these galaxies has an even bigger effect in reducing the significance of their color-magnitude slopes. We critically examine the evidence for nonlinearity versus bimodal metallicities as explanations for the characteristic double-peaked color histograms of giant ellipticals and conclude that the question remains open. We discuss the prospects for further theoretical and observational progress in constraining the models presented here and for uncovering the true metallicity distributions of extragalactic globular cluster systems.

The formation of close massive black-hole binaries is a challenge for binary evolutionary models, especially the intriguing system M33 X-7 which harbours one of the most massive stellar-mass black holes (16 solar masses) orbiting a 70 solar mass O-star every 3.5 days. In standard binary evolution theory an episode of mass transfer or common envelope is inevitable in a binary with such a small orbital period, which complicates the formation of a black hole with such a high mass. To explain this system, we discuss a new binary evolution channel (De Mink et al. 2009), in which rotational mixing plays an important role. In very massive close binaries, tides force the rotation rate of the stars to be so high that rotationally induced mixing becomes very efficient. Helium produced in the center is mixed throughout the envelope. Instead of expanding during their mainsequence evolution (with the inevitable consequence of mass transfer), these stars stay compact, and avoid filling their Roche lobe. They gradually evolve into massive helium stars. This scenario naturally leads to the formation of very massive black holes in a very close orbit with a less evolved massive companion such as M33 X-7.

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

Aims.We study the photometric and structural properties of the star cluster system in the late type Sc spiral NGC 3370. Methods. BVI observations from the Advanced Camera for Surveys on board of HST are used to analyse in detail the colours, magnitudes and spatial properties of cluster candidates. The final catalogue of sources used for the study is composed by 277 objects. Results. The colour distributions of cluster candidates appear multi-modal. Although firm age constraints need the use of more age sensitive indicators, the comparison of cluster candidate colours with the colours of Galactic and Magellanic Clouds star clusters, suggests an age difference between the various sub-peaks, with a red old sub-system, a rich population of intermediate age (~1 Gyr), and a blue tail of very young (below ~100 Myr) clusters. The luminosity functions appear normal for this type of galaxy, as for the distribution of cluster effective radii (Reff). Our analysis suggests the presence of a peak in the Reff distribution at ~3.0 pc, with blue (likely young) cluster candidates showing smaller radii respect to red (likely old) objects. Finally, inspecting the properties of candidate globular clusters, we find a colour distribution matching with Galactic Globulars, with a median [Fe/H] -1.5 dex, though a non negligible tail towards lower metallicities is also present.

We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\Lsun$, $10^{3.9}\Lsun$, and $10^{4.2}$\Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.

Models of rotating single stars can successfully account for a wide variety of observed stellar phenomena, such as the surface enhancements of N and He. However, recent observations have questioned the idea that rotational mixing is the main process responsible for the surface enhancements, emphasizing the need for a strong and conclusive test. We investigate the consequences of rotational mixing for massive main-sequence stars in short-period binaries. In these systems the tides spin up the stars to rapid rotation. We use a state-of-the-art stellar evolution code including the effect of rotational mixing, tides, and magnetic fields. We discuss the surface abundances expected in massive close binaries (M1~20 solar masses) and we propose using such systems to test the concept of rotational mixing. As these short-period binaries often show eclipses, their parameters can be determined with high accuracy, allowing for a direct comparison with binary evolution models. In more massive close systems (M1~50 solar masses, Porb<~2 days) we find that helium is efficiently mixed throughout the envelope. The star remains blue and compact during the main-sequence phase. It stays within its Roche lobe while it gradually becomes a helium star. It is the less massive star, in which the effects of rotational mixing are less pronounced, which fills its Roche lobe first. We propose that this evolution path provides an alternative channel for the formation of tight Wolf-Rayet binaries with a main-sequence companion and might explain massive black hole binaries such as the intriguing system M33 X-7.

We present some results of an observational and theoretical study on unresolved stellar systems based on the Surface Brightness Fluctuations (SBF) technique. It is shown that SBF magnitudes are a valuable tracer of stellar population properties, and a reliable distance indicator. SBF magnitudes, SBF-colors, and SBF-gradients can help to constrain within relatively narrow limits the metallicity and age of the dominant stellar component in distant stellar systems, especially if coupled with other spectro-photometric indicators.

We measure F814W Surface Brightness Fluctuations (SBFs) for a sample of distant shell galaxies observed with the Advanced Camera for Survey (ACS) on board of HST. To evaluate the distance at galaxies, theoretical SBF magnitudes for the ACS@HST filters are computed for single burst stellar populations covering a wide range of ages (t=1.5-14 Gyr) and metallicities (Z=0.008-0.04). Using these stellar population models we provide the first M_SBF,F814W versus (F475W-F814W)0 calibration. The results suggest that optical SBFs can be measured at d>100 Mpc using high resolution spatial optical data.

We present some results on the study of stellar population properties and distances of galaxies using the SBF technique. The applications summarized here show that the Surface Brightness Fluctuations (SBF) method is able to i) provide accurate distances of resolved and unresolved stellar systems from ~10 Kpc to ~150 Mpc, and ii) to reliably constrain the physical properties (e.g. age and metallicity) of unresolved stellar systems.

One of the main uncertainties in evolutionary calculations of massive stars is the efficiency of internal mixing. It changes the chemical profile inside the star and can therefore affect the structure and further evolution. We demonstrate that eclipsing binaries, in which the tides synchronize the rotation period of the stars and the orbital period, constitute a potentially strong test for the efficiency of rotational mixing. We present detailed stellar evolutionary models of massive binaries assuming the composition of the Small Magellanic Cloud. In these models we find enhancements in the surface nitrogen abundance of up to 0.6 dex.

We present results from the first extensive study of convection zones in the envelopes of hot massive stars, which are caused by opacity peaks associated with iron and helium ionization. These convective regions can be located very close to the stellar surface. Recent observations of microturbulence in massive stars from the VLT-Flames survey are in good agreement with our predictions concerning the occurrence and the strength of sub-surface convection in hot stars. We argue further that convection close to the surface may trigger clumping at the base of the stellar wind of massive stars.

Thermohaline mixing has recently been proposed to occur in low mass red giants, with large consequences for the chemical yields of low mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We confirm that thermohaline mixing has the potential to destroy most of the ^3He which is produced earlier on the main sequence during the red giant stage. In our models we find that this process is working only in stars with initial mass M <~ 1.5Msun. Moreover, we report that thermohaline mixing can be present during core helium burning and beyond in stars which still have a ^3He reservoir. While rotational and magnetic mixing is negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with differential rotation and magnetic fields may be essential to establish the time scale of thermohaline mixing in red giants.

Rotational mixing is a very important but uncertain process in the evolution of massive stars. We propose to use close binaries to test its efficiency. Based on rotating single stellar models we predict nitrogen surface enhancements for tidally locked binaries. Furthermore we demonstrate the possibility of a new evolutionary scenario for very massive (M > 40 solar mass) close (P < 3 days) binaries: Case M, in which mixing is so efficient that the stars evolve quasi-chemically homogeneously, stay compact and avoid any Roche-lobe overflow, leading to very close (double) WR binaries.

We review the role of rotation in massive close binary systems. Rotation has been advocated as an essential ingredient in massive single star models. However, rotation clearly is most important in massive binaries where one star accretes matter from a close companion, as the resulting spin-up drives the accretor towards critical rotation. Here, we explore our understanding of this process, and its observable consequences. When accounting for these consequences, the question remains whether rotational effects in massive single stars are still needed to explain the observations.

We measure F814W Surface Brightness Fluctuations (SBF) for a sample of distant shell galaxies with radial velocities ranging from 4000 to 8000 km/s. The distance at galaxies is then evaluated by using the SBF method. For this purpose, theoretical SBF magnitudes for the ACS@HST filters are computed for single burst stellar populations covering a wide range of ages (t=1.5-14 Gyr) and metallicities (Z=0.008-0.04). Using these stellar population models we provide the first $\bar{M}_{F814W}$ versus $(F475W-F814W)_0$ calibration and we extend the previous I-band versus $(B-I)_0$ color relation to colors $(B-I)_{0}\leq 2.0$ mag. Coupling our SBF measurements with the theoretical calibration we derive distances with a statistical uncertainty of $\sim 8%$, and systematic error of $\sim 6 %$. The procedure developed to analyze data ensures that the indetermination due to possible unmasked residual shells is well below $\sim 12 %$. The results suggest that \emph{optical} SBFs can be measured at $d \geq 100 Mpc$ with ACS@HST imaging. SBF-based distances coupled with recession velocities corrected for peculiar motion, allow us obtain $H_{0} = 76 \pm 6$ (statistical) $\pm 5$ (systematic) km/s/Mpc.

We discuss recent models on the evolution of massive stars at very low metallicity including the effects of rotation, magnetic fields and binarity. Very metal poor stars lose very little mass and angular momentum during the main sequence evolution, and rotation plays a dominant role in their evolution. In rapidly rotating massive stars, the rotationally induced mixing time scale can be even shorter than the nuclear time scale throughout the main sequence. The consequent quasi-chemically homogeneous evolution greatly differs from the standard massive star evolution that leads to formation of red giants with strong chemical stratification. Interesting outcomes of such a new mode of evolution include the formation of rapidly rotating massive Wolf-Rayet stars that emit large amounts of ionizing photons, the formation of a long gamma-ray bursts and a hypernovae, and the production of large amounts of primary nitrogen. We show that binary interactions can further enhance the effects of rotation, as mass accretion in a close binary spins up the secondary.

We discuss how rotation and binary interactions may be related to the diversity of type Ibc supernovae and long gamma-ray bursts. After presenting recent evolutionary models of massive single and binary stars including rotation, the Tayler-Spruit dynamo and binary interactions, we argue that the nature of SNe Ibc progenitors from binary systems may not significantly differ from that of single star progenitors in terms of rotation, and that most long GRB progenitors may be produced via the quasi-chemically homogeneous evolution at sub-solar metallicity. We also briefly discuss the possible role of magnetic fields generated in the convective core of a massive star for the transport of angular momentum, which is potentially important for future stellar evolution models of supernova and GRB progenitors.

We perform a series of numerical experiments to study how the nonlinear metallicity--color relations predicted by different stellar population models affect the color distributions observed in extragalactic globular cluster systems. % We present simulations in the $UBVRIJHK$ bandpasses based on five different sets of simple stellar population (SSP) models. The presence of photometric scatter in the colors is included as well. % We find that unimodal metallicity distributions frequently ``project'' into bimodal color distributions. The likelihood of this effect depends on both the mean and dispersion of the metallicity distribution, as well as of course on the SSP model used for the transformation. % Adopting the Teramo-SPoT SSP models for reference, we find that optical--to--near-IR colors should be favored with respect to other colors to avoid the bias effect in globular cluster color distributions discussed by \citet{yoon06}. In particular, colors such as \vh\ or \vk are more robust against nonlinearity of the metallicity--color relation, and an observed bimodal distribution in such colors is more likely to indicate a true underlying bimodal metallicity distribution. Similar conclusions come from the simulations based on different SSP models, although we also identify exceptions to this result.

The collapsar model for long gamma-ray bursts requires a rapidly rotating Wolf-Rayet star as progenitor. We test the idea of producing rapidly rotating Wolf-Rayet stars in massive close binaries through mass accretion and consecutive quasi-chemically homogeneous evolution; the latter had previously been shown to provide collapsars below a certain metallicity threshold for single stars. The binary channel presented here may provide a means for massive stars to obtain the high rotation rates required to evolve quasi-chemically homogeneous and fulfill the collapsar scenario. Moreover, it suggests that a possibly large fraction of long gamma-ray bursts occurs in runaway stars.

Thermohaline mixing has recently been proposed to occur in low mass red giants, with large consequence for the chemical yields of low mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1 Msun and 3 Msun. We use a stellar evolution code which includes rotational mixing and internal magnetic fields. We confirm that thermohaline mixing has the potential to destroy most of the helium 3 which is produced earlier on the main sequence during the red giant stage, in stars below 1.5Msun. We find this process to continue during core helium burning and beyond. We find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, even if the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants.

Pair creation supernovae (PCSN) are thought to be produced from very massive low metallicity stars. The spectacularly bright SN 2006gy does show signatures expected from PCSNe. Here, we investigate the metallicity threshold below which PCSN can form and estimate their occurrence rate. We perform stellar evolution calculations for stars of 150$\mso$ and 250$\mso$ of low metallicity (Z$_{\odot}$/5 and Z$_{\odot}$/20), and analyze their mass loss rates. We find that the bifurcation between quasi-chemically homogeneous evolution for fast rotation and conventional evolution for slower rotation, which has been found earlier for massive low metallicity stars, persists in the mass range considered here. Consequently, there are two separate PCSN progenitor types: (I) Fast rotators produce PCSNe from very massive Wolf-Rayet stars, and (II) Slower rotators that generate PCSNe in hydrogen-rich massive yellow hypergiants. We find that hydrogen-rich PCSNe could occur at metallicities as high as Z$_{\odot}$/3, which -- assuming standard IMFs are still valid to estimate their birth rates -- results in a rate of about one PCSN per 1000 supernovae in the local universe, and one PCSN per 100 supernovae at a redshift of $z=5$. PCSNe from WC-type Wolf-Rayet stars are restricted to much lower metallicity.

We derive Surface Brightness Fluctuations (SBF) and integrated magnitudes in the V- and I-bands using Advanced Camera for Surveys (ACS) archival data. The sample includes 14 galaxies covering a wide range of physical properties: morphology, total absolute magnitude, integrated color. We take advantage of the latter characteristic of the sample to check existing empirical calibrations of absolute SBF magnitudes both in the I- and V-passbands. Additionally, by comparing our SBF and color data with the Teramo-SPoT simple stellar population models, and other recent sets of population synthesis models, we discuss the feasibility of stellar population studies based on fluctuation magnitudes analysis. The main result of this study is that multiband optical SBF data and integrated colors can be used to significantly constrain the chemical composition of the dominant stellar system in the galaxy, but not the age in the case of systems older than 3 Gyr. SBF color gradients are also detected and analyzed. These SBF gradient data, together with other available data, point to the existence of mass dependent metallicity gradients in galaxies, with the more massive objects showing a non--negligible SBF versus color gradient. The comparison with models suggests that such gradients imply more metal rich stellar populations in the galaxies' inner regions with respect to the outer ones.

Recent applications have proved that the Surface Brightness Fluctuations (SBF) technique is a reliable distance indicator in a wide range of distances, and a promising tool to analyze the physical and chemical properties of unresolved stellar systems, in terms of their metallicity and age. We present the preliminary results of a project aimed at studying the evolutionary properties and distance of the stellar populations in external galaxies based on the SBF method. On the observational side, we have succeeded in detecting I-band SBF gradients in six bright ellipticals imaged with the ACS, for these same objects we are now presenting also B-band SBF data. These B-band data are the first fluctuations magnitude measurements for galaxies beyond 10 Mpc. To analyze the properties of stellar populations from the data, accurate SBF models are essential. As a part of this project, we have evaluated SBF magnitudes from Simple Stellar Population (SSP) models specifically optimized for the purpose. A wide range of chemical compositions and ages, as well as different choices of the photometric system have been investigated. All models are available at the Teramo-Stellar Populations Tools web site: www.oa-teramo.inaf.it/SPoT.

We introduce a procedure developed by the ``Teramo Stellar Populations Tools'' group (Teramo-SPoT), specifically optimized to obtain realistic simulations of CCD images of elliptical galaxies. Particular attention is devoted to include the Surface Brightness Fluctuation (SBF) signal observed in ellipticals and to simulate the Globular Cluster (GC) system in the galaxy, and the distribution of background galaxies present in real CCD frames. In addition to the physical properties of the simulated objects - galaxy distance and brightness profile, luminosity function of GC and background galaxies, etc. - the tool presented allows the user to set some of the main instrumental properties - FoV, zero point magnitude, exposure time, etc.

Taking advantage of the exceptional capabilities of ACS on board of HST, we derive Surface Brightness Fluctuation (SBF) measurements in the B and I bands from images of six elliptical galaxies with $1500 \leq cz \leq 3500$. Given the low S/N ratio of the SBF signal in the blue band images, the reliability of the measurements is verified both with numerical simulations and experimental data tests. This paper presents the first published B- and I-band SBF measurements for distant ($\geq$ 20 Mpc) galaxies, essential for the comparisons of the models to observations of normal ellipticals. By comparing I-band data with our new Simple Stellar Population (SSP) models we find an excellent agreement and we confirm that I-band SBF magnitudes are mainly sensitive to the metallicity of the dominant stellar component in the galaxy, and are not strongly affected by the contribution of possible secondary stellar components. As a consequence I-band fluctuations magnitudes are ideal for distance studies. On the other hand, we show that standard SSP models do not reproduce the B-band SBF magnitudes of red ($(B-I)_0 \gsim 2.1$) galaxies in our sample. We explore the capability of two non--canonical models in properly reproducing the high sensitivity of B SBF to the presence of even small fractions of bright, hot stars (metal poor stars, hot evolved stars, etc.). The disagreement is solved both by taking into account hot (Post--AGB) stars in SSP models and/or by adopting Composite Stellar Population models. Finally, we suggest a limit value of the S/N for the B-band SBF signal required to carry out a detailed study of stellar population properties based on this technique.

We present a binary channel for the progenitors of long gamma-ray bursts. We test the idea of producing rapidly rotating Wolf-Rayet stars in massive close binaries through mass accretion and consecutive quasi-chemically homogeneous evolution. The binary channel presented here may provide a means for massive stars to obtain the required high rotation rates. Moreover, it suggests that a possibly large fraction of long gamma-ray bursts occurs in runaway stars. This can have important observational consequences for both the positions of GRBs, and their afterglow properties.

We present theoretical surface-brightness fluctuations (SBF) amplitudes for single-burst stellar populations of young and intermediate age ($25 Myr \leq t \leq 5 Gyr$), and $Z=0.0003$, 0.001, 0.004, 0.008, 0.01, 0.02, and 0.04. We pay attention to the contribution of thermal-pulses asymptotic giant branch (TP-AGB) stars. The sensitivity of the predicted SBF to changes in the mass-loss scenario along the TP-AGB phase is examined. Below 0.6-$1 Gyr$ both optical and NIR SBF models exhibit a strong dependence on age and mass loss. We also evaluate SBF amplitudes by using Monte Carlo techniques to reproduce the random variation in the number of stars experiencing bright and fast evolutionary phases (Red Giant Branch, AGB, TP-AGB). We provide constraints upon the faintest integrated flux of real stellar populations required to derive reliable and meaningful SBF measurements. We analyze a technique for deriving SBF amplitudes of star clusters from the photometry of individual stars, and estimate the indetermination due to statistical effects, which may impinge on the procedure. The first optical SBF measurements for 11 Large Magellanic Cloud (LMC) star rich clusters - with age ranging from few $Myr$ to several $Gyr$ - are derived by using Hubble Space Telescope observations. The measurements are successfully compared to our SBF predictions, thus providing a good agreement with models of metallicity $Z = 0.0001$-0.01. A mass loss as a power-law function of the TP-AGB star luminosity is required in order to properly reproduce the optical SBF data of the LMC clusters. We suggest how to overcome the general problem of SBF models in reproducing the details of the near-infrared SBF measurements of the Magellanic Cloud star clusters.

We study surface brightness fluctuations (SBF) in a sample of 8 elliptical galaxies using Advanced Camera for Surveys (ACS) Wide Field Channel (WFC) data drawn from the Hubble Space Telescope (HST) archive. SBF magnitudes in the F814W bandpass, and galaxy colors from F814W, F435W, and F606W images -- when available -- are presented. Galaxy surface brightness profiles are determined as well. We present the first SBF--broadband color calibration for the ACS/WFC F814W bandpass, and (relative) distance moduli estimates for 7 of our galaxies. We detect and study in detail the SBF variations within individual galaxies as a probe of possible changes in the underlying stellar populations. Inspecting both the SBF and color gradients in comparison to model predictions, we argue that SBF, and SBF-gradients, can in principle be used for unraveling the different evolutionary paths taken by galaxies, though a more comprehensive study of this issue would be required. We confirm that the radial variation of galaxy stellar population properties should be mainly connected to the presence of radial chemical abundance gradients, with the outer galaxy regions being more metal poor than the inner ones.

It has been suggested that a monitoring program employing the pixel lensing method to search for microlensing events towards galactic globular clusters may increase the statistics and discriminate among different halo models. Stimulated by this proposal, we evaluate an upper limit to the pixel lensing event rate for such a survey. Four different dark halo models have been considered changing both the flattening and the slope of the mass density profile. The lenses mass function has been modelled as a homogenous power - law for $\mu \in (\mu_l, \mu_u)$ and both the mass limits and the slope of the mass function have been varied to investigate their effect on the rate. The target globular clusters have been selected in order to minimize the disk contribution to the event rate. We find that a pixel lensing survey towards globular clusters is unable to discriminate among different halo models since the number of detectable events is too small to allow any reliable statistical analysis.

We present new theoretical models for Surface Brightness Fluctuations (SBF) both for optical and near-infrared bands in standard ground-based and Hubble Space Telescope filter systems. Simple Stellar Population simulations are adopted. Models cover the age and metallicity ranges from $t=5$ to $15~Gyr$ and from $Z=0.0001$ to 0.04 respectively. Effects due to the variation of the Initial Mass Function and the stellar color-temperature relations are explored. Particular attention is devoted to very bright stars in the color-magnitude diagram and to investigate the effects of mass loss along the Red Giant Branch (RGB) and the Asymptotic Giant Branch (AGB). It is found that $U$ and $B$ bands SBF amplitudes are powerful diagnostics for the morphology of the Horizontal Branch and the Post-AGB stars population. We point out that a careful treatment of mass loss process along the RGB and AGB is fundamental in determining reliable SBF evaluations. The SBF measurements are used to give robust constraints on the evolution of AGB stars, suggesting that mass loss activity on AGB stars should be twice more efficient than on the RGB stars. Our models are able to reproduce the absolute SBF magnitudes of the Galactic Globular Clusters and of galaxies, and their integrated colors. New calibrations of absolute SBF magnitude in $V$, $R$, $I$, and $K$ photometric filters are provided, which appear reliable enough to directly gauge distances bypassing other distance indicators. The SBF technique is also used as stellar population tracer to derive age and metallicity of a selected sample of galaxies of known distances. Finally, {\it SBF color} versus {\it integrated color} diagrams are proposed as particularly useful in removing the well known {\it age-metallicity degeneracy} affecting our knowledge of remote stellar systems.

We present new theoretical evaluations of optical and near-IR Surface Brightness Fluctuations (SBF) magnitudes for single-burst stellar populations in the age range t=5-15 Gyr and metallicity from $Z_{\sun}/200$ to $2Z_{\sun}$. Our theoretical predictions can be successfully used to derive reliable distance evaluations. They also appear to be a new and valuable tool to trace the properties of unresolved stellar populations.