next up previous
Next: About this document

Radiation gasdynamics of planetary nebulae -- VI. The evolution of aspherical planetary nebulae

Garrelt Mellema

Dept. of Mathematics, UMIST, P.O. Box 88, Manchester M60 1QD, UK

This paper reports the results of the numerical study of the formation of aspherical Planetary Nebulae through the Generalised Interacting Winds Model, taking into account the effects caused by the evolving central star and fast wind. The results show for the first time that aspherical nebulae do form within the required time scale. Consideration of the development of the nebula shows that in the early stages it is the ionization of the aspherical AGB wind that contributes considerably to the shaping of the nebula. Furthermore the passing through of the ionization front may modify the density distribution in the slow wind, leading to the formation of a surrounding envelope, and sometimes a different morphology for the nebula than was to be expected from the initial conditions. I consider how the different phases of ionization fronts and wind swept bubbles can be observationally distinguished.

Accepted for publication in MNRAS Preprints can be obtained by

  1. contacting gam@ast.ma.umist.ac.uk
  2. via WWW on http://saturn.ma.umist.ac.uk:8000/Preprints2.html
  3. via anonnymous ftp on ftp://axp0.ast.man.ac.uk/outgoing/gm/mn_pnVI.v2f.ps.Z (1.4Mb)

The Outflow Structure of NGC 6537

L. Cuesta, J. P. Phillips and A. Mampaso

Instituto de Astrofísica de Canarias, E-38200 La Laguna, Tenerife, Spain

We present high and low resolution spectroscopy and narrow-band imaging of the unusual bipolar outflow source NGC 6537. As a consequence, we are able for the first time to construct large scale velocity maps of this source; note the presence of strongly variable core outflow densities, rising to a peak value ; find evidence for substantial core ionization stratification; and demonstrate the presence of a spatially extended, intermediate velocity wind extending over a range , together with a more tightly constrained component having . This latter wind appears co-spatial with the slower, higher intensity core emission, and may be responsible for driving denser clumps of material to a typical velocity v.

Temperatures appear to vary significantly over the core, from a peak value close to the the nucleus, to for projected radial offsets of , whilst extinction is more or less constant, yielding a mean estimate .

Application of a shock outflow model suggests that most of these features are readily explained providing wind mass-loss rates are of order --, and the flow is strongly collimated by an interior disk. In particular, we are able to provide excellent simulations for both the complex kinematics, and unusual bipolar.

Accepted for publication in Astronomy and Astrophysics For preprints, contact lcc@iac.es

Polarization of Astronomical Maser Radiation. III. Arbitrary Zeeman Splitting and Anisotropic Pumping

Moshe Elitzur

Dept. of Physics & Astronomy, University of Kentucky, Lexington, KY 40506, USA

=

= =

General solutions of the maser polarization problem are presented for arbitrary absorption coefficients. The results are used to calculate polarization for masers permeated by magnetic fields with arbitrary values of , the ratio of Zeeman splitting to Doppler linewidth, and for anisotropic (m-dependent) pumping. In the case of magnetic fields, one solution describes the polarization for overlapping Zeeman components, . The limit of this solution reproduces the linear polarization derived in previous studies, which were always conducted at this unphysical limit. Terms of higher order in have a negligible effect on the magnitude of q. However, these terms produce some major new results: (1) The solution is realized only when the Zeeman splitting is sufficiently large that , where is the source function and is the saturation intensity (pumping schemes typically have ). When this condition is met, the linear polarization requires , where J is the angle-averaged intensity. This condition generally requires considerable amplification, but is met long before saturation (). (2) The linear polarization is accompanied by circular polarization, proportional to . Because is proportional to the transition wavelength, the circular polarization of SiO masers should decrease with rotation quantum number, as observed. In the absence of theory for , previous estimates of magnetic fields from detected maser circular polarization had to rely on conjectures in this case and generally need to be revised downward. The fields in SiO masers are 2--10 G and were overestimated by a factor of 8. The OH maser regions around supergiants have fields of 0.1--0.5 mG, which were overestimated by factors of 10--100. The fields were properly estimated for OH/IR masers ( 0.1 mG) and HO masers in star-forming regions ( 15--50 mG). (3) Spurious solutions that required stability analysis for their removal in all previous studies are never reproduced here; in particular, there are no stationary physical solutions for propagation at , where is the angle from the direction of the magnetic field, so such radiation is unpolarized. These spurious solutions can be identified as the = 0 limits of non-physical solutions and they never arise at finite values of , however small. (4) Allowed values of are limited by bounds that depend both on Zeeman splitting and frequency shift from line center. At , the allowed phase space region encompasses essentially all frequencies and . As the field strength increases, the allowed angular region shrinks at a frequency-dependent rate, leading to contraction of the allowed spectral region. This can result in narrow maser features with linewidths smaller than the Doppler width and substantial circular polarization in sources with . When , all frequencies and directions are prohibited for the stationary solution and the radiation is unpolarized.

Another solution describes the polarization when the Zeeman components separate. This occurs at line center when and at one Doppler width when . The solution is identical to that previously identified in the limit, and applies to OH masers around H II regions. A significant new result involves the substantial differences between the - and -components for most propagation directions, differences that persist into the saturated domain. Overall, H II /OH regions should display a preponderance of -components. The absence of any -components in W3(OH) finds a simple explanation as maser action in a magnetic field aligned within to the line of sight.

Accepted by The Astrophysical Journal (January 20, 1996). For preprints, contact moshe@pa.uky.edu

SiO maser survey of the galactic bulge IRAS sources:
IV. Observationnal properties of SiO masers

B. W. Jiang, S. Deguchi, H. Izumiura, Y. Nakada, I. Yamamura

Department of Astronomical Sciences, The Graduate University for Advanced Studies,
Nobeyama Radio Observatory, National Astronomical Observatory
Minamimaki, Minamisaku, Nagano 384-13, Japan
On leave from Beijing Astronomical Observatory, China
Nobeyama Radio Observatory, National Astronomical Observatory
Minamimaki, Minamisaku, Nagano 384-13, Japan
Department of Astronomy and Earth Sciences, Tokyo Gakugei University
4-1-1 Nukui-kita, Koganei, Tokyo 184, Japan
Kiso Observatory, Faculty of Science, University of Tokyo Mitake, Kiso, Nagano 397-01, Japan
Kiso Observatory, Faculty of Science, University of Tokyo Mitake, Kiso, Nagano 397-01, Japan

Statistical properties of SiO masers are studied based on the observations of color selected IRAS sources in the direction of the Galactic bulge. Among 313 observed late-type stellar objects in the region of and 3, 193 were detected in the SiO v=1, J=1-0 line and/or the v=2, J=1-0 line. We discuss the relation between SiO detection rate and infrared properties for these sources. The sources are divided into bulge and disk members in terms of the luminosity distance. The SiO detection rate is found to be appreciably higher for the bulge sources than for the disk sources, and that the detection rate is higher for variable IRAS sources (with IRAS variability index above 90) than that for the nonvariable sources (with the index below 90). We also find that the center velocities of the v=1 and v=2, J=1-0 lines are almost equal to the average velocity of OH 1612 MHz double peaks. This implies that the SiO velocity is a good indicator of the stellar radial velocity, approving the analysis of bulge stellar kinematics using SiO radial velocities. On average, the line-intensity ratio of the v=1 to v=2 lines is approximately unity and there seems no systematic difference on the ratio between disk and bulge sources. The SiO line intensities are very weakly correlated with IRAS fluxes for the bulge sources. No clear correlation is found between the SiO line width and the stellar expansion velocity. This indicates that the SiO line width cannot be used as a measure of luminosity of the infrared sources.

Accepted for publication in Publ. Astron. Soc. Japan., Dec. 25, 1995 issue For preprints, contact deguchi@nro.nao.ac.jp

Modelling of the Dust and Gas Outflows from OH26.5+0.6
: The Superwind

K. Justtanont, C.J. Skinner, A.G.G.M. Tielens, M. Meixner, and F. Baas

NASA Ames Research Center, MS 245-3, Moffett Field, CA 94035

Institut d'Astrophysique de Paris, 98 bis Bd. Arago, 75014 Paris, France

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218

Institute of Geophysics and Planetary Physics, Lawrence Livermore National Laboratory, PO Box 808, CA 94581-9900

Astronomy Dept, MC 221, 1002 W. Green St., University of Illinois, Urban, IL 61801

Leiden Observatory, PO Box 9513, NL-2300, RA Leiden, The Netherlands

Joint Astronomy Center, 660 N. A'ohoku Place, University Park, Hilo, HI 96720

We have observed the extreme OH/IR star, OH26.5+0.6, in the infrared dust continuum and in the submillimeter rotational lines of CO. Mid-infrared images reveal the compact nature of the circumstellar shell (). A deep 9.7m absorption feature and an absorption at 18m show that the dust mass loss rate is very high. However, the low antenna temperatures of CO J=1-0 and 2-1 lines suggest that the outer part of the circumstellar shell is much more tenuous. In order to resolve this discrepancy, we have observed the J=3-2 and 4-3 CO rotational transitions.

We have developed a model for the circumstellar shell for OH26.5+0.6 which is consistent with the infrared and submillimeter observations. The dust and gas data are well fitted by a two-shell model, consisting of a dense shell surrounded by a more tenuous shell. The former, we identify with the superwind ( M/yr) and the latter with mass loss on the AGB ( M/yr). The transition between the two mass loss phases is shown to be rather abrupt ( years). Depending on the mass of the progenitor, this superwind phase may be the last thermal pulse (for M 1.5 M), or the first of a series of the superwind phases (for up to 8M), punctuated by a period of low mass loss rates, before the star evolves off the AGB.

Accepted by Astrophysical Journal For preprints, contact kj@dusty.arc.nasa.gov

Hubble Space Telescope Observations of Planetary Nebulae in the Magellanic Clouds IV: [O III] Images and Evolutionary Ages

M.A. Dopita, E. Vassiliadis, S.J. Meatheringham, R.C. Bohlin, H.C. Ford, J.P. Harrington, P.R. Wood, T.P. Stecher and S.P Maran Mt. Stromlo and Siding Spring Observatories, Institute of Advanced Studies, The Australian National University Private Bag, Weston Creek P.O., ACT 2611, Australia.
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA.
Johns Hopkins University, Department of Physics & Astronomy, Charles & 34th Sts., Baltimore, MD 21218, USA.
Dept. of Astronomy, The University of Maryland, College Park, MD 20742-2421, USA.
NASA Goddard Space Flight Center, Code 680, Greenbelt, MD 20771, USA.

The Planetary Camera of the Hubble Space Telescope has been used to obtain images in the [O III] 5007 Å line for a sample of 15 planetary nebulae in the Magellanic Clouds selected to eliminate any selection bias in either excitation class or in flux (for erg cm s). These images are used to derive the physical dimensions, the spatial structure and the kinematic ages of the nebulae. The raw images were deconvolved using the Richardson-Lucy Image Restoration Algorithm, and based on extensive tests of model images, a limit of 100 iterations of the algorithm was adopted. The restored images show clear evidence for size evolution across the H-R Diagram. The younger, low excitation, compact planetary nebulae tend to be systematically smaller than photoionisation models based on ground-based data would predict, suggesting these planetary nebulae have a central reservoir of dense atomic and molecular gas. This gas lies close to the central star and is undergoing ionisation and being accelerated into outflow. Planetary Nebulae previously classified as nitrogen-rich objects with massive central stars (Peimbert Type I), show the bipolar "butterfly" symmetry that is also a characteristic of their Galactic counterparts. The derived kinematic ages range from less than 1000 years up to almost 5000 years, but show little sign of systematic increase along the evolutionary tracks. The true ages of the larger objects are systematically underestimated because of acceleration of the nebular shell during its lifetime. Using both the empirical fit that we had previously derived for the expansion velocity as a function of the position on the H-R Diagram, and the theoretical evolutionary tracks of the central star, we have derived two semi-empirical estimates for the evolutionary timescales based upon the nebular size, and the measured dynamical age. If these evolutionary timescales are to be consistent with the evolutionary age derived from theory, then He-burners outnumber H-burners in the approximate ratio 2:1.

Accepted by Astrophysical Journal. For preprints, contact Michael.Dopita@anu.edu.au

The Formation of Barium and CH Stars and related objects

Zhanwen Han, Peter P. Eggleton, Philipp Podsiadlowski and Christopher A. Tout

Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK

We investigate the formation of barium (Ba) stars and CH stars via binary interactions. We consider four evolutionary channels for their formation: wind accretion, wind exposure, stable Roche lobe overflow (RLOF), and common-envelope (CE) ejection, and carry out Monte-Carlo simulations. We also explore the effects of model parameters, such as a tidal-enhancement parameter for stellar wind, a maximum stellar mass for s-processing, and a minimum core mass for thermal pulsation, on our results. We also explore the effects of various assumptions about age, mass-ratio distribution and wind velocity.

Our results show that binary interaction is successful in explaining the formation of Ba and CH stars. We successfully explain the distribution of orbital periods, the distribution of mass functions and the number of observed Ba stars. We also support the views (a) that a tidally enhanced stellar wind exists, (b) that thermal pulsation may begin at an earlier phase during the asymptotic giant branch (AGB) than is usually assumed, and (c) that the maximum stellar mass for s-processing is about . We also find that the degree of Ba pollution is strongly correlated with orbital period. The average mass of strong Ba stars is estimated to be , and of all Ba stars ; the average mass for the white dwarf (WD) companions in the Ba-star binaries is estimated to be . The average mass of strong CH stars is estimated to be , and the average mass for their WD companions is estimated to be . The total number of Ba stars in the Galaxy brighter than 10th apparent magnitude is estimated to be 6000, while that of CH stars is much less.

We also investigate the formation of pre-Ba/CH stars, CVs, Algols, DDs and symbiotics. Some encouraging results are found, although much work needs to be done in order to understand fully their formation, especially of pre-Ba/CH stars and of symbiotics.

Accepted by MNRAS For preprints, contact zhanwen@ast.cam.ac.uk Manuscript is available by anonymous ftp to ftp.ast.cam.ac.uk or 131.111.69.186, in directory pub/zhanwen.

Post Doctoral Research Associate, University of Washington

Applications are invited for an observational or theoretical postdoctoral position at the University of Washington in Seattle in the research areas of nebular astrophysics, nebular gas dynamics (with special emphasis on the effects of stellar winds), stellar mass-loss processes in the nuclei of planetary nebulae and their progenitors, or closely related fields. The start date is flexible. The starting salary is $3,031 per month.

Postdocs at UW have access to the 3.5m telescope at Apache Point Observatory (equipped with a double mid-resolution spectrograph, an Echelle, and a 2-5-micron IR imager/grism system) and the 0.8m telescope at Manastash Ridge Observatory. The two-year position may be extended to three years if additional funding is secured.

Applicants must have a PhD in astronomy. Send a 2-page statement of relevant experience, a 1-page statement of professional goals for the position, a curriculum vitae and bibliography, copies of two recent publications, and a letter stating optimum range of start dates and the names and addresses (with phone and email) of three people who can provide knowledgeable professional recommendations, all in one mailing, before 15 October 1995.

The University of Washington is building a culturally diverse faculty and strongly encourages applications from women and minority candidates. AAE/EOE.

UNIVERSITY OF WASHINGTON
Attention: Prof. Bruce Balick
Astronomy Dept, Box 351580
Seattle WA 98195-1580
email: balick@astro.washington.edu





next up previous
Next: About this document



Peter te Lintel Hekkert
Sat Nov 25 20:55:00 EST 1995
Staff space
Public