Showing 1–50 of 54 results for author: Otsuka, M

Deep JWST/NIRCam imaging of Supernova 1987A

Authors: Mikako Matsuura, M. Boyer, Richard G. Arendt, J. Larsson, C. Fransson, A. Rest, A. P. Ravi, S. Park, P. Cigan, T. Temim, E. Dwek, M. J. Barlow, P. Bouchet, G. Clayton, R. Chevalier, J. Danziger, J. De Buizer, I. De Looze, G. De Marchi, O. Fox, C. Gall, R. D. Gehrz, H. L. Gomez, R. Indebetouw, T. Kangas , et al. (24 additional authors not shown)

Abstract: JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the… ▽ More JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS; 18 pages

Seimei KOOLS-IFU mapping of the gas and dust distributions in Galactic PNe: Unveiling the origin and evolution of Galactic halo PN H4-1

Authors: Masaaki Otsuka, Toshiya Ueta, Akito Tajitsu

Abstract: H4-1 is a planetary nebula (PN) located in the Galactic halo, and is notably carbon-rich and one of the most metal-deficient PNe in the Milky Way. To unveil its progenitor evolution through the accurate measurement of the gas mass, we conducted a comprehensive investigation of H4-1, using the newly obtained Seimei/KOOLS-IFU spectra and multiwavelength spectro-photometry data. The emission line ima… ▽ More H4-1 is a planetary nebula (PN) located in the Galactic halo, and is notably carbon-rich and one of the most metal-deficient PNe in the Milky Way. To unveil its progenitor evolution through the accurate measurement of the gas mass, we conducted a comprehensive investigation of H4-1, using the newly obtained Seimei/KOOLS-IFU spectra and multiwavelength spectro-photometry data. The emission line images generated from the KOOLS-IFU datacube successfully resolve the ellipsoidal nebula and the equatorial flattened disk that are frequently seen in bipolar PNe evolved from massive progenitors. By a fully data-driven method, we directly derived the seven elemental abundances, the gas-to-dust mass ratio, and the gas and dust masses based on our own distance scale. By comparing the observed quantities with both the photoionization model and the binary nucleosynthesis model, we conclude that the progenitors of an initial mass of 1.87 Msun and 0.82 Msun are second generation stars formed ~4 Gyrs after the Big Bang, and underwent mass-transfers, binary merger, and ultimately evolved into a PN showing unique chemical abundances. Our binary model successfully reproduces the observed abundances and also explains evolutionary time scale of H4-1. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 19 pages, 5 figures, 11 tables, accepted for publication in PASJ

Probing the origin of the two-component structure of broad line region by reverberation mapping of an extremely variable quasar

Authors: Shumpei Nagoshi, Fumihide Iwamuro, Satoshi Yamada, Yoshihiro Ueda, Yuto Oikawa, Masaaki Otsuka, Keisuke Isogai, Shin Mineshige

Abstract: The physical origins of quasar components, such as the broad line region (BLR) and dust torus, remain under debate. To gain insights into them, we focused on Changing-State Quasars (CSQs) which provide a unique perspective through structural changes associated with accretion disk state transitions. We targeted SDSS J125809.31+351943.0, an extremely variable CSQ, to study its central core structure… ▽ More The physical origins of quasar components, such as the broad line region (BLR) and dust torus, remain under debate. To gain insights into them, we focused on Changing-State Quasars (CSQs) which provide a unique perspective through structural changes associated with accretion disk state transitions. We targeted SDSS J125809.31+351943.0, an extremely variable CSQ, to study its central core structure and kinematics. We conducted reverberation mapping with optical spectroscopy to explore the structure of the BLR and estimate the black hole mass. The results from H$β$ reverberation mapping indicated a black hole mass of $10^{9.64^{+0.11}_{-0.20}}\rm{M_\odot}$. Additionally, we analyzed variations in the optical to X-ray spectral indices, $α_{\rm{ox}}$, before and after the state transition, to investigate the accretion disk. These variations in $α_{\rm{ox}}$ and the Eddington ratio (from 0.4 \% to 2.4 \%) exhibitied behavior similar to state transitions observed in X-ray binary systems. Spectral analysis of H$β$ revealed a predominantly double-peaked profile during dim periods, transitioning to include a single-peaked component as the quasar brightened, suggesting that H$β$ contains a mixture of two components. Each of these components has its distinct characteristics: the first is a double-peaked profile that remains stable despite changes in the accretion rate, while the second is a variable single-peaked profile. Using time lags from reverberation mapping, we estimated the spatial relationships between these BLR components, the accretion disk, and the dust torus. Our results suggest that the BLR consists of two distinct components, each differing in location and origin. △ Less

Submitted 5 February, 2024; v1 submitted 24 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS. Comments are welcome

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Authors: Orsola De Marco, Muhammad Akashi, Stavros Akras, Javier Alcolea, Isabel Aleman, Philippe Amram, Bruce Balick, Elvire De Beck, Eric G. Blackman, Henri M. J. Boffin, Panos Boumis, Jesse Bublitz, Beatrice Bucciarelli, Valentin Bujarrabal, Jan Cami, Nicholas Chornay, You-Hua Chu, Romano L. M. Corradi, Adam Frank, Guillermo Garcia-Segura, D. A. Garcia-Hernandez, Jorge Garcia-Rojas, Veronica Gomez-Llanos, Denise R. Goncalves, Martin A. Guerrero , et al. (44 additional authors not shown)

Abstract: Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures,… ▽ More Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 32 pages, 5 figures for the main article. 12 pages 8 figures for the supplementary material

Journal ref: Nature Astronomy, 2022, Vol. 6, p. 1421

The [S II] Spectral Images of the Planetary Nebula NGC 7009: II. Major Axis

Authors: Siek Hyung, Seong-Jae Lee, Masaaki Otsuka

Abstract: We derived position-velocity density distribution diagrams along the major (PA = 77$^{\circ}$) axis of the elliptical planetary nebula NGC 7009 with the Keck HIRES [SII] 6716/6731 \A, doublet spectral images. The average densities of the main shell and knots of NGC 7009 derived from the [SII] 6716/673 \A, fluxes integrated over the line of sight indicate a density range of $N_{\rm e}$ =… ▽ More We derived position-velocity density distribution diagrams along the major (PA = 77$^{\circ}$) axis of the elliptical planetary nebula NGC 7009 with the Keck HIRES [SII] 6716/6731 \A, doublet spectral images. The average densities of the main shell and knots of NGC 7009 derived from the [SII] 6716/673 \A, fluxes integrated over the line of sight indicate a density range of $N_{\rm e}$ = $10^{3.4}$ to 10$^{3.9}$ $cm^{-3}$, while the local densities from the volume fraction resolved in radial velocities along the line of sight show a considerably large range of about 10$^{2.8}$ -- 10$^{4.7}$ $cm^{-3}$. The derived projection angle of the major axis of the main shell is about $ψ\sim$18.3($\pm$2)$^{\circ}$. Assuming that the main shell is an ellipsoidal shell with $a \simeq$16$''$ and $b \simeq 6''$, we found the range of expansion velocity, radius, and latitude of four knots and a hot bubble. The four knots at the points in symmetrical positions are roughly aligned with the same axis of expansion of latitudes $φ\sim \pm 34.5(\pm 0.6)^{\circ}$: One pair expands at about 35 $kms^{-1}$ close to the main ellipsoidal shell, and the other expands rapidly at about 60 $kms^{-1}$ at a distance of $r \sim 16''$. In the latitude range $φ= $65 -- 75$^{\circ}$, the hot bubble of a relatively large structure expands rapidly with a velocity of 130 -- 150 $kms^{-1}$. Four knots and hot bubble points that expand faster than the main shell appear to have been formed by two to three eruptions at a different epoch than the primary structure formation. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 18 pages, 7 figures Accepted for publication in ApJS

Time-Variable Jet Ejections from RW Aur A, RY Tau and DG Tau

Authors: Michihiro Takami, Hans Moritz Guenther, P. Christian Schneider, Tracy L. Beck, Jennifer L. Karr, Youichi Ohyama, Roberto Galvan-Madrid, Taichi Uyama, Marc White, Konstantin Grankin, Deirdre Coffey, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activ… ▽ More We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activities of mass accretion and/or the stellar magnetic fields. All of these jets consist of several moving knots with tangential velocities of 70-240 km s-1, ejected from the star with different velocities and at irregular time intervals. Via comparison with literature, we identify significant differences in tangential velocities for the DG Tau jet between 1985-2008 and 2008-2021. The sizes of the individual knots appear to increase with time, and in turn, their peak brightnesses in the 1.644-micron emission decreased up to a factor of ~30 during the epochs of our observations. A variety of the decay timescales measured in the [Fe II] 1.644 micron emission can be attributed to different pre-shock conditions if the moving knots are unresolved shocks. However, our data do not exclude the possibility that these knots are due to non-uniform density/temperature distributions with another heating mechanism, or in some cases due to stationary shocks without proper motions. Spatially resolved observations of these knots with significantly higher angular resolutions are necessary to better understand their physical nature. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 24 pages, 11 figures, 8 tables, accepted for ApJS

On the N/O abundance ratio and the progenitor mass for the most luminous planetary nebulae of M 31

Authors: Toshiya Ueta, Masaaki Otsuka

Abstract: Plasma diagnostics are the bases of investigation into the physical and chemical properties of line-emitting gaseous systems. To perform plasma diagnostics properly, it is essential to correct the input spectrum for extinction properly. This is simply because determining the degree of extinction is dependent on the physical properties of the line-emitting gas. Hence, both extinction correction and… ▽ More Plasma diagnostics are the bases of investigation into the physical and chemical properties of line-emitting gaseous systems. To perform plasma diagnostics properly, it is essential to correct the input spectrum for extinction properly. This is simply because determining the degree of extinction is dependent on the physical properties of the line-emitting gas. Hence, both extinction correction and plasma diagnostics have to be performed simultaneously and self-consistently. By comparing the results of analyses performed for a sample of nine bright planetary nebulae in M 31 with and without the proper extinction correction and plasma diagnostics, we demonstrate how initial assumptions for the physical conditions of the line-emitting gas in extinction correction would compromise the results of the entire analyses. While the electron density/temperature are relatively immune to the imposed inconsistent assumptions, the compromised extinction would cause systematic offsets in the extinction-corrected line strengths, which consequently would impose adverse effects on the resulting ionic and elemental abundances, and other inferences made from the incorrect results. We find that this M 31 PN sample simply represents those around the high-mass end of the mass range for low-mass planetary nebula progenitor stars as expected from the existing theoretical models. It appears that the suspicion raised in the previous study - these PNe being anomalously nitrogen overabundant - is simply caused by the apparent underestimate in extinction that originates from the imposed inconsistent assumptions in extinction correction. In a larger context, the results of plasma diagnostics in the literature without seeking simultaneous self-consistency with extinction correction have to be handled cautiously. Ideally, such previous results should be re-evaluated by seeking simultaneous self-consistency. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 5 pages (main text with 2 tables) plus 5 supplemental pages (5 tables in Appendix); to be published in Astronomy and Astrophysics Letter to the Editor

Journal ref: A&A 667, L8 (2022)

Double Shells of the Planetary Nebula NGC 7009 Minor Axis

Authors: Seong-Jae Lee, Siek Hyung, Masaaki Otsuka

Abstract: We analyzed the minor axis spectra of the elliptical planetary nebula (PN) NGC 7009 observed with the Keck HIRES with a 0.862$"$ $\times$ 10$"$ slit placed at about $\sim$7.5$"$ and 10$"$ away from the center and a 0.862$"$ $\times$ 14$"$ slit at the center. The mean densities derived from the integrated [SII] 6716/6731Å, fluxes along the Keck HIRES slit length indicate a density range of 10… ▽ More We analyzed the minor axis spectra of the elliptical planetary nebula (PN) NGC 7009 observed with the Keck HIRES with a 0.862$"$ $\times$ 10$"$ slit placed at about $\sim$7.5$"$ and 10$"$ away from the center and a 0.862$"$ $\times$ 14$"$ slit at the center. The mean densities derived from the integrated [SII] 6716/6731Å, fluxes along the Keck HIRES slit length indicate a density range of 10$^{3.7}$ to 10$^{4.1}$ cm$^{-3}$, while the local densities derived from the slit spectral images show a large local density variation of about 10$^{2.8}$ - 10$^{4.6}$ cm$^{-3}$: local densities vary substantially more than values integrated over the line of sight. The expansion rates of the main and outer shells obtained by [SII] are about 21.7 and 30.0 kms$^{-1}$, respectively. The kinematic results of the [SII] spectral lines correspond to the outermost regions of the two shells and are not representative of the whole PN but are closely related to the other emission lines observed in the shell gas. We conclude that the density contrast leads to the formation of the inner shell, while the change in ionization state leads to the formation of the outer shell. We suggest that the inner main and outer shells result from two successive major ejections. The physical conditions of the central star must have been different when these shells first formed. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 2022,

CO line observations of OH/IR stars in the inner Galactic Bulge: Characteristics of stars at the tip of the AGB

Authors: H. Olofsson, T. Khouri, B. A. Sargent, A. Winnberg, J. A. D. L. Blommaert, M. A. T. Groenewegen, S. Muller, J. H. Kastner, M. Meixner, M. Otsuka, N. Patel, N. Ryde, S. Srinivasan

Abstract: 12CO and 13CO lines, as well as a mm-wave continuum, have been observed for a sample of 22 OH/IR stars in directions within 2 degrees of the Galactic Centre. Photometry data have been gathered from the literature to construct SEDs and to determine pulsational variability. Radiative transfer models have been used to interpret the data. All stars in the sample were detected in at least one CO line,… ▽ More 12CO and 13CO lines, as well as a mm-wave continuum, have been observed for a sample of 22 OH/IR stars in directions within 2 degrees of the Galactic Centre. Photometry data have been gathered from the literature to construct SEDs and to determine pulsational variability. Radiative transfer models have been used to interpret the data. All stars in the sample were detected in at least one CO line, and 8 objects were detected in 324 GHz continuum. Based on luminosity criteria, the sample is divided into 17 objects that most likely lie within the inner Galactic Bulge, and 5 objects that are most likely foreground objects. The median luminosity of the inner-Galactic-Bulge objects, 5600 Lsun, corresponds to an initial mass in the range 1.2-1.6 Msun, indicating that these OH/IR stars descend from solar-type stars. The objects in this sub-sample are further divided into two classes based on their SED characteristics: 11 objects have SEDs that are well matched by models invoking dust envelopes extending from a few stellar radii and outwards, while 6 objects are better modelled as having detached dust envelopes with inner radii in the range 200-600 au and warmer central stars. The former objects have periodic variability, while the latter objects are predominantly non-periodic. The median gas-mass-loss rate, gas terminal expansion velocity, gas-to-dust mass ratio, and circumstellar 12CO/13CO abundance ratio have been estimated to be 2x10{-5} Msun/yr, 18 km/s, 200 (excluding the sources with detached dust envelopes, which show markedly lower gas-to-dust ratios), and 5, respectively, for the inner-Galactic-Bulge objects. The inner-Galactic-Bulge OH/IR stars studied here constitute an excellent sample of equidistant objects for the purpose of understanding the evolution of the mass-loss-rate characteristics at the tip of the AGB. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: 32 pages, 17 figures

Seimei KOOLS-IFU Mapping of the Gas and Dust Distributions in Galactic Planetary Nebulae: the Case of IC2165

Authors: Masaaki Otsuka

Abstract: We investigated the physical and chemical properties of the gas and dust components in a carbon-rich planetary nebula (PN) IC2165 using two-dimensional (2-D) emission-line maps with superior resolution. The extinction map is generated in a self-consistent and assumption-free manner. The circumstellar gas-to-dust mass ratio (GDR) map ranges radially from 1210 in the central nebula filled with hot g… ▽ More We investigated the physical and chemical properties of the gas and dust components in a carbon-rich planetary nebula (PN) IC2165 using two-dimensional (2-D) emission-line maps with superior resolution. The extinction map is generated in a self-consistent and assumption-free manner. The circumstellar gas-to-dust mass ratio (GDR) map ranges radially from 1210 in the central nebula filled with hot gas plasma to 120 near the ionisation front. The determined GDR is comparable to ~400, which is commonly adopted for carbon-rich asymptotic giant branch (AGB) stars, and ~100 for ISM. Except for the inner regions, the GDR in IC2165 is nearly the same as in such AGB stars, indicating that most dust grains withstand the harsh radiation field without being destroyed. The gas and dust mass distributions concentrated in the equatorial plane may be related to the nonisotropic mass loss during the AGB phase and nebula shaping. The spatial distributions of electron densities/temperatures and ionic/elemental abundances were investigated herein. We determined 13 elemental abundances using PSF-matched spatially integrated multiwavelength spectra extracted from the same aperture. Their values are consistent with values predicted by a theoretical model for stars of initially 1.75 Msun and Z = 0.003. Finally, we constructed the photoionisation model using our distance measurement to be consistent with all derived quantities, including the GDR and gas and dust masses and post-AGB evolution. Thus, we demonstrate the capability of Seimei/KOOLS-IFU and how the spatial variation of the gas and dust components in PNe derived from IFU observations can help understand the evolution of the circumstellar/interstellar medium. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 25 pages, 16 figures, 19 tables, accepted for publication in MNRAS

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

Authors: P. Scicluna, F. Kemper, I. McDonald, S. Srinivasan, A. Trejo, S. H. J. Wallström, J. G. A. Wouterloot, J. Cami, J. Greaves, Jinhua He, D. T. Hoai, Hyosun Kim, O. C. Jones, H. Shinnaga, C. J. R. Clark, T. Dharmawardena, W. Holland, H. Imai, J. Th. van Loon, K. M. Menten, R. Wesson, H. Chawner, S. Feng, S. Goldman, F. C. Liu , et al. (67 additional authors not shown)

Abstract: The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a n… ▽ More The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

Proper Plasma Analysis Practice (PPAP), an Integrated Procedure of the Extinction Correction and Plasma Diagnostics: a Demo with an HST/WFC3 Image Set of NGC6720

Authors: Toshiya Ueta, Masaaki Otsuka

Abstract: In this work, we propose a proper plasma analysis practice (PPAP), an updated procedure of plasma diagnostics in the era of spatially-resolved spectroscopy. In particular, we emphasize the importance of performing both of the extinction correction and the direct method of plasma diagnostics simultaneously as an integrated process. This approach is motivated by the reciprocal dependence between cri… ▽ More In this work, we propose a proper plasma analysis practice (PPAP), an updated procedure of plasma diagnostics in the era of spatially-resolved spectroscopy. In particular, we emphasize the importance of performing both of the extinction correction and the direct method of plasma diagnostics simultaneously as an integrated process. This approach is motivated by the reciprocal dependence between critical parameters in these analyses, which can be resolved by iteratively seeking a converged solution. The use of PPAP allows us to eliminate unnecessary assumptions that prevent us from obtaining an exact solution at each element of the spectral imaging data. Using a suite of HST/WFC3 narrowband images of the planetary nebula, NGC 6720, we validate PPAP by (1) simultaneously and self-consistently deriving the extinction, c(Hb), and electron density/temperature distribution, (n_e, T_e), maps that are consistent with each other, and (2) obtaining identical metal abundance distribution maps, (n(N^+)/n(H^+), n(S^+)/n(H^+)), from multiple emission line maps at different wavelengths/transition energies. We also determine that the derived c(Hb) consists both of the ISM and circumsource components and that the ionized gas-to-dust mass ratio in the main ring is at least 437 and as high as about 1600. We find that, unless we deliberately seek self-consistency, uncertainties at tens of per cent can easily arise in outcomes, making it impossible to discern actual spatial variations that occurs at the same level, defeating the purpose of conducting spatially resolved spectroscopic observations. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: 27 pages, 18 figures, to be published as a Tutorial in the Publication of the Astronomical Society of the Pacific

Magnetism of Fullerene C60 Compared with Graphene Molecule by DFT Calculation, Laboratory Experiment and Astronomical Observation

Authors: Norio Ota, Aigen Li, Laszlo Nemes, Masaaki Otsuka

Abstract: Magnetism of fullerene C60 was studied by three methods of the density functional theory (DFT) calculation, laboratory experiment and astronomical observation. DFT revealed that the most stable spin state was non-magnetic one of Sz=0/2. This is contrary to our recent study on void induced graphene molecules of C23 and C53 to be magnetic one of Sz=2/2. Two graphene molecules combined model suggeste… ▽ More Magnetism of fullerene C60 was studied by three methods of the density functional theory (DFT) calculation, laboratory experiment and astronomical observation. DFT revealed that the most stable spin state was non-magnetic one of Sz=0/2. This is contrary to our recent study on void induced graphene molecules of C23 and C53 to be magnetic one of Sz=2/2. Two graphene molecules combined model suggested that two up-spin at every carbon pentagon ring may cancel each other to bring Sz=0/2. Similar cancelation may occur on C60. Molecular vibrational infrared spectrum of C60 show four major bands, which coincide with gas-phase laboratory experiment, also with astronomically observed one of carbon rich planetary nebula Tc1 and Lin49. However, there remain many unidentified bands on astronomical one. We supposed multiple voids on graphene sheet, which may create both C60 and complex graphene molecules. It was revealed that spectrum of two voids induced graphene molecule coincident well with major astronomical bands. Simple sum of C60 and graphene molecules could successfully reproduce astronomical bands in detail. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 10 pages, 11 figures, 2 tables. arXiv admin note: text overlap with arXiv:2104.09745

Journal ref: Journal of the Magnetics Society of Japan, Volume 45 Issue 6 Pages 142-148 (2021)

Intermediate Luminosity Type Iax SN 2019muj With Narrow Absorption Lines: Long-Lasting Radiation Associated With a Possible Bound Remnant Predicted by the Weak Deflagration Model

Authors: Miho Kawabata, Keiichi Maeda, Masayuki Yamanaka, Tatsuya Nakaoka, Koji S. Kawabata, Kentaro Aoki, G. C. Anupama, Umut Burgaz, Anirban Dutta, Keisuke Isogai, Masaru Kino, Naoto Kojiguchi, Iida Kota, Brajesh Kumar, Daisuke Kuroda, Hiroyuki Maehara, Kazuya Matsubayashi, Kumiko Morihana, Katsuhiro L. Murata, Tomohito Ohshima, Masaaki Otsuka, Devendra K. Sahu, Avinash Singh, Koji Sugitani, Jun Takahashi , et al. (1 additional authors not shown)

Abstract: We present comprehensive spectroscopic and photometric analyses of the intermediate luminosity Type Iax supernova (SN Iax) 2019muj based on multi-band datasets observed through the framework of the OISTER target-of-opportunity program. SN 2019muj exhibits almost identical characteristics with the subluminous SNe Iax 2008ha and 2010ae in terms of the observed spectral features and the light curve e… ▽ More We present comprehensive spectroscopic and photometric analyses of the intermediate luminosity Type Iax supernova (SN Iax) 2019muj based on multi-band datasets observed through the framework of the OISTER target-of-opportunity program. SN 2019muj exhibits almost identical characteristics with the subluminous SNe Iax 2008ha and 2010ae in terms of the observed spectral features and the light curve evolution at the early phase, except for the peak luminosity. The long-term observations unveil the flattening light curves at the late time as seen in a luminous SN Iax 2014dt. This can be explained by the existence of an inner dense and optically-thick component possibly associated with a bound white dwarf remnant left behind the explosion. We demonstrate that the weak deflagration model with a wide range of the explosion parameters can reproduce the late-phase light curves of other SNe Iax. Therefore, we conclude that a common explosion mechanism operates for different subclass SNe Iax. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: 37 pages, 22 figures, 6 tables, accepted for publication in PASJ

Void defect induced magnetism and structure change of carbon material-2, Graphene molecules

Authors: Norio Ota, Aigen Li, Laszlo Nemes, Masaaki Otsuka

Abstract: Void-defect is a possible origin of ferromagnetic feature on pure carbon materials. In our previous paper, void-defect on graphene-nanoribbon show highly polarized spin configuration. In this paper, we studied cases for graphene molecules by quantum theory, by astronomical observation and by laboratory experiment. Model molecules for the density functional theory are graphene molecules of C23 and… ▽ More Void-defect is a possible origin of ferromagnetic feature on pure carbon materials. In our previous paper, void-defect on graphene-nanoribbon show highly polarized spin configuration. In this paper, we studied cases for graphene molecules by quantum theory, by astronomical observation and by laboratory experiment. Model molecules for the density functional theory are graphene molecules of C23 and C53 induced by a void-defect. They have carbon pentagon ring within a hexagon network. Single void has three radical carbons, holding six spins. Those spins make several spin-states, which affects to molecular structure and molecular vibration, finally to infrared spectrum. The stable spin state was triplet, not singlet. This suggests magnetic pure carbon molecule. It was a surprise that those molecules show close infrared spectrum with astronomically observed one, especially observed on carbon rich planetary nebulae. We could assign major band at 18.9 micrometer, and sub-bands at 6.6, 7.0, 7.6, 8.1, 8.5, 9.0 and 17.4 micrometer. Also, calculated spectrum roughly coincides with that of laboratory experiment by the laser-induced carbon plasma, which is an analogy of cosmic carbon creation in interstellar space. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 9 pages, 6 figures, 1 table. arXiv admin note: text overlap with arXiv:2007.03862

Journal ref: J. Magn. Soc. Jpn. 45, 41 (2021)

Spectroscopic and Photometric Observations of Dwarf Nova Superoutbursts by the 3.8 m Telescope Seimei and the Variable Star Network

Authors: Yusuke Tampo, Keisuke Isogai, Naoto Kojiguchi, Hiroyuki Maehara, Kenta Taguchi, Taichi Kato, Mariko Kimura, Yasuyuki Wakamatsu, Masaaki Shibata, Daisaku Nogami, Miho Kawabata, Keiichi Maeda, Kosuke Namekata, Soshi Okamoto, Masaaki Otsuka, Burgaz Umut, Shumpei Nagoshi, Hiroshi Itoh, Tonny Vanmunster, Tamas Tordai, Geoffrey Stone, Katsura Matsumoto, Daiti Fujii, Ryota Matsumura, Momoka Nakagawa , et al. (37 additional authors not shown)

Abstract: We present spectroscopic and photometric observations of 17 dwarf-nova superoutbursts obtained by KOOLS-IFU mounted on the 3.8 m telescope Seimei at Okayama Observatory of Kyoto University and through VSNET collaboration. Our spectroscopic observations for six outbursts were performed within 1 d from their optical peak. 11 objects (TCP J00590972+3438357. ASASSN-19ado, TCP J06073081-0101501, ZTF20a… ▽ More We present spectroscopic and photometric observations of 17 dwarf-nova superoutbursts obtained by KOOLS-IFU mounted on the 3.8 m telescope Seimei at Okayama Observatory of Kyoto University and through VSNET collaboration. Our spectroscopic observations for six outbursts were performed within 1 d from their optical peak. 11 objects (TCP J00590972+3438357. ASASSN-19ado, TCP J06073081-0101501, ZTF20aavnpug, ASASSN-19ady, MASTER OT J061642.05+435617.9, TCP J20034647+1335125, ASASSN-20kv, ASASSN-20kw, MASTER OT J213908.79+161240.2, and ASASSN-20mf) were previously unknown systems, and our observations enabled quick classification of their transient type. These results illustrate that Seimei telescope has the capability to conduct quick follow-up observations of unknown transients. Our photometric observations yielded that 11 objects are WZ Sge-type dwarf novae and their candidates, and the other six objects are SU UMa-type dwarf novae and their candidates. The He II 4686Å~ emission line was clearly detected among ASASSN-19ado, TCP J06073081-0101501 and MASTER OT J213908.79+161240.2, whose association with a spiral arm structure in an accretion disk has been suggested in the previous studies. Our result suggests that a higher-inclination system shows a stronger emission line of He II 4686Å, as well as larger-amplitude early superhumps. △ Less

Submitted 11 April, 2021; originally announced April 2021.

Comments: 21 pages, 14 figures, 3 tables, accepted for publication in PASJ

K-band High-Resolution Spectroscopy of Embedded Massive Protostars

Authors: Tien-Hao Hsieh, Michihiro Takami, Michael S. Connelley, Sheng-Yuan Liu, Yu-Nung Su, Naomi Hirano, Motohide Tamura, Masaaki Otsuka, Jennifer L. Karr, Tae-Soo Pyo

Abstract: A classical paradox in high-mass star formation is that powerful radiation pressure can halt accretion, preventing further growth of a central star. Disk accretion has been proposed to solve this problem, but the disks and the accretion process in high-mass star formation are poorly understood. We executed high-resolution ($R$=35,000-70,000) iSHELL spectroscopy in $K$-band for eleven high-mass pro… ▽ More A classical paradox in high-mass star formation is that powerful radiation pressure can halt accretion, preventing further growth of a central star. Disk accretion has been proposed to solve this problem, but the disks and the accretion process in high-mass star formation are poorly understood. We executed high-resolution ($R$=35,000-70,000) iSHELL spectroscopy in $K$-band for eleven high-mass protostars. Br-$γ$ emission was observed toward eight sources, and the line profiles for most of these sources are similar to those of low-mass PMS stars. Using an empirical relationship between the Br-$γ$ and accretion luminosities, we tentatively estimate disk accretion rates ranging from $\lesssim$10$^{-8}$ and $\sim$10$^{-4}$ $M_\odot$ yr$^{-1}$. These low-mass-accretion rates suggest that high-mass protostars gain more mass via episodic accretion as proposed for low-mass protostars. Given the detection limits, CO overtone emission ($v$=2-0 and 3-1), likely associated with the inner disk region ($r \ll 100$ au), was found towards two sources. This low-detection rate compared with Br-$γ$ emission is consistent with previous observations. Ten out of the eleven sources show absorption at the $v$=0-2 ${\rm R(7)-R(14)}$ CO R-branch. Most of them are either blueshifted or redshifted, indicating that the absorption is associated with an outflow or an inflow with a velocity of up to $\sim50$ km s$^{-1}$. Our analysis indicates that the absorption layer is well thermalized (and therefore $n_{\mathrm H_2} \gtrsim 10^6$ cm$^{-3}$) at a single temperature of typically 100-200 K, and located within 200-600 au of the star. △ Less

Submitted 14 March, 2021; originally announced March 2021.

Comments: 18 page, 7 figures, accepted to ApJ

Possible Time Correlation Between Jet Ejection and Mass Accretion for RW Aur A

Authors: Michihiro Takami, Tracy L. Beck, P. Christian Schneider, Hans Moritz Guenther, Marc White, Konstantin Grankin, Jennifer L. Karr, Youichi Ohyama, Deirdre Coffey, Hauyu Baobab Liu, Roberto Galvan-Madrid, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities be… ▽ More For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities between (A) the Ca II 8542 A and O I 7772 A line profiles associated with magnetospheric accretion, and (B) optical continuum fluxes. The jet images and their proper motions show that four knot ejections occurred at the star over the past ~15 years with an irregular interval of 2-6 years. The time scale and irregularity of these intervals are similar to those of the dimming events seen in the optical photometry data. Our observations show a possible link between remarkable (Delta_V < -1 mag.) photometric rises and jet knot ejections. Observations over another few years may confirm or reject this trend. If confirmed, this would imply that the location of the jet launching region is very close to the star (r <<0.1 au) as predicted by some jet launching models. Such a conclusion would be crucial for understanding disk evolution within a few au of the star, and therefore possible ongoing planet formation at these radii. △ Less

Submitted 31 July, 2020; originally announced July 2020.

Comments: 14 pages, 5 figures, Accepted by Astrophysical Journal

Graphene Molecules Contributing to the Infrared Bands of Carbon Rich Planetary Nebulae

Authors: Norio Ota, Aigen Li, Laszlo Nemes, Masaaki Otsuka

Abstract: 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 f… ▽ More 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 astronomical assumption, that is, single void in 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 are classified as different species. Single void is surrounded by 3 radical carbons, holding 6 spins. Spin state affects molecular configuration and vibrational spectrum. It was a surprise that the triplet state is stable than the singlet. Most of charge neutral and triplet spin state species show closely resembling spectra with observed one of carbon rich planetary nebulae 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 plasma, 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. △ Less

Submitted 30 July, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: 21 pages, 8 figures, 1 table, 1 appendix

Calcium-rich Transient SN 2019ehk in A Star-Forming Environment: Yet Another Candidate for An Ultra-Stripped Envelope Supernova

Authors: Tatsuya Nakaoka, Keiichi Maeda, Masayuki Yamanaka, Masaomi Tanaka, Miho Kawabata, Takashi J. Moriya, Koji S. Kawabata, Nozomu Tominaga, Kengo Takagi, Fumiya Imazato, Tomoki Morokuma, Shigeyuki Sako, Ryou Ohsawa, Takashi Nagao, Ji-an Jiang, Umut Burgaz, Kenta Taguchi, Makoto Uemura, Hiroshi Akitaya, Mahito Sasada, Keisuke Isogai, Masaaki Otsuka, Hiroyuki Maehara

Abstract: 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,… ▽ More 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. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: submitted to ApJ on Mar 27, 2020

Physical properties of the fluorine and neutron-capture element-rich PN Jonckheere 900

Authors: Masaaki Otsuka, Siek Hyung

Abstract: 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 exhib… ▽ More 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. △ Less

Submitted 23 December, 2019; v1 submitted 4 November, 2019; originally announced November 2019.

Comments: 22 pages, 11 figures, accepted for publication by MNRAS. The title is corrected for accordance with the paper published by MNRAS

Type Ia SN 2019ein: New Insights into the Similarities and diversities among High-Velocity SNe Ia

Authors: Miho Kawabata, Keiichi Maeda, Masayuki Yamanaka, Tatsuya Nakaoka, Koji Kawabata, Ryo Adachi, Hiroshi Akitaya, Umut Burgaz, Hidekazu Hanayama, Takashi Horiuchi, Ryohei Hosokawa, Kota Iida, Fumiya Imazato, Keisuke Isogai, Ji-an Jiang, Noriyuki Katoh, Hiroki Kimura, Masaru Kino, Daisuke Kuroda, Hiroyuki Maehara, Kazuya Matsubayashi, Kumiko Morihana, Katsuhiro Murata, Takashi Nagao, Masafumi Niwano , et al. (16 additional authors not shown)

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 16 March, 2020; v1 submitted 8 August, 2019; originally announced August 2019.

Comments: 15 pages, 12 figres, Accepted to ApJ

An Imaging Spectroscopic Survey of the Planetary Nebula NGC 7009 with MUSE

Authors: J. R. Walsh, A. Monreal-Ibero, M. J. Barlow, T. Ueta, R. Wesson, A. A. Zijlstra, S. Kimeswenger, M. L. Leal-Ferreira, M. Otsuka

Abstract: 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 ove… ▽ More 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) △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: Comments: 27 pages, 25 figures, 4 tables. Accepted for publication in A&A

Journal ref: A&A 620, A169 (2018)

Physical properties of the fullerene C60-containing planetary nebula SaSt 2-3

Authors: Masaaki Otsuka

Abstract: 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 va… ▽ More 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. △ Less

Submitted 5 October, 2018; v1 submitted 3 October, 2018; originally announced October 2018.

Comments: 21 pages, 13 figures, 16 tables, Accepted for publication by MNRAS

Morpho-Kinematics of the Circumstellar Environments around Post-AGB Stars

Authors: Toshiya Ueta, Hideyuki Izumiura, Issei Yamamura, Masaaki Otsuka

Abstract: 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. 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. △ Less

Submitted 15 September, 2018; originally announced September 2018.

Comments: 2 pages, 2 figures; contribution to Proceedings of IAU Symp 343

Herschel Planetary Nebula Survey (HerPlaNS): Hydrogen Recombination Laser Lines in Mz 3

Authors: Isabel Aleman, Katrina Exter, Toshiya Ueta, Samuel Walton, A. G. G. M. Tielens, Albert Zijlstra, Rodolfo Montez Jr., Zulema Abraham, Masaaki Otsuka, Pedro P. B. Beaklini, Peter A. M. van Hoof, Eva Villaver, Marcelo L. Leal-Ferreira, Edgar Mendoza, Jacques D. R. Lepine

Abstract: 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 spe… ▽ More 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 $μ$m range covered by the PACS and SPIRE instruments (H11$α$ to H21$α$ and H14$β$). 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 $α$ 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. △ Less

Submitted 16 May, 2018; originally announced May 2018.

Comments: Published in MNRAS. ESA Portal Story. 13 Pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018 sty966, https://doi.org/10.1093/mnras/sty966

The Herschel Planetary Nebula Survey (HerPlaNS) - a comprehensive dusty photoionization model of NGC6781

Authors: Masaaki Otsuka, Toshiya Ueta, Peter A. M. van Hoof, Raghvendra Sahai, Isabel Aleman, Albert A. Zijlstra, You-Hua Chu, Eva Villaver, Marcelo L. Leal-Ferreira, Joel Kastner, Ryszard Szczerba, Katrina M. Exter

Abstract: 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 gian… ▽ More 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. △ Less

Submitted 20 July, 2017; originally announced July 2017.

Comments: 34 pages, 13 Figures and 16 Tables, accepted for publication in the Astrophysical Journal Supplement Series

Physical properties of the very young PN Hen3-1357 (Stingray Nebula) based on multiwavelength observations

Authors: Masaaki Otsuka, M. Parthasarathy, A. Tajitsu, S. Hubrig

Abstract: 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 supp… ▽ More 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. △ Less

Submitted 6 March, 2017; originally announced March 2017.

Comments: 24 pages, 7 figures, 12 tables, Accepted for publication by the Astrophysical Journal

XSHOOTER spectroscopy of the enigmatic planetary nebula Lin49 in the Small Magellanic Cloud

Authors: Masaaki Otsuka, F. Kemper, M. L. Leal-Ferreira, I. Aleman, J. Bernard-Salas, J. Cami, B. Ochsendorf, E. Peeters, P. Scicluna

Abstract: 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… ▽ More 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. △ Less

Submitted 4 July, 2016; originally announced July 2016.

Comments: 24 pages, 17 figures, 17 tables, Accepted for publication by MNRAS

First-generation Science Cases for Ground-based Terahertz Telescopes

Authors: Hiroyuki Hirashita, Patrick M. Koch, Satoki Matsushita, Shigehisa Takakuwa, Masanori Nakamura, Keiichi Asada, Hauyu Baobab Liu, Yuji Urata, Ming-Jye Wang, Wei-Hao Wang, Satoko Takahashi, Ya-Wen Tang, Hsian-Hong Chang, Kuiyun Huang, Oscar Morata, Masaaki Otsuka, Kai-Yang Lin, An-Li Tsai, Yen-Ting Lin, Sundar Srinivasan, Pierre Martin-Cocher, Hung-Yi Pu, Francisca Kemper, Nimesh Patel, Paul Grimes , et al. (11 additional authors not shown)

Abstract: 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 the… ▽ More 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. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

Chemical Abundances in the PN Wray16-423 in the Sagittarius Dwarf Spheroidal Galaxy: Constraining the Dust Composition

Authors: Masaaki Otsuka

Abstract: 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 ext… ▽ More 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. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 25 pages, 17 figures, 16 tables, Accepted for publication by MNRAS

Chemical Abundances and Dust in the Halo Planetary Nebula K648 in M15: Its Origin and Evolution based on an Analysis of Multiwavelength Data

Authors: Masaaki Otsuka, Siek Hyung, Akito Tajitsu

Abstract: 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 dete… ▽ More 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. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 23 pages, 14 Figures and 20 Tables, accepted for publication in the Astrophysical Journal Supplement Series

CO J = 2 - 1 Emission from Evolved Stars in the Galactic Bulge

Authors: Benjamin A. Sargent, N. A. Patel, M. Meixner, M. Otsuka, D. Riebel, S. Srinivasan

Abstract: 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] A… ▽ More 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. △ Less

Submitted 29 July, 2014; originally announced July 2014.

Comments: 40 pages, 16 figures, appeared in the 1 March 2013 issue of the Astrophysical Journal

MSC Class: 85A04

Journal ref: Ap. J. 765 (2013) 20

Molecular hydrogen emission in the interstellar medium of the Large Magellanic Cloud

Authors: Naslim Neelamkodan, Francisca Kemper, Suzanne Madden, Sacha Hony, You-Hua Chu, Frederic Galliano, Caroline Bot, Yao-Lun Yang, Ji Yeon Seok, Joana M. Oliveira, Jacco Th. van Loon, Margaret Meixner, Aigen Li, Annie Hughes, Karl D. Gordon, Masaaki Otsuka, Hiroyuki Hirashita, Oscar Morata, Vianney Lebouteiller, Remy Indebetouw, Sundar Srinivasan, Jean-philippe Bernard, William T. Reach

Abstract: 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 μm}$ are detected in the IRS spectra for ten regions. The higher level tran… ▽ More 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 μ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$σ$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm μ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. △ Less

Submitted 3 November, 2014; v1 submitted 29 July, 2014; originally announced July 2014.

Comments: Accepted for publication in MNRAS

Interstellar and circumstellar fullerenes

Authors: J. Bernard-Salas, J. Cami, A. P. Jones, E. Peeters, E. R. Micelotta, M. Otsuka, G. C. Sloan, F. Kemper, M. Groenewegen

Abstract: 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 rece… ▽ More 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. △ Less

Submitted 3 July, 2014; originally announced July 2014.

Comments: 6 pages, 2 figures. Proceedings of "The Life Cycle of Dust in the Universe: Observations, Theory, and laboratory Experiments". Taipei, Taiwan 18-22 November 2013. Editors: A. Andersen, M. Baes, H. Gomez, C. Kemper, D. Watson

Journal ref: Proceedings of Science (LDCU 2013)032

The Herschel Planetary Nebula Survey (HerPlaNS) I. Data Overview and Analysis Demonstration with NGC 6781

Authors: T. Ueta, D. Ladjal, K. M. Exter, M. Otsuka, R. Szczerba, N. Siódmiak, I. Aleman, P. A. M. van Hoof, J. H. Kastner, R. Montez, I. McDonald, M. Wittkowski, C. Sandin, S. Ramstedt, O. De Marco, E. Villaver, Y. -H. Chu, W. Vlemmings, H. Izumiura, R. Sahai, J. A. Lopez, B. Balick, A. Zijlstra, A. G. G. M. Tielens, R. E. Rattray , et al. (8 additional authors not shown)

Abstract: 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, dat… ▽ More 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. △ Less

Submitted 11 March, 2014; originally announced March 2014.

Comments: 27 pages, 16 figures, accepted for publication in A&A

Spitzer Space Telescope spectra of post-AGB stars in the Large Magellanic Cloud ---polycyclic aromatic hydrocarbons at low metallicities

Authors: Mikako Matsuura, Jeronimo Bernard-Salas, T. Lloyd Evans, Kevin M. Volk, Bruce J. Hrivnak, G. C. Sloan, You-Hua Chu, Robert Gruendl, Kathleen E. Kraemer, Els Peeters, R. Szczerba, P. R. Wood, Albert A. Zijlstra, S. Hony, Yoshifusa Ita, Devika Kamath, Eric Lagadec, Quentin A Parker, Warren A. Reid, Takashi Shimonishi, H. Van Winckel, Paul M. Woods, F. Kemper, Margaret Meixner, M. Otsuka , et al. (4 additional authors not shown)

Abstract: 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 pap… ▽ More 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. △ Less

Submitted 3 January, 2014; originally announced January 2014.

Comments: Accepted by MNRAS

Physical Properties of Fullerene-containing Galactic Planetary Nebulae

Authors: Masaaki Otsuka, F. Kemper, J. Cami, E. Peeters, J. Bernard-Salas

Abstract: 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… ▽ More 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. △ Less

Submitted 29 October, 2013; originally announced October 2013.

Comments: 18 pages, 13 figures, 13 tables, Accepted for publication by MNRAS

Chemical abundances in the extremely Carbon and Xenon-rich halo planetary nebula H4-1

Authors: Masaaki Otsuka, Akito Tajitsu

Abstract: 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 d… ▽ More 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. △ Less

Submitted 3 October, 2013; originally announced October 2013.

Comments: Accepted for publication by the Astrophysical Journal

The detection of C60 in the well-characterized planetary nebula M1-11

Authors: Masaaki Otsuka, F. Kemper, S. Hyung, B. A. Sargent, M. Meixner, A. Tajitsu, K. Yanagisawa

Abstract: 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 amo… ▽ More 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. △ Less

Submitted 29 January, 2013; originally announced January 2013.

Comments: Accepted for publication by the Astrophysical Journal

Optical and Infrared Analysis of Type II SN 2006BC

Authors: Joseph S. Gallagher, B. E. K. Sugerman, Geoffrey C. Clayton, J. E. Andrews, J. Clem, M. J. Barlow, B. Ercolano, J. Fabbri, M. Otsuka, R. Wesson, M. Meixner

Abstract: 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$α$ line profile to be red-shifted with respect to the host galaxy H$α$ 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}$N… ▽ More 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$α$ line profile to be red-shifted with respect to the host galaxy H$α$ 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$α$ 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. △ Less

Submitted 24 May, 2012; originally announced May 2012.

Comments: 6 pages, 10 figures, accepted for publication in ApJ

Thirty years of SN 1980K: Evidence for light echoes

Authors: Ben E. K. Sugerman, Jennifer E. Andrews, Michael J. Barlow, Geoffrey C. Clayton, Barbara Ercolano, Parviz Ghavamian, Robert C. Kennicutt Jr., Oliver Krause, Margaret Meixner, Masaaki Otsuka

Abstract: 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- ana… ▽ More 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. △ Less

Submitted 15 February, 2012; v1 submitted 14 February, 2012; originally announced February 2012.

Comments: Accepted 2/14/12 to be published in ApJ. 15 pages, 10 figures

Late-time Light Curves of Type II Supernovae: Physical Properties of SNe and Their Environment

Authors: Masaaki Otsuka, Margaret Meixner, Nino Panagia, Joanna Fabbri, Michael J. Barlow, Geoffrey C. Clayton, Joseph S. Gallagher, Ben E. K. Sugerman, Roger Wesson, Jennifer E. Andrews, Barbara Ercolano, Douglas Welch

Abstract: 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… ▽ More 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. △ Less

Submitted 23 November, 2011; v1 submitted 21 November, 2011; originally announced November 2011.

Comments: We corrected the 56Ni mass of SN2005cs and Figures 8 (a) and 8 (c)

The effects of dust on the optical and infrared evolution of SN 2004et

Authors: J. Fabbri, M. Otsuka, M. J. Barlow, Joseph S. Gallagher, R. Wesson, B. E. K. Sugerman, Geoffrey C. Clayton, M. Meixner, J. E. Andrews, D. L. Welch, B. Ercolano

Abstract: 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, alon… ▽ More 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. △ Less

Submitted 14 October, 2011; originally announced October 2011.

Comments: 24 pages, 12 figures, 9 tables, published in MNRAS

Herschel Detects a Massive Dust Reservoir in Supernova 1987A

Authors: M. Matsuura, E. Dwek, M. Meixner, M. Otsuka, B. Babler, M. J. Barlow, J. Roman-Duval, C. Engelbracht, K. Sandstrom, M. Lakicevic, J. Th. van Loon, G. Sonneborn, G. C. Clayton, K. S. Long, P. Lundqvist, T. Nozawa, K. D. Gordon, S. Hony, P. Panuzzo, K. Okumura, K. A. Misselt, E. Montiel, M. Sauvage

Abstract: 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… ▽ More 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. △ Less

Submitted 14 July, 2011; v1 submitted 7 July, 2011; originally announced July 2011.

Comments: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. To appear in Science; available in Science Express on 7th July, DOI:10.1126/science.1205983 at http://www.sciencemag.org/lookup/doi/10.1126/science.1205983 ; 14 pages, 2 figures, 4 tables. High resolution images are available from the authors

Evidence for Pre-Existing Dust in the Bright Type IIn SN 2010jl

Authors: J. E. Andrews, Geoffrey C. Clayton, R. Wesson, B. E. K. Sugerman, M. J. Barlow, J. Clem, B. Ercolano, J. Fabbri, J. S. Gallagher, A. Landolt, M. Meixner, M. Otsuka, D. Riebel, D. L. Welch

Abstract: 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 l… ▽ More 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 η Carinae. △ Less

Submitted 2 June, 2011; originally announced June 2011.

Comments: Accepted in AJ

Photometric and Spectroscopic Evolution of the Type IIP SN 2007it to Day 944

Authors: J. E. Andrews, B. E. K Sugerman, Geoffrey Clayton, J. S. Gallagher, M. J. Barlow, J. Clem, B. Ercolano, J. Fabbri, M. Meixner, M. Otsuka, D. L. Welch, R. Wesson

Abstract: 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 emissi… ▽ More 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. △ Less

Submitted 11 February, 2011; originally announced February 2011.

Comments: Accepted for publication in ApJ

Dust and Chemical Abundances of the Sagittarius dwarf Galaxy Planetary Nebula Hen2-436

Authors: Masaaki Otsuka, Margaret Meixner, David Riebel, Siek Hyung, Akito Tajitsu, Hideyuki Izumiura

Abstract: 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 pres… ▽ More 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. △ Less

Submitted 6 December, 2010; originally announced December 2010.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

The Origin and Evolution of the Halo PN BoBn 1: From a Viewpoint of Chemical Abundances Based on Multiwavelength Spectra

Authors: M. Otsuka, A. Tajitsu, S. Hyung, H. Izumiura

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2010; originally announced August 2010.

Comments: 28 pages, 14 figures, accepted for publication in ApJ

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

Authors: M. Meixner, F. Galliano, S. Hony, J. Roman-Duval, T. Robitaille, P. Panuzzo, M. Sauvage, K. Gordon, C. Engelbracht, K. Misselt, K. Okumura, T. Beck, J. -P. Bernard, A. Bolatto, C. Bot, M. Boyer, S. Bracker, L. R. Carlson, G. C. Clayton, C. -H. R. Chen, E. Churchwell, Y. Fukui, M. Galametz, J. L. Hora, A. Hughes , et al. (30 additional authors not shown)

Abstract: 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 scie… ▽ More 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 μm and SPIRE 250, 350 and 500 μ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 μ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. △ Less

Submitted 4 June, 2010; originally announced June 2010.

Comments: To be published in Astronomy & Astrophysics, Herschel First Results Issue