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

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

Aims. As a new approach to the study of resolved stellar populations in nearby galaxies, our goal is to demonstrate in NGC300 that integral field spectroscopy with high spatial resolution and excellent seeing conditions reaches an unprecedented depth in severely crowded fields. Methods. MUSE observations with seven pointings in NGC300 have resulted in datacubes that are analyzed in four ways: (1) PSF-fitting 3D spectroscopy with PampelMUSE yields deblended spectra of individually distinguishable stars. The technique also provides samples of planetary nebulae that are complete down to m5007=28. (2) pseudo-monochromatic images, created at the wavelengths of the most important emission lines and corrected for continuum light by using the P3D visualization tool, provide maps of HII regions, SNR, and the diffuse ISM at a high level of sensitivity, allowing for the discovery of planetary nebulae, WR stars etc. (3) The use of the P3D line-fitting tool yields emission line fluxes, surface brightness, and kinematic information for gaseous objects, corrected for absorption line profiles of the underlying stellar population. (4) Visual inspection of the datacubes is demonstrated to be effcient for data mining and the discovery of background galaxies and unusual objects. Results. We present a catalogue of luminous stars, rare stars such as WR and other emission line stars, carbon stars, symbiotic star candidates, planetary nebulae, HII regions, supernova remnants, giant shells, peculiar diffuse and filamentary emission line objects, and background galaxies, along with their spectra. Conclusions. The technique of crowded-field 3D spectroscopy is capable of deblending individual bright stars, the unresolved background of faint stars, gaseous nebulae, and the diffuse component of the interstellar medium, resulting in unprecedented legacy value for observations of nearby galaxies with MUSE.

The "Antennae Galaxy" (NGC 4038/39) is the closest major interacting galaxy system and therefore often taken as merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star-formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect HII regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of HII regions, showing fainter luminosities. By comparing HII region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each HII region against ionization-parameter sensitive emission line ratios. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that just from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking HII regions for the diffuse ionized gas in the Antennae.

ABRIGED: Quantifying the number, type and distribution of W-R stars is a key component in the context of galaxy evolution, since they put constraints on the age of the star formation bursts. Nearby galaxies (d<5 Mpc) are particularly relevant in this context since they fill the gap between studies in the Local Group, where individual stars can be resolved, and galaxies in the Local Volume and beyond. We intend to characterize the W-R star population in NGC625, a low-metallicity dwarf galaxy suffering a currently declining burst of star formation. Optical IFS data have been obtained with the VIMOS-IFU covering the starburst region. We estimate the number of W-R stars using a linear combination of 3 W-R templates: 1 early-type nitrogen (WN) star, 1 late-type WN star and 1 carbon-type (WC) star (or oxygen-type (WO) star). Fits using several ensembles of templates were tested. Results were confronted with: i) high spatial resolution HST photometry; ii) numbers of W-R stars in nearby galaxies; iii) model predictions. The W-R star population is spread over the main body of the galaxy, not necessarily coincident with the overall stellar distribution. Our best estimation for the number of W-R stars yields a total of 28 W-R stars in the galaxy, out of which 17 are early- type WN, 6 are late-type WN and 5 are WC stars. The width of the stellar features nicely correlates with the dominant W-R type found in each aperture. The distribution of the different types of WR in the galaxy is roughly compatible with the way star formation has propagated in the galaxy, according to previous findings using HST images. Fits using templates at the metallicity of the LMC yield more reasonable number of W-R than those using templates at the metallicity of the SMC. Given the metallicity of NGC 625, this suggests a non-linear relation between the metallicity and the luminosity of the W-R spectral features.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

Heat conduction has been found a plausible solution to explain discrepancies between expected and measured temperatures in hot bubbles of planetary nebulae (PNe). While the heat conduction process depends on the chemical composition, to date it has been exclusively studied for pure hydrogen plasmas in PNe. A smaller population of PNe show hydrogen-deficient and helium- and carbon-enriched surfaces surrounded by bubbles of the same composition; considerable differences are expected in physical properties of these objects in comparison to the pure hydrogen case. The aim of this study is to explore how a chemistry-dependent formulation of the heat conduction affects physical properties and how it affects the X-ray emission from PN bubbles of hydrogen-deficient stars. We extend the description of heat conduction in our radiation hydrodynamics code to work with any chemical composition. We then compare the bubble-formation process with a representative PN model using both the new and the old descriptions. We also compare differences in the resulting X-ray temperature and luminosity observables of the two descriptions. The improved equations show that the heat conduction in our representative model of a hydrogen-deficient PN is nearly as efficient with the chemistry-dependent description; a lower value on the diffusion coefficient is compensated by a slightly steeper temperature gradient. The bubble becomes somewhat hotter with the improved equations, but differences are otherwise minute. The observable properties of the bubble in terms of the X-ray temperature and luminosity are seemingly unaffected.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

Diffuse Interstellar Bands (DIBs) are non-stellar weak absorption features of unknown origin found in the spectra of stars viewed through one or several clouds of Interstellar Medium (ISM). Research of DIBs outside the Milky Way is currently very limited. Specifically spatially resolved investigations of DIBs outside of the Local Group is, to our knowledge, inexistent. Here, we explore the capability of the high sensitivity Integral Field Spectrograph, MUSE, as a tool to map diffuse interstellar bands at distances larger than 100 Mpc. We use MUSE commissioning data for AM 1353-272 B, the member with highest extinction of the "The Dentist's Chair", an interacting system of two spiral galaxies. High signal-to-noise spectra were created by co-adding the signal of many spatial elements distributed in a geometry of concentric elliptical half-rings. We derived decreasing radial profiles for the equivalent width of the $\lambda$5780.5 DIB both in the receding and approaching side of the companion galaxy up to distances of $\sim$4.6 kpc from the center of the galaxy. Likewise, interstellar extinction, as derived from the Halpha/Hbeta line ratio displays a similar trend, with decreasing values towards the external parts. This translates into an intrinsic correlation between the strength of the DIB and the extinction within AM 1353-272 B consistent with the current existing global trend between these quantities when using measurements for both Galactic and extragalactic sight lines. Mapping of DIB strength in the Local Universe as up to now only done for the Milky Way seems feasible. This offers a new approach to study the relationship between DIBs and other characteristics and species of the ISM in different conditions as those found in our Galaxy to the use of galaxies in the Local Group and/or single sightlines towards supernovae, quasars and galaxies outside the Local Group.

Studies of deep photometry of galaxies have presented discoveries of excess light in surface-brightness and colour profiles at large radii in the form of diffuse faint haloes and thick discs. In a majority of the cases, it has seemed necessary to use exotic stellar populations or alternative physical solutions to explain the excess. Few studies have carefully scrutinized the role of scattered light in this context. I explore the influence of scattered light on ground-based observations of haloes and thick discs around edge-on galaxies, haloes around face-on disc galaxies, host galaxies around blue compact galaxies (BCGs), and haloes around elliptical galaxies. Surface-brightness structures of all considered types of galaxies are modelled and analysed to compare scattered-light haloes and thick discs with measurements. I simulate the influence of scattered light and accurate sky subtraction on simplified S\'ersic-type and face-on disc galaxy models. All galaxy models are convolved with both lower-limit and brighter point spread functions (PSFs); for a few galaxies it was possible to use dedicated PSFs. The results show bright scattered-light haloes and high amounts of red excess at large radii and faint surface brightnesses for nearly all types of galaxies; exceptions are the largest elliptical-type galaxies where the influence of scattered light is smaller. Studies have underestimated the role of scattered light to explain their surface-brightness profiles. My analysis shows surface-brightness profiles that include scattered light that are very similar to and overlap measurements at all radii. The derivation of physical properties of haloes, thick discs, and BCG hosts from diffuse data is misleading since accurate and radially extended PSFs are non-existent. Significantly improved analyses that include new measurements of PSFs are required to study diffuse haloes further.

We present X-ray spectral analysis of 20 point-like X-ray sources detected in Chandra Planetary Nebula Survey (ChanPlaNS) observations of 59 planetary nebulae (PNe) in the solar neighborhood. Most of these 20 detections are associated with luminous central stars within relatively young, compact nebulae. The vast majority of these point-like X-ray-emitting sources at PN cores display relatively "hard" ($\geq0.5$~keV) X-ray emission components that are unlikely to be due to photospheric emission from the hot central stars (CSPN). Instead, we demonstrate that these sources are well modeled by optically-thin thermal plasmas. From the plasma properties, we identify two classes of CSPN X-ray emission: (1) high-temperature plasmas with X-ray luminosities, $L_{\rm X}$, that appear uncorrelated with the CSPN bolometric luminosity, $L_{\rm bol}$; and (2) lower-temperature plasmas with $L_{\rm X}/L_{\rm bol}\sim10^{-7}$. We suggest these two classes correspond to the physical processes of magnetically active binary companions and self-shocking stellar winds, respectively. In many cases this conclusion is supported by corroborative multiwavelength evidence for the wind and binary properties of the PN central stars. By thus honing in on the origins of X-ray emission from PN central stars, we enhance the ability of CSPN X-ray sources to constrain models of PN shaping that invoke wind interactions and binarity.

It is not yet clear whether magnetic fields play an essential role in shaping planetary nebulae (PNe), or whether stellar rotation alone and/or a close binary companion can account for the variety of the observed nebular morphologies. In a quest for empirical evidence verifying or disproving the role of magnetic fields in shaping PNe, we follow up on previous attempts to measure the magnetic field in a representative sample of PN central stars. We obtained low-resolution polarimetric spectra with FORS 2 at VLT for a sample of twelve bright central stars of PNe with different morphology, including two round nebulae, seven elliptical nebulae, and three bipolar nebulae. Two targets are Wolf-Rayet type central stars. For the majority of the observed central stars, we do not find any significant evidence for the existence of surface magnetic fields. However, our measurements may indicate the presence of weak mean longitudinal magnetic fields of the order of 100 Gauss in the central star of the young elliptical planetary nebula IC 418, as well as in the Wolf-Rayet type central star of the bipolar nebula Hen2-113 and the weak emission line central star of the elliptical nebula Hen2-131. A clear detection of a 250 G mean longitudinal field is achieved for the A-type companion of the central star of NGC 1514. Some of the central stars show a moderate night-to-night spectrum variability, which may be the signature of a variable stellar wind and/or rotational modulation due to magnetic features. We conclude that strong magnetic fields of the order of kG are not widespread among PNe central stars. Nevertheless, simple estimates based on a theoretical model of magnetized wind bubbles suggest that even weak magnetic fields below the current detection limit of the order of 100 G may well be sufficient to contribute to the shaping of PNe throughout their evolution.

We present results from the most recent set of observations obtained as part of the Chandra X-ray observatory Planetary Nebula Survey (ChanPlaNS), the first comprehensive X-ray survey of planetary nebulae (PNe) in the solar neighborhood (i.e., within ~1.5 kpc of the Sun). The survey is designed to place constraints on the frequency of appearance and range of X-ray spectral characteristics of X-ray-emitting PN central stars and the evolutionary timescales of wind-shock-heated bubbles within PNe. ChanPlaNS began with a combined Cycle 12 and archive Chandra survey of 35 PNe. ChanPlaNS continued via a Chandra Cycle 14 Large Program which targeted all (24) remaining known compact (R_neb <~ 0.4 pc), young PNe that lie within ~1.5 kpc. Results from these Cycle 14 observations include first-time X-ray detections of hot bubbles within NGC 1501, 3918, 6153, and 6369, and point sources in HbDs 1, NGC 6337, and Sp 1. The addition of the Cycle 14 results brings the overall ChanPlaNS diffuse X-ray detection rate to ~27% and the point source detection rate to ~36%. It has become clearer that diffuse X-ray emission is associated with young (<~5x10^3 yr), and likewise compact (R_neb<~0.15 pc), PNe with closed structures and high central electron densities (n_e>~1000 cm^-3), and rarely associated with PNe that show H_2 emission and/or pronounced butterfly structures. Hb 5 is one such exception of a PN with a butterfly structure that hosts diffuse X-ray emission. Additionally, of the five new diffuse X-ray detections, two host [WR]-type CSPNe, NGC 1501 and NGC 6369, supporting the hypothesis that PNe with central stars of [WR]-type are likely to display diffuse X-ray emission.

All telescopes and instruments are to some degree affected by scattered light. It is possible to estimate the amount of such scattered light, and even correct for it, with a radially extended point spread function (PSF). The outer parts of the PSF have only rarely been determined, since they are faint and therefore difficult to measure. A mostly complete overview of existing properties and measurements of radially extended PSFs is presented, to both show their similarities and to indicate how bright extended objects can be used to measure the faintest regions. The importance of the far wings of the PSF and their possible temporal variations are demonstrated in three edge-on galaxy models. The same study is applied to the first edge-on galaxy where earlier observations reveal a halo, NGC 5907. All PSFs were collected in two diagrams, after they were offset or normalized, when that was possible. Surface-brightness structures of edge-on galaxies were modelled and analysed to study scattered-light haloes that result with an exponential disc. The models were convolved with both a lower-limit PSF and a more average PSF. The PSF of the observed data could be used in the case of NGC 5907. The comparison of the PSFs demonstrates a lower-limit $r^{-2}$ power-law decline at larger radii. The analysis of the galaxy models shows that also the outer parts of the PSF are important to correctly model and analyse observations and, in particular, fainter regions. The reassessed analysis of the earlier measurements of NGC 5907 reveals an explanation for the faint halo in scattered light, within the quoted level of accuracy.

We used the PMAS integral field spectrograph to obtain large sets of radial velocities in the central regions of three northern Galactic globular clusters: M3, M13, and M92. By applying the novel technique of crowded field 3D spectroscopy, we measured radial velocities for about 80 stars within the central ~ 10 arcsec of each cluster. These are by far the largest spectroscopic datasets obtained in the innermost parts of these clusters up to now. To obtain kinematical data across the whole extent of the clusters, we complement our data with measurements available in the literature. We combine our velocity measurements with surface brightness profiles to analyse the internal dynamics of each cluster using spherical Jeans models, and investigate whether our data provide evidence for an intermediate-mass black hole in any of the clusters. The surface brightness profiles reveal that all three clusters are consistent with a core profile, although shallow cusps cannot be excluded. We find that spherical Jeans models with a constant mass-to-light ratio provide a good overall representation of the kinematical data. A massive black hole is required in none of the three clusters to explain the observed kinematics. Our 1sigma (3sigma) upper limits are 5300 M_sun (12000 M_sun) for M3, 8600 M_sun (13000 M_sun) for M13, and 980 M_sun (2700 M_sun) for M92. A puzzling circumstance is the existence of several potential high velocity stars in M3 and M13, as their presence can account for the majority of the discrepancies that we find in our mass limits compared to M92.

We present the result of a study on the expansion properties and internal kinematics of round/elliptical planetary nebulae of the Milky Way disk, the halo, and of the globular cluster M15. The purpose of this study is to considerably enlarge the small sample of nebulae with precisely determined expansion properties. To this aim, we selected a representative sample of objects with different evolutionary stages and metallicities and conducted high-resolution echelle spectroscopy. In most cases, we succeeded in detecting the weak signals from the outer nebular shell which are attached to the main line emission from the bright nebular rim. Next to the measurement of the motion of the rim gas by decomposition of the main line components into Gaussians, we were able to measure separately, for most objects for the first time, the gas velocity immediately behind the leading shock of the shell, i.e. the post-shock velocity. We more than doubled the number of objects for which the velocities of both rim and shell are known and confirm that the overall expansion of planetary nebulae is accelerating with time. There are, however, differences between the expansion behaviour of the shell and the rim. This observed distinct velocity evolution of both rim and shell is explained by radiation-hydrodynamics simulations, at least qualitatively. Because of the time-dependent boundary conditions, a planetary nebula will never evolve into a simple self-similar expansion. Also the metal-poor objects behave as theory predicts: The post-shock velocities are higher and the rim flow velocities are equal or even lower compared to disk objects at similar evolutionary stage. We detected, for the first time, in some objects an asymmetric expansion behaviour: The relative expansions between rim and shell appear to be different for the receding and approaching parts of the nebular envelope.

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.

We present and discuss new visual wavelength-range observations of the inner regions of the supernova remnant SNR 0540-69.3 that is located in the Large Magellanic Cloud (LMC). These observations provide us with more spatial and spectral information than were previously available for this object. We use these data to create a detailed three-dimensional model of the remnant, assuming linear expansion of the ejecta. With the observations and the model we study the general three-dimensional structure of the remnant, and the influence of an active region in the remnant - a "blob" - that we address in previous papers. We used the fibre-fed integral-field spectrograph VIMOS at the Very Large Telescope of the European Southern Observatory. The observations provide us with three-dimensional data in [OIII]5007 and [SII]6717,6731 at an 0.33x0.33 spatial sampling and a velocity resolution of about 35 km/s. We decomposed the two, partially overlapping, sulphur lines and used them to calculate electron densities across the remnant at high signal-to-noise ratio. Our analysis reveals a structure that stretches from the position of the "blob", and into the plane of the sky at a position angle of about 60 degrees. We speculate that the pulsar is positioned along this activity axis, where it has a velocity along the line of sight of a few hundred km/s. The "blob" is most likely a region of shock activity, as it is mainly bright in [SII]; future observations of [OII]3729 would be useful to test whether the S/O abundance ratio is higher than average for that location in the remnant. The striking resemblance in X-rays between the pulsar wind nebula (PWN) of SNR 0540-69.3 and the Crab, in combination with our findings in this paper, suggests that the symmetry axis is part of a torus in the PWN. (abridged)

Progress is being made in code discoverability and preservation, but as discussed at ADASS XXI, many codes still remain hidden from public view. With the Astrophysics Source Code Library (ASCL) now indexed by the SAO/NASA Astrophysics Data System (ADS), the introduction of a new journal, Astronomy & Computing, focused on astrophysics software, and the increasing success of education efforts such as Software Carpentry and SciCoder, the community has the opportunity to set a higher standard for its science by encouraging the release of software for examination and possible reuse. We assembled representatives of the community to present issues inhibiting code release and sought suggestions for tackling these factors. The session began with brief statements by panelists; the floor was then opened for discussion and ideas. Comments covered a diverse range of related topics and points of view, with apparent support for the propositions that algorithms should be readily available, code used to produce published scientific results should be made available, and there should be discovery mechanisms to allow these to be found easily. With increased use of resources such as GitHub (for code availability), ASCL (for code discovery), and a stated strong preference from the new journal Astronomy & Computing for code release, we expect to see additional progress over the next few years.

[Abridged] Detecting cosmic ray hits (cosmics) in fiber-fed IFS data of single exposures is a challenging task, because of the complex signal recorded by IFS instruments. Existing detection algorithms are commonly found to be unreliable in the case of IFS data and the optimal parameter settings are usually unknown a-priori for a given dataset. The CALIFA survey generates hundreds of IFS datasets for which a reliable and robust detection algorithm for cosmics is required as an important part of the fully automatic CALIFA data reduction pipeline. We developed a novel algorithm, PyCosmic, which combines the edge-detection algorithm of L.A.Cosmic with a point-spread function convolution scheme. We generated mock data to compute the efficiency of different algorithms for a wide range of characteristic fibre-fed IFS datasets using the PMAS and VIMOS IFS instruments as representative cases. PyCosmic is the only algorithm that achieves an acceptable detection performance for CALIFA data. We find that PyCosmic is the most robust tool with a detection rate of >~90% and a false detection rate <5% for any of the tested IFS data. It has one less free parameter than the L.A.Cosmic algorithm. Only for strongly undersampled IFS data does L.A.Cosmic exceed the performance of PyCosmic by a few per cent. DCR never reaches the efficiency of the other two algorithms and should only be used if computational speed is a concern. Thus, PyCosmic appears to be the most versatile cosmics detection algorithm for IFS data. It is implemented in the new CALIFA data reduction pipeline as well as in recent versions of the multi-instrument IFS pipeline P3D.

Astronomical instrumentation is most of the time faced with challenging requirements in terms of sensitivity, stability, complexity, etc., and therefore leads to high performance developments that at first sight appear to be suitable only for the specific design application at the telescope. However, their usefulness in other disciplines and for other applications is not excluded. The ERA2 facility is a lab demonstrator, based on a high-performance astronomical spectrograph, which is intended to explore the innovation potential of fiber-coupled multi-channel spectroscopy for spatially resolved spectroscopy in life science, material sciences, and other areas of research.

We present an overview of the initial results from the Chandra Planetary Nebula Survey (ChanPlaNS), the first systematic (volume-limited) Chandra X-ray Observatory survey of planetary nebulae (PNe) in the solar neighborhood. The first phase of ChanPlaNS targeted 21 mostly high-excitation PNe within ~1.5 kpc of Earth, yielding 4 detections of diffuse X-ray emission and 9 detections of X-ray-luminous point sources at the central stars (CSPNe) of these objects. Combining these results with those obtained from Chandra archival data for all (14) other PNe within ~1.5 kpc that have been observed to date, we find an overall X-ray detection rate of ~70%. Roughly 50% of the PNe observed by Chandra harbor X-ray-luminous CSPNe, while soft, diffuse X-ray emission tracing shocks formed by energetic wind collisions is detected in ~30%; five objects display both diffuse and point-like emission components. The presence of X-ray sources appears correlated with PN density structure, in that molecule-poor, elliptical nebulae are more likely to display X-ray emission (either point-like or diffuse) than molecule-rich, bipolar or Ring-like nebulae. All but one of the X-ray point sources detected at CSPNe display X-ray spectra that are harder than expected from hot (~100 kK) central star photospheres, possibly indicating a high frequency of binary companions to CSPNe. Other potential explanations include self-shocking winds or PN mass fallback. Most PNe detected as diffuse X-ray sources are elliptical nebulae that display a nested shell/halo structure and bright ansae; the diffuse X-ray emission regions are confined within inner, sharp-rimmed shells. All sample PNe that display diffuse X-ray emission have inner shell dynamical ages <~5x10^3 yr, placing firm constraints on the timescale for strong shocks due to wind interactions in PNe.

Based on time-dependent radiation-hydrodynamics simulations of the evolution of Planetary Nebulae (PNe), we have carried out a systematic parameter study to address the non-trivial question of how the diffuse X-ray emission of PNe with closed central cavities is expected to depend on the evolutionary state of the nebula, the mass of the central star, and the metallicity of stellar wind and circumstellar matter. We have also investigated how the model predictions depend on the treatment of thermal conduction at the interface between the central `hot bubble' and the `cool' inner nebula, and compare the results with recent X-ray observations. Our study includes models whose properties resemble the extreme case of PNe with Wolf-Rayet type central stars. Indeed, such models are found to produce the highest X-ray luminosities.

Integral Field Spectroscopy is a powerful observing technique for Astronomy that is becoming available at most ground-based observatories as well as in space. The complex data obtained with this technique require new approaches for visualization. Typical requirements and the p3d tool, as an example, are discussed.

By means of hydrodynamical models we do the first investigations of how the properties of planetary nebulae are affected by their metal content and what can be learned from spatially unresolved spectrograms of planetary nebulae in distant stellar systems. We computed a new series of 1D radiation-hydrodynamics planetary nebulae model sequences with central stars of 0.595 M_sun surrounded by initial envelope structures that differ only by their metal content. At selected phases along the evolutionary path, the hydrodynamic terms were switched off, allowing the models to relax for fixed radial structure and radiation field into their equilibrium state with respect to energy and ionisation. The analyses of the line spectra emitted from both the dynamical and static models enabled us to systematically study the influence of hydrodynamics as a function of metallicity and evolution. We also recomputed selected sequences already used in previous publications, but now with different metal abundances. These sequences were used to study the expansion properties of planetary nebulae close to the bright cut-off of the planetary nebula luminosity function. Our simulations show that the metal content strongly influences the expansion of planetary nebulae: the lower the metal content, the weaker the pressure of the stellar wind bubble, but the faster the expansion of the outer shell because of the higher electron temperature. This is in variance with the predictions of the interacting-stellar-winds model (or its variants) according to which only the central-star wind is thought to be responsible for driving the expansion of a planetary nebula. Metal-poor objects around slowly evolving central stars become very dilute and are prone to depart from thermal equilibrium because then adiabatic expansion contributes to gas cooling. ...abridged abstract.

The reduction of integral-field spectrograph (IFS) data is demanding work. Many repetitive operations are required in order to convert raw data into, typically a large number of, spectra. This effort can be markedly simplified through the use of a tool or pipeline, which is designed to complete many of the repetitive operations without human interaction. Here we present our semi-automatic data-reduction tool p3d that is designed to be used with fiber-fed IFSs. Important components of p3d include a novel algorithm for automatic finding and tracing of spectra on the detector, and two methods of optimal spectrum extraction in addition to standard aperture extraction. p3d also provides tools to combine several images, perform wavelength calibration and flat field data. p3d is at the moment configured for four IFSs. In order to evaluate its performance we have tested the different components of the tool. For these tests we used both simulated and observational data. We demonstrate that for three of the IFSs a correction for so-called cross-talk due to overlapping spectra on the detector is required. Without such a correction spectra will be inaccurate, in particular if there is a significant intensity gradient across the object. Our tests showed that p3d is able to produce accurate results. p3d is a highly general and freely available tool. It is easily extended to include improved algorithms, new visualization tools and support for additional instruments. The program code can be downloaded from the p3d-project web site http://p3d.sourceforge.net

The actual value of the oxygen abundance of the metal-poor planetary nebula PN G135.9+55.9 has frequently been debated in the literature. We wanted to clarify the situation by making an improved abundance determination based on a study that includes both new accurate observations and new models. We made observations using the method of integral field spectroscopy with the PMAS instrument, and also used ultraviolet observations that were measured with HST-STIS. In our interpretation of the reduced and calibrated spectrum we used for the first time, recent radiation hydrodynamic models, which were calculated with several setups of scaled values of mean Galactic disk planetary nebula metallicities. For evolved planetary nebulae, such as PN G135.9+55.9, it turns out that departures from thermal equilibrium can be significant, leading to much lower electron temperatures, hence weaker emission in collisionally excited lines. Based on our time-dependent hydrodynamic models and the observed emission line [OIII]5007, we found a very low oxygen content of about 1/80 of the mean Galactic disk value. This result is consistent with emission line measurements in the ultraviolet wavelength range. The C/O and Ne/O ratios are unusually high and similar to those of another halo object, BoBn-1.

In the final stages of stellar evolution low- to intermediate-mass stars lose their envelope in increasingly massive stellar winds. Such winds affect the interstellar medium and the galactic chemical evolution as well as the circumstellar envelope where planetary nebulae form subsequently. Characteristics of this mass loss depend on both stellar properties and properties of gas and dust in the wind formation region. In this paper we present an approach towards studies of mass loss using both observations and models, focusing on the stage where the stellar envelope is nearly empty of mass. In a recent study we measure the mass-loss evolution, and other properties, of four planetary nebulae in the Galactic Disk. Specifically we use the method of integral field spectroscopy on faint halos, which are found outside the much brighter central parts of a planetary nebula. We begin with a brief comparison between our and other observational methods to determine mass-loss rates in order to illustrate how they differ and complement each other. An advantage of our method is that it measures the gas component directly requiring no assumptions of properties of dust in the wind. Thereafter we present our observational approach in more detail in terms of its validity and its assumptions. In the second part of this paper we discuss capabilities and assumptions of current models of stellar winds. We propose and discuss improvements to such models that will allow meaningful comparisons with our observations. Currently the physically most complete models include too little mass in the model domain to permit a formation of winds with as high mass-loss rates as our observations show.

Strong mass loss off stars at the tip of the asymptotic giant branch (AGB) profoundly affects properties of these stars and their surroundings, including the subsequent planetary nebula (PN) stage. With this study we wanted to determine physical properties of mass loss by studying weakly emitting halos, focusing on objects in the galactic disk. Halos surround the, up to several thousand times, brighter central regions of PNe. Young halos, specifically, still contain information of the preceeding final mass loss stage on the AGB. In the observations we used the method of integral field spectroscopy with the PMAS instrument. This is the first committed study of halos of PNe that uses this technique. We improved our data analysis by a number of steps. In a study of the influence of scattered light we found that a moderate fraction of intensities in the inner halo originate in adjacent regions. As we combine line intensities of distant wavelengths, and because radial intensity gradients are steep, we corrected for effects of differential atmospheric refraction. In order to increase the signal-to-noise of weak emission lines we introduced a dedicated method to bin spectra of individual spatial elements. We also developed a general technique to subtract telluric lines - without using separate sky exposures. By these steps we avoided introducing errors of several thousand Kelvin to our temperature measurements in the halo. For IC3568 we detected a halo. For M2-2 we found a halo radius that is 2.5 times larger... (abridged)

Models describing dust-driven winds are important for understanding the physical mechanism and properties of mass loss on the asymptotic giant branch. These models are becoming increasingly realistic with more detailed physics included, but also more computationally demanding. The purpose of this study is to clarify to what extent the applied numerical approach affects resulting physical structures of modelled winds, and to discuss resulting changes. Following the previously developed radiation hydrodynamic model - which includes descriptions for time-dependent dust formation and gas-dust drift - and using its physical assumptions and parameters, numerical improvements are introduced. Impacts of the so-called adaptive grid equation and advection schemes are assessed from models calculated with different numerical setups. Results show that wind models are strongly influenced by numerical imprecision, displaying differences in calculated physical properties of up to one hundred per cent. Using a non-adaptive grid, models become periodic (in multiples of stellar pulsation periods), instead of irregular, as obtained previously. Furthermore, the numerical improvements reveal changes in physical structures. The influence of gas-dust drift is confirmed to be highly important, in particular for the dust component. Gas and dust are less tightly coupled than previously, and drastically larger amounts of dust form assuming drift.

The nebular evolution is followed from the vicinity of the asymptotic-giant branch across the Hertzsprung-Russell diagram until the white-dwarf domain is reached, using various central-star models coupled to different initial envelope configurations. Along each sequence the relevant line emissions of the nebulae are computed and analysed. Maximum line luminosities in Hbeta and [OIII] 5007A are achieved at stellar effective temperatures of about 65000K and 95000-100000K, respectively, provided the nebula remains optically thick for ionising photons. In the optically thin case, the maximum line emission occurs at or shortly after the thick/thin transition. Our models suggest that most planetary nebulae with hotter (>~ 45000K) central stars are optically thin in the Lyman continuum, and that their [OIII] 5007A emission fails to explain the bright end of the observed planetary nebulae luminosity function. However, sequences with central stars of >~ 0.6 Msun and rather dense initial envelopes remain virtually optically thick and are able to populate the bright end of the luminosity function. Individual luminosity functions depend strongly on the central-star mass and on the variation of the nebular optical depth with time. Hydrodynamical simulations of planetary nebulae are essential for any understanding of the basic physics behind their observed luminosity function. In particular, our models do not support the claim of Marigo et.al (2004) according to which the maximum 5007A luminosity occurs during the recombination phase well beyond 100 000K when the stellar luminosity declines and the nebular models become, at least partially, optically thick. Consequently, there is no need to invoke relatively massive central stars of, say > 0.7 Msun, to account for the bright end of the luminosity function.

We report on our studies of the physical structure of the planetary nebula (PN) NGC7662. Using (3D) Integral Field Spectroscopy we have been able to measure the electron temperature more accurately and at a larger number of radial locations than before. Here we briefly present our method by which we find a strong positive temperature gradient with increasing radius. According to hydrodynamic models a hot halo, when compared to the central star, can be the product of the passage of an ionization front (e.g. Marten 1993). Such a gradient is not found in equilibrium models, and this finding - when confirmed for other objects - strongly advocates the use of hydrodynamic models when modeling PN halos.

In addition to study extragalactic stellar populations in their integrated light, the detailed analysis of individual resolved objects has become feasible, mainly for luminous giant stars and for extragalactic planetary nebulae (XPNe) in nearby galaxies. A recently started project at the Astrophysical Institute Potsdam (AIP), called ``XPN--Physics'', aims to verify if XPNe are useful probes to measure the chemical abundances of their parent stellar population. The project involves theoretical and observational work packages.

Three closely related stumbling blocks of solar mean field dynamo theory are discussed: how dominant are the small scale fields, how is the alpha effect quenched, and whether magnetic and current helicity fluxes alleviate the quenching? It is shown that even at the largest currently available resolution there is no clear evidence of power law scaling of the magnetic and kinetic energy spectra in turbulence. However, using subgrid scale modeling, some indications of asymptotic equipartition can be found. The frequently used first order smoothing approach to calculate the alpha effect and other transport coefficients is contrasted with the superior minimal tau approximation. The possibility of catastrophic alpha quenching is discussed as a result of magnetic helicity conservation. Magnetic and current helicity fluxes are shown to alleviate catastrophic quenching in the presence of shear. Evidence for strong large scale dynamo action, even in the absence of helicity in the forcing, is presented.

The standard theory of the solar cycle in terms of an alpha-Omega dynamo hinges on a proper understanding of the nonlinear alpha effect. Boundary conditions play a surprisingly important role in determining the magnitude of alpha. For closed boundaries, the total magnetic helicity is conserved, and since the alpha effect produces magnetic helicity of one sign in the large scale field, it must simultaneously produce magnetic helicity of the opposite sign. It is this secondary magnetic helicity that suppresses the dynamo in a potentially catastrophic fashion. Open boundaries allow magnetic helicity to be lost. Simulations are presented that allow an estimate of alpha in the presence of open or closed boundaries, either with or without solar-like differential rotation. In all cases the sign of the magnetic helicity agrees with that observed at the solar surface (negative in the north, positive in the south), where significant amounts of magnetic helicity can be ejected via coronal mass ejections. It is shown that open boundaries tend to alleviate catastrophic alpha quenching. The importance of looking at current helicity instead of magnetic helicity is emphasized and the conceptual advantages are discussed.

A new simulation set-up is proposed for studying mean field dynamo action. The model combines the computational advantages of local cartesian geometry with the ability to include a shear profile that resembles the sun's differential rotation at low latitudes. It is shown that in a two-dimensional mean field model this geometry produces cyclic solutions with dynamo waves traveling away from the equator -- as expected for a positive alpha effect in the northern hemisphere. In three dimensions with turbulence driven by a helical forcing function, an alpha effect is self-consistently generated in the presence of a finite imposed toroidal magnetic field. The results suggest that, due to a finite flux of current helicity out of the domain, alpha quenching appears to be non-catastrophic -- at least for intermediate values of the magnetic Reynolds number. For larger values of the magnetic Reynolds number, however, there is evidence for a reversal of the trend and that $\alpha$ may decrease with increasing magnetic Reynolds number. Control experiments with closed boundaries confirm that in the absence of a current helicity flux, but with shear as before, alpha quenching is always catastrophic and alpha decreases inversely proportional to the magnetic Reynolds number. For solar parameters, our results suggest a current helicity flux of about 0.001 G^2/s. This corresponds to a magnetic helicity flux, integrated over the northern hemisphere and over the 11 year solar cycle, of about 10^{46}Mx^2.

A proper treatment of the non-equilibrium dust formation process is crucial in models of AGB star winds. In this paper the micro-physics of this process is treated in detail, with an emphasis on the effects of drift (drift models). We summarize the description of the dust formation process and make a few additions to previous work. A detailed study shows that different growth species dominate the grain growth rates at different drift velocities. The new models show that the net effect of drift is to significantly increase the amounts of dust, seemingly without affecting the mean wind properties, such as e.g., the mass loss rate. In some cases there is several times more dust in drift models, compared to the values in the corresponding non-drift models. We study the formation of a dust shell in the inner parts of the wind and find that drift plays an active role in accumulating dust to certain narrow regions. In view of the results presented here it is questionable if drift -- under the current assumptions -- can be ignored in the grain growth rates.

We present three-component wind models for carbon rich pulsating AGB stars. In particular we study the effects of drift in models of long-period variables, meaning that the dust is allowed to move relative to the gas (drift models). In addition we investigate the importance of the degree of variability of the wind structures. The wind model contains separate conservation laws for each of the three components of gas, dust and the radiation field. We use two different representations for the gas opacity, resulting in models with different gas densities in the wind. The effects which we investigate here are important for the understanding of the wind mechanism and mass loss of AGB stars. This study is hereby a necessary step towards more reliable interpretations of observations. We find that the effects of drift generally are significant. They cannot be predicted from models calculated without drift. Moreover, the non-drift models showing the lowest mass loss rates, outflow velocities, and the smallest variability in the degree of condensation do not form drift model winds. The wind formation in drift models is, except for a few cases, generally less efficient and the mass loss consequently lower than in the corresponding non-drift models. The effects of drift are generally larger in the more realistic models using that representation of the gas opacity which results in lower densities. The outflow properties of these models are also -- for all cases we have studied -- sensitive to the period of the stellar pulsations. A check of the mass loss rates against a (recent) fit formula shows systematically lower values, in particular in the more realistic models with a low density. The fit is in its current form inapplicable to the new models presented here.

Radiative acceleration of newly-formed dust grains and transfer of momentum from the dust to the gas plays an important role for driving winds of AGB stars. Therefore a detailed description of the interaction of gas and dust is a prerequisite for realistic models of such winds. In this paper we present the method and first results of a three-component time-dependent model of dust-driven AGB star winds. With the model we plan to study the role and effects of the gas-dust interaction on the mass loss and wind formation. The wind model includes separate conservation laws for each of the three components of gas, dust and the radiation field and is developed from an existing model which assumes position coupling between the gas and the dust. As a new feature we introduce a separate equation of motion for the dust component in order to fully separate the dust phase from the gas phase. The transfer of mass, energy and momentum between the phases is treated by interaction terms. We also carry out a detailed study of the physical form and influence of the momentum transfer term (the drag force) and three approximations to it. In the present study we are interested mainly in the effect of the new treatment of the dust velocity on dust-induced instabilities in the wind. As we want to study the consequences of the additional freedom of the dust velocity on the model we calculate winds both with and without the separate dust equation of motion. The wind models are calculated for several sets of stellar parameters. We find that there is a higher threshold in the carbon/oxygen abundance ratio at which winds form in the new model. The winds of the new models, which include drift, differ from the previously stationary winds, and the winds with the lowest mass loss rates no longer form.