Spatio-kinematical Model Research Articles

A dynamic view of V Hydrae

Context. The well studied carbon star V Hydrae is known to exhibit a complex asymmetric environment made of a dense equatorial wind and high-velocity outflows, hinting at its transition from the AGB phase to the asymmetric planetary nebula phase. In addition, V Hydrae also exhibits a long secondary period of 17 yr in its light curve, suggesting the presence of a binary companion that could shape the circumstellar environment. Aims. In this paper, we aim to confirm the binary nature of V Hydrae by deriving its orbital parameters and investigating the effect of the orbital motion on the circumbinary environment. Methods. In a first step, we used a radial-velocity monitoring performed with the HERMES spectrograph to disentangle the pulsation signal of the AGB from its orbital motion and to obtain the spectroscopic orbit. We combined the spectroscopic results with astrometric information to get the complete set of orbital parameters, including the system inclination. Next, we reported the time variations of the sodium and potassium resonance doublets. Finally, following the methods used for post-AGB stars, we carried out spatio-kinematic modelling of a conical jet to reproduce the observed spectral-line modulation. Results. We found the orbital solution of V Hydrae for a period of 17 yr. We correlated the companion passage across the line of sight with the obscuration event and the blue-shifted absorption of alkaline resonant lines. Those variations were modelled by a conical jet emitted from the companion, whose opening angle is wide and whose sky-projected orientation is found to be consistent with the axis of the large-scale bipolar outflow previously detected in the radio-emission lines of CO. Conclusions. We show that the periodic variation seen for V Hydrae is likely to be due to orbital motion. The presence of a conical jet offers a coherent model to explain the various features of V Hydrae environment.

The structure of jets launched from post-AGB binary systems

Context. In this paper, we focus on post-asymptotic giant branch (post-AGB) binaries and study the interaction between the different components of these complex systems. These components comprise the post-AGB primary, a main sequence secondary, a circumbinary disk, as well as a fast bipolar outflow (jet) launched by the companion. We obtained well-sampled time series of high resolution optical spectra over the last decade and these spectra provide the basis of our study. Aims. We aim to use the time-series data to quantify the velocity and density structure of the jets in nine of these post-AGB binaries. This complements our earlier work and this amounts to the analyses of 16 jet-launching systems in total. Methods. The jet is detected in absorption, at superior conjunction, when the line of sight towards the primary goes through the bipolar cone. Our spectral time series scan the jets during orbital motion. Our spatio-kinematic model is constrained by these dynamical spectra. We complement this with a radiative-transfer model in which the Balmer series are used to derive total mass-loss rates in the jets. Results. The jets are found to be wide (> 30°) and display an angle-dependent density structure with a dense and slower outer region near the jet cone and a fast inner part along the jet symmetry axes. The deprojected outflow velocities confirm that the companions are main sequence companions. The total mass-loss rates are large (10−8 − 10−5 M⊙ yr−1), from which we can infer that the mass-accretion rates onto the companion star must be high as well. The circumbinary disk is likely the main source for the accretion disk around the companion. All systems with full disks that start near the sublimation radius show jets, whereas for systems with evolved transition disks this lowers to a detection rate of 50%. Objects without an infrared excess do not show jets. Conclusions. We conclude that jet creation in post-AGB binaries is a mainstream process. Our geometric spatio-kinematic model is versatile enough to model the variety of spectral time series. The interaction between the circumbinary disks and the central binary provide the needed accretion flow, but the presence of a circumbinary disk does not seem to be the only prerequisite to launch a jet.

Spatiokinematic models of five nova remnants: correlations between nova shell axial ratio, expansion velocity, and speed class

ABSTRACT We present long-slit intermediate-dispersion spectroscopic observations and narrow-band direct imaging of four classical nova shells, namely T Aur, HR Del, DQ Her, and QU Vul, and the nova-like source CK Vul. These are used to construct models of their nebular remnants using the morphokinematic modelling tool shape to reveal their 3D shape. All these nova remnants but CK Vul can be described by prolate ellipsoidal shells with different eccentricity degree, from the spherical QU Vul to the highly elongated shell with an equatorial component HR Del. On the other hand, CK Vul shows a more complex structure, with two pairs of nested bipolar lobes. The spatiokinematic properties of the ellipsoidal nova shells derived from our models include their true axial ratios. This parameter is expected to correlate with the expansion velocity and decline time t3 (i.e. their speed class) of a nova as the result the interaction of the ejecta with the circumstellar material and rotation speed and magnetic field of the white dwarf. We have compared these three parameters including data available in the literature for another two nova shells, V533 Her and FH Ser. There is an anticorrelation between the expansion velocity and the axial ratio and decline time t3 for nova remnants with ellipsoidal morphology, and a correlation between their axial ratios and decline times t3, confirming theoretical expectations that the fastest expanding novae have the smallest axial ratios. We note that the high expansion velocity of the nova shell HR Del of 615 km s−1 is inconsistent with its long decline time t3 of 250 d.

3D spatio-kinematic modelling of Abell 48, a planetary nebula around a Wolf–Rayet [WN] star

ABSTRACT The planetary nebula (PN) Abell 48 (PN G029.0+00.4) is around a rare Wolf–Rayet [WN5] star whose stellar history is as yet unknown. Using the integral field observations of the H αλ6563 and [N ii] λ6584 line emissions, we conducted a comprehensive spatio-kinematic analysis of this PN. A three-dimensional spatio-kinematic ionization model was developed with the kinematic modelling tool shape to replicate the observed spatially resolved velocity channels and position–velocity diagrams. According to our kinematic analysis of the H α emission, this object possesses a deformed elliptic toroidal shell with an outer radius of 23 arcsec and a thickness of 15 arcsec associated with an integrated H α emission-line expansion of ∼35 ± 5 km s−1, a maximum poloidal expansion of around 70 ± 20 km s−1 at an inclination angle of ∼30° with respect to the line of sight, and a position angle of ∼130° measured from east towards north in the equatorial coordinate system. Furthermore, [N ii] kinematic modelling reveals the presence of narrow (∼3 arcsec) exterior low-ionization structures surrounding the main elliptical shell, which could have formed as a result of shock collisions with the interstellar medium. The torus-shaped morphology of this PN could be related to its unusual hydrogen-deficient [WN] nucleus that needs to be inspected further.

The ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES). III. A Young Molecular Outflow Driven by a Decelerating Jet

Abstract We present a spatio-kinematical analysis of the CO (J = 2 → 1) line emission, observed with the Atacama Large Millimeter/submillimeter Array (ALMA), of the outflow associated with the most massive core, ALMA1, in the 70 μm dark clump G010.991–00.082. The position–velocity (PV) diagram of the molecular outflow exhibits a peculiar S-shaped morphology that has not been seen in any other star-forming region. We propose a spatio-kinematical model for the bipolar molecular outflow that consists of a decelerating high-velocity component surrounded by a slower component whose velocity increases with distance from the central source. The physical interpretation of the model is in terms of a jet that decelerates as it entrains material from the ambient medium, which has been predicted by calculations and numerical simulations of molecular outflows in the past. One side of the outflow is shorter and shows a stronger deceleration, suggesting that the medium through which the jet moves is significantly inhomogeneous. The age of the outflow is estimated to be τ ≈ 1300 yr, after correction for a mean inclination of the system of ≈57°.

The post-common-envelope binary central star of the planetary nebula ETHOS 1

ABSTRACT We present a detailed study of the binary central star of the planetary nebula ETHOS 1 (PN G068.1+11.0). Simultaneous modelling of light and radial velocity curves reveals the binary to comprise a hot and massive pre-white dwarf with an M-type main-sequence companion. A good fit to the observations was found with a companion that follows expected mass–temperature–radius relationships for low-mass stars, indicating that despite being highly irradiated, it is consistent with not being significantly hotter or larger than a typical star of the same mass. Previous modelling indicated that ETHOS 1 may comprise the first case where the orbital plane of the central binary does not lie perpendicular to the nebular symmetry axis, at odds with the expectation that the common envelope is ejected in the orbital plane. We find no evidence for such a discrepancy, deriving a binary inclination in agreement with that of the nebula as determined by spatio-kinematic modelling. This makes ETHOS 1 the ninth post-common-envelope planetary nebula in which the binary orbital and nebular symmetry axes have been shown to be aligned, with as yet no known counter-examples. The probability of finding such a correlation by chance is now less than 0.000 02 per cent.

Determining mass-accretion and jet mass-loss rates in post-asymptotic giant branch binary systems

Aims. In this study we determine the morphology and mass-loss rate of jets emanating from the companion in post-asymptotic giant branch (post-AGB) binary stars with a circumbinary disc. In doing so we also determine the mass-accretion rates onto the companion, and investigate the source feeding the circum-companion accretion disc. Methods. We perform a spatio-kinematic modelling of the jet of two well-sampled post-AGB binaries, BD+46°442 and IRAS 19135+3937, by fitting the orbital phased time series of Hα spectra. Once the jet geometry, velocity, and scaled density structure are computed, we carry out radiative transfer modelling of the jet for the first four Balmer lines to determine the jet densities, thus allowing us to compute the jet mass-loss rates and mass-accretion rates. We distinguish the origin of the accretion by comparing the computed mass-accretion rates with theoretically estimated mass-loss rates, both from the post-AGB star and from the circumbinary disc. Results. The spatio-kinematic model of the jet reproduces the observed absorption feature in the Hα lines. The jets have an inner region with extremely low density in both objects. The jet model for BD+46°442 is tilted by 15° with respect to the orbital axis of the binary system. IRAS 19135+3937 has a smaller tilt of 6°. Using our radiative transfer model, we find the full 3D density structure of both jets. By combining these results, we can compute the mass-loss rates of the jets, which are of the order of 10−7 − 10−5 M⊙ yr−1. From this we estimate mass-accretion rates onto the companion of 10−6 − 10−4 M⊙ yr−1. Conclusions. Based on the mass-accretion rates found for these two objects, we conclude that the circumbinary disc is most likely the source feeding the circum-companion accretion disc. This is in agreement with the observed depletion patterns in post-AGB binaries, which is caused by re-accretion of gas from the circumbinary disc that is under-abundant in refractory elements. The high accretion rates from the circumbinary disc imply that the lifetime of the disc will be short. Mass transfer from the post-AGB star cannot be excluded in these systems, but it is unlikely to provide a sufficient mass-transfer rate to sustain the observed jet mass-loss rates.

A spatio-kinematic model for jets in post-AGB stars,

Aims. We aim to determine the geometry, density gradient, and velocity structure of jets in post-asymptotic giant branch (post-AGB) binaries. Methods. Our high cadence time series of high-resolution optical spectra of jet-creating post-AGB binary systems provide us with a unique tomography of the jet. We determine the spatio-kinematic structure of the jets based on these data by fitting the synthetic spectral line profiles created by our model to the observed, orbital phase-resolved, Hα-line profiles of these systems. The fitting routine is provided with an initial spectrum and is allowed to test three configurations, derived from three specific jet launching models: a stellar jet launched by the star, an X-wind, and a disk wind configuration. We apply a Markov-chain Monte Carlo routine in order to fit our model to the observations. Our fitting code is tested on the post-AGB binary IRAS 19135+3937. Results. We find that a model using the stellar jet configuration gives a marginally better fit to our observations. The jet has a wide half-opening angle of about 76° and reaches velocities up to 870 km s−1. Conclusions. Our methodology is successful in determining some parameters for jets in post-AGB binaries. The model for IRAS 19135+3937 includes a transparent, low density inner region (for a half-opening angle < 40°). The source feeding the accretion disk around the companion is most likely the circumbinary disk. We will apply this jet fitting routine to other jet-creating post-AGB stars in order to provide a more complete description of these objects.

Spatio-kinematical model of the collimated molecular outflow in the water-fountain nebula IRAS 16342–3814

Context. Water-fountain nebulae are asymptotic giant branch (AGB) and post-AGB objects that exhibit high-velocity outflows traced by water-maser emission. Their study is important for understanding the interaction between collimated jets and the circumstellar material that leads to the formation of bipolar and/or multi-polar morphologies in evolved stars. Aims. The aim of this paper is to describe the three-dimensional morphology and kinematics of the molecular gas of the water-fountain nebula IRAS 16342−3814. Methods. Data was retrieved from the ALMA archive for analysis using a simple spatio-kinematical model. The software SHAPE was employed to construct a three-dimensional, spatio-kinematical model of the molecular gas in IRAS 16342−3814, and to then reproduce the intensity distribution and position-velocity diagram of the CO emission from the ALMA observations to derive the morphology and velocity field of the gas. Data from CO(J = 1 → 0) supported the physical interpretation of the model. Results. A spatio-kinematical model that includes a high-velocity collimated outflow embedded within material expanding at relatively lower velocity reproduces the images and position-velocity diagrams from the observations. The derived morphology is in good agreement with previous results from IR and water-maser emission observations. The high-velocity collimated outflow exhibits deceleration across its length, while the velocity of the surrounding component increases with distance. The morphology of the emitting region, the velocity field, and the mass of the gas as function of velocity are in excellent agreement with the properties predicted for a molecular outflow driven by a jet. The timescale of the molecular outflow is estimated to be ~70–100 yr. The scalar momentum carried by the outflow is much larger than it can be provided by the radiation of the central star. An oscillating pattern was found associated with the high-velocity collimated outflow. The oscillation period of the pattern is T ≈ 60–90 yr and its opening angle is θop ≈ 2°. Conclusions. The CO (J = 3 → 2) emission in IRAS 16342−3814 is interpreted in terms of a jet-driven molecular outflow expanding along an elongated region. The position-velocity diagram and the mass spectrum reveal a feature due to entrained material that is associated with the driving jet. This feature is not seen in other more evolved objects that exhibit more developed bipolar morphologies. It is likely that the jet in those objects has already disappeared since it is expected to last only for a couple hundred years. This strengthens the idea that water fountain nebulae are undergoing a very short transition during which they develop the collimated outflows that shape the circumstellar envelopes. The oscillating pattern seen in the CO high-velocity outflow is interpreted as due to precession with a relatively small opening angle. The precession period is compatible with the period of the corkscrew pattern seen at IR wavelengths. We propose that the high-velocity molecular outflow traces the underlying primary jet that produces such a pattern.

A submillimeter polarization analysis of Frosty Leo

We present a polarimetric investigation of the protoplanetary nebula Frosty Leo performed with the Submillimeter Array. We were able to detect, in the low continuum level (peak at 14.4 mJy beam$^{-1}$), a marginal polarization at $\sim2.6\sigma$. The molecular line investigation based on the CO $J=3\rightarrow2$ emission shows a peak emission of 68.1 Jy beam$^{-1}$ km s$^{-1}$ and the polarization detection in this CO line is also marginal, with a peak at $\sim3.8\sigma$. In both cases, it was therefore not possible to use the electric vector maps (B-field) to accurately trace the magnetic field (B-field) within the PPN. The spatio-kinematic modelling realised with the different velocity channel maps indicates three main structures: a distorted torus accompanied by a bipolar outflow or jet aligned with its axis and a flattened spherical "cap". The comparison of the CO polarization segments with our model suggests that the polarized emission probably arises in the first two components.

High-velocity Bullets from V Hydrae, an Asymptotic Giant Branch Star in Transition: Ejection History and Spatio-kinematic Modeling

Abstract The carbon star V Hydrae (V Hya) provides new insight into the nature of the launching mechanism of jet-like outflows that are believed to be the cause of the poorly understood transition phase of asymptotic giant branch stars into aspherical planetary nebulae. V Hya has been shown to periodically eject collimated gas blobs at high velocities (“bullets”). By analyzing data from Hubble Space Telescope/Space Telescope Imaging Spectrograph 2D spectra, obtained at six epochs spaced over a decade that show four successively ejected bullets with a spacing of 8.5 yr, we have created kinematic models of the dynamical evolution of a specific bullet (#1) for the first three observed epochs (2002, 2003, 2004) using a 3D spatio-kinematic code, SHAPE. Using these models, we fit the observed morphology, line-of-sight velocity, proper motion, and intensity for the extended, gaseous bullet as a function of time over a period of 2 yr, in order to constrain its 3D movement and the evolution of its physical properties over this period. Our results suggest that although bullet #1’s motion is predominantly ballistic, there are small but significant changes in the position angle and inclination angle of the long (symmetry) axis of the bullet that tilt it progressively toward the symmetry axis of the bipolar molecular nebula around V Hya. In contrast, bullet #3 shows strong acceleration soon after ejection. We discuss the possibilities that bullet acceleration is caused by a nonradial magnetic field and/or by hydrodynamic interaction with the ambient gas through which the bullet is traveling.

The interaction of the halo around the butterfly planetary nebula NGC 650-1 with the interstellar medium

With its bright and wide equatorial waist seen almost edge-on (the butterfly body) and the faint and broad bipolar extensions (the butterfly wings), NGC650-1 is the archetypical example of bipolar planetary nebula (PN) with butterfly morphology. We present here deep high-resolution broad- and narrow-band optical images that expose the rich and intricate fine structure of this bipolar PN, with small-scale bubble-like features and collimated outflows. A SHAPE spatio-kinematical model indicates that NGC650-1 has a broad central torus with an inclination angle of 75 degrees with respect to the line of sight, whereas that of the bipolar lobes, which are clearly seen in the position-velocity maps, is 85 degrees. Large field of view deep images show, for first time, an arc-like diffuse envelope in low- and high-excitation emission lines located up to 180 arcsec towards the East-Southeast of the central star, well outside the main nebula. This morphological component is confirmed by Spitzer MIPS and WISE infrared imaging, as well as by long-slit low- and high-dispersion optical spectroscopic observations. HST images of NGC650-1 obtained at two different epochs (14 yrs apart), reveal the proper motion of the central star along this direction. We propose that this motion of the star through the interstellar medium compresses the remnant material of a slow Asymptotic Giant Branch wind, producing this bow-shock-like feature.

SALT HRS discovery of a long-period double-degenerate binary in the planetary nebula NGC 1360

Whether planetary nebulae (PNe) are predominantly the product of binary stellar evolution as some population synthesis models (PSM) suggest remains an open question. Around 50 short period binary central stars ($P\sim1$ d) are known, but with only four with measured orbital periods over 10 d, our knowledge is severely incomplete. Here we report on the first discovery from a systematic SALT HRS survey for long period binary central stars. We find a 142 d orbital period from radial velocities of the central star of NGC~1360, HIP~16566. NGC~1360 appears to be the product of common-envelope (CE) evolution, with nebula features similar to post-CE PNe, albeit with an orbital period considerably longer than expected to be typical of post-CE PSM. The most striking feature is a newly-identified ring of candidate low-ionisation structures (LIS). Previous spatio-kinematic modelling of the nebula gives a nebula inclination of $30\pm10$ deg, and assuming the binary nucleus is coplanar with the nebula, multi-wavelength observations best fit a more massive, evolved WD companion. A WD companion in a 142 d orbit is not the focus of many PSM, making NGC~1360 a valuable system with which to improve future PSM work. HIP~16566 is amongst many central stars in which large radial velocity variability was found by low-resolution surveys. The discovery of its binary nature may indicate long period binaries may be more common than PSM models predict.

ALMA high spatial resolution observations of the dense molecular region of NGC 6302.

The mechanism behind the shaping of bipolar planetary nebulae is still poorly understood. It is becoming increasingly clear that the main agents must operate at their innermost regions, where a significant equatorial density enhancement should be present and related to the collimation of light and jet launching from the central star preferentially towards the polar directions. Most of the material in this equatorial condensation must be lost during the asymptotic giant branch as stellar wind and later released from the surface of dust grains to the gas phase in molecular form. Accurately tracing the molecule-rich regions of these objects can give valuable insight into the ejection mechanisms themselves. We investigate the physical conditions, structure and velocity field of the dense molecular region of the planetary nebula NGC 6302 by means of ALMA band 7 interferometric maps. The high spatial resolution of the 12CO and 13CO J=3-2 ALMA data allows for an analysis of the geometry of the ejecta in unprecedented detail. We built a spatio-kinematical model of the molecular region with the software SHAPE and performed detailed non-LTE calculations of excitation and radiative transfer with the shapemol plug-in. We find that the molecular region consists of a massive ring out of which a system of fragments of lobe walls emerge and enclose the base of the lobes visible in the optical. The general properties of this region are in agreement with previous works, although the much greater spatial resolution of the data allows for a very detailed description. We confirm that the mass of the molecular region is 0.1 M⊙. Additionally, we report a previously undetected component at the nebular equator, an inner, younger ring inclined ~60° with respect to the main ring, showing a characteristic radius of 7.5×1016 cm, a mass of 2.7×10-3 M⊙, and a counterpart in optical images of the nebula. This inner ring has the same kinematical age as the northwest optical lobes, implying it was ejected approximately at the same time, hundreds of years after the ejection of the bulk of the molecular ring-like region. We discuss a sequence of events leading to the formation of the molecular and optical nebulae, and briefly speculate on the origin of this intriguing inner ring.

OBSERVATIONAL CONFIRMATION OF A LINK BETWEEN COMMON ENVELOPE BINARY INTERACTION AND PLANETARY NEBULA SHAPING

ABSTRACT A current issue in the study of planetary nebulae with close binary central stars (CSs) is the extent to which the binaries affect the shaping of the nebulae. Recent studies have begun to show a high coincidence rate between nebulae with large-scale axial or point symmetries and close binary stars. In addition, combined binary-star and spatiokinematic modeling of the nebulae have demonstrated that all of the systems studied to date appear to have their central binary axis aligned with the primary axis of the nebula. Here we add two more systems to the list, the CSs and nebulae of NGC 6337 and Sp 1. We show both systems to be low inclination, with their binary axis nearly aligned with our line of sight. Their inclinations match published values for the inclinations of their surrounding nebulae. Including these two systems with the existing sample statistically demonstrates a direct link between the central binary and the nebular morphology. In addition to the systems’ inclinations we give ranges for other orbital parameters from binary modeling, including updated orbital periods for the binary CSs of NGC 6337 and Sp 1.

SHAPEMOL: a 3D code for calculating CO line emission in planetary and protoplanetary nebulae

Modern instrumentation in radioastronomy constitutes a valuable tool for studying the Universe: ALMA has reached unprecedented sensitivities and spatial resolution, while Herschel/HIFI has opened a new window for probing molecular warm gas (~50-1000 K). On the other hand, the software SHAPE has emerged in the past few years as a standard tool for determining the morphology and velocity field of different kinds of gaseous emission nebulae via spatio-kinematical modelling. SHAPE implements radiative transfer solving, but it is only available for atomic species and not for molecules. Being aware of the growing importance of the development of tools for simplifying the analyses of molecular data, we introduce shapemol, a complement to SHAPE, with which we intend to fill the so-far under-developed molecular niche. shapemol enables user-friendly, spatio-kinematic modelling with accurate non-LTE calculations of excitation and radiative transfer in CO lines. It allows radiative transfer solving in the 12CO and 13CO J=1-0 to J=17-16 lines, but its implementation permits easily extending the code to different molecular species. shapemol allows easily generating synthetic maps and line profiles to match against interferometric or single-dish observations. We fully describe shapemol and discuss its limitations and the sources of uncertainty to be expected in the final synthetic profiles or maps. As an example of the power and versatility of shapemol, we build a model of the molecular envelope of the planetary nebula NGC 6302 and compare it with 12CO and 13CO J=2-1 interferometric maps from SMA and high-J transitions from Herschel/HIFI. We find the molecular envelope to have a complex, broken ring-like structure with an inner, hotter region and several 'fingers' and high-velocity blobs, emerging outwards from the plane of the ring. We derive a mass of 0.11 Msun for the molecular envelope.

The molecular envelope of CRL 618: A new model based onHerschel/HIFI observations

We study the physical properties and molecular excitation of the different warm gas components found in the protoplanetary nebula CRL 618. We revise our previous Herschel/HIFI observations, which consist of several 12CO and 13CO lines in the far-infrared/sub-mm band. These data have been re-analyzed in detail by improving calibration, the signal-to-noise-ratio, and baseline substraction. We identify the contributions of the different nebular components to the line profiles. We have used a spatio-kinematical model to better constrain the temperature, density, and kinematics of the molecular components probed by the improved CO observations. The 12CO and 13CO J=16-15, J=10-9, and J=6-5 transitions are detected in this source. The line profiles present a composite structure showing spectacular wings in some cases, which become dominant as the energy level increases. Our analysis of the high-energy CO emission with the already known low-energy J=2-1 and J=1-0 lines confirms that the high-velocity component, or fast bipolar outflow, is hotter than previously estimated with a typical temperature of ~300 K. This component may then be an example of a very recent acceleration of the gas by shocks that has not yet cooled down. We also find that the dense central core is characterized by a very low expansion velocity, ~5 km/s, and a strong velocity gradient. We conclude that this component is very likely to be the unaltered circumstellar layers that are lost in the last AGB phase, where the ejection velocity is particularly low. The physical properties of the other two nebular components, the diffuse halo and the double empty shell, more or less agrees with the estimations derived in previous models.

Spatio-kinematic modelling of Abell 65, a double-shelled planetary nebula with a binary central star★

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell 65, which is known to host a post-common envelope, binary, central star system. As such, this object is of great interest in studying the link between nebular morphology and central star binarity. [OIII]5007A and H-alpha+[NII]6584A longslit spectra and imagery of Abell 65 were obtained with the Manchester Echelle Spectrometer on the 2.1-m telescope at the San Pedro Martir Observatory (MES-SPM). Further [OIII]5007A longslit spectra were obtained with the Ultraviolet and Visual Echelle Spectrograph on the Very Large Telescope (VLT-UVES). These data were used to develop a spatio-kinematical model for the [OIII]5007A emission from Abell 65. A "best-fit" model was found by comparing synthetic spectra and images rendered from the model to the data. The model comprises an outer shell and an inner shell, with kinematical ages of 15000 +/- 5000 yr kpc^-1 and 8000 +/- 3000 yr kpc^-1, respectively. Both shells have peanut-shaped bipolar structures with symmetry axes at inclinations of 55 +/- 10 deg (to the line-of-sight) for the outer shell and 68 +/- 10 deg for the inner shell. The near-alignment between the nebular shells and the binary orbital inclination (of 68 +/- 2 deg) is strongly indicative that the binary is responsible for shaping the nebula. Abell 65 is one of a growing number of planetary nebulae (seven to date, including Abell 65 itself) for which observations and modelling support the shaping influence of a central binary.

MOLECULAR LINE OBSERVATIONS OF THE CARBON-RICH CIRCUMSTELLAR ENVELOPE CIT 6 AT 7 mm WAVELENGTHS

We present a {lambda}7 mm spectral line survey of the carbon-rich circumstellar envelope (CSE) CIT 6 using the 45 m telescope at the Nobeyama Radio Observatory. A total of 25 spectral features belonging to five molecular species (HC{sub 3}N, HC{sub 5}N, HC{sub 7}N, SiO, and CS) are detected, enabling us to investigate the chemistry of cyanopolyyne chains. The line strengths are compared with those of the proto-typical carbon-rich CSE IRC+10216. The results show that the cyanopolyyne molecules are enhanced in CIT 6, suggesting that it is more evolved than IRC+10216. In order to investigate the structure of CIT 6, we have constructed a three-dimensional spatiokinematic model. By comparing the observed line profiles with the models, we conclude that this envelope is asymmetric and is composed of several incomplete shells.

Modeling the physical and excitation conditions of the molecular envelope of NGC 7027

The link between the shaping of bipolar planetary nebulae and the mass ejection activity of their central stars is still poorly understood. Appropriately characterizing the evolution of the shells ejected during the late stages of evolution is vital to gain insight into the mechanism of nebular shaping. Herschel/HIFI provides an invaluable tool by opening a new window from which to probe warm molecular gas (~50-1000 K). We present a radiative-transfer, spatio-kinematic modeling of the molecular envelope of the young planetary nebula NGC 7027 in high- and low-J 12CO and 13CO transitions observed by Herschel/HIFI and IRAM 30-m, and discuss the structure and dynamics of the molecular envelope. We have developed a code which, used along with the existing SHAPE software, implements spatio-kinematic modeling with accurate non-LTE calculations of line excitation and radiative transfer in molecular species. We have used this code to build a relatively simple "russian doll" model of the molecular envelope of NGC 7027. The model nebula consists of four nested, mildly bipolar shells plus a pair of high-velocity blobs. The innermost shell is the thinnest and shows a significant jump in physical conditions (temperature, density, abundance and velocity) with respect to the adjacent shell. This is a clear indication of a shock front in the system. Each of the high-velocity blobs is divided into two sections with different physical conditions. The presence of H2O in NGC 7027, a C-rich nebula, is likely due to photo-induced chemistry from the hot central star. The computed molecular mass of the nebula is 1.3 Msun, compatible with estimates from previous works.