-
Three-component modelling of O-rich AGB star winds I. Effects of drift using forsterite
Authors:
C. Sandin,
L. Mattsson,
K. L. Chubb,
M. Ergon,
P. M. Weilbacher
Abstract:
Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic mod…
▽ More
Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic model code T-800 that includes frequency-dependent radiative transfer, dust extinction based on Mie scattering, grain growth and ablation, gas-to-dust drift using one mean grain size, a piston that simulates stellar pulsations, and an accurate high spatial resolution numerical scheme. To enable this study, we calculated new gas opacities based on the ExoMol database, and we extended the model code to handle the formation of minerals that may form in M-type stars. We determine the effects of drift by comparing drift models to our new and extant non-drift models. Three out of four new drift models show high drift velocities, 87-310 km/s. Our new drift model mass-loss rates are 1.7-13 per cent of the corresponding values of our non-drift models, but compared to the results of two extant non-drift models that use the same stellar parameters, these same values are 0.33-1.5 per cent. Meanwhile, a comparison of other properties such as the expansion velocity and grain size show similar values. Our results, which are based on single-component forsterite particles, show that the inclusion of gas-to-drift is of fundamental importance in stellar wind models driven by such transparent grains. Assuming that the drift velocity is insignificant, properties such as the mass-loss rate may be off from more realistic values by a factor of 50 or more.
△ Less
Submitted 13 July, 2023; v1 submitted 3 January, 2023;
originally announced January 2023.
-
Spectral characterisation of inertial particle clustering in turbulence
Authors:
N. E. L. Haugen,
A. Brandenburg,
C. Sandin,
L. Mattsson
Abstract:
Clustering of inertial particles is important for many types of astrophysical and geophysical turbulence, but it has been studied predominately for incompressible flows. Here we study compressible flows and compare clustering in both compressively (irrotationally) and vortically (solenoidally) forced turbulence. Vortically and compressively forced flows are driven stochastically either by solenoid…
▽ More
Clustering of inertial particles is important for many types of astrophysical and geophysical turbulence, but it has been studied predominately for incompressible flows. Here we study compressible flows and compare clustering in both compressively (irrotationally) and vortically (solenoidally) forced turbulence. Vortically and compressively forced flows are driven stochastically either by solenoidal waves or by circular expansion waves, respectively. For compressively forced flows, the power spectrum of the density of inertial particles is a useful tool for displaying particle clustering relative to the fluid density enhancement. Power spectra are shown to be particularly sensitive for studying large-scale particle clustering, while conventional tools such as radial distribution functions are more suitable for studying small-scale clustering. Our primary finding is that particle clustering through shock interaction is particularly prominent in turbulence driven by spherical expansion waves. It manifests itself through a double-peaked distribution of spectral power as a function of Stokes number. The two peaks are associated with two distinct clustering mechanisms; shock interaction for smaller Stokes numbers and the centrifugal sling effect for larger values. The clustering of inertial particles is associated with the formation of caustics. Such caustics can only be captured in the Lagrangian description, which allows us to assess the relative importance of caustics in vortically and compressively forced turbulence. We show that the statistical noise resulting from the limited number of particles in the Lagrangian description can be removed from the particle power spectra, allowing us a more detailed comparison of the residual spectra. We focus on the Epstein drag law relevant for rarefied gases, but show that our findings apply also to the usual Stokes drag.
△ Less
Submitted 16 December, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Wide Field Spectral Imaging with Shifted Excitation Raman Difference Spectroscopy Using the Nod and Shuffle Technique
Authors:
Florian Korinth,
Elmar Schmälzlin,
Clara Stiebing,
Tanya Urrutia,
Genoveva Micheva,
Christer Sandin,
André Müller,
Martin Maiwald,
Bernd Sumpf,
Christoph Krafft,
Günther Tränkle,
Martin M. Roth,
Jürgen Popp
Abstract:
Wide field Raman imaging using the integral field spectroscopy approach was used as a fast, one shot imaging method for the simultaneous collection of all spectra composing a Raman image. For the suppression of autofluorescence and background signals such as room light, shifted excitation Raman difference spectroscopy (SERDS) was applied to remove background artifacts in Raman spectra. To reduce a…
▽ More
Wide field Raman imaging using the integral field spectroscopy approach was used as a fast, one shot imaging method for the simultaneous collection of all spectra composing a Raman image. For the suppression of autofluorescence and background signals such as room light, shifted excitation Raman difference spectroscopy (SERDS) was applied to remove background artifacts in Raman spectra. To reduce acquisition times in wide field SERDS imaging, we adapted the nod and shuffle technique from astrophysics and implemented it into a wide field SERDS imaging setup. In our adapted version, the nod corresponds to the change in excitation wavelength, whereas the shuffle corresponds to the shifting of charges up and down on a Charge-Coupled Device (CCD) chip synchronous to the change in excitation wavelength. We coupled this improved wide field SERDS imaging setup to diode lasers with 784.4/785.5 and 457.7/458.9 nm excitation and applied it to samples such as paracetamol and aspirin tablets, polystyrene and polymethyl methacrylate beads, as well as pork meat using multiple accumulations with acquisition times in the range of 50 to 200 ms. The results tackle two main challenges of SERDS imaging: gradual photobleaching changes the autofluorescence background, and multiple readouts of CCD detector prolong the acquisition time.
△ Less
Submitted 27 November, 2020;
originally announced November 2020.
-
The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
Authors:
A. Brandenburg,
A. Johansen,
P. A. Bourdin,
W. Dobler,
W. Lyra,
M. Rheinhardt,
S. Bingert,
N. E. L. Haugen,
A. Mee,
F. Gent,
N. Babkovskaia,
C. -C. Yang,
T. Heinemann,
B. Dintrans,
D. Mitra,
S. Candelaresi,
J. Warnecke,
P. J. Käpylä,
A. Schreiber,
P. Chatterjee,
M. J. Käpylä,
X. -Y. Li,
J. Krüger,
J. R. Aarnes,
G. R. Sarson
, et al. (12 additional authors not shown)
Abstract:
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the…
▽ More
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Three-component modelling of C-rich AGB star winds V. Effects of frequency-dependent radiative transfer including drift
Authors:
Christer Sandin,
Lars Mattsson
Abstract:
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 extinc…
▽ More
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.
△ Less
Submitted 4 September, 2020; v1 submitted 19 June, 2020;
originally announced June 2020.
-
Origins Space Telescope Mission Concept Study Report
Authors:
M. Meixner,
A. Cooray,
D. Leisawitz,
J. Staguhn,
L. Armus,
C. Battersby,
J. Bauer,
E. Bergin,
C. M. Bradford,
K. Ennico-Smith,
J. Fortney,
T. Kataria,
G. Melnick,
S. Milam,
D. Narayanan,
D. Padgett,
K. Pontoppidan,
A. Pope,
T. Roellig,
K. Sandstrom,
K. Stevenson,
K. Su,
J. Vieira,
E. Wright,
J. Zmuidzinas
, et al. (44 additional authors not shown)
Abstract:
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 lar…
▽ More
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.
△ Less
Submitted 23 December, 2019; v1 submitted 12 December, 2019;
originally announced December 2019.
-
MUSE crowded field 3D spectroscopy in NGC300 I. First results from central fields
Authors:
M. M. Roth,
C. Sandin,
S. Kamann,
T. -O. Husser,
P. M. Weilbacher,
A. Monreal-Ibero,
R. Bacon,
M. den Brok,
S. Dreizler,
A. Kelz,
R. A. Marino,
M. Steinmetz
Abstract:
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…
▽ More
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.
△ Less
Submitted 11 June, 2018;
originally announced June 2018.
-
On the Origin of Diffuse Ionized Gas in the Antennae Galaxy
Authors:
Peter M. Weilbacher,
Ana Monreal-Ibero,
Anne Verhamme,
Christer Sandin,
Matthias Steinmetz,
Wolfram Kollatschny,
Davor Krajnović,
Sebastian Kamann,
Martin M. Roth,
Santiago Erroz-Ferrer,
Raffaella Anna Marino,
Michael V. Maseda,
Martin Wendt,
Roland Bacon,
Stefan Dreizler,
Johan Richard,
Lutz Wisotzki
Abstract:
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 ti…
▽ More
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.
△ Less
Submitted 14 December, 2017; v1 submitted 12 December, 2017;
originally announced December 2017.
-
The Wolf-Rayet star population in the dwarf galaxy NGC 625
Authors:
Ana Monreal-Ibero,
Jeremy R. Walsh,
Jorge Iglesias-Paramo,
Christer Sandin,
M. Relano,
E. Perez-Montero,
J. Vilchez
Abstract:
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 beyo…
▽ More
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.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
Hot bubbles of planetary nebulae with hydrogen-deficient winds I. Heat conduction in a chemically stratified plasma
Authors:
Christer Sandin,
Matthias Steffen,
Detlef Schönberner,
Ute Rühling
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 December, 2015;
originally announced December 2015.
-
Extended Lyman alpha haloes around individual high-redshift galaxies revealed by MUSE
Authors:
L. Wisotzki,
R. Bacon,
J. Blaizot,
J. Brinchmann,
E. C. Herenz,
J. Schaye,
N. Bouché,
S. Cantalupo,
T. Contini,
C. M. Carollo,
J. Caruana,
J. -B. Courbot,
E. Emsellem,
S. Kamann,
J. Kerutt,
F. Leclercq,
S. J. Lilly,
V. Patrício,
C. Sandin,
M. Steinmetz,
L. A. Straka,
T. Urrutia,
A. Verhamme,
P. M. Weilbacher,
M. Wendt
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 December, 2015; v1 submitted 17 September, 2015;
originally announced September 2015.
-
A MUSE map of the central Orion Nebula (M 42)
Authors:
Peter M. Weilbacher,
Ana Monreal-Ibero,
Wolfram Kollatschny,
Adam Ginsburg,
Anna F. McLeod,
Sebastian Kamann,
Christer Sandin,
Ralf Palsa,
Lutz Wisotzki,
Roland Bacon,
Fernando Selman,
Jarle Brinchmann,
Joseph Caruana,
Andreas Kelz,
Thomas Martinsson,
Arlette Pécontal-Rousset,
Johan Richard,
Martin Wendt
Abstract:
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…
▽ More
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.
△ Less
Submitted 30 June, 2015;
originally announced July 2015.
-
The influence of diffuse scattered light II. Observations of galaxy haloes and thick discs and hosts of BCGs
Authors:
Christer Sandin
Abstract:
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 ligh…
▽ More
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érsic-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.
△ Less
Submitted 25 February, 2015;
originally announced February 2015.
-
Towards DIB mapping in galaxies beyond 100 Mpc. A radial profile of the $λ$5780.5 diffuse interstellar band in AM 1353-272 B
Authors:
Ana Monreal-Ibero,
Peter M. Weilbacher,
Martin Wendt,
Fernando Selman,
Rosine Lallement,
Jarle Brinchmann,
Sebastian Kamann,
Christer Sandin
Abstract:
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 capabil…
▽ More
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 $λ$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.
△ Less
Submitted 26 February, 2015; v1 submitted 23 February, 2015;
originally announced February 2015.
-
The Chandra Planetary Nebulae Survey (ChanPlaNS): III. X-ray Emission from the Central Stars of Planetary Nebulae
Authors:
R. Montez Jr.,
J. H. Kastner,
B. Balick,
E. Behar,
E. Blackman,
V. Bujarrabal,
Y. -H. Chu,
R. L. M. Corradi,
O. De Marco,
A. Frank,
M. Freeman,
D. J. Frew,
M. A. Guerrero,
D. Jones,
J. A. Lopez,
B. Miszalski,
J. Nordhaus,
Q. A. Parker,
R. Sahai,
C. Sandin,
D. Schonberner,
N. Soker,
J. L. Sokoloski,
M. Steffen,
J. A. Toalá
, et al. (3 additional authors not shown)
Abstract:
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 "ha…
▽ More
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.
△ Less
Submitted 8 December, 2014;
originally announced December 2014.
-
Weak magnetic fields in central stars of planetary nebulae?
Authors:
M. Steffen,
S. Hubrig,
H. Todt,
M. Schöller,
W. -R. Hamann,
C. Sandin,
D. Schönberner
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 September, 2014; v1 submitted 13 August, 2014;
originally announced August 2014.
-
The Chandra Planetary Nebula Survey (ChanPlaNS). II. X-ray Emission from Compact Planetary Nebulae
Authors:
M. Freeman,
R. Montez Jr.,
J. H. Kastner,
B. Balick,
D. J. Frew,
D. Jones,
B. Miszalski,
R. Sahai,
E. Blackman,
Y. -H. Chu,
O. De Marco,
A. Frank,
M. A. Guerrero,
J. A. Lopez,
A. Zijlstra,
V. Bujarrabal,
R. L. M. Corradi,
J. Nordhaus,
Q. A. Parker,
C. Sandin,
D. Schönberner,
N. Soker,
J. L. Sokoloski,
M. Steffen,
J. A. Toalá
, et al. (2 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 August, 2014; v1 submitted 15 July, 2014;
originally announced July 2014.
-
The influence of diffuse scattered light I. The PSF and its role to observations of the edge-on galaxy NGC 5907
Authors:
Christer Sandin
Abstract:
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 measure…
▽ More
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.
△ Less
Submitted 20 June, 2014;
originally announced June 2014.
-
The central dynamics of M3, M13, and M92: Stringent limits on the masses of intermediate-mass black holes
Authors:
Sebastian Kamann,
Lutz Wisotzki,
Martin M. Roth,
Joris Gerssen,
Tim-Oliver Husser,
Christer Sandin,
Peter Weilbacher
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 April, 2014;
originally announced April 2014.
-
A hydrodynamical study of multiple-shell planetary nebulae. III. Expansion properties and internal kinematics: Theory versus observation
Authors:
D. Schoenberner,
R. Jacob,
H. Lehmann,
G. Hildebrandt,
M. Steffen,
A. Zwanzig,
C. Sandin,
R. L. M. Corradi
Abstract:
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 evoluti…
▽ More
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.
△ Less
Submitted 1 April, 2014;
originally announced April 2014.
-
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.
-
Properties of the three-dimensional structure in the central region of the supernova remnant SNR 0540-69.3
Authors:
C. Sandin,
P. Lundqvist,
N. Lundqvist,
C. -I. Björnsson,
G. Olofsson,
Yu. A. Shibanov
Abstract:
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 expansi…
▽ More
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)
△ Less
Submitted 15 April, 2013;
originally announced April 2013.
-
Bring out your codes! Bring out your codes! (Increasing Software Visibility and Re-use)
Authors:
Alice Allen,
Bruce Berriman,
Robert Brunner,
Dan Burger,
Kimberly DuPrie,
Robert J. Hanisch,
Robert Mann,
Jessica Mink,
Christer Sandin,
Keith Shortridge,
Peter Teuben
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 9 December, 2012;
originally announced December 2012.
-
PyCosmic: a robust method to detect cosmics in CALIFA and other fiber-fed integral-field spectroscopy datasets
Authors:
B. Husemann,
S. Kamann,
C. Sandin,
S. F. Sánchez,
R. García-Benito,
D. Mast
Abstract:
[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 f…
▽ More
[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.
△ Less
Submitted 1 November, 2012; v1 submitted 8 August, 2012;
originally announced August 2012.
-
The ERA2 facility: towards application of a fiber-based astronomical spectrograph for imaging spectroscopy in life sciences
Authors:
Martin M. Roth,
Karl Zenichowski,
Nicolae Tarcea,
Jürgen Popp,
Silvia Adelhelm,
Marvin Stolz,
Andreas Kelz,
Christer Sandin,
Svend-Marian Bauer,
Thomas Fechner,
Thomas Jahn,
Emil Popow,
Bernhard Roth,
Paul Singh,
Mudit Srivastava,
Dieter Wolter
Abstract:
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 facilit…
▽ More
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.
△ Less
Submitted 17 July, 2012;
originally announced July 2012.
-
The Chandra X-ray Survey of Planetary Nebulae (ChanPlaNS): Probing Binarity, Magnetic Fields, and Wind Collisions
Authors:
J. H. Kastner,
R. Montez Jr,
B. Balick,
D. J. Frew,
B. Miszalski,
R. Sahai,
E. Blackman,
Y. -H. Chu,
O. De Marco,
A. Frank,
M. A. Guerrero,
J. A. Lopez,
V. Rapson,
A. Zijlstra,
E. Behar,
V. Bujarrabal,
R. L. M. Corradi,
J. Nordhaus,
Q. Parker,
C. Sandin,
D. Schönberner,
N. Soker,
J. L. Sokoloski,
M. Steffen,
T. Ueta
, et al. (1 additional authors not shown)
Abstract:
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-luminou…
▽ More
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.
△ Less
Submitted 31 May, 2012; v1 submitted 25 April, 2012;
originally announced April 2012.
-
Modeling the diffuse X-ray emission of Planetary Nebulae with different chemical composition
Authors:
M. Steffen,
C. Sandin,
R. Jacob,
D. Schoenberner
Abstract:
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 circums…
▽ More
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.
△ Less
Submitted 28 September, 2011;
originally announced September 2011.
-
Visualization of Data from Integral Field Spectroscopy and the P3d Tool
Authors:
Martin M. Roth,
Christer Sandin
Abstract:
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.
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.
△ Less
Submitted 29 November, 2010;
originally announced November 2010.
-
The evolution of planetary nebulae VII. Modelling planetary nebulae of distant stellar systems
Authors:
D. Schönberner,
R. Jacob,
C. Sandin,
M. Steffen
Abstract:
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…
▽ More
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.
△ Less
Submitted 6 September, 2010;
originally announced September 2010.
-
p3d: a general data-reduction tool for fiber-fed integral-field spectrographs
Authors:
C. Sandin,
T. Becker,
M. M. Roth,
J. Gerssen,
A. Monreal-Ibero,
P. Böhm,
P. Weilbacher
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-r…
▽ More
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
△ Less
Submitted 23 February, 2010;
originally announced February 2010.
-
The evolution of planetary nebulae VI. On the chemical composition of the metal-poor PN G135.9+55.9
Authors:
C. Sandin,
R. Jacob,
D. Schönberner,
M. Steffen,
M. M. Roth
Abstract:
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 al…
▽ More
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.
△ Less
Submitted 30 December, 2009;
originally announced December 2009.
-
Measuring The Mass Loss Evolution at The Tip of The Asymptotic Giant Branch
Authors:
C. Sandin,
M. M. Roth,
D. Schönberner
Abstract:
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 w…
▽ More
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.
△ Less
Submitted 3 December, 2009; v1 submitted 17 August, 2009;
originally announced August 2009.
-
Spatially resolved spectroscopy of planetary nebulae and their halos I. Five galactic disk objects
Authors:
Christer Sandin,
Detlef Schönberner,
Martin M. Roth,
Matthias Steffen,
Petra Böhm,
Ana Monreal-Ibero
Abstract:
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, b…
▽ More
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)
△ Less
Submitted 23 May, 2008; v1 submitted 26 February, 2008;
originally announced February 2008.
-
Three-component modeling of C-rich AGB star winds IV. Revised interpretation with improved numerical descriptions
Authors:
Christer Sandin
Abstract:
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 structu…
▽ More
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.
△ Less
Submitted 2 November, 2007;
originally announced November 2007.
-
The evolution of planetary nebulae IV. On the physics of the luminosity function
Authors:
D. Schoenberner,
R. Jacob,
M. Steffen,
C. Sandin
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 31 August, 2007;
originally announced August 2007.
-
New observations of the halo radial temperature structure in NGC7662
Authors:
C. Sandin,
D. Schoenberner,
M. M. Roth,
M. Steffen,
A. Monreal-Ibero,
P. Boehm,
U. Tripphahn
Abstract:
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 mod…
▽ More
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.
△ Less
Submitted 11 September, 2006;
originally announced September 2006.
-
3D-Spectroscopy of extragalactic planetary nebulae as diagnostic probes for galaxy evolution
Authors:
Andreas Kelz,
Ana Monreal-Ibero,
Martin M. Roth,
Christer Sandin,
Detlef Schoenberner,
Matthias Steffen
Abstract:
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 m…
▽ More
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.
△ Less
Submitted 21 July, 2006;
originally announced July 2006.
-
The problem of small and large scale fields in the solar dynamo
Authors:
A. Brandenburg,
N. E. L. Haugen,
P. J. Käpylä,
C. Sandin
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 December, 2004;
originally announced December 2004.
-
Helical coronal ejections and their role in the solar cycle
Authors:
Axel Brandenburg,
Christer Sandin,
Petri J. Käpylä
Abstract:
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 si…
▽ More
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.
△ Less
Submitted 29 July, 2004;
originally announced July 2004.
-
Catastrophic alpha quenching alleviated by helicity flux and shear
Authors:
Axel Brandenburg,
Christer Sandin
Abstract:
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 e…
▽ More
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 $α$ 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.
△ Less
Submitted 29 July, 2004; v1 submitted 14 January, 2004;
originally announced January 2004.
-
Three-component modeling of C-rich AGB star winds III. Micro-physics of drift-dependent dust formation
Authors:
C. Sandin,
S. Hoefner
Abstract:
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 th…
▽ More
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.
△ Less
Submitted 30 September, 2003;
originally announced September 2003.
-
Three-component modeling of C-rich AGB star winds. II. The effects of drift in long-period variables
Authors:
C. Sandin,
S. Hoefner
Abstract:
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 thr…
▽ More
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.
△ Less
Submitted 15 April, 2003;
originally announced April 2003.
-
Three-component modeling of C-rich AGB star winds I. Method and first results
Authors:
C. Sandin,
S. Hoefner
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 November, 2002;
originally announced November 2002.