It is well known since 2010 that fullerene C60 is widespread through the interstellar space. Also, it is well known that graphene is a source material for synthesizing fullerene. Here, we simply assume the occurrence of graphene in space. Infrared spectra of graphene molecules are calculated to compare both to astronomical observational spectra and to laboratory experimental one. Model molecules for DFT calculation are selected by one simple astronomical assumption, that is, single void creation on charge neutral graphene of C13, C24 and C54, resulting C12, C23 and C53. They have a carbon pentagon ring within a hexagon network. Different void positions should be classified as different species. In every species, single void has 3 radical carbons holding 6 spins. These spins select stable state by recombination of carbon atoms to reduce molecular energy, and to give infrared spectrum due to molecular vibration. It was found that stable spin state was the triplet, not singlet. Most of species show prominent infrared spectral features closely resembling the astronomically observed infrared emission seen in the carbon rich planetary nebulae of Tc1 and Lin49. We could assign major bands at 18.9 micrometer, and sub-bands at 6.6, 7.0, 7.6, 8.1, 8.5, 9.0 and 17.4 micrometer. It is interesting that those graphene species were also assigned in the laboratory experiments on laser-induced carbon plasmas, which are analogies of carbon cluster creation in space. The conclusion is that graphene molecules could potentially contribute to the infrared emission bands of carbon-rich planetary nebulae.

We present optical and near-infrared observations of SN~Ib~2019ehk. We show that it evolved to a Ca-rich transient according to its spectral properties and evolution in late phases. It, however, shows a few distinguishable properties from the canonical Ca-rich transients: a short-duration first peak in the light curve, high peak luminosity, and association with a star-forming environment. Indeed, some of these features are shared with iPTF14gqr and iPTF16hgs, which are candidates for a special class of core-collapse SNe (CCSNe): the so-called ultra-stripped envelope SNe, i.e., a relatively low-mass He (or C+O) star explosion in a binary as a precursor of double neutron star binaries. The estimated ejecta mass ($0.43 M_\odot$) and explosion energy ($1.7 \times 10^{50} $~erg) are consistent with this scenario. The analysis of the first peak suggests existence of dense circumstellar material in the vicinity of the progenitor, implying a CCSN origin. Based on these analyses, we suggest SN 2019ehk is another candidate for an ultra-stripped envelope SN. These ultra-stripped envelope SN candidates seem to form a subpopulation among Ca-rich transients, associated with young population. We propose that the key to distinguishing this population is the early first peak in their light curves.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

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

The spatial structure of the emission lines and continuum over the 50 arcsecond extent of the nearby, O-rich, planetary nebula NGC 7009 (Saturn Nebula) have been observed with the MUSE integral field spectrograph on the ESO Very Large Telescope. Science Verification data, in <0.6 arcsecond seeing, have been reduced and analysed as images over the wavelength range 4750-9350A. Emission line maps over the bright shells are presented, from neutral to the highest ionization available (He II and [Mn V]). For collisionally excited lines (CELs), maps of electron temperature (T_e from [N II] and [S III]) and electron density (N_e from [S II] and [Cl III]) are available and for optical recombination lines (ORLs) temperature (from the Paschen jump and ratio of He I lines) and density (from high Paschen lines). These estimates are compared: for the first time, maps of the differences in CEL and ORL T_e's have been derived, and correspondingly a map of t^2 between a CEL and ORL temperature, showing considerable detail. Total abundances of He and O were formed, the latter using three ionization correction factors. However the map of He/H is not flat, departing by ~2% from a constant value, with remnants corresponding to ionization structures. Ionization correction factor methods are compared for O abundance, but none delivers a flat map. An integrated spectrum over an area of 2340 square arcseconds was also formed and compared to 1D photoionization models. The spatial variation of a range of nebular parameters illustrates the complexity of the ionized media in NGC 7009. These MUSE data are very rich with detections of many lines over areas of hundreds of square arcseconds and follow-on studies are indicated. (Abridged)

We perform a detailed analysis of the fullerene C60-containing planetary nebula (PN) SaSt2-3 to investigate the physical properties of the central star (B0-1II) and nebula based on our own Subaru/HDS spectra and multiwavelength archival data. By assessing the stellar absorption, we derive the effective temperature, surface gravity, and photospheric abundances. For the first time, we report time variability of the central star's radial velocity, strongly indicating a binary central star. Comparison between the derived elemental abundances and those predicted values by asymptotic giant branch (AGB) star nucleosynthesis models indicates that the progenitor is a star with initial mass of ~1.25 Msun and metallicity Z = 0.001/alpha-element/Cl-rich ([alpha,Cl/Fe] ~ +0.3-0.4). We determine the distance (11.33 kpc) to be consistent with the post-AGB evolution of 1.25 Msun initial mass stars with Z = 0.001. Using the photoionisation model, we fully reproduce the derived quantities by adopting a cylindrically shaped nebula. We derive the mass fraction of the C-atoms present in atomic gas, graphite grain, and C60. The highest mass fraction of C60 (~0.19%) indicates that SaSt2-3 is the C60-richest PN amongst Galactic PNe. From comparison of stellar/nebular properties with other C60 PNe, we conclude that the C60 formation depends on the central star's properties and its surrounding environment (e.g., binary disc), rather than the amount of C-atoms produced during the AGB phase.

We observed two proto-planetary nebulae, HD 56126 representing a source with an elliptical circumstellar shell and IRAS 16594-4656 representing a source with a bipolar circumstellar shell, with ALMA in the 12CO and 13CO J=3-2 lines and neighboring continuum to see how the morpho-kinematics of CO gas and dust emission properties in their circumstellar environments differ.

The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.

We perform a comprehensive analysis of the planetary nebula (PN) NGC6781 to investigate the physical conditions of each of its ionized, atomic, and molecular gas and dust components and the object's evolution, based on panchromatic observational data ranging from UV to radio. Empirical nebular elemental abundances, compared with theoretical predictions via nucleosynthesis models of asymptotic giant branch (AGB) stars, indicate that the progenitor is a solar-metallicity, 2.25-3.0 Msun initial-mass star. We derive the best-fit distance of 0.46 kpc by fitting the stellar luminosity (as a function of the distance and effective temperature of the central star) with the adopted post-AGB evolutionary tracks. Our excitation energy diagram analysis indicate high excitation temperatures in the photodissociation region (PDR) beyond the ionized part of the nebula, suggesting extra heating by shock interactions between the slow AGB wind and the fast PN wind. Through iterative fitting using the Cloudy code with empirically-derived constraints, we find the best-fit dusty photoionization model of the object that would inclusively reproduce all of the adopted panchromatic observational data. The estimated total gas mass (0.41 Msun) corresponds to the mass ejected during the last AGB thermal pulse event predicted for a 2.5 Msun initial-mass star. A significant fraction of the total mass (about 70 percent) is found to exist in the PDR, demonstrating the critical importance of the PDR in PNe that are generally recognized as the hallmark of ionized/H+ regions.

We have carried out a detailed analysis of the interesting and important very young planetary nebula (PN) Hen3-1357 (Stingray Nebula) based on a unique dataset of optical to far-IR spectra and photometric images. We calculated the abundances of nine elements using collisionally excited lines (CELs) and recombination lines (RLs). The RL C/O ratio indicates that this PN is O-rich, which is also supported by the detection of the broad 9/18 um bands from amorphous silicate grain. The observed elemental abundances can be explained by asymptotic giant branch (AGB) nucleosynthesis models for initially 1-1.5 Msun stars with Z = 0.008. The Ne overabundance might be due to the enhancement of 22Ne isotope in the He-rich intershell. By using the spectrum of the central star synthesized by Tlusty as the ionization/heating source of the PN, we constructed the self-consistent photoionization model with Cloudy to the observed quantities, and we derived the gas and dust masses, dust-to-gas mass ratio, and core-mass of the central star. About 80 % of the total dust mass is from warm-cold dust component beyond ionization front. Comparison with other Galactic PNe indicates that Hen3-1357 is an ordinary amorphous silicate rich and O-rich gas PN. Among other studied PNe, IC4846 shows many similarities in properties of the PN to Hen3-1357, although their post-AGB evolution is quite different from each other. Further monitoring observations and comparisons with other PNe such as IC4846 are necessary to understand the evolution of Hen3-1357.

We performed a detailed spectroscopic analysis of the fullerene C60-containing planetary nebula (PN) Lin49 in the Small Magellanic Cloud using XSHOOTER at the ESO VLT and the Spitzer/IRS instruments. We derived nebular abundances for nine elements. We used TLUSTY to derive photospheric parameters for the central star. Lin49 is C-rich and metal-deficient PN (Z~0.0006). The nebular abundances are in good agreement with Asymptotic Giant Branch nucleosynthesis models for stars with initial mass 1.25 Msun and metallicity Z = 0.001. Using the TLUSTY synthetic spectrum of the central star to define the heating and ionising source, we constructed the photoionisation model with CLOUDY that matches the observed spectral energy distribution (SED) and the line fluxes in the UV to far-IR wavelength ranges simultaneously. We could not fit the ~1-5 um SED using a model with 0.005-0.1 um-sized graphite grains and a constant hydrogen density shell owing to the prominent near-IR excess, while at other wavelengths the model fits the observed values reasonably well. We argue that the near-IR excess might indicate either (1) the presence of very small particles in the form of small carbon clusters, small graphite sheets, or fullerene precursors, or (2) the presence of a high-density structure surrounding the central star. We found that SMC C60 PNe show a near-IR excess component to lesser or greater degree. This suggests that these C60 PNe might maintain a structure nearby their central star.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

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

We report an investigation of the extremely metal-poor and C-rich planetary nebula (PN) K648 in the globular cluster M15 using the UV to far-IR data obtained using the Subaru, HST, FUSE, Spitzer, and Herschel. We determined the nebular abundances of ten elements. The enhancement of F ([F/H]=+0.96) is comparable to that of the halo PN BoBn1. The central stellar abundances of seven elements are determined. The stellar C/O ratio is similar to the nebular C/O ratios from recombination line and from collisionally excited line (CEL) within error, and the stellar Ne/O ratio is also close to the nebular CEL Ne/O ratio. We found evidence of carbonaceous dust grains and molecules including Class B 6-9 um and 11.3 um polycyclic aromatic hydrocarbons and the broad 11 um feature. The profiles of these bands are similar to those of the C-rich halo PNe H4-1 and BoBn1. Based on the theoretical model, we determined the physical conditions of the gas and dust and their masses, i.e., 0.048 Msun and 4.95x10^{-7} Msun, respectively. The observed chemical abundances and gas mass are in good agreement with an asymptotic giant branch nucleosynthesis model prediction for stars with an initial 1.25 Msun plus a 2.0x10^{-3} Msun partial mixing zone (PMZ) and stars with an initial mass of 1.5 Msun without a PMZ. The core-mass of the central star is approximately 0.61-0.63 Msun. K648 is therefore likely to have evolved from a progenitor that experienced coalescence or tidal disruption during the early stages of evolution, and became a ~1.25-1.5 Msun blue straggler.

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm \mu m}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu m}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

We observe a sample of 8 evolved stars in the Galactic Bulge in the CO J = 2 - 1 line using the Submillimeter Array (SMA) with angular resolution of 1 - 4 arcseconds. These stars have been detected previously at infrared wavelengths, and several of them have OH maser emission. We detect CO J = 2 - 1 emission from three of the sources in the sample: OH 359.943 +0.260, [SLO2003] A12, and [SLO2003] A51. We do not detect the remaining 5 stars in the sample because of heavy contamination from the galactic foreground CO emission. Combining CO data with observations at infrared wavelengths constraining dust mass loss from these stars, we determine the gas-to-dust ratios of the Galactic Bulge stars for which CO emission is detected. For OH 359.943 +0.260, we determine a gas mass-loss rate of 7.9 (+/- 2.2) x 10^-5 M_Sun/year and a gas-to-dust ratio of 310 (+/- 89). For [SLO2003] A12, we find a gas mass-loss rate of 5.4 (+/- 2.8) x 10^-5 M_Sun/year and a gas-to-dust ratio of 220 (+/- 110). For [SLO2003] A51, we find a gas mass-loss rate of 3.4 (+/- 3.0) x 10^-5 M_Sun/year and a gas-to-dust ratio of 160 (+/- 140), reflecting the low quality of our tentative detection of the CO J = 2 - 1 emission from A51. We find the CO J = 2 - 1 detections of OH/IR stars in the Galactic Bulge require lower average CO J = 2 - 1 backgrounds.

Fullerenes are a particularly stable class of carbon molecules in the shape of a hollow sphere or ellipsoid that might be formed in the outflows of carbon stars. Once injected into the interstellar medium (ISM), these stable species survive and are thus likely to be widespread in the Galaxy where they contribute to interstellar extinction, heating processes, and complex chemical reactions. In recent years, the fullerene species C60 (and to a lesser extent C70) have been detected in a wide variety of circumstellar and interstellar environments showing that when conditions are favourable, fullerenes are formed efficiently. Fullerenes are the first and only large aromatics firmly identified in space. The detection of fullerenes is thus crucial to provide clues as to the key chemical pathways leading to the formation of large complex organic molecules in space, and offers a great diagnostic tool to describe the environment in which they reside. Since fullerenes share many physical properties with PAHs, understanding how fullerenes form, evolve and respond to their physical environment will yield important insights into one of the largest reservoirs of organic material in space. In spite of all these detections, many questions remain about precisely which members of the fullerene family are present in space, how they form and evolve, and what their excitation mechanism is. We present here an overview of what we know from astronomical observations of fullerenes in these different environments, and discuss current thinking about the excitation process. We highlight the various formation mechanisms that have been proposed, discuss the physical conditions conducive to the formation and/or detection of fullerenes in carbon stars, and their possible connection to PAHs, HACs and other dust features.

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme. We found that five post-AGB stars showed a broad feature with a peak at 7.7 micron, that had not been classified before. Further, the 10--13 micron PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 micron rather than two distinct sharp peaks at 11.3 and 12.7 micron, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to PNe, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars.

We searched the Spitzer Space Telescope data archive for Galactic planetary nebulae (PNe), that show the characteristic 17.4 and um features due to C60, also known as buckminsterfullerene. Out of 338 objects with Spitzer/IRS data, we found eleven C60-containing PNe, six of which (Hen2-68, IC2501, K3-62, M1-6, M1-9, and SaSt2-3) are new detections. We analyzed the spectra, along with ancillary data, using the photo-ionization code CLOUDY to establish the atomic line fluxes, and determine the properties of the radiation field, as set by the effective temperature of the central star. In addition, we measured the infrared spectral features due to dust grains. We find that the Polycyclic Aromatic Hydrocarbon (PAH) profile over 6-9 um in these C60-bearing carbon-rich PNe is of the more chemically-processed class A. The intensity ratio of 3.3 um to 11.3 um PAH indicates that the number of C-atoms per PAH in C60-containing PNe is small compared to that in non-C60 PNe. The Spitzer spectra also show broad dust features around 11 and 30 um. Analysis of the 30-um feature shows that it is strongly correlated with the continuum, and we propose that a single carbon-based carrier is responsible for both the continuum and the feature. The strength of the 11-um feature is correlated to the temperature of the dust, suggesting that it is at least partially due to a solid-state carrier. The chemical abundances of C60-containing PNe can be explained by AGB nucleosynthesis models for initially 1.5-2.5 Msun stars with Z=0.004. We plotted the locations of C60-containing PNe on a face-on map of the Milky Way and we found that most of these PNe are outside the solar circle, consistent with low metallicity values. Their metallicity suggests that the progenitors are an older population.

We performed detailed chemical abundance analysis of the extremely metal-poor ([Ar/H]-2) halo planetary nebula H4-1 based on the multi-wavelength spectra from Subaru/HDS, GALEX, SDSS, and Spitzer/IRS and determined the abundances of 10 elements. The C and O abundances were derived from collisionally excited lines (CELs) and are almost consistent with abundances from recombination lines (RLs). We demonstrated that the large discrepancy in the C abundance between CEL and RL in H4-1 can be solved using the temperature fluctuation model. We reported the first detection of the [Xe III]5846 A line in H4-1 and determination of its elemental abundance ([Xe/H]>+0.48). H4-1 is the most Xe-rich PN among the Xe-detected PNe. The observed abundances are close to the theoretical prediction by a ~2.0 Msun single star model with initially r-process element rich ([r/Fe]=+2.0 dex). The observed Xe abundance would be a product of the r-process in primordial SNe. The [C/O]-[Ba/(Eu or Xe)] diagram suggests that the progenitor of H4-1 shares the evolution with two types of carbon-enhanced metal-poor stars (CEMP), CEMP-r/s and CEMP-no stars. The progenitor of H4-1 is a presumably binary formed in an r-process rich environment.

We performed multiwavelength observations of the young planetary nebula (PN) M1-11 and obtained its elemental abundances, dust mass, and the evolutionary status of the central star. The AKARI/IRC, VLT/VISIR, and Spitzer/IRS spectra show features due to carbon-rich dust, such as the 3.3, 8.6, and 11.3 um features due to polycyclic aromatic hydrocarbons (PAHs), a smooth continuum attributable to amorphous carbon, and the broad 11.5 and 30 um features often ascribed to SiC and MgS, respectively. We also report the presence of an unidentified broad feature at 16-22 um, similar to the feature found in Magellanic Cloud PNe with either C-rich or O-rich gas-phase compositions. We identify for the first time in M1-11 spectral lines at 8.5 (blended with PAH), 17.3, and 18.9 um that we attribute to the C60 fullerene. This identification is strengthened by the fact that other Galactic PNe in which fullerenes are detected, have similar central stars, similar gas-phase abundances, and a similar dust composition to M1-11. The weak radiation field due to the relatively cool central stars in these PNe may provide favorable conditions for fullerenes to survive in the circumstellar medium. Using the photo-ionization code CLOUDY, combined with a modified blackbody, we have fitted the ~0.1-90 um spectral energy distribution and determined the dust mass in the nebula to be ~3.5x10^{-4} Msun$. Our chemical abundance analysis and SED model suggest that M1-11 is perhaps a C-rich PN with C/O ratio in the gas-phase of +0.19 dex, and that it evolved from a 1-1.5 Msun star.

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type II SN 2006bc. Observations reveal the central wavelength of the symmetric H$\alpha$ line profile to be red-shifted with respect to the host galaxy H$\alpha$ emission by day 325. Such an phenomenon has been argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of $^{56}$Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual blue-shifting of this H$\alpha$ peak which is indicative of dust formation in the ejecta. Although showing a normal peak brightness, V $\sim$ -17.2, for a core-collapse SN, 2006bc fades by $\sim$6 mag during the first 400 days suggesting either a relatively low $^{56}$Ni yield, an increase in extinction due to new dust, or both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting an optical light echo. Based on the narrow time window of this echo, we discuss implications on the location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper limit of 2 x 10$^{-3}$ M$_{\odot}$ of dust around SN 2006bc. In the event that all of this dust were formed during the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to explain the large quantities of dust observed in the early universe.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies.

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts.

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

We have estimated elemental abundances of the planetary nebula (PN) Hen2-436 in the Sagittarius (Sgr) spheroidal dwarf galaxy using ESO/VLT FORS2, Magellan/MMIRS, and Spitzer/IRS spectra. We have detected candidates of [F II] 4790A, [Kr III] 6826A, and [P II] 7875A lines and successfully estimated the abundances of these elements ([F/H]=+1.23, [Kr/H]=+0.26, [P/H]=+0.26) for the first time. We present a relation between C, F, P, and Kr abundances among PNe and C-rich stars. The detections of F and Kr support the idea that F and Kr together with C are synthesized in the same layer and brought to the surface by the third dredge-up. We have estimated the N^2+ and O^2+ abundances using optical recombination lines (ORLs) and collisionally excited lines (CELs). The discrepancy between the abundance derived from the O ORL and that derived from the O CEL is >1 dex. To investigate the status of the central star of the PN, nebula condition, and dust properties, we construct a theoretical SED model with CLOUDY. By comparing the derived luminosity and temperature of the central star with theoretical evolutionary tracks, we conclude that the initial mass of the progenitor is likely to be ~1.5-2.0 Msun and the age is ~3000 yr after the AGB phase. The observed elemental abundances can be explained by a theoretical nucleosynthesis model with a star of initial mass 2.25 Msun, Z=0.008 and LMC compositions. We have estimated the dust mass to be 2.9x10^-4 Msun (amorphous carbon only) or 4.0x10^-4 Msun (amorphous carbon and PAH). Based on the assumption that most of the observed dust is formed during the last two thermal pulses and the dust-to-gas mass ratio is 5.58x10^-3, the dust mass-loss rate and the total mass-loss rate are <3.1x10^-8 Msun / yr and <5.5x10^-6 Msun / yr, respectively. Our estimated dust mass-loss rate is comparable to a Sgr dwarf galaxy AGB star with similar metallicity and luminosity.

We have performed a comprehensive chemical abundance analysis of the extremely metal-poor ([Ar/H]<-2) halo planetary nebula (PN) BoBn 1 based on IUE archive data, Subaru/HDS spectra, VLT/UVES archive data, and Spitzer/IRS spectra. We have detected over 600 lines in total and calculated ionic and elemental abundances of 13 elements using detected optical recombination lines (ORLs) and collisionally excited lines (CELs). The estimations of C, N, O, and Ne abundances from the ORLs and Kr, Xe, and Ba from the CELs are done the first for this nebula, empirically and theoretically. The C, N, O, and Ne abundances from ORLs are systematically larger than those from CELs. The abundance discrepancies apart from O could be explained by a temperature fluctuation model, and that of O might be by a hydrogen deficient cold component model. We have detected 5 fluorine and several s-process elements. The amounts of [F/H], [Kr/H], and [Xe/H] suggest that BoBn 1 is the most F-rich among F detected PNe and is a heavy s-process element rich PN. We have confirmed dust in the nebula that is composed of amorphous carbon and PAHs with a total mass of 5.8 x 10^-6 Msun. The photo-ionization models built with non-LTE theoretical stellar atmospheres indicate that the progenitor was a 1-1.5 Msun star that would evolve into a white dwarf with an ~0.62 Msun core mass and ~0.09 Msun ionized nebula. The derived elemental abundances have been reviewed from the standpoint of theoretical nucleosynthesis models. It is likely that the elemental abundances except N could be explained either by a 1.5 Msun single star model or by a binary model composed of 0.75 Msun + 1.5 Msun stars. Careful examination implies that BoBn 1 has evolved from a 0.75 Msun + 1.5 Msun binary and experienced coalescence during the evolution to become a visible PN, similar to the other extremely metal-poor halo PN, K 648 in M 15.

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

We investigate the dust associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC) as observed with the Herschel Space Observatory. N49 is unusually bright because of an interaction with a molecular cloud along its eastern edge. We have used PACS and SPIRE to measure the far IR flux densities of the entire SNR and of a bright region on the eastern edge of the SNR where the SNR shock is encountering the molecular cloud. Using these fluxes supplemented with archival data at shorter wavelengths, we estimate the dust mass associated with N49 to be about 10 Msun. The bulk of the dust in our simple two-component model has a temperature of 20-30 K, similar to that of nearby molecular clouds. Unfortunately, as a result of the limited angular resolution of Herschel at the wavelengths sampled with SPIRE, the uncertainties are fairly large. Assuming this estimate of the dust mass associated with the SNR is approximately correct, it is probable that most of the dust in the SNR arises from regions where the shock speed is too low to produce significant X-ray emission. The total amount of warm 50-60 K dust is ~0.1 or 0.4 Msun, depending on whether the dust is modeled in terms of carbonaceous or silicate grains. This provides a firm lower limit to the amount of shock heated dust in N49.

We present development of the collimated bipolar jets from the symbiotic prototype Z And that appeared and disappeared during its 2006 outburst. In 2006 July Z And reached its historical maximum at U ~ 8.0. During this period, rapid photometric variations with Dm ~ 0.06 mag on the timescale of hours developed. Simultaneously, high-velocity satellite components appeared on both sides of the H-alpha and H-beta emission line profiles. They were launched asymmetrically with the red/blue velocity ratio of 1.2 - 1.3. From about mid-August they became symmetric. Their spectral properties indicated ejection of bipolar jets collimated within an average opening angle of 6.1 degrees. We estimated average outflow rate via jets to dM(jet)/dt ~ 2xE10-6(R(jet)/1AU)**(1/2) M(Sun)/year, during their August - September maximum, which corresponds to the emitting mass in jets, M(jet, emitting) ~ 6xE-10(Rjet)/1AU)^{3/2} M(Sun). During their lifetime, the jets released the total mass of M(jet, total) approx 7.4x1E-7 M(Sun). Evolution in the rapid photometric variability and asymmetric ejection of jets around the optical maximum can be explained by a disruption of the inner parts of the disk caused by radiation-induced warping of the disk.

We have found the fluorine lines [F IV] 3996.92A,4059.90A in the extremely metal-poor ([Ar/H] = -2.10+/-0.21) halo planetary nebula (PN) BoBn 1 in high-dispersion spectra from the 8.2-m VLT/UVES archive. Chemical abundance analysis shows that the fluorine abundance is [F/H] = +1.06+/-0.08, making BoBn 1 the most fluorine-enhanced and metal-poor PN among fluorine-detected PNe and providing new evidence that fluorine is enhanced by nucleosynthesis in low mass metal-poor stars. A comparison with the abundances of carbon-enhanced metal-poor (CEMP) stars suggests that BoBn 1 shares their origin and evolution with CEMP-s stars such as HE1305+0132. BoBn 1 might have evolved from a binary consisting of ~2 Msun primary and ~0.8 Msun secondary stars.

Aims: To investigate structure of the hot object in the symbiotic prototype Z And during its major 2000-03 active phase. Methods: Analysis of the far ultraviolet, optical low- and high-resolution spectroscopy and UBVR photometry. Reconstruction of the spectral energy distribution (SED) during the outburst. The Raman scattering process. Results: At the initial stages of the outburst the hot object was characterized by the two-temperature spectrum (a warm stellar radiation and a strong nebular emission) with signatures of a mass-outflow at moderate (100-200 km/s) and very high (1000-2000 km/s) velocities. The corresponding structure of the hot object consists of an optically thick, slowly-expanding disk-like material encompassing the accretor at the orbital plane and a fast optically thin wind over the remainder of the star. A striking similarity of [FeVII]6087 and Raman 6825 profiles at/after the dilution of the disk suggested their origin within the interaction zone where the winds from the binary components collide.