We present the results of high-resolution (R $\ge$ 30,000) optical and near-infrared spectroscopic monitoring observations of a FU Orionis-type object, V960 Mon, which underwent an outburst in 2014 November. We have monitored this object with the Bohyunsan Optical Echelle Spectrograph (BOES) and the Immersion GRating INfrared Spectrograph (IGRINS) since 2014 December. Various features produced by a wind, disk, and outflow/jet were detected. The wind features varied over time and continually weakened after the outburst. We detected double-peaked line profiles in the optical and near-infrared, and the line widths tend to decrease with increasing wavelength, indicative of Keplerian disk rotation. The disk features in the optical and near-infrared spectra fit well with G-type and K-type stellar spectra convolved with a kernel to account for the maximum projected disk rotation velocity of about 40.3$\pm$3.8 km s$^{-1}$ and 36.3$\pm$3.9 km s$^{-1}$, respectively. We also report the detection of [S II] and H$_{2}$ emission lines, which are jet/outflow tracers and rarely found in FUors.

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescales of the order of a Hubble time, are small enough to display some collisional behaviour and big enough to be impossible to simulate with direct $N$-body codes and current hardware. We introduce a new method to simulate collisional stellar systems, and validate it by comparison with direct N-body codes on small-N simulations. The Multi-Particle collision for Dense stellar systems Code (MPCDSS) is a new code for evolving stellar systems with the Multi-Particle Collision method. Such method amounts to a stochastic collision rule that allows to conserve exactly the energy and momentum over a cluster of particles experiencing the collision. The code complexity scales with $N \log N$ in the number of particles. Unlike Monte-Carlo codes, MPCDSS can easily model asymmetric, non-homogeneous, unrelaxed and rotating systems, while allowing us to follow the orbits of individual stars. We evolve small ($N = 3.2 \times 10^4$) star clusters with MPCDSS and with the direct-summation code NBODY6, finding a similar evolution of key indicators. We then simulate different initial conditions in the $10^4 - 10^6$ star range. MPCDSS bridges the gap between small, collisional systems that can be simulated with direct $N$-body codes and large noncollisional systems. MPCDSS in principle allows us to simulate globular clusters such as Omega Cen and M54 and even the nuclear star cluster, beyond the limits of current direct N-body codes in terms of the number of particles.

We present that the spin$-$orbit alignment (SOA; i.e., the angular alignment between the spin vector of a halo and the orbital angular momentum vector of its neighbor) provides an important clue to how galactic angular momenta develop. In particular, we identify virial-radius-wise contact halo pairs with mass ratios from 1/3 to 3 in a set of cosmological $N$-body simulations, and divide them into merger and flyby subsamples according to their total (kinetic+potential) energy. In the spin$-$orbit angle distribution, we find a significant SOA in that $75.0\pm0.6$ % of merging neighbors and $58.7\pm0.6$ % of flybying neighbors are on the prograde orbit. The overall SOA of our sample is mainly driven by fast-rotating halos, corroborating that a well-aligned interaction spins a halo faster. More interestingly, we find for the first time a strong number excess of nearly perpendicular but still prograde interactions ($\sim75^{\circ}$) in the spin$-$orbit angle distribution for both the merger and flyby cases. Such prograde-polar interactions predominate for slow-rotating halos, testifying that misaligned interactions reduce the halos' spin. The frequency of the prograde-polar interactions correlates with the halo mass, yet anticorrelates with the large-scale density. This instantly invokes the spin-flip phenomenon that is conditional on the mass and environment. The prograde-polar interaction will soon flip the spin of a slow-rotator to align with its neighbor's orbital angular momentum. Finally, we propose a scenario that connects the SOA to the ambient large-scale structure based on the spin-flip argument.

We performed simultaneous observations of the H2O 6(1,6) - 5(2,3) (22.235080 GHz) and SiO v= 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1, 3 - 2 (43.122080, 42.820587, 86.243442, and 129.363359 GHz) masers towards the suspected D-type symbiotic star, V627 Cas, using the Korean VLBI Network. Here, we present astrometrically registered maps of the H2O and SiO v = 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1 masers for five epochs from January 2016 to June 2018. Distributions of the SiO maser spots do not show clear ring-like structures, and those of the H2O maser are biased towards the north-north-west to west with respect to the SiO maser features according to observational epochs. These asymmetric distributions of H2O and SiO masers are discussed based on two scenarios of a bipolar outflow and the presence of the hot companion, a white dwarf, in V627 Cas. We carried out ring fitting of SiO v = 1, and v = 2 masers and estimated the expected position of the cool red giant. The ring radii of the SiO v = 1 maser are slightly larger than those of the SiO v = 2 maser, as previously known. Our assumption for the physical size of the SiO maser ring of V627 Cas to be the typical size of a SiO maser ring radius (\sim4 au) of red giants yields the distance of V627 Cas to be \sim1 kpc.

We present an observational study of the interacting pair of dwarf galaxies, MCG+07-20-052, in the vicinity of Milky Way mass spiral galaxy NGC 2998. MCG+07-20-052 is located at a sky-projected distance of 105 kpc from NGC 2998 and the two have a relative line-of-sight velocity of 60 kms. We observed tidal tail-like extensions on both members (D1 and D2) of the interacting pair MCG+07-20-052. The interacting dwarf galaxies, D1 and D2, have B-band absolute magnitudes of $-$17.17 and $-$17.14 mag, respectively, and D2 is significantly bluer than D1. We obtained HI 21 cm line data of the NGC 2998 system using the Giant Metrewave Radio Telescope (GMRT) to get a more detailed view of the neutral hydrogen (HI) emission in the interacting dwarf galaxies and in the galaxy members of the NGC 2998 group. Evidence of a merger between the dwarf galaxies in the MCG+07-20-052 is also present in the HI kinematics and morphology where we find that the HI is mostly concentrated around D2, which also shows a higher level of star-forming activity and bluer $g-r$ color index compared to D1. In addition, we detect extended tenuous HI emission around another member galaxy, NGC 3006, located close to the MCG+07-20-052-pair at a sky-projected distance of 41 kpc. We compare here our results from the MCG+07-20-052 pair-NGC 2998 system with other known LMC-SMC-Milky Way type systems and discuss the possible origin of the dwarf-dwarf interaction.

A few dozen dwarf satellite galaxies of the Milky Way have been discovered, which are often viewed as the remaining building blocks of our Galaxy. The follow-up spectroscopy showed that dwarf galaxies have a sizeable spread in their metallicities. Several scenarios were suggested to explain the metallicity spread, which can be tested by the structural patterns of stellar subpopulations with distinct metallicities. However, such chemical plus structural examination, to which we refer to as "chemostructural study" is hindered by the lack of stars with spectroscopic metallicity. Here we propose the Ca$-$$by$ photometry as an alternative way to secure metallicities for a 2$-$3 orders of magnitude larger stellar sample than the spectroscopic sample and thus enable us to perform a chemostructural study on dwarf galaxies. In particular, we use the $hk$ index [$\equiv($Ca$-b)-(b-y)$], whose validity as a photometric metallicity indicator (and crass insensitivity to age) for red-giant-branch stars was upheld via Galactic globular clusters, and observe three dwarf spheroidal galaxies$-$Draco, Sextans, and Canes Venatici I (CVnI)$-$with Subaru/Suprime-Cam. We find that in all the galaxies the metal-rich stellar populations are more centrally concentrated than the metal-poor counterparts, suggesting that the central regions of the galaxies underwent extended star formation. Such a negative radial metallicity gradient for Sextans and CVnI opposes to the traditional spectroscopic results. We also find that their metallicity distribution functions (MDFs) can be characterized by a unimodal, skewed Gaussian shape with a metal-rich peak and a metal-poor tail. We discuss their features in the chemostructure and MDF in terms of dwarf galaxy formation theories.

Cassiopeia A is a nearby young supernova remnant that provides a unique laboratory for the study of core-collapse supernova explosions. Cassiopeia A is known to be a Type IIb supernova from the optical spectrum of its light echo, but the immediate progenitor of the supernova remains uncertain. Here we report results of near-infrared, high-resolution spectroscopic observations of Cassiopeia A where we detected the pristine circumstellar material of the supernova progenitor. Our observations revealed a strong emission line of iron (Fe) from a circumstellar clump that has not yet been processed by the supernova shock wave. A comprehensive analysis of the observed spectra, together with an HST image, indicates that the majority of Fe in this unprocessed circumstellar material is in the gas phase, not depleted onto dust grains as in the general interstellar medium. This result is consistent with a theoretical model of dust condensation in material that is heavily enriched with CNO-cycle products, supporting the idea that the clump originated near the He core of the progenitor. It has been recently found that Type IIb supernovae can result from the explosion of a blue supergiant with a thin hydrogen envelope, and our results support such a scenario for Cassiopeia A.

We present a statistical analysis of flybys of dark matter halos compared to mergers using cosmological $N$-body simulations. We mainly focus on gravitationally interacting target halos with mass of $10^{10.8}-10^{13.0}h^{-1}M_{\odot}$, and their neighbors are counted only when the mass ratio is 1:3$-$3:1 and the distance is less than the sum of the virial radii of target and neighbor. The neighbors are divided into the flyby or merger samples if the pair's total energy is greater or smaller, respectively, than the capture criterion with consideration of dynamical friction. The main results are as follows: (a) The flyby fraction increases by up to a factor of 50 with decreasing halo mass and by up to a factor of 400 with increasing large-scale density, while the merger fraction does not show any significant dependencies on these two parameters; (b) The redshift evolution of the flyby fraction is twofold, increasing with redshift at $0<z<1$ and remaining constant at $z>1$, while the merger fraction increases monotonically with redshift at $z=0\sim4$; (c) The multiple interactions with two or more neighbors are on average flyby-dominated, and their fraction has a mass and environment dependence similar to that for the flyby fraction; (d) Given that flybys substantially outnumber mergers toward $z=0$ (by a factor of five) and the multiple interactions are flyby-dominated, the flyby's contribution to galactic evolution is stronger than ever at the present epoch, especially for less massive halos and in the higher density environment. We propose a scenario that connects the evolution of the flyby and merger fractions to the hierarchical structure formation process.

Recent interpretation of the color$-$magnitude diagrams of the Milky Way (MW) bulge has suggested that the observed double red-clump feature can be a natural consequence of He-enhanced stellar populations in the MW bulge. This implies that globular clusters (GCs), where the He-enhanced second-generation (SG) stars can be efficiently created, are the most likely candidate contributors of He-rich stars to the MW bulge. We extend this idea to the Galactic inner halo and investigate the fraction of the SG stars as a function of the Galactocentric distance. We use bluer blue-horizontal branch (bBHB) stars, which are assumed to be originated from He-rich SG populations, as proxies of SG stars, and find that the fraction of bBHB stars increases with decreasing Galactocentric distance. Simulations of the GC evolution in the MW tidal field qualitatively support the observed trend of bBHB enhancement in the inner halo. In these simulations, the increasing tidal force with decreasing Galactocentric distance leads to stripping of stars not only from the outskirts but also from the central regions of GCs, where SG stars are more abundant. We discuss the implication and prospect of our findings concerning the formation history of the bulge and inner halo of the MW.

Galaxies in pairs show enhanced star formation (SF) compared to their counterparts in isolation, which is often explained by the tidal effect of neighboring galaxies. Recent observations, however, reported that galaxies paired with early-type neighbors do not undergo the SF enhancement. Here we revisit the influence of neighbors using a large sample of paired galaxies from the Sloan Digital Sky Survey and a carefully constructed control sample of isolated counterparts. We find that star-forming neighbors enhance SF, and even more so for more star-forming (and closer) neighbors, which can be attributed to collisions of interstellar medium (ISM) leading to SF. We further find that, contrary to the anticipated tidal effect, quiescent neighbors quench SF, and even more so for more quiescent (and closer) neighbors. This seems to be due to removal of gas reservoirs via ram pressure stripping and gas accretion cutoff by hot gas halos of quiescent neighbors, on top of their paucity of ISM to collide to form stars. Our findings, especially the intimate connection of SF to the status and strength of neighbors' SF, imply that the hydrodynamic mechanisms, along with the tidal effect, play a crucial role during the early phase of galactic interactions.

Helium accretion induced explosions in CO white dwarfs (WDs) are considered promising candidates for a number of observed types of stellar transients, including supernovae (SNe) of Type Ia and Type Iax. However, a clear favorite outcome has not yet emerged. We explore the conditions of helium ignition in the white dwarf and the final fates of helium star-WD binaries as function of their initial orbital periods and component masses. We compute 274 model binary systems with the Binary Evolution Code (BEC), where both components are fully resolved. Stellar and orbital evolution is computed simultaneously, including mass and angular momentum transfer, tides, and gravitational wave emission, as well as differential rotation and internal hydrodynamic and magnetic angular momentum transport. We find that helium detonations are expected only in systems with the shortest initial orbital periods, and for initially massive white dwarfs (MWD > 1.0 MSun ) and lower mass donors (Mdonor < 0.8 MSun), with accumulated helium layers mostly exceeding 0.1 MSun. Upon detonation, these systems would release the donor as a hypervelocity pre-WD runaway star, for which we predict the expected range of kinematic and stellar properties. Systems with more massive donors or initial periods exceeding 1.5 h will likely undergo helium deflagrations after accumulating 0.1 - 0.001 MSun of helium. Helium ignition in the white dwarf is avoided in systems with helium donor stars below - 0.6 MSun, and lead to three distinctly different groups of double white dwarf systems. The size of the parameter space open to helium detonation corresponds to only about 3 % of the galactic SN Ia rate, and to 10 % of the SN Iax rate, while the predicted large amounts of helium (>0.1 MSun) in progenitors cannot easily be reconciled with observations of archetypical SN Ia. ...

The Large Synoptic Survey Telescope (LSST) will open a discovery frontier for faint and fast transients with its ability to detect variable flux components down to $\sim$24.5 mag in a $\sim$30 second exposure. Spectroscopic follow-up of such phenomena - which are necessary for understanding the physics of stellar explosions - can require a rapid response and several hours with a 8-10m telescope, making it both expensive and difficult to acquire. The future Extremely Large Telescopes (ELTs) would be able to provide not only spectroscopy but capabilities such as spectropolarimetry and high-resolution diffraction-limited imaging that would contribute to future advances in our physical understanding of stellar explosions. In this white paper we focus on several specific scientific impacts in the field of explosive transient astrophysics that will be generated by the combination of LSST's discovery abilities and ELTs' follow-up capacities. First, we map the uncharted frontier of discovery phase-space in terms of intrinsic luminosity and timescales for explosive transients, where we expect the unexpected. We then focus on six areas with open science questions for known transients: the progenitors of thermonuclear supernovae (SNe), mass loss prior to core collapse, asymmetry in stellar explosions, light echoes, high-$z$ transients, and strongly lensed SNe. We conclude with a brief discussion of the practical aspects of ELT & LSST synergy.

The large aperture and sensitive optical and near infrared imager spectrographs will enable an ELT system to observe some supernovae at large distances, deep into cosmological history when supernovae first began to occur.

An ELT system with its large aperture and sensitive optical and near infrared imager spectrographs will make major contributions to the study of stars ripped apart by supermassive black holes.

Complex organic molecules (COMs), which are the seeds of prebiotic material and precursors of amino acids and sugars, form in the icy mantles of circumstellar dust grains but cannot be detected remotely unless they are heated and released to the gas phase. Around solar-mass stars, water and COMs only sublimate in the inner few au of the disk, making them extremely difficult to spatially resolve and study. Sudden increases in the luminosity of the central star will quickly expand the sublimation front (so-called snow line) to larger radii, as seen previously in the FU Ori outburst of the young star V883 Ori. In this paper, we take advantage of the rapid increase in disk temperature of V883 Ori to detect and analyze five different COMs, methanol, acetone, acetonitrile, acetaldehyde, and methyl formate, in spatially-resolved submillimeter observations. The COMs abundances in V883 Ori is in reasonable agreement with cometary values. This result suggests that outbursting young stars can provide a special opportunity to study the ice composition of material directly related to planet formation.

Intensive Monitoring Survey of Nearby Galaxies (IMSNG) is a high cadence observation program monitoring nearby galaxies with high probabilities of hosting supernovae (SNe). IMSNG aims to constrain the SN explosion mechanism by inferring sizes of SN progenitor systems through the detection of the shock-heated emission that lasts less than a few days after the SN explosion. To catch the signal, IMSNG utilizes a network of 0.5-m to 1-m class telescopes around the world and monitors the images of 60 nearby galaxies at distances D < 50 Mpc to a cadence as short as a few hours. The target galaxies are bright in near-ultraviolet (NUV) with M_NUV < -18.4 AB mag and have high probabilities of hosting SNe (0.06 SN/yr per galaxy). With this strategy, we expect to detect the early light curves of 3.4 SNe per year to a depth of R ~ 19.5 mag, enabling us to detect the shock-heated emission from a progenitor star with a radius as small as 0.1 R_sun. The accumulated data will be also useful for studying faint features around the target galaxies and other science projects. So far, 18 SNe have occurred in our target fields (16 in IMSNG galaxies) over 5 years, confirming our SN rate estimate of 0.06 SN/yr per galaxy.

We report a possible discovery of a relic star cluster in the Sextans dwarf spheroidal galaxy. Using the \textit{hk} index ($\equiv$($Ca-b$)$-$($b-y$)) as a photometric metallicity indicator, we discriminate the metal-poor and metal-rich stars in the galaxy and find unexpected number density excess of metal-poor stars located 7'.7 ($\sim$190 pc in projected distance) away from the known galactic center. The $V-I$ color$-$magnitude diagram (CMD) for stars around the density excess reveals that both the main sequence and the giant branch are considerably narrower and redder than the bulk of field stars in Sextans. Our stellar population models show (a) that the narrow CMD is best reproduced by a simple stellar population with an age of $\sim$13 Gyr and [Fe/H] of $\sim$$-$2.3 dex, and (b) that the redder $V-I$ color of the $hk$-weak population is explained $only$ if it is $\sim$2 Gyr older than the field stars. The results lead us to conclude that the off-centered density peak is likely associated with an old, metal-poor globular cluster. The larger spatial extent ($>$ 80 pc in radius) and the smaller number of stars ($\sim$1000) than typical globular clusters point to a star cluster that is in the process of dissolution. The finding serves as the first detection of a surviving star cluster in Sextans, supporting previous suggestions of the presence of star cluster remnants in the galaxy. If confirmed, the survival of a relic star cluster until now implies a $cored$ dark matter halo profile for this dwarf galaxy.

The color distributions of globular clusters (GCs) in individual early-type galaxies show great diversity in their morphology. Based on the conventional "linear" relationship between the colors and metallicities of GCs, the GC metallicity distributions inferred from colors and in turn their formation histories, should be as diverse as they appear. In contrast, here we show that an alternative scenario rooted in the "nonlinear" nature of the metallicity-to-color transformation points to a simpler and more coherent picture. Our simulations of the color distributions for $\sim$80 GC systems in early-type galaxies from the ACS Virgo and Fornax Cluster Surveys suggest that the majority ($\sim$70$\%$) of early-type galaxies have old ($\sim$13 Gyr) and coeval GCs. Their variety in the color distribution morphology stems mainly from one parameter, the mean metallicity of a GC system. Furthermore, the color distributions of the remaining ($\sim$30$\%$) GC systems are also explained by the nonlinearity scenario, assuming additional young or intermediate-age GCs with a number fraction of $\sim$20$\%$ of underlying old GCs. Our results reinforce the nonlinearity explanation for the GC color bimodality and provide a new perspective on early-type galaxy formation in the cluster environment, such as the Virgo and Fornax galaxy clusters.

The shape of the radial distribution of Blue Straggler Stars (BSS), when normalized to a reference population of Horizontal Branch (HB) stars, has been found to be a powerful indicator of the dynamical evolution reached by a Globular Cluster (GC). In particular, observations suggest that the BSS distribution bimodality is modulated by the dynamical age of the host GC, with dynamically unrelaxed GCs showing a flat BSS distribution, and more relaxed GCs showing a minimum at a radius that increases for increasing dynamical age, resulting in a natural dynamical clock. While direct N-body simulations are able to reproduce the general trend, thus supporting its dynamical origin, the migration of the minimum of the distribution appears to be noisy and not well defined. Here we show that a simple unidimensional model based on dynamical friction (drift) and Brownian motion (diffusion) correctly reproduces the qualitative motion of the minimum, without adjustable parameters except for the BSS to HB stars mass-ratio. Differential dynamical friction effects combine with diffusion in creating a bimodality in the BSS distribution and determining its evolution, driving the migration of the minimum to larger radii over time. The diffusion coefficient is strongly constrained by the need to reproduce the migratory behaviour of the minimum, and the radial dependence of diffusion set by fundamental physical arguments automatically satisfies this constraint. Therefore, our model appears to capture the fluctuation-dissipation dynamics that underpins the dynamical clock.

We obtained, for the first time, astrometrically registered maps of the 22.2 GHz H2O and 42.8, 43.1, and 86.2 GHz SiO maser emission toward the semiregular b-type variable (SRb) R Crateris, at three epochs (2015 May 21, and 2016 January 7 and 26) using the Korean Very-long-baseline Interferometry Network. The SiO masers show a ring-like spatial structure, while the H2O maser shows a very asymmetric one-side outflow structure, which is located at the southern part of the ring-like SiO maser feature. We also found that the 86.2 GHz SiO maser spots are distributed in an inner region, compared to those of the 43.1 GHz SiO maser, which is different from all previously known distributions of the 86.2 GHz SiO masers in variable stars. The different distribution of the 86.2 GHz SiO maser seems to be related to the complex dynamics caused by the overtone pulsation mode of the SRb R Crateris. Furthermore, we estimated the position of the central star based on the ring fitting of the SiO masers, which is essential for interpreting the morphology and kinematics of a circumstellar envelope. The estimated stellar coordinate corresponds well to the position measured by Gaia.

We investigate the effect of mixing of radioactive nickel ($^{56}$Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of $^{56}$Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak $^{56}$Ni mixing is characterized by three different phases: initial rapid reddening, blueward evolution due to the delayed effect of $^{56}$Ni heating, and redward evolution thereafter until the transition to the nebular phase. With a strong $^{56}$Ni mixing, the second phase disappears. We compare our models with the early-time color evolution of several SNe Ib/Ic (SN1999ex, SN 2008D, SN 2009jf, iPTF13bvn, SN 1994I, SN 2007gr, SN 2013ge, and 2017ein) and find signatures of relatively weak and strong $^{56}$Ni mixing for SNe Ib and SNe Ic, respectively. This suggests that SNe Ib progenitors are distinct from SN Ic progenitors in terms of helium content and that $^{56}$Ni mixing is generally stronger in the carbon-oxygen core and weaker in the helium-rich envelope. We conclude that the early-time color evolution is a powerful probe of $^{56}$Ni mixing in SNe Ib/Ic.

A large fraction of massive stars evolve in interacting binary systems, which dramatically modifies the outcome of stellar evolution. We investigated the properties of blue supergiants in binary systems and whether they are suitable for extragalactic distance determinations using the flux-weighted gravity luminosity relationship (FGLR). This is a relationship between the absolute bolometric magnitude $M_{\rm bol}$ and the spectroscopically determined flux-weighted gravity $g_{\rm F} = g/T^4_{\rm eff}$. We computed a grid of binary stellar evolution models with MESA and use the v2.1 BPASS models to examine whether they are compatible with the relatively small scatter shown by the observed relationship. Our models have initial primary masses of 9 - 30 $M_\odot$, initial orbital periods of 10 - 2511 days, mass ratio $q$ = 0.9, and metallicity $Z$ = 0.02. We find that the majority of primary stars that produce blue supergiant stages are consistent with the observed FGLR, with a small offset towards brighter bolometric magnitudes. In between 1\% - 24\% of cases, binary evolution may produce blue supergiants after a mass transfer episode, that lie below the observed FGLR. A small number of such stars have been found in extragalactic FGLR studies, suggesting that they may have evolved through binary interaction. Some models with shorter periods could resemble blue hypergiants and LBVs. We used CMFGEN radiative transfer models to investigate the effects of unresolved secondaries on diagnostics for $T_{\rm eff}$ and $g$, and the biases on the determination of interstellar reddening and $M_{\rm bol}$. We find that the effects are small and within the observed scatter, but could lead to a small overestimate of the luminosity, of $T_{\rm eff}$ and of $g$ for extreme cases. We conclude that the observed FGLR can, in principle, be well reproduced by close binary evolution models.

We present a long-exposure (~10 hr) image of the supernova (SN) remnant Cassiopeia A (Cas A) obtained with the UKIRT 3.8-m telescope using a narrow band filter centered at 1.644 um emission. The passband contains [Fe II] 1.644 um and [Si I] 1.645 um lines, and our `deep [Fe II]+[Si I] image' provides an unprecedented panoramic view of Cas A, showing both shocked and unshocked SN ejecta together with shocked circumstellar medium at subarcsec (~0.7 arcsec or 0.012 pc) resolution. The diffuse emission from the unshocked SN ejecta has a form of clumps, filaments, and arcs, and their spatial distribution correlates well with that of the Spitzer [Si II] infrared emission, suggesting that the emission is likely due to [Si I] line not [Fe II] line as in shocked material. The structure of the optically-invisible western area of Cas A is clearly seen for the first time. The area is filled with many Quasi-Stationary Flocculi (QSFs) and fragments of the disrupted ejecta shell. We suggest that the anomalous radio properties in this area could be due to the increased number of such dense clumps. We identified 309 knots in the deep [Fe II]+[Si I] image and classified them into QSFs and fast-moving knots (FMKs). The total H+He mass of QSFs is ~0.23 Msun, implying that the mass fraction of dense clumps in the progenitor's red-supergiant wind is 4--13%. The spatial distribution of QSFs suggests that there had been a highly asymmetric mass loss $10^4$--$10^5$ yr before the SN explosion. The mass of the [Fe II] line-emitting, shocked dense Fe ejecta is ~3x$10^{-5}$ Msun. The comparison with the ionic S-line dominated Hubble Space Telescope WFC3/IR image suggests that the outermost FMKs in the southeastern area are Fe-rich.

Type II supernovae (SNe) originate from the explosion of hydrogen-rich supergiant massive stars. Their first electromagnetic signature is the shock breakout, a short-lived phenomenon which can last from hours to days depending on the density at shock emergence. We present 26 rising optical light curves of SN II candidates discovered shortly after explosion by the High cadence Transient Survey (HiTS) and derive physical parameters based on hydrodynamical models using a Bayesian approach. We observe a steep rise of a few days in 24 out of 26 SN II candidates, indicating the systematic detection of shock breakouts in a dense circumstellar matter consistent with a mass loss rate $\dot{M} > 10^{-4} M_\odot yr^{-1}$ or a dense atmosphere. This implies that the characteristic hour timescale signature of stellar envelope SBOs may be rare in nature and could be delayed into longer-lived circumstellar material shock breakouts in most Type II SNe.

We study a gas rich merging dwarf system KUG 0200-096. Deep optical image reveals an optically faint tail with a length of 20 kpc, giving a visual impression of tidal antenna similar to NGC 4038/39. The interacting dwarf galaxies have B-band absolute magnitudes -18.06 and -16.63 mag. We identify a young stellar clump of stellar mass of 2$\times$10$^{7}$ M$_{\sun}$ at the tip of the antenna, possibly a Tidal Dwarf Galaxy (TDG). The putative TDG candidate is quite blue with $g-r$ color index of -0.07 mag, whereas the interacting dwarf galaxies have $g-r$ color indices 0.29 and 0.19 mag. The TDG is currently forming star at the rate of 0.02 M$_{\sun}$/yr. We obtained HI 21 cm line data of KUG 0200-096 using the GMRT to get a more detailed view of neutral hydrogen (HI) emission in interacting dwarf galaxies and its TDG. Evidence of merger between the dwarf galaxy pair is also presence in HI kinematics and morphology where we find the HI contents of interacting pair is disturbed, forming a extended tail toward the TDG. The HI velocity field shows strong gradient along the HI tidal tail extension. We present a comparative study between the Antennae galaxy, NGC 4038/39, and KUG 0200-096 in both optical and HI gas properties and discuss possible origin of KUG 0200-096 TDG.

We introduce the first synthetic light-curve model set of Type IIP supernovae exploded within circumstellar media in which the acceleration of the red supergiant winds is taken into account. Because wind acceleration makes the wind velocities near the progenitors low, the density of the immediate vicinity of the red supergiant supernova progenitors can be higher than that extrapolated by using a constant terminal wind velocity. Therefore, even if the mass-loss rate of the progenitor is relatively low, it can have a dense circumstellar medium at the immediate stellar vicinity and the early light curves of Type IIP supernovae are significantly affected by it. We adopt a simple beta velocity law to formulate the wind acceleration. We provide bolometric and multicolor light curves of Type IIP supernovae exploding within such accelerated winds from the combinations of three progenitors, 12 - 16 Msun; five beta, 1-5; seven mass-loss rates, 1e-5 - 1e-2 Msun/yr; and four explosion energies, 0.5e51 - 2e51 erg. All the light curve models are available at https://goo.gl/o5phYb. When the circumstellar density is sufficiently high, our models do not show a classical shock breakout as a consequence of the interaction with the dense and optically-thick circumstellar media. Instead, they show a delayed 'wind breakout', substantially affecting early light curves of Type IIP supernovae. We find that the mass-loss rates of the progenitors need to be 1e-3 - 1e-2 Msun/yr to explain typical rise times of 5 - 10 days in Type IIP supernovae assuming a dense circumstellar radius of 1e15 cm.

Massive stars with a core-halo structure are interesting objects for stellar physics and hydrodynamics. Using simulations for stellar evolution, radiation hydrodynamics, and radiative transfer, we study the explosion of stars with an extended and tenuous envelope (i.e., stars in which 95% of the mass is contained within 10% of the surface radius or less). We consider both H-rich supergiant and He-giant progenitors resulting from close-binary evolution and dying with a final mass of 2.8-5Msun. An extended envelope causes the supernova (SN) shock to brake and a reverse shock to form, sweeping core material into a dense shell. The shock deposited energy, which suffers little degradation from expansion, is trapped in ejecta layers of moderate optical depth, thereby enhancing the SN luminosity at early times. With the delayed 56Ni heating, we find that the resulting optical and near-IR light curves all exhibit a double-peak morphology. We show how an extended progenitor can explain the blue and featureless optical spectra of some Type IIb and Ib SNe. The dense shell formed by the reverse shock leads to line profiles with a smaller and near-constant width. This ejecta property can explain the statistically narrower profiles of Type IIb compared to Type Ib SNe, as well as the peculiar Halpha profile seen in SN1993J. At early times, our He-giant star explosion model shows a high luminosity, a blue color, and featureless spectra reminiscent of the Type Ib SN2008D, suggesting a low-mass progenitor.

Recent studies on the temperatures of red supergiants (RSGs) in the local universe provide us with an excellent observational constraint on RSG models. We calibrate the mixing length parameter by comparing model predictions with the empirical RSG temperatures in Small and Large Magellanic Clouds, Milky Way, and M31, which are inferred from the TiO band and the spectral energy distribution (SED). Although our RSG models are computed with the MESA code, our result may be applied to other stellar evolution codes, including the BEC and TWIN codes. We find evidence that the mixing length increases with increasing metallicity for both cases where the TiO and SED temperatures of RSGs are used for the calibration. Together with the recent finding of a similar correlation in low-mass red giants by Tayar et al, this implies that the metallicity dependence of the mixing length is a universal feature in post-main sequence stars of both low and high masses. Our result implies that typical Type IIP supernova (SN IIP) progenitors with initial masses of $\sim 10 - 16~M_\odot$ have a radius range of $400 R_\odot \lesssim R \lesssim 800 R_\odot$ regardless of metallicity. As an auxiliary result of this study, we find that the hydrogen-rich envelope mass of SN IIP progenitors for a given initial mass is predicted to be largely independent of metallicity if the Ledoux criterion with slow semiconvection is adopted, while the Schwarzschild models predict systematically more massive hydrogen-rich envelopes for lower metallicity.

The discovery of multiple stellar populations in Milky Way globular clusters (GCs) has stimulated various follow-up studies on helium-enhanced stellar populations. Here we present the evolutionary population synthesis models for the spectro-photometric evolution of simple stellar populations (SSPs) with varying initial helium abundance ($Y_{\rm ini}$). We show that $Y_{\rm ini}$ brings about {dramatic} changes in spectro-photometric properties of SSPs. Like the normal-helium SSPs, the integrated spectro-photometric evolution of helium-enhanced SSPs is also dependent on metallicity and age for a given $Y_{\rm ini}$. {We discuss the implications and prospects for the helium-enhanced populations in relation to the second-generation populations found in the Milky Way GCs.} All of the models are available at \url{http://web.yonsei.ac.kr/cosmic/data/YEPS.htm}.

Hydrogen-deficient Wolf-Rayet (WR) stars are potential candidates of Type Ib/Ic supernova (SN Ib/Ic) progenitors and their evolution is governed by mass loss. Stellar evolution models with the most popular prescription for WR mass-loss rates given by Nugis & Lamers have difficulties in explaining the luminosity distribution of WR stars of WC and WO types and the SN Ic progenitor properties. Here we suggest some improvements in the WR mass-loss rate prescription and discuss its implications for the evolution of WR stars and SN Ib/Ic progenitors. Recent studies on Galactic WR stars clearly indicate that the mass-loss rates of WC stars are systematically higher than those of WNE stars for a given luminosity. The luminosity and initial metallicity dependencies of WNE mass-loss rates are also significantly different from those of WC stars. These factors have not been adequately considered together in previous stellar evolution models. We also find that an overall increase of WR mass loss rates by about 60 per cent compared to the empirical values obtained with a clumping factor of 10 is needed to explain the most faint WC/WO stars. This moderate increase with our new WR mas-loss rate prescription results in SN Ib/Ic progenitor models more consistent with observations than those given by the Nugis & Lamers prescription. In particular, our new models predict that the properties of SN Ib and SN Ic progenitors are distinctively different, rather than they form a continuous sequence.

Type IIP supernova progenitors are often surrounded by dense circumstellar media that may result from mass-loss enhancement of the progenitors shortly before their explosions. Previous light-curve studies suggest that the mass-loss rates are enhanced up to ~ 0.1 Msun/yr, assuming a constant wind velocity. However, density of circumstellar media at the immediate stellar vicinity can be much higher than previously inferred for a given mass-loss rate if wind acceleration is taken into account. We show that the wind acceleration has a huge impact when we estimate mass-loss rates from early light curves of Type IIP supernovae by taking SN 2013fs as an example. We perform numerical calculations of the interaction between supernova ejecta and circumstellar media with a constant mass-loss rate but with a beta-law wind velocity profile. We find that the mass-loss rate of the progenitor of SN 2013fs shortly before the explosion, which was inferred to be ~ 0.1 Msun/yr with a constant wind velocity of 10 km/s by a previous light-curve modeling, can be as low as ~ 1e-3 Msun/yr with the same terminal wind velocity of 10 km/s but with a wind velocity profile with beta ~ 5. In both cases, the mass of the circumstellar medium is similar (~ 0.5 Msun). Therefore, the beginning of the progenitor's mass-loss enhancement in our interpretation is ~ 100 years before the explosion, not several years. Our result indicates that the immediate dense environment of Type II supernova progenitors may be significantly influenced by wind acceleration.

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.

Based mostly on stellar models which do not include rotation, CO white dwarfs which accrete helium at rates of about $\sim 10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$ have been put forward as candidate progenitors for a number of transient astrophysical phenomena, including supernovae of Type Ia, and the peculiar and fainter Type Iax supernovae. Here we study the impact of accretion-induced spin-up including the subsequent magnetic field generation, angular momentum transport, and viscous heating on the white dwarf evolution up to the point of helium ignition. We resolve the structure of the helium accreting white dwarf models with a one dimensional Langrangian hydrodynamic code, modified to include rotational and magnetic effects. We find magnetic angular momentum transport, which leads to quasi solid-body rotation, profoundly impacts the evolution of the white dwarf models. Our rotating lower mass ($0.54$ and $0.82~\mathrm{M}_\odot$) models accrete up to 50\% more mass up to ignition compared to the non-rotating case and the opposite for our more massive models. Furthermore, we find that rotation leads to up to 10-times smaller helium ignition densities, except for the lowest adopted initial white dwarf mass. Ignition densities of order $10^6\,\mathrm{g/cm}^3$ are only found for the lowest accretion rates and for large amounts of accreted helium $\gtrsim 0.4\,\mathrm{M}_\odot$. However, correspondingly massive donor stars would transfer mass at much higher rates. We expect explosive He-shell burning to mostly occur as deflagrations and at $\dot{M}>2\cdot10^{-8}~\mathrm{M}_\odot/\mathrm{yr}$, regardless of white dwarf mass. Our results imply that helium accretion onto CO white dwarfs at the considered rates is unlikely to lead to explosions of the CO core or to classical Type,Ia supernovae, but may instead produce events like faint and fast hydrogen-free supernovae.

To investigate AGN outflows as a tracer of AGN feedback on star-formation, we perform integral-field spectroscopy of 20 type 2 AGNs at z<0.1, which are luminous AGNs with the [O III] luminosity >10$^{41.5}$ erg/s, and exhibit strong outflow signatures in the [O III] kinematics. By decomposing the emission-line profile, we obtain the maps of the narrow and broad components of [O III] and H$\alpha$ lines, respectively. The broad components in both [O III] and H$\alpha$ represent the non-gravitational kinematics, i.e., gas outflows, while the narrow components, especially in H$\alpha$, represent the gravitational kinematics, i.e., rotational disk. By using the integrated spectra within the flux-weighted size of the narrow-line region, we estimate the energetics of the gas outflows. The ionized gas mass is 1.0-38.5$\times 10^5$ $M_{\odot}$, and the mean mass outflow rate is 4.6$\pm$4.3 $M_{\odot}$/yr, which is a factor of ~260 higher than the mean mass accretion rate 0.02$\pm$0.01 $M_{\odot}$/yr. The mean energy injection rate of the sample is 0.8$\pm$0.6% of the AGN bolometric luminosity, while the momentum flux is (5.4$\pm$3.6)$\times$ $L_{bol}$/c on average, except for two most kinematically energetic AGNs with low $L_{bol}$, which are possibly due to the dynamical timescale of the outflows. The estimated outflow energetics are consistent with the theoretical expectations for energy-conserving outflows from AGNs, yet we find no supporting evidence of instantaneous quenching of star formation due to the outflows.

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, $15^{+9}_{-8}$%, of core-collapse supernovae are `late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass ($4-8M_{\odot}$). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by $14^{+15}_{-14}$% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that $\phi$ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.

We explore properties of Type Ib and IIb SN progenitors that are produced by stable mass transfer in binary systems, using a new grid of stellar evolution models from an initial primary mass in the range of 10 - 18 $\mathrm{M_\odot}$ at solar and Large Magellanic Cloud metallicities. We find that blue and yellow supergiant SN IIb progenitors (e.g., of SN 2008ax, SN 2011dh, SN 2016gkg) have a hydrogen envelope mass less than about 0.15 $\mathrm{M_\odot}$, mostly resulting from early Case B mass transfer with relatively low initial masses and/or low metallicity. Red supergiant (RSG) SN IIb progenitors (e.g., of SN 1993J, SN 2013df) are produced via late Case B mass transfer and have a more massive hydrogen envelope ($M_\mathrm{H,env} > 0.15~\mathrm{M_\odot}$). SN Ib progenitors are predominantly produced by early Case B mass transfer. Our models predict that SN IIb progenitors are systematically more luminous in the optical ($-8.0 \lesssim M_\mathrm{V} \lesssim -5.0$) than the majority of SN Ib progenitors ($M_\mathrm{V} \gtrsim -5.0$) for our considered initial mass range. However, the optically bright progenitor of SN Ib iPTF13bvn (i.e., $M_\mathrm{V} \simeq -6.5$) can be well explained by a relatively low-mass progenitor with a final mass of $\sim 3.0~\mathrm{M_\odot}$. The event rate of blue and yellow SN IIb progenitors would increase as metallicity decreases, while the event rate of SN Ib progenitors would decrease instead. By contrast, the population of RSG SN IIb progenitors would not be significantly affected by metallicity.

Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.

The missing historical record of the Cassiopeia A (Cas A) supernova (SN) event implies a large extinction to the SN, possibly greater than the interstellar extinction to the current SN remnant. Here we investigate the possibility that the guest star that appeared near Cas A in 1592--1593 in Korean history books could have been an `impostor' of the Cas A SN, i.e., a luminous transient that appeared to be a SN but did not destroy the progenitor star, with strong mass loss to have provided extra circumstellar extinction. We first review the Korean records and show that a spatial coincidence between the guest star and Cas A cannot be ruled out, as opposed to previous studies. Based on modern astrophysical findings on core-collapse SN, we argue that Cas A could have had an impostor and derive its anticipated properties. It turned out that the Cas A SN impostor must have been bright ($M_V =-14.7 \pm 2.2$ mag) and an amount of dust with visual extinction of $\ge 2.8\pm 2.2 $ mag should have formed in the ejected envelope and/or in a strong wind afterwards. The mass loss needs to have been spherically asymmetric in order to see the light echo from the SN event but not the one from the impostor event.

We review our current understanding of the progenitors of both long and short duration gamma-ray bursts (GRBs). Constraints can be derived from multiple directions, and we use three distinct strands; i) direct observations of GRBs and their host galaxies, ii) parameters derived from modeling, both via population synthesis and direct numerical simulation and iii) our understanding of plausible analog progenitor systems observed in the local Universe. From these joint constraints, we describe the likely routes that can drive massive stars to the creation of long GRBs, and our best estimates of the scenarios that can create compact object binaries which will ultimately form short GRBs, as well as the associated rates of both long and short GRBs. We further discuss how different the progenitors may be in the case of black hole engine or millisecond-magnetar models for the production of GRBs, and how central engines may provide a unifying theme between many classes of extremely luminous transient, from luminous and super-luminous supernovae to long and short GRBs.

We perform integrated spectroscopy of 24 Galactic globular clusters. Spectra are observed from one core radius for each cluster with a high wavelength resolution of ~2.0 A FWHM. In combination with two existing data sets from Puzia et al. (2002) and Schiavon et al. (2005), we construct a large database of Lick spectral indices for a total of 53 Galactic globular clusters with a wide range of metallicities, -2.4 < [Fe/H] < 0.1, and various horizontal-branch morphologies. The empirical index-to-metallicity conversion relationships are provided for the 20 Lick indices for the use of deriving metallicities for remote, unresolved stellar systems.

We present new binary stellar evolution models that include the effects of tidal forces, rotation, and magnetic torques with the goal of testing Planetary Nebulae (PNe) shaping via binary interaction. We explore whether tidal interaction with a companion can spin up the AGB envelope. To do so we have selected binary systems with main sequence masses of 2.5 \Mo and of 0.8 \Mo and evolve them allowing initial separations of 5, 6, 7, and 8 AU. The binary stellar evolution models have been computed all the way to the PNe formation phase or until Roche lobe overflow (RLOF) is reached, whatever happens first. We show that with initial separations of 7 and 8 AU, the binary avoids entering into RLOF, and the AGB star reaches moderate rotational velocities at the surface ($\sim 3.5 $ and $\sim 2 $ \kms respectively) during the inter-pulse phases, but after the thermal pulses it drops to a final rotational velocity of only $\sim 0.03 $ \kms. For the closest binary separations explored, 5 and 6 AU, the AGB star reaches rotational velocities of $\sim 6 $ and $\sim 4 $ \kms respectively when the RLOF is initiated. We conclude that the detached binary models that avoid entering the RLOF phase during the AGB will not shape bipolar PNe, since the acquired angular momentum is lost via the wind during the last two thermal pulses. This study rules out tidal spin-up in non-contact binaries as a sufficient condition to form bipolar PNe.

We explore the role of the modified Eddington limit due to rapid rotation (the so-called $\Omega\Gamma-$limit) in the formation of Population III stars. We performed one-dimensional stellar evolution simulations of mass-accreting zero-metallicity protostars at a very high rate ($\dot{M} \sim 10^{-3}~\mathrm{M_\odot~yr^{-1}}$) and dealt with stellar rotation as a separate post-process. The protostar would reach the Keplerian rotation very soon after the onset of mass accretion, but mass accretion would continue as stellar angular momentum is transferred outward to the accretion disk by viscous stress. The protostar envelope expands rapidly when the stellar mass reaches $5\sim7~\mathrm{M_\odot}$ and the Eddington factor sharply increases. This makes the protostar rotate critically at a rate that is significantly below the Keplerian value (i.e., the $\Omega\Gamma-$limit). The resultant positive gradient of the angluar velocity in the boundary layer between the protostar and the Keplerian disk prohibits angular momentum transport from the star to the disk, and consequently further rapid mass accretion. This would prevent the protostar from growing significantly beyond $20 - 40~\mathrm{M_\odot}$. Another important consequence of the $\Omega\Gamma-$limit is that the protostar can remain fairly compact ($R \lesssim 50~\mathrm{R_\odot}$) and avoid a fluffy structure ($R \gtrsim 500~\mathrm{R_\odot}$) that is usually found with a very high mass accretion rate. This effect would make the protostar less prone to binary interactions during the protostar phase. Although our analysis is based on Pop III protostar models, this role of the $\Omega\Gamma-$limit would be universal in the formation process of massive stars, regardless of metallicity.

Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRBs. These include signatures related to the prompt emission mechanism, as well as to the afterglow radiation, where the surrounding intergalactic medium might imprint a telltale absorption spectrum. We emphasize important remaining uncertainties in our emerging theoretical framework.

Type Ia supernovae (SNe Ia) have been an important tool for astronomy for quite some time; however, the nature of their progenitors remains somewhat mysterious. Recent theoretical studies indicated the possibility of producing thermonuclear detonations of carbon-oxygen white dwarfs (CO WDs) at masses less than the Chandrasekhar mass through accretion of helium-rich matter, which would, depending on mass accretion rate, mass, and initial temperature of the WD, spectrally resemble either a normal SN Ia or a peculiar one. This study aims to further resolve the state of binary systems comprised of a sub-Chandrasekhar-mass CO WD and a helium star at the point where an accretion-induced detonation occurs and constrains the part of the initial parameter space where this kind of phenomenon is possible. Preexisting data obtained through simulations of single, constantly accreting CO WDs is used as an indicator for the behavior of new binary models in which the WD is treated as a point mass and which include the non-degenerate partner as a fully resolved stellar model. We parameterize the ignition of the accumulated helium layer, changes in the WD temperature, and changes in the CO core mass depending on the mass transfer rate. The initial conditions allowing for detonation do not form a single contiguous area in the parameter space, whose shape is profoundly influenced by the behavior of the donor star. Mass loss due to Nova outbursts acts in favor of detonation. According to our criteria, about 10% of the detonations in this study can be expected to show spectra consistent with ordinary SNe Ia; the rest exhibit peculiar features.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

General-relativistic radiative transfer (GRRT) calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter (mm/sub-mm) Very Long Baseline Interferometry (VLBI) observations of the supermassive black holes at the centres of Sgr A^{*} and M87. To this end we introduce Odyssey, a Graphics Processing Unit(GPU)-based code for ray tracing and radiative transfer in the Kerr spacetime. On a single GPU, the performance of Odyssey can exceed 1 nanosecond per photon, per Runge-Kutta integration step. Odyssey is publicly available, fast, accurate, and flexible enough to be modified to suit the specific needs of new users. Along with a Graphical User Interface (GUI) powered by a video-accelerated display architecture, we also present an educational software tool, Odyssey_Edu, for showing in real time how null geodesics around a Kerr black hole vary as a function of black hole spin and angle of incidence onto the black hole.

We study multiple period states of a two-component unpolarized superfluid Fermi gas in an optical lattice along the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate (BEC) crossover. The existence of states whose period is a multiple of the lattice spacing is a direct consequence of the non-linear behavior of the gas, which is due to the presence of the order parameter associated with superfluidity. By solving Bogoliubov-de Gennes equations for a superfluid flow with finite quasimomentum, we find that, in the BCS side of the crossover, the multiple period states can be energetically favorable compared to the normal Bloch states and their survival time against dynamical instability drastically increases, suggesting that these states can be accessible in current experiments, in sharp contrast to the situation in BECs.

The metallicity distribution function of globular clusters (GCs) in galaxies is a key to understanding galactic formation and evolution. The calcium II triplet (CaT) index has recently become a popular metal abundance indicator thanks to its sensitivity to GC metallicity. Here we revisit and assess the reliability of CaT as a metallicity indicator using our new stellar population synthesis simulations based on empirical, high-resolution fluxes. The model shows that the CaT strength of old ($>$ 10 Gyr) GCs is proportional to ${\rm [Fe/H]}$ below $-0.5$. In the modest metal-rich regime, however, CaT does not increase anymore with ${\rm [Fe/H]}$ due to the little contribution from coolest red giant stars to the CaT absorption. The nonlinear nature of the color--$CaT$ relation is confirmed by the observations of GCs in nearby early-type galaxies. This indicates that the CaT should be used carefully when deriving metallicities of metal-rich stellar populations. Our results offer an explanation for the observed sharp difference between the color and $CaT$ distributions of GCs in the same galaxies. We take this as an analogy to the view that metallicity--color and metallicity--Lick index nonlinearity of GCs is primarily responsible for their observed "bimodal" distributions of colors and absorption indices.

The main progenitor candidate of Type Ia supernovae (SNe Ia) is white dwarfs in binary systems where the companion star is another white dwarf (double degenerate system) or a less evolved non-degenerate star with R* >~ 0.1 Rsun (single degenerate system), but no direct observational evidence exists that tells which progenitor system is more common. Recent studies suggest that the light curve of a supernova shortly after its explosion can be used to set a limit on the progenitor size, R*. Here, we report a high cadence monitoring observation of SN 2015F, a normal SN Ia, in the galaxy NGC 2442 starting about 84 days before the first light time. With our daily cadence data, we catch the emergence of the radioactively powered light curve, but more importantly detect with a > 97.4% confidence a possible dim precursor emission that appears at roughly 1.5 days before the rise of the radioactively powered emission. The signal is consistent with theoretical expectations for a progenitor system involving a companion star with R* = ~0.1 -- 1 Rsun or a prompt explosion of a double degenerate system, but inconsistent with a typically invoked size of white dwarf progenitor of R* ~ 0.01 Rsun. Upper limits on the precursor emission also constrain the progenitor size to be R* < 0.1 Rsun, and a companion star size of R* < ~1.0 Rsun, excluding a very large companion star in the progenitor system. Additionally, we find that the distance to SN 2015F is 23.9 +-0.4 Mpc.

Gravitational wave bursts in the formation of neutron stars and black holes in energetic core-collapse supernovae (CC-SNe) are of potential interest to LIGO-Virgo and KAGRA. Events nearby are readily discovered using moderately sized telescopes. CC-SNe are competitive with mergers of neutron stars and black holes, if the fraction producing an energetic output in gravitational waves exceeds about 1\%. This opportunity motivates the design of a novel Sejong University Core-CollapsE Supernova Survey (SUCCESS), to provide triggers for follow-up searches for gravitational waves. It is based on the 76 cm Sejong University Telescope (SUT) for weekly monitoring of nearby star-forming galaxies, i.e., M51, M81-M82 and Blue Dwarf Galaxies from the Unified Nearby Galaxy Catalog with an expected yield of a few hundred per year. Optical light curves will be resolved for the true time-of-onset for probes of gravitational waves by broadband time-sliced matched filtering.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

Massive rotating single stars with an initial metal composition appropriate for the dwarf galaxy I Zw 18 ([Fe/H]=$-$1.7) are modelled during hydrogen burning for initial masses of 9-300 M$_{\odot}$ and rotational velocities of 0-900 km s$^{-1}$. Internal mixing processes in these models were calibrated based on an observed sample of OB-type stars in the Magellanic Clouds. Even moderately fast rotators, which may be abundant at this metallicity, are found to undergo efficient mixing induced by rotation resulting in quasi chemically-homogeneous evolution. These homogeneously-evolving models reach effective temperatures of up to 90 kK during core hydrogen burning. This, together with their moderate mass-loss rates, make them Transparent Wind Ultraviolet INtense stars (TWUIN star), and their expected numbers might explain the observed HeII ionizing photon flux in I Zw 18 and other low-metallicity HeII galaxies. Our slowly rotating stars above $\sim$80 M$_{\odot}$ evolve into late B- to M-type supergiants during core hydrogen burning, with visual magnitudes up to 19$^{\mathrm{m}}$ at the distance of I Zw 18. Both types of stars, TWUIN stars and luminous late-type supergiants, are only predicted at low metallicity. Massive star evolution at low metallicity is shown to differ qualitatively from that in metal-rich environments. Our grid can be used to interpret observations of local star-forming dwarf galaxies and high-redshift galaxies, as well as the metal-poor components of our Milky Way and its globular clusters.

In several recent observational studies on Type Ib/c supernovae (SNe Ib/c), the inferred ejecta masses have a peak value of 2.0 -- 4.0 $M_\odot$, in favor of the binary scenario for their progenitors rather than the Wolf-Rayet star scenario. To investigate the observational properties of relatively low-mass helium stars in binary systems as SN Ib/c progenitors, we constructed atmospheric models with the non-LTE radiative transfer code CMFGEN, using binary star evolution models. We find that these helium stars can be characterized by relatively narrow helium emission lines if the mass-loss rate during the final evolutionary phase is significantly enhanced as implied by many SN Ib/c observations. The optical brightness of helium star progenitors can be meaningfully enhanced with a strong wind for $M \gtrsim 4.4 M_\odot$, but hardly affected or slightly weakened for relatively low-mass of $\sim 3.0 M_\odot$, compared to the simple estimate using blackbody approximation. We further confirm the previous suggestion that the optical brightness would be generally higher for a less massive SN Ib/c progenitor. In good agreement with previous studies, our results indicate that the optical magnitudes and colors of the recently detected progenitor of the SN Ib iPTF13bvn can be well explained by a binary progenitor with a final helium star mass of about 3.0 -- 4.4 $M_\odot$.

Type Ib/c supernovae (SNe Ib/c) mark the deaths of hydrogen-deficient massive stars. The evolutionary scenarios for SNe Ib/c progenitors involve many important physical processes including mass loss by winds and its metallicity dependence, stellar rotation, and binary interactions. This makes SNe Ib/c an excellent test bed for stellar evolution theory. We review the main results of evolutionary models for SN Ib/c progenitors available in the literature and their confrontation with recent observations. We argue that the nature of SN Ib/c progenitors can be significantly different for single and binary systems, and that binary evolution models can explain the ejecta masses derived from SN Ib/c light curves, the distribution of SN Ib/c sites in their host galaxies, and the optical magnitudes of the tentative progenitor candidate of iPTF13bvn. We emphasize the importance of early-time observations of light curves and spectra, accurate measurements of helium mass in SN Ib/c ejecta, and systematic studies about the metallicity dependence of SN Ib/c properties, to better constrain theories.

A fraction of the first generation of stars in the early Universe may be very massive ($\gtrsim 300~\mathrm{M_\odot}$) as they form in metal-free environments. Formation of black holes from these stars can be accompanied by supermassive collapsars to produce long gamma-ray bursts of a unique type having a very high total energy ($\sim 10^{54}~\mathrm{erg}$) as recently suggested by several authors. We present new stellar evolution models of very massive Population III stars including the effect of rotation to provide theoretical constraints on super-collapsar progenitors. We find that the angular momentum condition for super-collapsar can be fulfilled if magnetic torques are ignored, in which case Eddington-Sweet circulations play the dominant role for the transport of angular momentum. We further find that the initial mass range for super-collapsar progenitors would be limited to $300~\mathrm{M_\odot} \lesssim M \lesssim 700~\mathrm{M_\odot}$. However, all of our very massive star models of this mass range end their lives as red supergiants rather than blue supergiants, in good agreement with most of the previous studies. The predicted final fate of these stars is either a jet-powered type IIP supernova or an ultra-long, relatively faint gamma-ray transient, depending on the initial amount of angular momentum.

We investigate the power spectra of outflow-driven turbulence through high-resolution three-dimensional isothermal numerical simulations where the turbulence is driven locally in real-space by a simple spherical outflow model. The resulting turbulent flow saturates at an average Mach number of ~2.5 and is analysed through density and velocity power spectra, including an investigation of the evolution of the solenoidal and compressional components. We obtain a shallow density power spectrum with a slope of ~-1.2 attributed to the presence of a network of localised dense filamentary structures formed by strong shock interactions. The total velocity power spectrum slope is found to be ~-2.0, representative of Burgers shock dominated turbulence model. The density weighted velocity power spectrum slope is measured as ~-1.6, slightly less than the expected Kolmogorov scaling value (slope of -5/3) found in previous works. The discrepancy may be caused by the nature of our real space driving model and we suggest there is no universal scaling law for supersonic compressible turbulence. We find that on average, solenoidal modes slightly dominate in our turbulence model as the interaction between strong curved compressible shocks generates solenoidal modes, and compressible modes decay faster.

Intensity interferometry, based on the Hanbury Brown-Twiss effect, is a simple and inexpensive method for optical interferometry at microarcsecond angular resolutions; its use in astronomy was abandoned in the 1970s because of low sensitivity. Motivated by recent technical developments, we argue that the sensitivity of large modern intensity interferometers can be improved by factors up to approximately 25,000, corresponding to 11 photometric magnitudes, compared to the pioneering Narrabri Stellar Interferometer. This is made possible by (i) using avalanche photodiodes (APD) as light detectors, (ii) distributing the light received from the source over multiple independent spectral channels, and (iii) use of arrays composed of multiple large light collectors. Our approach permits the construction of large (with baselines ranging from few kilometers to intercontinental distances) optical interferometers at the cost of (very) long-baseline radio interferometers. Realistic intensity interferometer designs are able to achieve limiting R-band magnitudes as good as ~14, sufficient for spatially resolved observations of main-sequence O-type stars in the Magellanic Clouds. Multi-channel intensity interferometers can address a wide variety of science cases: (i) linear radii, effective temperatures, and luminosities of stars; (ii) mass-radius relationships of compact stellar remnants; (iii) stellar rotation; (iv) stellar convection and the interaction of stellar photospheres and magnetic fields; (v) the structure and evolution of multiple stars; (vi) direct measurements of interstellar distances; (vii) the physics of gas accretion onto supermassive black holes; and (viii) calibration of amplitude interferometers by providing a sample of calibrator stars.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

Warped disks are almost ubiquitous among spiral galaxies. Here we revisit and test the `fly-by scenario' of warp formation, in which impulsive encounters between galaxies are responsible for warped disks. Based on N-body simulations, we investigate the morphological and kinematical evolution of the stellar component of disks when galaxies undergo fly-by interactions with adjacent dark matter halos. We find that the so-called `S'-shaped warps can be excited by fly-bys and sustained for even up to a few billion years, and that this scenario provides a cohesive explanation for several key observations. We show that disk warp properties are governed primarily by the following three parameters; (1) the impact parameter, i.e., the minimum distance between two halos, (2) the mass ratio between two halos, and (3) the incident angle of the fly-by perturber. The warp angle is tied up with all three parameters, yet the warp lifetime is particularly sensitive to the incident angle of the perturber. Interestingly, the modeled S-shaped warps are often non-symmetric depending on the incident angle. We speculate that the puzzling U- and L-shaped warps are geometrically superimposed S-types produced by successive fly-bys with different incident angles, including multiple interactions with a satellite on a highly elongated orbit.

The United Kingdom Infrared Telescope (UKIRT) Widefield Infrared Survey for Fe$^+$ (UWIFE) is a 180 deg$^2$ imaging survey of the first Galactic quadrant (7$^{\circ}$ < l < 62$^{\circ}$; |b| < 1.5$^{\circ}$) using a narrow-band filter centered on the [Fe II] 1.644 {\mu}m emission line. The [Fe II] 1.644 {\mu}m emission is a good tracer of dense, shock-excited gas, and the survey will probe violent environments around stars: star-forming regions, evolved stars, and supernova remnants, among others. The UWIFE survey is designed to complement the existing UKIRT Widefield Infrared Survey for H2 (UWISH2; Froebrich et al. 2011). The survey will also complement existing broad-band surveys. The observed images have a nominal 5{\sigma} detection limit of 18.7 mag for point sources, with the median seeing of 0.83". For extended sources, we estimate surface brightness limit of 8.1 x 10$^{-20}$ W m$^{-2}$ arcsec$^{-2}$ . In this paper, we present the overview and preliminary results of this survey.

Both recent observations and stellar evolution models suggest that pair-instability supernovae (PISNe) could occur in the local Universe, at metallicities below Z_Sun/3. Previous PISN models were mostly produced at very low metallicities in the context of the early Universe. We present new PISNe models at a metallicity of Z=0.001, which are relevant for the local Universe. We take the self-consistent stellar evolutionary models of pair-instability progenitors with initial masses of 150 and 250 solar masses at metallicity of Z=0.001 by Langer et al. (2007) and follow the evolution of these models through the supernova explosions, using a hydrodynamics stellar evolution code with an extensive nuclear network including 200 isotopes. Both models explode as PISNe without leaving a compact stellar remnant. Our models produce a nucleosynthetic pattern that is generally similar to that of Population III PISN models, which is mainly characterized by the production of large amounts of alpha-elements and a strong deficiency of the odd-charged elements. However, the odd-even effect in our models is significantly weaker than that found in Population III models. The comparison with the nucleosynthetic yields from core-collapse supernovae at a similar metallicity (Z=0.002) indicates that PISNe could have strongly influenced the chemical evolution below Z=0.002, assuming a standard initial mass function. The odd-even effect is predicted to be most prominent for the intermediate mass elements between silicon and calcium. With future observations of chemical abundances in Population II stars, our result can be used to constrain the number of PISNe that occurred during the past evolution of our Galaxy.

We investigate the formation of dust in a stellar wind during the red-supergiant (RSG) phase of a very massive Population III star with the zero-age main sequence mass of 500 M_sun. We show that, in a carbon-rich wind with a constant velocity, carbon grains can form with a lognormal-like size distribution, and that all of the carbon available for dust formation finally condense into dust for wide ranges of the mass-loss rate ((0.1-3)x10^{-3} M_sun yr^{-1}) and wind velocity (1-100 km s^{-1}). We also find that the acceleration of the wind driven by newly formed dust suppresses the grain growth but still allows more than half of gas-phase carbon to be finally locked up in dust grains. These results indicate that at most 1.7 M_sun of carbon grains can form in total during the RSG phase of 500 M_sun Population III stars. Such a high dust yield could place very massive primordial stars as important sources of dust at the very early epoch of the universe if the initial mass function of Population III stars was top-heavy. We also briefly discuss a new formation scenario of carbon-rich ultra-metal-poor stars considering the feedback from very massive Population III stars.

So called superluminous supernovae have been recently discovered in the local Universe. It appears possible that some of them originate from stellar explosions induced by the pair instability mechanism. Recent stellar evolution models also predict pair instability supernovae (PISNe) from very massive stars at fairly high metallicities (i.e. Z~0.004). We provide supernova (SN) models and synthetic light curves for two progenitor models, a 150 Msun red-supergiant and a 250 Msun yellow-supergiant at a metallicity of Z=0.001, for which the evolution from the main sequence to collapse, and the initiation of the PISN itself, has been previously computed in a realistic and self-consistent way. We are using the radiation hydrodynamics code STELLA to describe the SN evolution of both models over a time frame of about 500 days. We describe the shock-breakout phases of both SNe which are characterized by a higher luminosity, a longer duration and a lower effective temperature than those of ordinary SNeIIP. We derive the bolometric as well as the U, B, V, R and I light curves of our PISN models, which show a long-lasting plateau phase with maxima at Mbol=-19.3 mag and -21.3 mag for our lower and higher mass model, respectively. We also describe the photospheric composition and velocity as function of time. We conclude that the light curve of the explosion of our initially 150 Msun star resembles those of relatively bright SNeIIP, whereas its photospheric velocity at early times is smaller. Its 56Ni mass of 0.04 Msun also falls well into the range found in ordinary core collapse SNe. The light curve and photospheric velocity of our 250 Msun models has a striking resemblance with that of the superluminous SN2007bi, strengthening its interpretation as PISN. We conclude that PISNe may occur more frequently in the local universe than previously assumed.

We have computed new stellar evolution models that include the effects of rotation and magnetic torques under different hypothesis. The goal is to test if a single star can sustain in the envelope the rotational velocities needed for the magneto hydrodynamical (MHD) simulations to shape bipolar Planetary Nebulae (PNe) when the high mass-loss rates take place. Stellar evolution models with main sequence masses of 2.5 and 5 Mo, and initial rotational velocities of 250 km/s have been followed all the way to the PNe formation phase. We find that stellar cores have to be spun down using magnetic torques in order to reproduce the rotation rates observed for white dwarfs. During the asymptotic giant branch phase and beyond, the magnetic braking of the core has a practically null effect in increasing the rotational velocity of the envelope since the stellar angular momentum is removed efficiently by the wind. We have, as well, tested best possible case scenarios in rather non-physical contexts to give enough angular momentum to the envelope. We find that we cannot get the envelope of a single star rotating at the speeds needed by the MHD simulations to form bipolar PNe. We conclude that single stellar rotators are unlikely to be the progenitors of bipolar PNe under the current MHD model paradigm.

Phosphorus ($^{31}$P), which is essential for life, is thought to be synthesized in massive stars and dispersed into interstellar space when these stars explode as supernovae (SNe). Here we report on near-infrared spectroscopic observations of the young SN remnant Cassiopeia A, which show that the abundance ratio of phosphorus to the major nucleosynthetic product iron ($^{56}$Fe) in SN material is up to 100 times the average ratio of the Milky Way, confirming that phosphorus is produced in SNe. The observed range is compatible with predictions from SN nucleosynthetic models but not with the scenario in which the chemical elements in the inner SN layers are completely mixed by hydrodynamic instabilities during the explosion.

Recent discoveries of weak and fast optical transients raise the question of their origin. We investigate the minimum ejecta mass associated with core-collapse supernovae (SNe) of Type Ic. We show that mass transfer from a helium star to a compact companion can produce an ultra-stripped core which undergoes iron core collapse and leads to an extremely fast and faint SN Ic. In this Letter, a detailed example is presented in which the pre-SN stellar mass is barely above the Chandrasekhar limit, resulting in the ejection of only ~0.05-0.20 M_sun of material and the formation of a low-mass neutron star. We compute synthetic light curves of this case and demonstrate that SN 2005ek could be explained by our model. We estimate that the fraction of such ultra-stripped to all SNe could be as high as 0.001-0.01. Finally, we argue that the second explosion in some double neutron star systems (for example, the double pulsar PSR J0737-3039B) was likely associated with an ultra-stripped SN Ic.

Recent spectroscopic observations of the outer halo globular cluster (GC) NGC 2419 show that it is unique among GCs, in terms of chemical abundance patterns, and some suggest that it was originated in the nucleus of a dwarf galaxy. Here we show, from the Subaru narrow-band photometry employing a calcium filter, that the red giant-branch (RGB) of this GC is split into two distinct subpopulations. Comparison with spectroscopy has confirmed that the redder RGB stars in the $hk$[=(Ca$-b)-(b-y)$] index are enhanced in [Ca/H] by $\sim$0.2 dex compared to the bluer RGB stars. Our population model further indicates that the calcium-rich second generation stars are also enhanced in helium abundance by a large amount ($\Delta$Y = 0.19). Our photometry, together with the results for other massive GCs (e.g., $\omega$ Cen, M22, and NGC 1851), suggests that the discrete distribution of RGB stars in the $hk$ index might be a universal characteristic of this growing group of peculiar GCs. The planned narrow-band calcium photometry for the Local Group dwarf galaxies would help to establish an empirical connection between these GCs and the primordial building blocks in the hierarchical merging paradigm of galaxy formation.

Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-AIC evolution of close binaries. We investigated the mass transfer onto a massive WD in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. We demonstrate that the final properties of these MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. Finally, we discuss the link between AIC and young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. [Abridged]

Protostellar jets and outflows are signatures of star formation and promising mechanisms for driving supersonic turbulence in molecular clouds. We quantify outflow-driven turbulence through three-dimensional numerical simulations using an isothermal version of the robust total variation diminishing code. We drive turbulence in real-space using a simplified spherical outflow model, analyse the data through density probability distribution functions (PDF), and investigate the Core Formation Rate per free-fall time (CFR_ff). The real-space turbulence driving method produces a negatively skewed density PDF possessing an enhanced tail on the low-density side. It deviates from the log-normal distributions typically obtained from Fourier-space turbulence driving at low densities, but can provide a good fit at high-densities, particularly in terms of mass weighted rather than volume weighted density PDF. Due to this fact, we suggest that the CFR_ff determined from a Fourier-driven turbulence model could be comparable to that of our particular real-space driving model, which has a ratio of solenoidal to compressional components from the resulting turbulence velocity fields of ~0.6.

The mean color of globular clusters (GCs) in early-type galaxies is in general bluer than the integrated color of halo field stars in host galaxies. Metal-rich GCs often appear more associated with field stars than metal-poor GCs, yet show bluer colors than their host galaxy light. Motivated by the discovery of multiple stellar populations in Milky Way GCs, we present a new scenario in which the presence of second-generation (SG) stars in GCs is responsible for the color discrepancy between metal-rich GCs and field stars. The model assumes that the SG populations have enhanced helium abundance as evidenced by observations, and it gives a good explanation of the bluer optical colors of metal-rich GCs than field stars as well as strong Balmer lines and blue UV colors of metal-rich GCs. Ours may be complementary to the recent scenario suggesting the difference in stellar mass functions (MFs) as an origin for the GC-to-star color offset. A quantitative comparison is given between the SG and MF models.

Recent spectroscopy on the globular cluster (GC) system of M31 with unprecedented precision witnessed a clear bimodality in absorption-line index distributions of old GCs. Such division of extragalactic GCs, so far asserted mainly by photometric color bimodality, has been viewed as the presence of merely two distinct metallicity subgroups within individual galaxies and forms a critical backbone of various galaxy formation theories. Given that spectroscopy is a more detailed probe into stellar population than photometry, the discovery of index bimodality may point to the very existence of dual GC populations. However, here we show that the observed spectroscopic dichotomy of M31 GCs emerges due to the nonlinear nature of metallicity-to-index conversion and thus one does not necessarily have to invoke two separate GC subsystems. We take this as a close analogy to the recent view that metallicity-color nonlinearity is primarily responsible for observed GC color bimodality. We also demonstrate that the metallicity-sensitive magnesium line displays non-negligible metallicity-index nonlinearity and Balmer lines show rather strong nonlinearity. This gives rise to bimodal index distributions, which are routinely interpreted as bimodal metallicity distributions, not considering metallicity-index nonlinearity. Our findings give a new insight into the constitution of M31's GC system, which could change much of the current thought on the formation of GC systems and their host galaxies.

Color distributions of globular clusters (GCs) in most massive galaxies are bimodal. Assuming linear color-to-metallicity conversions, bimodality is viewed as the presence of merely two GC subsystems with distinct metallicities, which serves as a critical backbone of various galaxy formation theories. Recent studies, however, revealed that the color-metallicity relations (CMRs) often used to derive GC metallicities (e.g., CMRs of g-z, V-I and C-T1) are in fact inflected. Such inflection can create bimodal color distributions if the underlying GC metallicity spread is simply broad as expected from the hierarchical merging paradigm of galaxy formation. In order to test the nonlinear-CMR scenario for GC color bimodality, the u-band photometry is proposed because the u-related CMRs (e.g., CMRs of u-g and u-z) are theoretically predicted to be least inflected and most distinctive among commonly used optical CMRs. Here, we present Hubble Space Telescope (HST)/WFC3 F336W (u-band) photometry of the GC system in M84, a giant elliptical in the Virgo galaxy cluster. Combining the u data with the existing HST ACS/WFC g and z data, we find that the u-z and u-g color distributions are different from the g-z distribution in a very systematic manner and remarkably consistent with our model predictions based on the nonlinear-CMR hypothesis. The results lend further confidence to validity of the nonlinear-CMR scenario as an explanation for GC color bimodality. There are some GC systems showing bimodal spectroscopic metallicity, and in such systems the inflected CMRs often create stronger bimodality in the color domain.

We present wide-field multiband photometry of globular cluster (GC) systems in NGC 1399, NGC 1404, and NGC 1387 located at the central region of the Fornax galaxy cluster. Observation was carried out through U, B, V, and I bands, which marks one of the widest and deepest U-band studies on extragalactic GC systems. The present U-band photometry enables us to significantly reduce the contamination by a factor of two for faint sources (Vo~23.5). The main results based on some 2000 GC candidates around NGC 1399, NGC 1404, and NGC 1387 are as follows: (1) the GC system in each galaxy exhibits bimodal color distributions in all colors examined, but the shape of color histograms varies systematically depending on colors; (2) NGC 1399 shows that the mean colors of both blue and red GCs become bluer with increasing galactocentric radius; (3) NGC 1399 shows overabundance of GCs in the directions of NGC 1404 and NGC 1387, indicating their ongoing interactions; and (4) NGC 1399 also exhibits a ~0.5' offset between the center of the inner GC distribution and the galaxy's optical center, suggesting that NGC 1399 is not yet dynamically relaxed and may be undergoing merger events.

Despite the importance of their size evolution in understanding the dynamical evolution of globular clusters (GCs) of the Milky Way, studies are rare that focus specifically on this issue. Based on the advanced, realistic Fokker-Planck (FP) approach, we predict theoretically the initial size distribution (SD) of the Galactic GCs along with their initial mass function and radial distribution. Over one thousand FP calculations in a wide parameter space have pinpointed the best-fit initial conditions for the SD, mass function, and radial distribution. Our best-fit model shows that the initial SD of the Galactic GCs is of larger dispersion than today's SD, and that typical projected half-light radius of the initial GCs is ~4.6 pc, which is 1.8 times larger than that of the present-day GCs (~2.5 pc). Their large size signifies greater susceptibility to the Galactic tides: the total mass of destroyed GCs reaches 3-5x10^8 M_sun$, several times larger than the previous estimates. Our result challenges a recent view that the Milky Way GCs were born compact on the sub-pc scale, and rather implies that (1) the initial GCs are generally larger than the typical size of the present-day GCs, (2) the initially large GCs mostly shrink and/or disrupt as a result of the galactic tides, and (3) the initially small GCs expand by two-body relaxation, and later shrink by the galactic tides.

We present a series of papers on the year-2012 version of Yonsei Evolutionary Population Synthesis (YEPS) model which is constructed on over 20 years of heritage. This first paper delineates the spectroscopic aspect of integrated light from stellar populations older than 1 Gyr. The standard YEPS is based on the most up-to-date Yonsei-Yale stellar evolutionary tracks and BaSel 3.1 flux libraries, and provides absorption line indices of the Lick/IDS system and high-order Balmer lines for simple stellar populations as functions of stellar parameters, such as metallicity, age and {\alpha}-element mixture. Special care has been taken to incorporate systematic contribution from horizontal-branch stars which alters the temperature-sensitive Balmer lines significantly, resulting in up to 5 Gyr difference in age estimation of old, metal-poor stellar populations. We also find that the horizontal branches exert an appreciable effect not only on the Balmer lines but also on the metallicity-sensitive lines including the magnesium index. This is critical to explain the intriguing bimodality found in index distributions of globular clusters in massive galaxies and to derive spectroscopic metallicities accurately from various indices. A full set of the spectroscopic and photometric YEPS model data of the entire parameter space is currently downloadable at http://web.yonsei.ac.kr/cosmic/data/YEPS.htm

(shortened) The first couple of stellar generations may have been massive, of order 100 Msun, and to have played a dominant role in galaxy formation and the chemical enrichment of the early Universe. Some fraction of these objects may have died as pair-instability supernovae or gamma-ray bursts. The winds if these stars may have played an important role in determining these outcomes. As the winds are driven by radiation pressure on spectral lines, their strengths are expected to vary with metallicity. Until now, most mass-loss predictions for metal-poor O-type stars have assumed a scaled-down solar-abundance pattern. However, Population III evolutionary tracks show significant surface enrichment through rotational mixing of CNO-processed material, because even metal-poor stars switch to CNO-burning early on. We address the question of whether the CNO surface enhanced self-enrichment in the first few generations of stars could impact their mass-loss properties. For this, we employ Monte Carlo simulations to establish the local line-force and solve for the momentum equation of the stellar outflow, testing whether an outflow can actually be established by assessing the net acceleration at the sonic point of the flow. Stellar evolution models of rotating metal-poor stars are used to specify the surface chemical composition, focussing on the phases of early enrichment. We find that the mass-loss rates of CNO enhanced metal-poor stars are higher than those of non-enriched stars, but they are much lower than those rates where the CNO abundance is included in the total abundance Z. We present a heuristic formula that provides mass-loss estimates for CNO-dominated winds in relation to scaled-down solar abundances.

The progenitors of many Type II supernovae have been observationally identified but the search for Type Ibc supernova (SN Ibc) progenitors has thus far been unsuccessful, despite the expectation that they are luminous Wolf-Rayet (WR) stars. We investigate how the evolution of massive helium stars affects their visual appearances, and discuss the implications for the detectability of SN Ibc progenitors. Massive WR stars that rapidly lose their helium envelopes through stellar-wind mass-loss end their lives when their effective temperatures -- related to their hydrostatic surfaces -- exceed about 150kK.At their pre-supernova stage, their surface properties resemble those of hot Galactic WR stars of WO sub-type. These are visually faint with narrow-band visual magnitudes Mv = -1.5 ~ -2.5, despite their high bolometric luminosities (log L/Lsun = 5.6 ~ 5.7), compared to the bulk of Galactic WR stars (Mv < -4). In contrast, relatively low-mass helium stars that retain a thick helium envelope appear fairly bright in optical bands, depending on the final masses and the history of the envelope expansion during the late evolutionary stages. We conclude that SNe Ibc observations have so far not provided strong constraints on progenitor bolometric luminosities and masses, even with the deepest searches. We also argue that Ic progenitors are more challenging to identify than Ib progenitors in any optical images.

Using spatially-resolved spectra obtained with the Low Resolution Imaging Spectrometer at the Keck I telescope, we investigate the nature of ionizing sources and kinematic properties of emission-line gas in a LINER galaxy SDSS J091628.05+420818.7, which is a nearby (z = 0.0241) and bright (M_r = -20.2) early-type galaxy. After subtracting stellar absorption features using a combination of simple stellar population models, we measure the flux, line-of-sight velocity, and velocity dispersion of four emission lines, i.e., H{\alpha}, H{\beta}, [O III] {\lambda}5007, and [N II] {\lambda}6584, to study radial change of emission-line fluxes and velocities. Compared to the point-spread-function of the observation, the emission-line region is slightly extended but comparable to the seeing size. The central concentration of emission-line gas suggests that ionization is triggered by a nuclear source, excluding old stellar population as ionizing sources. We find that emission-line gas is counter-rotating with respect to stellar component and that the [O III] {\lambda}5007 line is blueshifted compared to other emission lines, possibly due to an outflow.

[Abridged] We present a new grid of massive population III star models including the effects of rotation on the stellar structure and chemical mixing, and magnetic torques for the transport of angular momentum. Based on the grid, we also present a phase diagram for the expected final fates of rotating massive Pop III stars. Our non-rotating models become redder than the previous models in the literature, given the larger overshooting parameter adopted in this study. In particular, convective dredge-up of the helium core material into the hydrogen envelope is observed in our non-rotating very massive star models (>~200 Msun), which is potentially important for the chemical yields. On the other hand, the stars become bluer and more luminous with a higher rotational velocity. With the Spruit-Tayler dynamo, our models with a sufficiently high initial rotational velocity can reach the critical rotation earlier and lose more mass as a result, compared to the previous models without magnetic fields. The most dramatic effect of rotation is found with the so-called chemically homogeneous evolution (CHE), which is observed for a limited mass and rotational velocity range. CHE has several important consequences: 1) Both primary nitrogen and ionizing photons are abundantly produced. 2) Conditions for gamma-ray burst progenitors are fulfilled for an initial mass range of 13 - 84 Msun. 3) Pair instability supernovae of type Ibc are expected for 84 -190 Msun and 4) Both a pulsational pair instability supernova and a GRB may occur from the same progenitor of about 56 - 84 Msun, which might significantly influence the consequent GRB afterglow. We find that CHE does not occur for very massive stars (> 190 Msun), in which case the hydrogen envelope expands to the red-supergiant phase and the final angular momentum is too low to make any explosive event powered by rotation.

Deep images of 10 early-type galaxies in low-density environments have been obtained with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The global properties of the globular cluster (GC) systems of the galaxies have been derived in order to investigate the role of the environment in galaxy formation and evolution. Using the ACS Virgo Cluster Survey (ACSVCS) as a high-density counterpart, the similarities and differences between the GC properties in high- and low-density environments are presented. We find a strong correlation of the GC mean colours and the degree of colour bimodality with the host galaxy luminosity in low-density environments, in good agreement with high-density environments. In contrast, the GC mean colours at a given host luminosity are somewhat bluer (\Delta(g-z) ~ 0.05) than those for cluster galaxies, indicating more metal-poor (\Delta[Fe/H] ~ 0.10-0.15) and/or younger (\Delta age > 2 Gyr) GC systems than those in dense environments. Furthermore, with decreasing host luminosity, the colour bimodality disappears faster, when compared to galaxies in cluster environments. Our results suggest that: (1) in both high- and low-density environments, the mass of the host galaxy has the dominant effect on GC system properties, (2) the local environment has only a secondary effect on the history of GC system formation, (3) GC formation must be governed by common physical processes across a range of environments.

We present GALEX data for 44 Galactic globular clusters obtained during 3 GALEX observing cycles between 2004 and 2008. This is the largest homogeneous data set on the UV photometric properties of Galactic globular clusters ever collected. The sample selection and photometric analysis are discussed, and color-magnitude diagrams are presented. The blue and intermediate-blue horizontal branch is the dominant feature of the UV color-magnitude diagrams of old Galactic globular clusters. Our sample is large enough to display the remarkable variety of horizontal branch shapes found in old stellar populations. Other stellar types that are obviously detected are blue stragglers and post core-He burning stars. The main features of UV color-magnitude diagrams of Galactic globular clusters are briefly discussed. We establish the locus of post-core He burning stars in the UV color-magnitude diagram and present a catalog of candidate AGB-manqu \'e, post early-AGB, and post-AGB stars within our cluster sample.

We investigate the energy band structure of the superfluid flow of ultracold dilute Fermi gases in a one-dimensional optical lattice along the BCS to BEC crossover within a mean-field approach. In each side of the crossover region, a loop structure (swallowtail) appears in the Bloch energy band of the superfluid above a critical value of the interaction strength. The width of the swallowtail is largest near unitarity. Across the critical value of the interaction strength, the profiles of density and pairing field change more drastically in the BCS side than in the BEC side. It is found that along with the appearance of the swallowtail, there exists a narrow band in the quasiparticle energy spectrum close to the chemical potential and the incompressibility of the Fermi gas consequently experiences a profound dip in the BCS side, unlike in the BEC side.

Massive stars lose a large fraction of their original mass over the course of their evolution. These stellar winds shape the surrounding medium according to parameters that are the result of the characteristics of the stars, varying over time as the stars evolve, leading to both permanent and temporary features that can be used to constrain the evolution of the progenitor star. Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the characteristics of the circumstellar medium (CSM) should be observable in the signal of GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the medium it collides with, and indirectly because the GRB can only be observed through the extended circumstellar bubble that surrounds each massive star. We use computer simulations to describe the circumstellar features that can be found in the vicinity of massive stars and discuss if, and how, they may appear in GRB observations. Specifically, we make hydrodynamical models of the circumstellar environment of a rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor candidate model. The simulations show that the star creates a large scale bubble of shocked wind material, which sweeps up the interstellar medium in an expanding shell. Within this bubble, temporary circumstellar shells, clumps and voids are created as a result of changes in the stellar wind. Most of these temporary features have disappeared by the time the star reaches the end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star in a high density environment simplifies the evolution of the CSM as the more confined bubble prohibits the formation of some of the temporary structures.

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly-peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm of galaxies and supports additionally the similar nature shared by GCs and field stars. Our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

The optical color distributions of globular clusters (GCs) in most large elliptical galaxies are bimodal. Based on the assumed linear relationship between GC colors and their metallicities, the bimodality has been taken as evidence of two GC subsystems with different metallicities in each galaxy and led to a number of theories in the context of galaxy formation. More recent observations and modeling of GCs, however, suggest that the color-metallicity relations (CMRs) are inflected, and thus colors likely trace metallicities in a nonlinear manner. The nonlinearity could produce bimodal color distributions from a broad underlying metallicity spread, even if it is unimodal. Despite the far-reaching implications, whether CMRs are nonlinear and whether the nonlinearity indeed causes the color bimodality are still open questions. Given that the spectroscopic refinement of CMRs is still very challenging, we here propose a new photometric technique to probe the possible nonlinear nature of CMRs. In essence, a color distribution of GCs is a "projected" distribution of their metallicities. Since the form of CMRs hinges on which color is used, the shape of color distributions varies depending significantly on the colors. Among other optical colors, the u-band related colors (e.g., u - g and u - z) are theoretically predicted to exhibit significantly less inflected CMRs than other preferred CMRs (e.g., for g - z). As a case study, we performed the HST/WFPC2 archival u-band photometry for the M87 (NGC 4486) GC system with confirmed color bimodality. We show that the u-band color distributions are significantly different from that of g - z, and consistent with our model predictions. With more u-band measurements, this method will support or rule out the nonlinear-CMR scenario for the origin of GC color bimodality with high confidence. [Abridged]

Recent observations and modeling of globular clusters with multiple populations strongly indicate the presence of super helium-rich subpopulations in old stellar systems. Motivated by this, we have constructed new population synthesis models with and without helium-enhanced subpopulations to investigate their impact on the UV-upturn phenomenon of quiescent early-type galaxies. We find that our models with helium- enhanced subpopulations can naturally reproduce the strong UV-upturns observed in giant elliptical galaxies assuming an age similar to that of old globular clusters in the Milky Way. The major source of far-UV (FUV) flux, in this model, is relatively metal- poor and helium-enhanced hot horizontal branch stars and their progeny. The Burstein et al. (1988) relation of the F U V - V color with metallicity is also explained either by the variation of the fraction of helium-enhanced subpopulations or by the spread in mean age of stellar populations in early-type galaxies.

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 non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

Recent studies find that some early-type galaxies host Type II or Ibc supernovae (SNe II, Ibc). This may imply recent star-formation activities in these SNe host galaxies, but a massive star origin of the SNe Ib so far observed in early-type galaxies has been questioned because of their intrinsic faintness and unusually strong Ca lines shown in the nebular phase. To address the issue, we investigate the properties of early-type SNe host galaxies using the data with Galaxy Evolution Explore(GALEX) ultraviolet photometry, and the Sloan Digital Sky Survey (SDSS) optical data. Our sample includes eight SNe II and one peculiar SN Ib (SN 2000ds) host galaxies as well as 32 SN Ia host galaxies. The host galaxy of SN 2005cz, another peculiar SN Ib, is also analysed using the GALEX data and the NASA/IPAC Extragalactic Database (NED) optical data. We find that the NUV-optical colors of SN II/Ib host galaxies are systematically bluer than those of SN Ia host galaxies, and some SN II/Ib host galaxies with NUV-r colors markedly bluer than the others exhibit strong radio emission. We perform a stellar population synthesis analysis and find a clear signature of recent star-formation activities in most of the SN II/Ib host galaxies. Our results generally support the association of the SNe II/Ib hosted in early-type galaxies with core-collapse of massive stars. We briefly discuss implications for the progenitors of the peculiar SNe Ib 2000ds and 2005cz.

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Msun for the primary components, and of single helium stars with initial masses of 2.8 - 20 Msun. We find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of angular momentum retained in the core at the presupernova stage in different binary model sequences converge to a value similar to those found in previous single star models. This amount is large enough to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary systems are predicted to have a wide range of final masses even up to 7 Msun, with helium envelopes of 0.16 - 1.5 Msun. Our results indicate that, if the lack of helium lines in the spectra of SNe Ic were due to small amounts of helium, the distribution of both initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin hydrogen layer (0.001 - 0.01 Msun) is expected to be present in many SN Ib progenitors at the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova explosion. This may have important consequences for the shock break-out and supernova light curve. We also argue that some SN progenitors with thin hydrogen layers produced via Case AB/B transfer might be related to Type IIb supernova progenitors with relatively small radii of about 10 Rsun.

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.

There is a growing body of evidence for the presence of multiple stellar populations in some globular clusters, including NGC 1851. For most of these peculiar globular clusters, however, the evidence for the multiple red giant-branches (RGBs) having different heavy elemental abundances as observed in Omega Centauri is hitherto lacking, although spreads in some lighter elements are reported. It is therefore not clear whether they also share the suggested dwarf galaxy origin of Omega Cen or not. Here we show from the CTIO 4m UVI photometry of the globular cluster NGC 1851 that its RGB is clearly split into two in the U - I color. The two distinct RGB populations are also clearly separated in the abundance of heavy elements as traced by Calcium, suggesting that the type II supernovae enrichment is also responsible, in addition to the pollutions of lighter elements by intermediate mass asymptotic giant branch stars or fast-rotating massive stars. The RGB split, however, is not shown in the V - I color, as indicated by previous observations. Our stellar population models show that this and the presence of bimodal horizontal-branch distribution in NGC 1851 can be naturally reproduced if the metal-rich second generation stars are also enhanced in helium.

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.

We follow the bright, highly energetic afterglow of Swift-discovered GRB 080721 at z=2.591 out to 36 days or 3e6 s since the trigger in the optical and X-ray bands. We do not detect a break in the late-time light curve inferring a limit on the opening angle of theta_j >= 7.3 deg and setting tight constraints on the total energy budget of the burst of E_gamma >= 9.9e51 erg within the fireball model. To obey the fireball model closure relations the GRB jet must be expanding into a homogeneous surrounding medium and likely lies behind a significant column of dust. The energy constraint we derive can be used as observational input for models of the progenitors of long gamma-ray bursts: we discuss how such high collimation-corrected energies could be accommodated with certain parameters of the standard massive star core-collapse models. We can, however, most probably rule out a magnetar progenitor for this GRB which would require 100% efficiency to reach the observed total energy.

We explore the age distribution of the globular cluster (GC) system of the nearby elliptical galaxy NGC 5128 using ultraviolet (UV) photometry from Galaxy Evolution Explorer (GALEX) observations, with UV - optical colors used as the age indicator. Most GCs in NGC 5128 follow the general trends of GCs in M31 and Milky Way in UV - optical color-color diagram, which indicates that the majority of GCs in NGC 5128 are old similar to the age range of old GCs in M31 and Milky Way. A large fraction of spectroscopically identified intermediate-age GC (IAGC) candidates with ~ 3-8 Gyr are not detected in the FUV passband. Considering the nature of intermediate-age populations being faint in the far-UV (FUV) passband, we suggest that many of the spectroscopically identified IAGCs may be truly intermediate in age. This is in contrast to the case of M31 where a large fraction of spectroscopically suggested IAGCs are detected in FUV and therefore may not be genuine IAGCs but rather older GCs with developed blue horizontal branch stars. Our UV photometry strengthens the results previously suggesting the presence of GC and stellar subpopulation with intermediate age in NGC 5128. The existence of IAGCs strongly indicates the occurrence of at least one more major star formation episode after a starburst at high redshift.

Strong Balmer absorption lines and the lack of Ha and [OII] emission lines signify that E+As are post-starburst systems. Recent studies suggest that E+As may undergo the transition from the `blue cloud' to the `red sequence' and eventually migrate to red sequence ETGs. An observational validation of this scenario is to identify the intervening galaxy population between E+As and the red-sequence. Motivated by recent findings with GALEX that a large fraction of ETGs exhibit UV-excess as a sign of RSF, we investigate the possible connection of the UV-excess galaxies to E+As. In particular, we examine the FP scaling relations of the largest sample of ~1,000 E+As selected from the SDSS and ~20,000 morphologically-selected SDSS ETGs with GALEX UV data. The FP parameters, combined with stellar population indicators, reveal a certain group of UV-excess ETGs that bridges between E+As and quiescent red galaxies. The newly identified galaxies are the post-starburst systems characterized by UV-excess but no Ha emission. This is a conceptual generalisation of "E+A", in that the Balmer absorption line in the "E+A" definition is replaced with UV-optical colours that are far more sensitive to RSF than the Balmer lines. We refer to these UV-excess galaxies as "E+a" galaxies, which stands for elliptical ("E") galaxies with a minority of A-type ("a") young stars. The species are either (1) galaxies that experienced starbursts weaker than those observed in E+As (1~10% of E+As, "mild E+As") or (2) the products of passively evolved E+As after quenching star formation quite a while ago (~1 Gyr, "old E+As"). We suggest that the latter type of E+a galaxies represents the most recent arrival to the red sequence in the final phase of the "E+A" to "red early-type" transition. (Abridged)

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 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.

Recent theoretical studies have revealed the possibly important role of the capture and annihilation process of weakly interacting massive particles (WIMPs) for the first stars. Using new evolutionary models of metal-free massive stars, we investigate the impact of such ``dark matter burning'' for the first stars in different environments of dark matter (DM) halos, in terms of the ambient WIMP density (rho_chi). We find that, in agreement with existing literature, stellar life times can be significantly prolonged for a certain range of rho_\chi (i.e., 10^{10} ~ < \rho_\chi GeV/cm3 ~ < 10^{11} with the current upper limit for the spin-dependent elastic scattering cross section sigma = 5*10^{-39} cm2}). This greatly enhances the role of rotationally induced chemical mixing in rotating stars, in favor of abundant production of primary nitrogen, massive helium stars and long gamma-ray bursts, from the first stars. We also find that stars with rho_chi > 2*10^{11} GeV/cm3 may not undergo nuclear burning stages, confirming the previous work, and that ionizing photon fluxes from such DM supported stars are very weak. Delayed metal enrichment and slow reionization in the early universe would have resulted if most of the first stars had been born in DM halos with such high rho_\chi, unless it had been lowered significantly below the threshold for efficient DM burning on a short time scale.

In recent years, increasing evidence for chemical complexity and multiple stellar populations in massive globular clusters (GCs) has emerged, including extreme horizontal branches (EHBs) and UV excess. Our goal is to improve our understanding of UV excess in the regime of both massive GCs and ultra-compact dwarf galaxies (UCDs). To this end, we use deep archival GALEX data of the central Fornax cluster to measure NUV and FUV magnitudes of UCDs and massive GCs. We obtain NUV photometry for a sample of 35 compact objects with -13.5<M_V<-10 mag. Of those, 21 objects also have FUV photometry. Roughly half of the sources fall into the UCD luminosity regime (M_V <=-11 mag). We find that seven out of 17 massive Fornax GCs exhibit a NUV excess with respect to expectations from stellar population models, even for models with enhanced Helium abundance. This suggests that not only He-enrichment has contributed to forming the EHB population of these GCs. The GCs extend to stronger UV excess than GCs in M31 and massive GCs in M87, at the 97% confidence level. Most of the UCDs with FUV photometry also show evidence for UV excess, but their UV colours can be matched by isochrones with enhanced Helium abundances and old ages 12-14 Gyrs. We find that Fornax compact objects with X-ray emission detected from Chandra images are almost disjunct in colour from compact objects with GALEX UV detection, with only one X-ray source among the 35 compact objects. However, since this source is one of the three most UV bright GCs, we cannot exclude that the physical processes causing X-ray emission also contribute to some of the observed UV excess.

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.

The prediction of the spins of the compact remnants is a fundamental goal of the theory of stellar evolution. Here, we confront the predictions for white dwarf spins from evolutionary models including rotation with observational constraints. We perform stellar evolution calculations for stars in the mass range 1... 3$\mso$, including the physics of rotation, from the zero age main sequence into the TP-AGB stage. We calculate two sets of model sequences, with and without inclusion of magnetic fields. From the final computed models of each sequence, we deduce the angular momenta and rotational velocities of the emerging white dwarfs. While models including magnetic torques predict white dwarf rotational velocities between 2 and 10 km s$^{-1}$, those from the non-magnetic sequences are found to be one to two orders of magnitude larger, well above empirical upper limits. We find the situation analogous to that in the neutron star progenitor mass range, and conclude that magnetic torques may be required in order to understand the slow rotation of compact stellar remnants in general.

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 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.

Despite the efforts of the past decade, the origin of the bimodal horizontal-branch (HB) found in some globular clusters (GCs) remains a conundrum. Inspired by the discovery of multiple stellar populations in the {\it most massive} Galactic GC, $\omega$ Centauri, we investigate the possibility that two distinct populations may coexist and are responsible for the bimodal HBs in the {\it third} and {\it fifth} brightest GCs, NGC 6388 and NGC 6441. Using the population synthesis technique, we examine two different chemical ``self-enrichment'' hypotheses in which a primordial GC was sufficiently massive to contain two or more distinct populations as suggested by the populations found in $\omega$ Cen: (1) the age-metallicity relation scenario in which two populations with different metallicity and age coexist, following an internal age-metallicity relation, and (2) the super-helium-rich scenario in which GCs contain a certain fraction of helium-enhanced stars, for instance, the second generation stars formed from the helium-enriched ejecta of the first. The comparative study indicates that the detailed color-magnitude diagram morphologies and the properties of the RR Lyrae variables in NGC 6388 and NGC 6441 support the latter scenario; i.e., the model which assumes a minor fraction ($\sim$ 15 %) of helium-excess (Y $\simeq$ 0.3) stars. The results suggest that helium content is the main driver behind the HB bimodality found most often in massive GCs. If confirmed, the GC-to-GC variation of helium abundance should be considered a {\it local} effect, further supporting the argument that age is the {\it global} second parameter of HB morphology.

Rapidly rotating, chemically homogeneously evolving massive stars are considered to be progenitors of long gamma-ray bursts. We present numerical simulations of the evolution of the circumstellar medium around a rapidly rotating 20 Msol star at a metallicity of Z=0.001. Its rotation is fast enough to produce quasi-chemically homogeneous evolution. While conventionally, a star of 20 Msol would not evolve into a Wolf-Rayet stage, the considered model evolves from the main sequence directly to the helium main sequence. We use the time-dependent wind parameters, such as mass loss rate, wind velocity and rotation-induced wind anisotropy from the evolution model as input for a 2D hydrodynamical simulation. While the outer edge of the pressure-driven circumstellar bubble is spherical, the circumstellar medium close to the star shows strong non-spherical features during and after the periods of near-critical rotation. We conclude that the circumstellar medium around rapidly rotating massive stars differs considerably from the surrounding material of non-rotating stars of similar mass. Multiple blue-shifted high velocity absorption components in gamma-ray burst afterglow spectra are predicted. As a consequence of near critical rotation and short stellar evolution time scales during the last few thousand years of the star's life, we find a strong deviation of the circumstellar density profile in the polar direction from the 1/R^2 density profile normally associated with stellar winds close to the star

We demonstrate that the combination of GALEX UV photometry in the FUV (~1530 angstroms) and NUV (~2310 angstroms) passbands with optical photometry in the standard U,B,V,R,I filters can efficiently break the age-metallicity degeneracy. We estimate well-constrained ages, metallicities and their associated errors for 42 GCs in M31, and show that the full set of FUV,NUV,U,B,V,R,I photometry produces age estimates that are ~90 percent more constrained and metallicity estimates that are ~60 percent more constrained than those produced by using optical filters alone. The quality of the age constraints is comparable or marginally better than those achieved using a large number of spectrscopic indices.

We have studied ~2100 early-type galaxies in the SDSS DR3 which have been detected by the GALEX Medium Imaging Survey (MIS), in the redshift range 0 < z < 0.11. Combining GALEX UV photometry with corollary optical data from the SDSS, we find that, at a 95 percent confidence level, at least ~30 percent of galaxies in this sample have UV to optical colours consistent with some recent star formation within the last Gyr. In particular, galaxies with a NUV - r colour less than 5.5 are very likely to have experienced such recent star formation, taking into account the possibility of a contribution to NUV flux from the UV upturn phenomenon. We find quantitative agreement between the observations and the predictions of a semi-analytical LCDM hierarchical merger model and deduce that early-type galaxies in the redshift range 0 < z < 0.11 have ~1 to 3 percent of their stellar mass in stars less than 1 Gyr old. The average age of this recently formed population is ~300 to 500 Myrs. We also find that monolithically evolving galaxies, where recent star formation can be driven solely by recycled gas from stellar mass loss, cannot exhibit the blue colours (NUV - r < 5.5) seen in a significant fraction (~30 percent) of our observed sample.

The evolution of single stars at low metallicity has attracted a large interest, while the effect of metallicity on binary evolution remains still relatively unexplored. We study the effect of metallicity on the number of binary systems that undergo different cases of mass transfer. We find that binaries at low metallicity are more likely to start transferring mass after the onset of central helium burning, often referred to as case C mass transfer. In other words, the donor star in a metal poor binary is more likely to have formed a massive CO core before the onset of mass transfer. At solar metallicity the range of initial binary separations that result in case C evolution is very small for massive stars, because they do not expand much after the ignition of helium and because mass loss from the system by stellar winds causes the orbit to widen, preventing the primary star to fill its Roche lobe. This effect is likely to have important consequences for the metallicity dependence of the formation rate of various objects through binary evolution channels, such as long GRBs, double neutron stars and double white dwarfs.

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.

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.

The observed association of Long Gamma-Ray Bursts (LGRBs) with peculiar Type Ic supernovae gives support to Woosley`s collapsar/hypernova model, in which the GRB is produced by the collapse of the rapidly rotating core of a massive star to a black hole. The association of LGRBs with small star-forming galaxies suggests low-metallicity to be a condition for a massive star to evolve to the collapsar stage. Both completely-mixed single star models and binary star models are possible. In binary models the progenitor of the GRB is a massive helium star with a close companion. We find that tidal synchronization during core-helium burning is reached on a short timescale (less than a few millennia). However, the strong core-envelope coupling in the subsequent evolutionary stages is likely to rule out helium stars with main-sequence companions as progenitors of hypernovae/GRBs. On the other hand, helium stars in close binaries with a neutron-star or black-hole companion can, despite the strong core-envelope coupling in the post-helium burning phase, retain sufficient core angular momentum to produce a hypernova/GRB.

We systematically explore the evolution of the merger of two carbon-oxygen (CO) white dwarfs. The dynamical evolution of a 0.9 Msun + 0.6 Msun CO white dwarf merger is followed by a three-dimensional SPH simulation. We use an elaborate prescription in which artificial viscosity is essentially absent, unless a shock is detected, and a much larger number of SPH particles than earlier calculations. Based on this simulation, we suggest that the central region of the merger remnant can, once it has reached quasi-static equilibrium, be approximated as a differentially rotating CO star, which consists of a slowly rotating cold core and a rapidly rotating hot envelope surrounded by a centrifugally supported disc. We construct a model of the CO remnant that mimics the results of the SPH simulation using a one-dimensional hydrodynamic stellar evolution code and then follow its secular evolution. The stellar evolution models indicate that the growth of the cold core is controlled by neutrino cooling at the interface between the core and the hot envelope, and that carbon ignition in the envelope can be avoided despite high effective accretion rates. This result suggests that the assumption of forced accretion of cold matter that was adopted in previous studies of the evolution of double CO white dwarf merger remnants may not be appropriate. Our results imply that at least some products of double CO white dwarfs merger may be considered good candidates for the progenitors of Type Ia supernovae. In this case, the characteristic time delay between the initial dynamical merger and the eventual explosion would be ~10^5 yr. (Abridged).

We present the GALEX UV photometry of the elliptical galaxies in Abell clusters at moderate redshifts (z < 0.2) for the study of the look-back time evolution of the UV upturn phenomenon. The brightest elliptical galaxies (M_r < -22) in 12 remote clusters are compared with the nearby giant elliptical galaxies of comparable optical luminosity in the Fornax and Virgo clusters. The sample galaxies presented here appear to be quiescent without signs of massive star formation or strong nuclear activity, and show smooth, extended profiles in their UV images indicating that the far-UV (FUV) light is mostly produced by hot stars in the underlying old stellar population. Compared to their counterparts in nearby clusters, the FUV flux of cluster giant elliptical galaxies at moderate redshifts fades rapidly with ~ 2 Gyrs of look-back time, and the observed pace in FUV - V color evolution agrees reasonably well with the prediction from the population synthesis models where the dominant FUV source is hot horizontal-branch stars and their progeny. A similar amount of color spread (~ 1 mag) in FUV - V exists among the brightest cluster elliptical galaxies at z ~ 0.1, as observed among the nearby giant elliptical galaxies of comparable optical luminosity.

We study the UV and optical properties of 38 massive GCs in M87, imaged using the HST/STIS and WFPC2 instruments. The majority of these GCs appear extremely bright in the FUV - roughly a magnitude brighter than their Galactic counterparts with similar metallicities. The observed FUV flux is several times larger than predictions of canonical old stellar population models. These canonical models, which assume a fiducial helium enrichment parameter, dY/dZ=2, are able to reproduce the observed FUV fluxes only if ages are 3--5 Gyr larger than the WMAP age of the Universe, although the same models fit the UV photometry of Galactic and M31 GCs for ages less than the WMAP age. A similar discrepancy is found between the mass-weighted and UV-luminosity weighted ages of the massive Galactic GC omega Cen, whose CMD (including peculiar features on its well-populated horizontal branch) can be accurately reproduced by invoking a small super-He-rich (dY/dZ ~ 90) stellar component. By comparison to omega Cen, we propose that the majority of M87 GCs in our sample contain strong signatures of similarly minor super-He-rich sub-components. Although we cannot prove or disprove the extreme helium scenario at the moment, we show that the same phenomenon that causes the extended horizontal branch of omega Cen explains the UV brightness of our sample. If this is indeed due to the extreme helium, this study would be the first to find its signatures in extragalactic objects.

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.

Although there is strong support for the collapsar engine as the power source of long-duration gamma-ray bursts (GRBs), we still do not definitively know the progenitor of these explosions. Here we review the current set of progenitor scenarios for long-duration GRBs and the observational constraints on these scenarios. Examining these, we find that single-star models cannot be the only progenitor for long-duration GRBs. Several binary progenitors can match the solid observational constraints and also have the potential to match the trends we are currently seeing in the observations. Type Ib/c supernovae are also likely to be produced primarily in binaries; we discuss the relationship between the progenitors of these explosions and those of the long-duration GRBs.

We present high resolution VLT UVES and low resolution FORS optical spectroscopy of supernova 2006aj and its host galaxy, associated with the nearby (z = 0.03342) gamma-ray burst GRB 060218. This host galaxy is a unique case, as it is one of the few nearby GRB host galaxies known, and it is only the second time high resolution spectra have been taken of a nearby GRB host galaxy (after GRB 980425). The resolution, wavelength range and S/N of the UVES spectrum combined with low resolution FORS spectra allow a detailed analysis of the circumburst and host galaxy environments. We analyse the emission and absorption lines in the spectrum, combining the high resolution UVES spectrum with low resolution FORS spectra and find the metallicity and chemical abundances in the host. We probe the geometry of the host by studying the emission line profiles. Our spectral analysis shows that the star forming region in the host is metal poor with 12 + log(O/H) = 7.54 (+0.17, -0.10) (~0.07 Z_sun), placing it among the most metal deficient subset of emission-line galaxies. It is also the lowest metallicity found so far for a GRB host from an emission-line analysis. Given the stellar mass of the galaxy of ~10^7 M_sun and the SFR (H alpha) = 0.065 +/- 0.005 M_sun/yr, the high specific star formation rate indicates an age for the galaxy of less than ~200 Myr. The brightest emission lines are clearly asymmetric and are well fit by two Gaussian components separated by ~22$ km/s. We detect two discrete Na I and Ca II absorption components at the same redshifts as the emission components. We tentatively interpret the two components as arising from two different starforming regions in the host, but high resolution imaging is necessary to confirm this.

We present ultraviolet (UV) photometry of M31 globular clusters (GCs) found in 23 Galaxy Evolution Explorer (GALEX) images covering the entirety of M31. We detect 485 and 273 GCs (and GC candidates) in the near-ultraviolet (NUV; 2267 A) and far-ultraviolet (FUV; 1516 A), respectively. Comparing M31 data with those of Galactic GCs in the UV with the aid of population models, we find that the age ranges of old GCs in M31 and the Galactic halo are similar. Three metal-rich ([Fe/H]>-1) GCs in M31 produce significant FUV flux making their FUV-V colors unusually blue for their metallicities. These are thought to be analogs of the two peculiar Galactic GCs NGC 6388 and NGC 6441 with extended blue HB stars. Based on the models incorporating helium enriched subpopulations in addition to the majority of the population that have a normal helium abundance, we suggest that even small fraction of super-helium-rich subpopulations in GCs can reproduce the observed UV bright metal-rich GCs. Young clusters in M31 show distinct UV and optical properties from GCs in Milky Way. Population models indicate that their typical age is less than ~ 2 Gyrs. A large fraction of young GCs have the kinematics of the thin, rapidly rotating disk component. However, a subset of the old GCs also shares the thin-disk kinematics of the younger clusters. The existence of young GCs on the outskirts of M31 disk suggests the occurrence of a significant recent star formation in the thin-disk of M31. Old thin-disk GCs may set constraints on the epoch of early formation of the M31 thin-disk. We detect 12 (10) intermediate-age GC candidates in NUV (FUV). We suggest that some of spectroscopically identified intermediate-age GCs may not be truly intermediate in age, but rather older GCs that possess developed HB.

We present grids of massive star evolution models at four different metallicities (Z=0.004, 0.002, 0.001, 0.00001). The effects of rotation on the stellar structure and the transport of angular momentum and chemical elements through the Spruit-Tayler dynamo and rotationally induced instabilities are considered. After discussing uncertainties involved with the adopted physics, we elaborate the final fate of massive stars as a function of initial mass and spin rate, at each considered metallicity. In particular, we investigate for which initial conditions long gamma-ray bursts (GRBs) are expected to be produced in the frame of the collapsar model. Then, using an empirical spin distribution of young massive metal-poor stars and a specified metallicity-dependent history of star-formation, we compute the expected GRB rate as function of metallicity and redshift based on our stellar evolution models. The GRB production in our models is limited to metallicities of Z \lsim 0.004, with the consequence that about 50 % of all GRBs are predicted to be found at redshifts above z = 4, with most supernovae occurring at redshifts below z\simeq 2.2. The average GRB/SN ratio predicted by our model is about 1/200 globally, and 1/1250 at low redshift. Future strategies for testing the considered GRB progenitor scenario are briefly discussed.

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method. Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The present day v_sini distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v_r) distribution that can be characterised by a mean velocity of about 160-190 km/s and an effective half width of 100-150 km/s. The v_r distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct, the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring v_r. The mass loss rates of the SMC objects having luminosities of log L/L_sun > 5.4 are in excellent agreement with predictions. However, for lower luminosity stars the winds are too weak to determine M_dot accurately from the optical spectrum. Two of three spectroscopically classified Vz stars from our sample are located close to the theoretical zero age main sequence, as expected.

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-Msun and a 1.92-Msun model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50-100ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a neutrino-driven wind that develops within the white dwarf, in a cone of ~40deg opening angle about the poles, with a mass loss rate of 5-8 x 10^{-3} Msun/yr. The ejecta have an entropy on the order of 20-50 k_B/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, at >600ms after bounce, the explosion energy has reached 3-4 x 10^{49}erg and the total ejecta mass has reached a few times 0.001Msun. We estimate the asymptotic explosion energies to be lower than 10^{50}erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion. We document in detail the numerous effects of the fast rotation of the progenitors: The neutron stars are aspherical; the ``nu_mu'' and anti-nu_e neutrino luminosities are reduced compared to the nu_e neutrino luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction depend strongly upon latitude (a la von Zeipel); and a quasi-Keplerian 0.1-0.5-Msun accretion disk is formed.

The UV upturn is the rising flux with decreasing wavelength between the Lyman limit and 2500\AA found virtually in all bright spheroidal galaxies. It has been a mystery ever since it was first detected by the OAO-2 space telescope (Code & Welch 1979) because such old metal-rich populations were not expected to contain any substantial number of hot stars. It was confirmed by following space missions, ANS (de Boer 1982), IUE (Bertola et al. 1982) and HUT (Brown et al. 1997). The positive correlation between the UV-to-optical colour (i.e., the strength of the UV upturn) and the Mg2 line strength found by Burstein et al. (1987) through IUE observations has urged theorists to construct novel scenarios in which metal-rich ($\gtrsim Z_{\odot}$) old ($\gtrsim$ a few Gyr) stars become UV bright (Greggio & Renzini 1990; Horch et al. 1992). Also interesting was to find using HUT that, regardless of the UV strength, the UV spectral slopes at 1000--2000\AA in the six UV bright galaxies were nearly identical suggesting a very small range of temperatures of the UV sources in these galaxies (Brown et al. 1997), which corresponds to $T_{\rm eff} \approx 20,000 \pm 3,000$ K. This, together with other evidence, effectively ruled out young stars as the main driver of the UV upturn. A good review on the observational side of the story is given in the next article by Tom Brown, as well as in the recent articles of Greggio & Renzini (1999) and O'Connell (1999).

The colors of globular clusters in most of large elliptical galaxies are bimodal. This is generally taken as evidence for the presence of two cluster subpopulations that have different geneses. However, here we find that, because of the non-linear nature of the metallicity-to-color transformation, a coeval group of old clusters with a unimodal metallicity spread can exhibit color bimodality. The models of cluster colors indicate that the horizontal-branch stars are the main drivers behind the empirical non-linearity. We show that the scenario gives simple and cohesive explanations for all the key observations, and could simplify theories of elliptical galaxy formation.

We use \textit{GALEX} (Galaxy Evolution Explorer) near-UV (NUV) photometry of a sample of early-type galaxies selected in \textit{SDSS} (Sloan Digital Sky Survey) to study the UV color-magnitude relation (CMR). $NUV-r$ color is an excellent tracer of even small amounts ($\sim 1$% mass fraction) of recent ($\la 1$ Gyr) star formation and so the $NUV-r$ CMR allows us to study the effect of environment on the recent star formation history. We analyze a volume-limited sample of 839 visually-inspected early-type galaxies in the redshift range $0.05 < z < 0.10$ brighter than $M_{r}$ of -21.5 with any possible emission-line or radio-selected AGN removed to avoid contamination. We find that contamination by AGN candidates and late-type interlopers highly bias any study of recent star formation in early-type galaxies and that, after removing those, our lower limit to the fraction of massive early-type galaxies showing signs of recent star formation is roughly $30 \pm 3%$ This suggests that residual star formation is common even amongst the present day early-type galaxy population. We find that the fraction of UV-bright early-type galaxies is 25% higher in low-density environments. However, the density effect is clear only in the lowest density bin. The blue galaxy fraction for the subsample of the brightest early-type galaxies however shows a very strong density dependence, in the sense that the blue galaxy fraction is lower in a higher density region.

We deduce the globular cluster formation history of the nearby elliptical NGC 5128 by using a phenomenological chemical enrichment model to reproduce its observed metallicity distribution function (MDF). We derive the observed MDF using recently obtained U and B photometry of the NGC 5128 GC system, with (U-B) used as the metallicity indicator. We test whether the observed MDF can be reproduced using a single starburst at high redshift followed by passive evolution (monolithic scenario). Our results strongly indicate that this scenario is incapable of reproducing the observed MDF of the NGC 5128 GC system. At least one more major starburst is required. A persistent feature of our best-fit models is a strong GC formation episode in the first 1-2 Gyrs, followed by a second formation episode typically after 10-12 Gyrs. Roughly one-third of the stellar mass in the GC system is created in this late starburst, producing clusters that are potentially less than 1 Gyr old. Although we do not have adequate resolution in the model or in the observed MDF to make a reliable claim about the possibility of additional GC formation episodes, our results suggest that the GCs in this galaxy did not form solely in a monolithic collapse and seem to support observational evidence which indicates that NGC 5128 is a merger remnant.

Despite the growing evidence that long Gamma-Ray Bursts (GRBs) are associated with deaths of Wolf-Rayet stars, the evolutionary path of massive stars to GRBs and the exact nature of GRB progenitors remained poorly known. However, recent massive star evolutionary models indicate that -- for sufficiently low metallicity -- initially very rapidly rotating stars can satisfy the conditions for collapsar formation. Even though magnetic torques are included in these models, a strong core spin-down is avoided through quasi-chemically homogeneous evolution induced by rotational mixing. Here, we explore for which initial mass and spin-range single stars of Z=Zsun/20 are expected to produce GRBs. We further find a dichotomy in the chemical structure of GRB progenitors, where lower initial masses end their lives with a massive helium envelope which still contains some amounts of hydrogen, while higher initial masses explode with C/O-dominated hydrogen-free atmospheres.

Recent models of rotating massive stars including magnetic fields prove it difficult for the cores of single stars to retain enough angular momentum to produce a collapsar and gamma-ray burst. At low metallicity, even very massive stars may retain a massive hydrogen envelope due to the weakness of the stellar winds, posing an additional obstacle to the collapsar model. Here, we consider the evolution of massive, magnetic stars where rapid rotation induces almost chemically homogeneous evolution. We find that in this case, the requirements of the collapsar model are rather easily fulfilled if the metallicity sufficiently small: 1) Rapidly rotating helium stars are formed without the need to remove the hydrogen envelope, avoiding mass-loss induced spin-down. 2) Angular momentum transport from the helium core to hydrogen envelope by magnetic torques is insignificant. We demonstrate this by calculating evolutionary models of massive stars with various metallicities, and derive an upper metallicity limit for this scenario based on currently proposed mass loss rates. Our models also suggest the existence of a lower CO-core mass limit of about 10 Msun -- which relates to an initial mass of only about 20 Msun within our scenario -- for GRB production. We argue that the relative importance of the considered GRB progenitor channel, compared to any channel related to binary stars, may increase with decreasing metallicity, and that this channel might be the major path to GRBs from first stars.

The collapsar model for gamma-ray bursts requires three essential ingredients: a massive core, removal of the hydrogen envelope, and enough angular momentum in the core. We study current massive star evolution models of solar metallicity to determine which massive star physics is capable of producing these ingredients. In particular, we investigate the role of hydrodynamic and magnetic internal angular momentum transport and binary mass and angular momentum transfer. We follow the evolution of rotating single stars and of binary systems that include rotational processes for both stars. Neglecting magnetic fields, we show that the cores of massive single stars can maintain a high specific angular momentum when evolved with the assumption that mean molecular weight gradients suppress rotational mixing processes. In binary systems that undergo mass transfer during core hydrogen burning the mass receiving star accretes large amounts of high angular momentum material, leading to a spin-up of the core. We find, however, that this merely compensates for the tidal angular momentum loss due to spin-orbit coupling, which leads to synchronous rotation before the mass transfer event. Therefore the resulting cores do not rotate faster than in single stars. We also present models that include magnetic fields generated by differential rotation and we consider the internal angular momentum transport by magnetic torques. We investigate the capability of magnetic torques to efficiently pump angular momentum into the cores of accreting stars. Despite our finding that this mechanism works, the magnetic coupling of core and envelope after the accreting star ends core hydrogen burning leads to slower rotation than in the non-magnetic case.

A recent study by Yoon & Langer (2004a) indicated that the inner cores of rapidly accreting (Mdot > 10^{-7} M_sun/yr) CO white dwarfs may rotate differentially, with a shear rate near the threshold value for the onset of the dynamical shear instability. Such differentially rotating white dwarfs obtain critical masses for thermonuclear explosion or electron-capture induced collapse which significantly exceed the canonical Chandrasekhar limit. Here, we construct two-dimensional differentially rotating white dwarf models with rotation laws resembling those of the one-dimensional models of Yoon & Langer (2004a). We derive analytic relations between the white dwarf mass, its angular momentum, and its rotational-, gravitational- and binding energy. We show that these relations are applicable for a wide range of angular velocity profiles, including solid body rotation. We demonstrate that pre-explosion and pre-collapse conditions of both, rigidly and differentially rotating white dwarfs are well established by the present work, which may facilitate future multi-dimensional simulations of Type Ia supernova explosions and studies of the formation of millisecond pulsars and gamma-ray bursts from collapsing white dwarfs.Our results lead us to suggest various possible evolutionary scenarios for progenitors of Type Ia supernovae, leading to a new paradigm of a variable mass of exploding white dwarfs, at values well above the classical Chandrasekhar mass. Based on our 2D-models, we argue for the supernova peak brightness being proportional to the white dwarf mass, which could explain various aspects of the diversity of Type Ia supernovae, such as their variation in brightness, the dependence of their mean luminosity on the host galaxy type, and the weak correlation between ejecta velocity and peak brightness.

Recent observations for the color-magnitude diagrams (CMDs) of the massive globular cluster Omega Centauri have shown that it has a striking double main sequence (MS), with a minority population of bluer and fainter MS well separated from a majority population of MS stars. Here we confirm, with the most up-to-date Y2 isochrones, that this special feature can only be reproduced by assuming a large variation (Delta Y = 0.15) of primordial helium abundance among several distinct populations in this cluster. We further show that the same helium enhancement required for this special feature on the MS can by itself reproduce the extreme horizontal-branch (HB) stars observed in Omega Cen, which are hotter than normal HB stars. Similarly, the complex features on the HBs of other globular clusters, such as NGC 2808, are explained by large internal variations of helium abundance. Supporting evidence for the helium-rich population is also provided by the far-UV (FUV) observations of extreme HB stars in these clusters, where the enhancement of helium can naturally explain the observed fainter FUV luminosity for these stars. The presence of super helium-rich populations in some globular clusters suggests that the third parameter, other than metallicity and age, also influences CMD morphology of these clusters.

The currently favored scenario for the progenitor evolution of Type Ia supernovae (SNe Ia) presumes that white dwarfs in close binary systems grow to the Chandrasekhar limit via mass accretion from their non-degenerate companions. However, the accreted hydrogen and/or helium usually participate thermally unstable or even violent nuclear reactions in a geometrically confined region, due to the compactness of the white dwarf. Since shell flashes induced by the thermal instability might induce significant loss of mass, efficient mass increase of white dwarfs by hydrogen and/or helium accretion has been seriously questioned. A good understanding of the stability of thermonuclear shell sources is therefore crucial in order to investigate the evolution of accreting white dwarfs as SNe Ia progenitors. Here, we present a quantitative criterion for the thermal stability of thermonuclear shell sources, and discuss the effects of rotation on the stability of helium shell burning in helium accreting CO white dwarfs with Mdot \approx 10^{-7} ... 10^{-6} Msun/yr. In particular, we show that, if the effects of rotation are properly considered, helium shell sources are significantly stabilized, which might increase the likelihood for accreting white dwarfs to grow to the Chandrasekhar limit.

We have used the GALEX UV photometric data to construct a first near-ultraviolet (NUV) color-magnitude relation (CMR) for the galaxies pre-classified as early-type by SDSS studies. The NUV CMR is a powerful tool for tracking the recent star formation history in early-type galaxies, owing to its high sensitivity to the presence of young stellar populations. Our NUV CMR for UV-weak galaxies shows a well-defined slope and thus will be useful for interpreting the restframe NUV data of distant galaxies and studying their star formation history. Compared to optical CMRs, the NUV CMR shows a substantially larger scatter, which we interpret as evidence of recent star formation activities. Roughly 15% of the recent epoch (z < 0.13) bright (M[r] < -22) early-type galaxies show a sign of recent (< 1Gyr) star formation at the 1-2% level (lower limit) in mass compared to the total stellar mass. This implies that low level residual star formation was common during the last few billion years even in bright early-type galaxies.

We present ultraviolet photometry for globular clusters (GCs) in M31 from 15 square deg of imaging using the Galaxy Evolution Explorer (GALEX). We detect 200 and 94 GCs with certainty in the near-ultraviolet (NUV; 1750 - 2750 Angstroms) and far-ultraviolet (FUV; 1350 - 1750 Angstroms) bandpasses, respectively. Our rate of detection is about 50% in the NUV and 23% in the FUV, to an approximate limiting V magnitude of 19. Out of six clusters with [Fe/H]>-1 seen in the NUV, none is detected in the FUV bandpass. Furthermore, we find no candidate metal-rich clusters with significant FUV flux, because of the contribution of blue horizontal-branch (HB) stars, such as NGC 6388 and NGC 6441, which are metal-rich Galactic GCs with hot HB stars. We show that our GALEX photometry follows the general color trends established in previous UV studies of GCs in M31 and the Galaxy. Comparing our data with Galactic GCs in the UV and with population synthesis models, we suggest that the age range of M31 and Galactic halo GCs are similar.

In order to investigate the origin of the far-UV (FUV) flux from the early-type galaxies, Galaxy Evolution Explorer (GALEX) is collecting the UV data for the elliptical-rich clusters at moderate redshifts (z < 0.2) where the dominant FUV source is predicted to be hot horizontal-branch (HB) stars and their post-HB progeny. Here we present our first result for the early-type galaxies in Abell 2670 at z = 0.076. Compared to NGC 1399, a nearby giant elliptical galaxy in the Fornax cluster, it appears that the rest-frame FUV - V color of the giant ellipticals gets redder by ~ 0.7 mag at the distance of Abell 2670 (z = 0.076; look-back time ~ 1.0 Gyr). Although a detailed comparison with the models is postponed until more cluster data are accumulated, it is interesting to note that this value is consistent with the variation predicted by the population synthesis models where the mean temperature of HB stars declines rapidly with increasing look-back time.

We present a quantitative criterion for the thermal stability of thermonuclear shell sources. We find the thermal stability of shell sources to depend on exactly three factors: they are more stable when they are geometrically thicker, less degenerate and hotter. This confirms and unifies previously obtained results in terms of the geometry, temperature and density of the shell source, by a simplified but quantitative approach to the physics of shell nuclear burning. We present instability diagrams in the temperature-density plane for hydrogen and helium shell burning, which allow a simple evaluation of the stability conditions of such shell sources in stellar models. The performance of our stability criterion is demonstrated in various numerical models: in a 3 Msun AGB star, in helium accreting CO white dwarfs, in a helium white dwarf which is covered by a thin hydrogen envelope, and in a 1.0 Msun giant.

We investigate the effects of rotation on the behavior of the helium burning shell source in accreting carbon-oxygen white dwarfs, in the context of the single degenerate Chandrasekhar mass progenitor scenario for Type Ia supernovae (SNe Ia). We model the evolution of helium accreting white dwarfs of initially 1 Msun, assuming four different constant accretion rates (2, 3, 5 and 10 times10^{-7} Msun/yr). In a one-dimensional approximation, we compute the mass accretion and subsequent nuclear fusion of helium into carbon and oxygen, as well as angular momentum accretion, angular momentum transport inside the white dwarf, and rotationally induced chemical mixing. Our models show two major effects of rotation: a) The helium burning nuclear shell source in the rotating models is much more stable than in corresponding non-rotating models -- which increases the likelihood of accreting white dwarfs to reach the stage of central carbon ignition. This effect is mainly due to rotationally induced mixing at the CO/He interface which widens the shell source, and due to the centrifugal force lowering the density and degeneracy at the shell source location. b) The C/O-ratio in the layers which experience helium shell burning -- which may affect the energy of a SN Ia explosion -- is strongly decreased by the rotationally induced mixing of alpha-particles into the carbon-rich layers. We discuss implications of our results for the evolution of SNe Ia progenitors.

Today, Type Ia supernovae are essential tools for cosmology, and recognized as major contributors to the chemical evolution of galaxies. The construction of detailed supernova progenitor models, however, was so far prevented by various physical and numerical difficulties in simulating binary systems with an accreting white dwarf component, e.g., unstable helium shell burning which may cause significant expansion and mass loss. Here, we present the first binary evolution calculation which models both stellar components and the binary interaction simultaneously, and where the white dwarf mass grows up to the Chandrasekhar limit by mass accretion. Our model starts with a 1.6 Msun helium star and a 1.0 Msun CO white dwarf in a 0.124 day orbit. Thermally unstable mass transfer starts when the CO core of the helium star reaches 0.53 Msun, with mass transfer rates of 1...8 times 10^{-6} Msun/yr. The white dwarf burns the accreted helium steadily until the white dwarf mass has reached ~ 1.3 Msun and weak thermal pulses follow until carbon ignites in the center when the white dwarf reaches 1.37 Msun. Although the supernova production rate through this channel is not well known, and this channel can not be the only one as its progenitor life time is rather short (~ 10^7 - 10^8 yr), our results indicate that helium star plus white dwarf systems form a reliable route for producing Type Ia supernovae.

We discuss the effects of rotation on the evolution of accreting carbon-oxygen white dwarfs, with the emphasis on possible consequences in Type Ia supernova (SN Ia) progenitors. Starting with a slowly rotating white dwarf, we simulate the accretion of matter and angular momentum from a quasi-Keplerian accretion disk. The role of the various rotationally induced hydrodynamic instabilities for the transport of angular momentum inside the white dwarf is investigated. We find that the dynamical shear instability is the most important one in the highly degenerate core. Our results imply that accreting white dwarfs rotate differentially throughout,with a shear rate close to the threshold value for the onset of the dynamical shear instability. As the latter depends on the temperature of the white dwarf, the thermal evolution of the white dwarf core is found to be relevant for the angular momentum redistribution. As found previously, significant rotation is shown to lead to carbon ignition masses well above 1.4 Msun. Our models suggest a wide range of white dwarf explosion masses, which could be responsible for some aspects of the diversity observed in SNe Ia. We analyze the potential role of the bar-mode and the r-mode instability in rapidly rotating white dwarfs, which may impose angular momentum loss by gravitational wave radiation. We discuss the consequences of the resulting spin-down for the fate of the white dwarf, and the possibility to detect the emitted gravitational waves at frequencies of 0.1 >... 1.0 Hz in nearby galaxies with LISA. Possible implications of fast and differentially rotating white dwarf cores for the flame propagation in exploding white dwarfs are also briefly discussed.

We present evolutionary models of helium accreting carbon-oxygen white dwarfs in which we include the effects of the spin-up of the accreting star induced by angular momentum accretion, rotationally induced chemical mixing and rotational energy dissipation. Initial masses of 0.6 Msun and 0.8 Msun and constant accretion rates of a few times 10^{-8} Msun/yr of helium rich matter have been considered, which is typical for the sub-Chandrasekhar mass progenitor scenario for Type Ia supernovae. It is found that the helium envelope in an accreting white dwarf is heated efficiently by friction in the differentially rotating spun-up layers. As a result, helium ignites much earlier and under much less degenerate conditions compared to the corresponding non-rotating case. Consequently, a helium detonation may be avoided, which questions the sub-Chandrasekhar mass progenitor scenario for Type Ia supernovae. We discuss implications of our results for the evolution of helium star plus white dwarf binary systems as possible progenitors of recurrent helium novae.

We present the results of our analysis of the photometric data of globular clusters in the elliptical galaxy NGC 5128 (Cen A). We show that the integrated colour U-B can be an effective metallicity indicator for simple stellar populations. This is because it is sensitive to metallicity via the opacity effect but relatively insensitive to the effective main sequence turn-off temperature of the population (and thus to age) when Teff ~ 7000-12000K, that is, when the Balmer Jump is strong. This flat U-B vs Teff relation is a result of the fact that the blueing effect of the optical continuum with increasing temperature is temporarily stopped by the Balmer Jump which becomes greater with temperature in this range. In this study we use U-B and B-V as metallicity and age indicators, respectively. We first show that the use of the U-B vs B-V two-colour diagram roughly yields the metallicities and ages of the Milky Way globular clusters independently determined, and then apply the technique to the clusters in NGC 5128. There is a large range in U-B, which corresponds to [Fe/H] of -2.0 through over +0.3. The large uncertainties from the data and the models currently prevent us from pinning down their ages and metallicities. Although a constant age for all these clusters cannot be ruled out, there is a hint of the metal-rich clusters being younger. Significance of these results and caveats of the analysis are discussed.

We discuss the first available binary evolution models which include up-to-date rotational physics for both components, as well as angular momentum accretion and spin-orbit coupling. These models allow a self-consistent computation of the mass transfer efficiency during Roche-lobe overflow phases, and a determination of the transition from quasi-conservative to non-conservative evolution. Applications to massive binary systems lead to predictions for the spin rates of compact objects in binaries, and for the occurrence of gamma-ray bursts from collapsars in binaries. Rotational effects in accreting white dwarfs are found to stabilise the shell burning and decrease the carbon abundance in progenitor models for Chandrasekhar-mass Type Ia supernovae, and to potentially avoid a detonation of the white dwarf within the sub-Chandrasekhar mass scenario.

One of the long-standing problems in modern astronomy is the curious division of Galactic globular clusters, the "Oosterhoff dichotomy," according to the properties of their RR Lyrae stars. Here, we find that most of the lowest-metallicity ([Fe/H]<-2.0) clusters, which are essential to an understanding of this phenomenon, display a planar alignment in the outer halo. This alignment, combined with evidence from kinematics and stellar population, indicates a captured origin from a satellite galaxy. We show that, together with the horizontal-branch evolutionary effect, the factor producing the dichotomy could be a small time-gap between the cluster-formation epochs in the Milky Way and the satellite. The results oppose the traditional view that the metal-poorest clusters represent the indigenous and oldest population of the Galaxy.

Our recent CCD photometry (Lee et al. 1999) has shown, for the first time, that omega Cen has several distinct stellar populations, which is reminiscent of the Sagittarius dwarf galaxy. Here we present more detailed analysis of the data along with the population models. We confirm the presence of several distinct red-giant-branches (RGBs) with a red metal-rich sequence well separated from other bluer metal-poor ones. Our population models suggest the red clump associated with the most metal-rich RGB is about 4 Gyr younger than the dominant metal-poor component, indicating that omega Cen was enriched over this timescale. These features, taken together with this cluster's other unusual characteristics, provide good evidence that omega Cen was once part of a more massive system that merged with the Milky Way, as the Sagittarius dwarf galaxy is in the process of doing now. Mergers probably were much more frequent in the early history of the Galaxy and omega Cen appears to be a relict of this era.

We have constructed synthetic color-magnitude diagrams (CMDs) for omega Cen and other massive globular clusters with apparently peculiar CMD morphology. Our population models, which adopt the most up-to-date input physics and parameters, show that the observed CMD of omega Cen can be reproduced by adopting (1) multimodal metallicity distribution function as derived from the observed color distribution of red-giant-branch (RGB) stars, and (2) an internal age-metallicty relation among the populations therein. Similar results were obtained for other massive globular clusters with bimodal horizontal-branches (HBs). In particular, we found that the peculiar CMD morphology (broad RGB, bimodal HB) and properties of RR Lyrae stars observed in NGC 6388 and NGC 6441 can be reproduced by the composite of two distinct populations with mild internal age-metallicity relations. This suggests that these clusters, as well as omega Cen, may represent the relicts of disrupted dwarf galaxies.

We propose a redshift estimator for the long ($T_{90} > 20$) gamma-ray bursts (GRBs) observed by the BATSE. It is based on an empirical relation between the redshift and the power-law index of power density spectra (PDSs) of the observed GRBs. This relation is constructed by using the fact that the power-law index is dependent upon a characteristic timescale of GRB light curves which are inevitably stretched by cosmological time dilation. We construct the empirical relation using both individual PDSs and averaged PDSs. An error estimates of $z$ are 1.09 and 1.11 for the empirical relation by individual PDS fits, 1.72 and 1.56 by averaged PDS fits, for the least squares fit and the maximum likelihood fit, respectively. We attempt to determine the spatial distribution of the GRBs observed by the BATSE as a function of redshifts on the basis of the resulting redshift estimator. We find that the obtained spatial distribution of the observed GRBs seems consistent with that of the GRBs whose redshifts are reported, even though the estimated errors are not very accurate. The GRBs observed by the BATSE seem distributed within $z \sim 5-6$. This result has implications on theoretical calculations of stellar formations at high redshifts and beaming geometry via a statistical study of the observed GRBs involving beaming-induced luminosity functions. We discuss such implications, and possible uncertainties of the suggested method.

We present high-precision V, B-V color-magnitude diagrams (CMDs) for the classic second parameter globular clusters M3 and M13 from wide-field deep CCD photometry. The data for the two clusters were obtained during the same photometric nights with the same instrument, allowing us to determine accurate relative ages. Based on a differential comparison of the CMDs using the Delta (B-V) method, an age difference of 1.7 +/- 0.7 Gyr is obtained between these two clusters. We compare this result with our updated horizontal-branch (HB) population models, which confirm that the observed age difference can produce the difference in HB morphology between the clusters. This provides further evidence that age is the dominant second parameter that influences HB morphology.

We present our recent revision of model constructions for the horizontal-branch (HB) morphology of globular clusters, which suggests the HB morphology is more sensitive to age compared to our earlier models. We also present our high precision CCD photometry for the classic second parameter pair M3 and M13. The relative age dating based on this photometry indicates that M13 is indeed older than M3 by 1.7 Gyr. This is consistent with the age difference predicted from our new models, which provides a further support that the HB morphology is a reliable age indicator in most population II stellar systems.

After taking into account, for the first time, the detailed systematic variation of horizontal-branch (HB) morphology with age and metallicity, our population synthesis models show that the integrated H beta index is significantly affected by the presence of blue HB stars. Our models indicate that the strength of the H beta index increases as much as 0.75 A due to blue HB stars. According to our models, a systematic difference between the globular cluster system in the Milky Way Galaxy and that in NGC 1399 in the H beta vs. Mg2 diagram is understood if globular cluster systems in giant elliptical galaxies are couple of billion years older, in the mean, than the Galactic counterpart.

The strength of the H$\beta$ index has been computed for the integrated spectra of model globular clusters from the evolutionary population synthesis. These models take into account, for the first time, the detailed systematic variation of horizontal-branch (HB) morphology with age and metallicity. Our models show that the H$\beta$ index is significantly affected by the presence of blue HB stars. Because of the contribution from blue HB stars, the H$\beta$ does not monotonically decrease as metallicity increases at a given age. Instead, it reaches a maximum strength when the distribution of HB stars is centered around 9500 K, the temperature where the H$\beta$ index becomes strongest. Our models indicate that the strength of the H$\beta$ index increases as much as 0.75 {\AA} due to the presence of blue HB stars. The comparison of the recent Keck observations of the globular cluster system in the Milky Way Galaxy with those in giant elliptical galaxies, NGC 1399 and M87, shows a systematic shift in the H$\beta$ against metallicity plane. Our models suggest that this systematic difference is understood if the globular cluster systems in giant elliptical galaxies are several billion years older, in the mean, than the Galactic counterpart. Further observations of globular cluster systems in the external galaxies from the large ground-based telescopes and space UV facilities will enable to clarify whether this difference is indeed due to the age difference or other explanations are also possible.

Surface temperature distribution of horizontal-branch (HB) stars is very sensitive to age in old stellar systems, which makes it an attractive age indicator. In this paper, we present the recent revision of our model calculations for the HB morphology of globular clusters. The result is more updated version of our earlier models (Lee, Demarque, & Zinn 1994), which suggests that the HB morphology is more sensitive (\sim 40%) to age. We also present our new model calculations on the effect of HB stars in dating old stellar systems using the H$\beta$ index. Our results indicate that the effect of HB stars is rather strong, and suggest a possibility for the systematic difference in age between the globular clusters in the Milky Way Galaxy and those in giant elliptical galaxies. Finally, we compare these results with the relative ages estimated from our far-UV dating technique for giant elliptical galaxies.