We present a new approach to study planetary nebulae using integral field spectroscopy. VLT@VIMOS datacube of the planetary nebula Abell 14 is analysed in three different ways by extracting: (i) the integrated spectrum, (ii) 1-dimensional simulated long slit spectra for different position angles and (iii) spaxel-by-spaxel spectra. These data are used to built emission-line diagnostic diagrams and explore the ionization structure and excitation mechanisms combining data from 1- and 3- dimensional photoionization models. The integrated and 1D simulated spectra are suitable for developing diagnostic diagrams, while the spaxel spectra can lead to misinterpretation of the observations. We find that the emission-line ratios of Abell 14 are consistent with UV photo-ionized emission, however there are some pieces of evidence of an additional thermal mechanism. The chemical abundances confirm its previous classification as a Type I planetary nebula, without spatial variation. We find, though, variation in the ionization correction factors (ICFs) as a function of the slit position angle. The star at the geometric centre of Abell 14 has an A5 spectral type with an effective temperature of Teff = 7909$\pm$135 K and surface gravity log(g) = 1.4$\pm$0.1 cm s$^{-2}$. Hence, this star cannot be responsible for the ionization state of the nebula. Gaia parallaxes of this star yield distances between 3.6 and 4.5 kpc in good agreement with the distance derived from a 3-dimensional photoionization modelling of Abell 14, indicating the presence of a binary system at the centre of the planetary nebula.

We present a detailed analysis of deep VLT/X-Shooter observations of the planetary nebula Tc 1. We calculate gas temperature, density, extinction, and abundances for several species from the empirical analysis of the total line fluxes. In addition, a spatially resolved analysis of the most intense lines provides the distribution of such quantities across the nebula. The new data reveal that several lines exhibit a double peak spectral profile consistent with the blue- and red-shifted components of an expanding spherical shell. The study of such components allowed us to construct for the first time a three-dimensional morphological model, which reveals that Tc 1 is a slightly elongated spheroid with an equatorial density enhancement seen almost pole on. A few bright lines present extended wings (with velocities up to a few hundred km/s), but the mechanism producing them is not clear. We constructed photoionization models for the main shell of Tc 1. The models predict the central star temperature and luminosity, as well as the nebular density and abundances similar to previous studies. Our models indicate that Tc 1 is located at a distance of approximately 2 kpc. We report the first detection of the [Kr III] 6825 A emission line, from which we determine the Krypton abundance. Our model indicates that the main shell of Tc 1 is matter bounded; leaking H ionizing photons may explain the ionization of its faint AGB-remnant halo.

In this work we revisit the issue of the rotation speed of the spiral arms and the location of the corotation radius of our Galaxy. This research was performed using homogeneous data set of young open clusters (age < 50 Myr) determined from Gaia DR2 data. The stellar astrometric membership were determined using proper motions and parallaxes, taking into account the full covariance matrix. The distance, age, reddening and metallicity of the clusters were determined by our non subjective multidimensional global optimization tool to fit theoretical isochrones to Gaia DR2 photometric data. The rotation speed of the arms is obtained from the relation between age and angular distance of the birthplace of the clusters to the present-day position of the arms. Using the clusters belonging to the Sagittarius-Carina, Local and Perseus arms, and adopting the Galactic parameters $R_0$ = 8.3 kpc and $V_0$ = 240 km\,s$^{-1}$, we determine a pattern speed of $28.2 \pm 2.1$ km\,s$^{-1}$\,kpc$^{-1}$, with no difference between the arms. This implies that the corotation radius is $R_c = 8.51 \pm 0.64$ kpc, close to the solar Galactic orbit ($R_c/R_0 = 1.02\pm0.07$).

We present a study of 7 southern open clusters based on UBVRI CCD photomety (Johnsons-Cousins system) and Gaia DR2 data. Dias 4, Dias 6 and four other clusters had UBVRI photometric observations determined for the first time. From the observational UBVRI data we obtained photometric membership probability estimates and, using the proper motions from the UCAC5 catalog, we also determined the kinematic membership. From Gaia DR2 astrometric data we determine the stellar membership using proper motions and parallaxes, taking into account the full covariance matrix. For both independent sets of data and membership we apply our non subjective multidimensional global optimization tool to fit theoretical isochrones to determine the distance, age, reddening, metallicity and binary fraction of the clusters. The results of the mean proper motions, distances and ages are in agreement, but the ones obtained from Gaia DR2 data are more precise in both membership selection and estimated parameters. In the case of NGC 6087, the Cepheid S Nor, member of the open cluster, was used to obtain an independent distance estimate, confirming the one determined by our fitting method. We also report a serendipitous discovery of two new clusters in the extended field near what was originally Dias 4.

We present mean proper motions and membership probabilities of individual stars for optically visible open clusters, which have been determined using data from the UCAC5 catalog. This follows our previous studies with the UCAC2 and UCAC4 catalogs, but now using improved proper motions in the GAIA reference frame. In the present study results were obtained for a sample of 1108 open clusters. For five clusters, this is the first determination of mean proper motion, and for the whole sample, we present results with a much larger number of identified astrometric member stars than on previous studies. It is the last update of our Open cluster Catalog based on proper motion data only. Future updates will count on astrometric, photometric and spectroscopic GAIA data as input for analyses.

The recent discovery of a roughly simultaneous periodic variability in the light curves of the BL Lac object PG 1553+113 at several electromagnetic bands represents the first case of such odd behavior reported in the literature. Motivated by this, we analyzed 15 GHz interferometric maps of the parsec-scale radio jet of PG 1553+113 to verify the presence of a possible counterpart of this periodic variability. We used the Cross-Entropy statistical technique to obtain the structural parameters of the Gaussian components present in the radio maps of this source. We kinematically identified seven jet components formed coincidentally with flare-like features seen in the $\gamma$-ray light curve. From the derived jet component positions in the sky plane and their kinematics (ejection epochs, proper motions, and sky position angles), we modeled their temporal changes in terms of a relativistic jet that is steadily precessing in time. Our results indicate a precession period in the observer's reference frame of $2.24\pm0.03$ years, compatible with the periodicity detected in the light curves of PG 1553+113. However, the maxima of the jet Doppler boosting factor are systematically delayed relative to the peaks of the main $\gamma$-ray flares. We propose two scenarios that could explain this delay, both based on the existence of a supermassive black hole binary system in PG 1553+113. We estimated the characteristics of this putative binary system that also would be responsible for driving the inferred jet precession.

The main objective of this work is to determine the cluster members of 1876 open clusters, using positions and proper motions of the astrometric catalogue UCAC4. For this purpose we apply three different methods, all them based on a Bayesian approach, but with different formulations: a purely parametric method, another completely non-parametric algorithm, and a third, recently developed by Sampedro & Alfaro, using both formulations at different steps of the whole process. The first and second statistical moments of the members phase-space subspace, obtained after applying the three methods, are compared for every cluster. Although, on average, the three methods yield similar results, specific differences between them, as well as for some particular clusters, are also present. The comparison with other published catalogues shows good agreement. We have also estimated for the first time the mean proper motion for a sample of 18 clusters. The results are organized in a single catalogue formed by two main files, one with the most relevant information for each cluster, partially including that in UCAC4, and the other showing the individual membership probabilities for each star in the cluster area. The final catalogue, with an interface design that enables an easy interaction with the user, is available in electronic format at SSG-IAA (http://ssg.iaa.es/en/content/sampedro-cluster-catalog) website.

A number of planetary nebulae show binary central stars and significant abundance discrepancies between values estimated from colisionally excited lines when compared to the same abundances estimated from recombination lines. One approach to investigate this yet unsolved problem is using spatially resolved images of emission lines in an attempt to detect a possibly distinct metal rich component in the nebula. In this work we present results of spatially resolved bundance analysis of NGC 6778 based on data gathered from VLT VIMOS-IFU. We discuss the spatial variations found as well as possible limitations of the method in answering questions about abundance variations.

The discrepancy between abundances computed using optical recombination lines (ORLs) and collisionally excited lines (CELs) is a major, unresolved problem with significant implications for the determination of chemical abundances throughout the Universe. In planetary nebulae (PNe), the most common explanation for the discrepancy is that two different gas phases coexist: a hot component with standard metallicity, and a much colder plasma enhanced in heavy elements. This dual nature is not predicted by mass loss theories, and direct observational support for it is still weak. In this work, we present our recent findings that demonstrate that the largest abundance discrepancies are associated with close binary central stars. OSIRIS-GTC tunable filter imaging of the faint O II ORLs and MUSE-VLT deep 2D spectrophotometry confirm that O II ORL emission is more centrally concentrated than that of [O III] CELs and, therefore, that the abundance discrepancy may be closely linked to binary evolution.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

In the second paper of the series we continue the investigation of open cluster fundamental parameters using a robust global optimization method to fit model isochrones to photometric data. We present optical UBVRI CCD photometry (Johnsons-Cousins system) observations for 24 neglected open clusters, of which 14 have high quality data in the visible obtained for the first time, as a part of our ongoing survey being carried out in the 0.6m telescope of the Pico dos Dias Observatory in Brazil. All objects were then analyzed with a global optimization tool developed by our group which estimates the membership likelihood of the observed stars and fits an isochrone from which a distance, age, reddening, total to selective extinction ratio $R_{V}$ (included in this work as a new free parameter) and metallicity are estimated. Based on those estimates and their associated errors we analyzed the status of each object as real clusters or not, finding that two are likely to be asterisms. We also identify important discrepancies between our results and previous ones obtained in the literature which were determined using 2MASS photometry.

We present the first direct image of the high-metallicity gas component in a planetary nebula (NGC 6778), taken with the OSIRIS Blue Tunable Filter centered on the O II 4649+50 angstroms optical recombination lines (ORLs) at the 10.4m Gran Telescopio Canarias. We show that the emission of these faint O II ORLs is concentrated in the central parts of the planetary nebula and is not spatially coincident either with emission coming from the bright [O III] 5007 angstroms collisionally excited line (CEL) or the bright H alpha recombination line. From monochromatic emission line maps taken with VIMOS at the 8.2m Very Large Telescope, we find that the spatial distribution of the emission from the auroral [O III] 4363 line resembles that of the O II ORLs but differs from nebular [O III] 5007 CEL distribution, implying a temperature gradient inside the planetary nebula. The centrally peaked distribution of the O II emission and the differences with the [O III] and H I emission profiles are consistent with the presence of an H-poor gas whose origin may be linked to the binarity of the central star. However, determination of the spatial distribution of the ORLs and CELs in other PNe, and a comparison of their dynamics is needed to further constrain the geometry and ejection mechanism of the metal-rich (H-poor) component and hence, understand the origin of the abundance discrepancy problem in PNe.

Abell 14 is a poorly studied object despite being considered a born again planetary nebula. We performed a detailed study of its 3D morphology and ionization structure using the SHAPE and MOCASSIN codes. We found that Abell 14 is a highly evolved, bipolar nebula with a kinematical age of $\sim$19,400 yr for a distance of 4 kpc. The high He abundance, and N/O ratio indicate a progenitor of 5 $M_{\odot}$ that has experienced the third dredge-up and hot bottom burning phases. The stellar parameters of the central source reveal a star at a highly evolved stage near to the white dwarf cooling track, being inconsistent with the born again scenario. The nebula shows unexpectedly strong [N I] $\lambda 5200$ and [O I] $\lambda 6300$ emission lines indicating possible shock interactions. Abell 14 appears to be a member of a small group of highly evolved, extreme Type-I PNe. The members of this group lie at the lower-left corner of the PNe regime on the [N II]/H$\alpha$ vs. [S II]/H$\alpha$ diagnostic diagram, where shock--excited regions/objects are also placed. The low luminosity of their central stars, in conjunction with the large physical size of the nebulae, result in a very low photo-ionization rate, which can make any contribution of shock interaction easily perceptible, even for small velocities.

We present BV and u'g'r'i' CCD photometry of the central region of NGC 2482. We also present BVu'g' CCD photometry of five clusters that have been poorly studied in the past: Ruprecht 42, Ruprecht 51, Ruprecht 153, Ruprecht 154; and AH03 J0748-26.9, which to our knowledge has not been studied before. Using a global optimization technique that eliminates much of the subjectivity previously inherent in main sequence fitting studies, we obtain values of the distances, ages, and metallicities of the clusters, with robust estimates of the uncertainties of these fundamental parameters. Four of our clusters are less than ~1.3 kpc beyond the Sun's distance from the Galactic Center and have essentially solar metallicity. The metallicities of those clusters more distant from the Galactic Center are consistent with a 0.3 dex step to lower [Fe/H] found in other studies.

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

We present 23 interferometric images of parsec-scale jet of the quasar PKS 1741--03 at 15, 24 and 43 GHz spanning about 13 yr. We model the images as a superposition of discrete two--dimensional elliptical Gaussian components, with parameters determined by the cross--entropy technique. All the images present a spatially unresolved component (core) and usually two or three components receding from it. The same components were found in simultaneous 24 and 43 GHz maps, showing the robustness of our model-fitting. The core-shift opacity effect between these frequencies is weak. We have identified seven components moving along straight lines at constant apparent superluminal speeds ($3.5\le\beta_\mathrm{obs}\le 6.1$), with different sky position angles ($-186\deg\le \eta\le -125\deg$). The core flux density tracks quite well the fluctuations seen in the historical single-dish light curve at 14.5 GHz, with no measurable delay. The total flux density from the moving jet components is delayed $\sim$2 yr in relation to the core light curve, roughly the same as the lag between the ejection epoch and the maximum flux density in the light curves of the jet components. We propose that there are three non-exclusive mechanisms for producing these delays. From the kinematics of the most robust jet components and the core brightness temperature, we determined the bulk Lorentz factor ($4.8\le\gamma\le 24.5$) and the jet viewing angle ($0.\deg 35 \le \theta \le 4.\deg 2$); these values agree with previous estimates from the spectral energy distribution of PKS 1741--03 and its radio variability.

For the past decade observations of the alignement of PNe symmetries with respect to the galactic disk have led to conflicting results. Recently observational evidence for alignment between PNe and local interstellar magnetic fields in the central part of the Galaxy ($b < 5^\circ$) has been found. We studied the role of the interstellar magnetic field on the dynamical evolution of a PN by means of an analytical model and from 3D MHD numerical simulations. We test under what conditions typical ejecta would have their dynamics severely modified by an interstellar magnetic field. We found that uniform fields of $> 100\mu$G are required in order to be dynamically dominant. This is found to occur only at later evolutionary stages, therefore being unable to change the general morphology of the nebula. However, the symmetry axis of bipolar and elliptical nebulae end up aligned to the external field. This result can explain why different samples of PNe result in different conclusions regarding the alignment of PNe. Objects located at high galactic latitudes, or at large radii, should present no preferential alignment with respect to the galactic plane. PNe located at the galactic centre and low latitudes would, on the other hand, be preferentiably aligned to the disk. Finally, we present synthetic polarization maps of the nebulae to show that the polarization vectors, as well as the field lines at the expanding shell, are not uniform even in the strongly magnetized case, indicating that polarization maps of nebulae are not adequate in probing the orientation, or intensity, of the dominant external field.

Optical integral-field spectroscopy was used to investigate the planetary nebula NGC 3242. We analysed the main morphological components of this source, including its knots, but not the halo. In addition to revealing the properties ofthe physical and chemical nature of this nebula, we also provided reliable spatially resolved constraints that can be used for future photoionisation modelling of the nebula. The latter is ultimately necessary to obtain a fully self-consistent 3D picture of the physical and chemical properties of the object. The observations were obtained with the VIMOS instrument attached to VLT-UT3. Maps and values for specific morphological zones for the detected emission-lines were obtained and analysed with routines developed by the authors to derive physical and chemical conditions of the ionised gas in a 2D fashion. We obtained spatially resolved maps and mean values of the electron densities, temperatures, and chemical abundances, for specific morphological structures in NGC 3242. These results show the pixel-to-pixel variations of the the small- and large-scale structures of the source. These diagnostic maps provide information free from the biases introduced by traditional single long-slit observations. In general, our results are consistent with a uniform abundance distribution for the object, whether we look at abundance maps or integrated fluxes from specified morphological structures. The results indicate that special care should be taken with the calibration of the data and that only data with extremely good signal-to-noise ratio and spectral coverage should be used to ensure the detection of possible spatial variations.

The metallicity is a critical parameter that affects the correct determination fundamental characteristics stellar cluster and has important implications in Galactic and Stellar evolution research. Fewer than 10 % of the 2174 currently catalog open clusters have their metallicity determined in the literature. In this work we present a method for estimating the metallicity of open clusters via non-subjective isochrone fitting using the cross-entropy global optimization algorithm applied to UBV photometric data. The free parameters distance, reddening, age, and metallicity simultaneously determined by the fitting method. The fitting procedure uses weights for the observational data based on the estimation of membership likelihood for each star, which considers the observational magnitude limit, the density profile of stars as a function of radius from the center of the cluster, and the density of stars in multi-dimensional magnitude space. We present results of [Fe/H] for nine well-studied open clusters based on 15 distinct UBV data sets. The [Fe/H] values obtained in the ten cases for which spectroscopic determinations were available in the literature agree, indicating that our method provides a good alternative to determining [Fe/H] by using an objective isochrone fitting. Our results show that the typical precision is about 0.1 dex.

BL Lacertae is the prototype of the BL Lac class of active galactic nuclei, exhibiting intensive activity on parsec (pc) scales, such as intense core variability and multiple ejections of jet components. In particular, in previous works the existence of precession motions in the pc-scale jet of BL Lacertae has been suggested. In this work we revisit this issue, investigating temporal changes of the observed right ascension and declination offsets of the jet knots in terms of our relativistic jet-precession model. The seven free parameters of our precession model were optimized via a heuristic cross-entropy method, comparing the projected precession helix with the positions of the jet components on the plane of the sky and imposing constraints on their maximum and minimum superluminal velocities. Our optimized best model is compatible with a jet having a bulk velocity of 0.9824c, which is precessing with a period of about 12.1 yr in the observer's reference frame and changing its orientation in relation to the line of sight between 4 and 5 degrees, approximately. Assuming that the jet precession has its origin in a supermassive binary black hole system, we show that the 2.3-yr periodic variation in the structural position angle of the very-long-baseline interferometry (VLBI) core of BL Lacertae reported by Stirling et al. is compatible with a nutation phenomenon if the secondary black hole has a mass higher than about six times that of the primary black hole.

We present a new technique for obtaining model fittings to VLBI images of astrophysical jets. The method minimizes a performance function proportional to the sum of the squared difference between the model and observed images. The model image is constructed by summing elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak intensity, eccentricity, amplitude and orientation angle of the major axis. We present results for the fitting of two main benchmark jets: the first, constructed from three individual Gaussian sources, the second formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio maps. We show that our technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the traditional AIPS task IMFIT when the image is relatively simple (e.g., few components). For more complex maps, our method displays superior performance in recovering the parameters of the jet components. Our methodology is also able to show quantitatively the number of individual components present in an image. An additional application of the cross-entropy technique to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy technique must be used in situations involving the analysis of complex emission regions having more than three sources, even though it is substantially slower than current model fitting tasks (at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). As in the case of any model fitting performed in the image plane, caution is required in analyzing images constructed from a poorly sampled (u,v) plane.

Continuing our series of papers on the 3-D structure and accurate distances of Planetary Nebulae (PNe), we present here the results obtained for the planetary nebula NGC\,40. Using data from different sources and wavelengths, we construct 3-D photoionization models and derive the physical quantitities of the ionizing source and nebular gas. The procedure, discussed in detail in the previous papers, consists of the use of 3-D photoionization codes constrained by observational data to derive the three-dimensional nebular structure, physical and chemical characteristics and ionizing star parameters of the objects by simultaneously fitting the integrated line intensities, the density map, the temperature map, and the observed morphologies in different emission lines. For this particular case we combined hydrodynamical simulations with the photoionization scheme in order to obtain self-consistent distributions of density and velocity of the nebular material. Combining the velocity field with the emission line cubes we also obtained the synthetic position-velocity plots that are compared to the observations. Finally, using theoretical evolutionary tracks of intermediate and low mass stars, we derive the mass and age of the central star of NGC\,40 as $(0.567 \pm 0.06)$M$_{\odot}$ and $(5810 \pm 600)$yrs, respectively. The distance obtained from the fitting procedure was $(1150 \pm 120)$pc.

In the last decades we witnessed an increase in studies of open clusters of the Galaxy, especially because of the good determination for a wide range of values of parameters such as age, distance, reddening, and proper motion. The reliable determination of the parameters strongly depends on the photometry available and especially on the U filter, which is used to obtain the color excess E(B-V) through the color-color diagram (U-B) by (B-V) by fitting a zero age main-sequence. Owing to the difficulty of performing photometry in the U band, many authors have tried to obtain E(B-V) without the filter. But because of the near linearity of the color-color diagrams that use the other bands, combined with the fact that most fitting procedures are highly subjective (many done "by eye") the reliability of those results has always been questioned. Our group has recently developed, a tool that performs isochrone fitting in open-cluster photometric data with a global optimization algorithm, which removes the need to visually perform the fits and thus removes most of the related subjectivity. Here we apply our method to a set of synthetic clusters and two observed open clusters (Trumpler 1 and Melotte 105) using only photometry for the BVRI bands. Our results show that, considering the cluster structural variance caused only by photometric and Poisson sampling errors, our method is able to recover the synthetic cluster parameters with errors of less than 10% for a wide range of ages, distances, and reddening, which clearly demonstrates its potential. The results obtained for Trumpler 1 and Melotte 105 also agree well with previous literature values.

In this paper we present an analysis of the physical and chemical conditions of the planetary nebula NGC 40 through spatially-resolved spectroscopic maps. We also introduce a new algorithm --2D_NEB-- based on the well-established IRAF nebular package, which was developed to enable the use of the spectroscopic maps to easily estimate the astrophysical quantities of ionised nebulae. The 2D_NEB was benchmarked, and we clearly show that it works properly, since it compares nicely with the IRAF nebular software. Using this software, we derive the maps of several physical parameters of NGC 40. From these maps, we conclude that Te[NII] shows only a slight temperature variation from region to region, with its values constrained between ~8,000 K and ~9,500 K. Electron densities, on the other hand, have a much more prominent spatial variation, as Ne[SII] values vary from ~1,000 cm^(-3) to ~3,000 cm^(-3). Maps of the chemical abundances also show significant variations. From the big picture of our work, we strongly suggest that analysis with spatial resolution be mandatory for more complete study of the physical and chemical properties of planetary nebulae.

We present a new technique to fit color-magnitude diagrams of open clusters based on the Cross-Entropy global optimization algorithm. The method uses theoretical isochrones available in the literature and maximizes a weighted likelihood function based on distances measured in the color-magnitude space. The weights are obtained through a non parametric technique that takes into account the star distance to the observed center of the cluster, observed magnitude uncertainties, the stellar density profile of the cluster among others. The parameters determined simultaneously are distance, reddening, age and metallicity. The method takes binary fraction into account and uses a Monte-Carlo approach to obtain uncertainties on the determined parameters for the cluster by running the fitting algorithm many times with a re-sampled data set through a bootstrapping procedure. We present results for 9 well studied open clusters, based on 15 distinct data sets, and show that the results are consistent with previous studies. The method is shown to be reliable and free of the subjectivity of most previous visual isochrone fitting techniques.

Evidence of jet precession in many galactic and extragalactic sources has been reported in the literature. Much of this evidence is based on studies of the kinematics of the jet knots, which depends on the correct identification of the components to determine their respective proper motions and position angles on the plane of the sky. Identification problems related to fitting procedures, as well as observations poorly sampled in time, may influence the follow up of the components in time, which consequently might contribute to a misinterpretation of the data. In order to deal with these limitations, we introduce a very powerful statistical tool to analyse jet precession: the cross-entropy method for continuous multi-extremal optimisation. Only based on the raw data of the jet components (right ascension and declination offsets from the core), the cross-entropy method searches for the precession model parameters that better represent the data. In this work we present a large number of tests to validate this technique, using synthetic precessing jets built from a given set of precession parameters. Aiming to recover these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated to the method. Our results have shown that even in the most challenging tests, the cross-entropy method was able to find the correct parameters within 1%-level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.

We show that the inclination to the line of sight of bipolar nebulae strongly affects some of their observed properties. We model these objects as having a spherically symmetric Planetary Nebula and a dusty equatorial density enhancement that produces extinction that varies with the viewing angle. Our sample of 29 nebulae taken from the literature shows a clear correlation between the inclination angle and the near-IR and optical photometric properties as well as the apparent luminosity of the objects. As the inclination angle increases --the viewing angle is closer to the equatorial plane-- the objects become redder, their average apparent luminosity decreases, and their average projected expansion velocity becomes smaller. We compute two-dimensional models of stars embedded in dusty disks of various shapes and compositions and show that the observed data can be reproduced by disk-star combinations with reasonable parameters. To compare with the observational data, we generate sets of model data by randomly varying the star and disk parameters within a physically meaningful range. We conclude that a only a smooth pole to equator density gradient agrees with the observed phenomena and that thin, equatorially concentrated disks can be discarded.

Continuing our series of papers on the three-dimensional (3-D) structures of and accurate distances to Planetary Nebulae (PNe), we present our study of the planetary nebula NGC6781. For this object we construct a 3-D photoionization model and, using the constraints provided by observational data from the literature we determine the detailed 3-D structure of the nebula, the physical parameters of the ionizing source and the first precise distance. The procedure consists in simultaneously fitting all the observed emission line morphologies, integrated intensities and the 2-D density map from the [SII] line ratios to the parameters generated by the model, and in an iterative way obtain the best fit for the central star parameters and the distance to NGC6781, obtaining values of 950+-143pc and 385 Lsun for the distance and luminosity of the central star respectively. Using theoretical evolutionary tracks of intermediate and low mass stars, we derive the mass of the central star of NGC6781 and its progenitor to be 0.60+-0.03 Msun and 1.5+-0.5 Msun respectively.

We provide the first age estimates for two recently discovered white dwarf-red subdwarf systems, LHS 193AB and LHS 300AB. These unusual systems provide a new opportunity for linking the reliable age estimates for the white dwarfs to the (measurable) metallicities of the red subdwarfs. We have obtained precise photometry in the $V_{J}R_{KC}I_{KC}JH$ bands and spectroscopy covering from 6000\AA to 9000\AA for the two new systems, as well as for a comparison white dwarf-main sequence red dwarf system, GJ 283 AB. Using model grids available in the literature, we estimate the cooling age as well as temperature, surface gravity, mass, progenitor mass and {\it total} lifetimes of the white dwarfs. The results indicate that the two new systems are probably ancient thick disk objects with ages of at least 6-9 Gyr. We also conduct searches of red dwarf and white dwarf compendia from SDSS data and the L{\'e}pine Shara Proper Motion (LSPM) catalog for additional common proper motion white dwarf-red subdwarf systems. Only seven new candidate systems are found, which indicates the rarity of these systems.

We present the results of a spatio-kinematic study of the planetary nebula Menzel 1 using spectro-photometric mapping and a 3-D photoionization code. We create several 2-D emission line images from our long-slit spectra, and use these to derive the line fluxes for 15 lines, the Halpha/Hbeta extinction map, and the [SII] line ratio density map of the nebula. We use our photoionization code constrained by these data to derive the three-dimensional nebular structure and ionizing star parameters of Menzel 1 by simultaneously fitting the integrated line intensities, the density map, and the observed morphologies in several lines, as well as the velocity structure. Using theoretical evolutionary tracks of intermediate and low mass stars, we derive a mass for the central star of 0.63+-0.05 Msolar. We also derive a distance of 1050+_150 pc to Menzel 1.

We present the results of mapping the planetary nebula NGC6369 using multiple long slit spectra taken with the CTIO 1.5m telescope. We create two dimensional emission line images from our spectra, and use these to derive fluxes for 17 lines, the Halpha/Hbeta extinction map, the [SII] line ratio density map, and the [NII] temperature map of the nebula. We use our photoionization code constrained by these data to determine the distance, the ionizing star characteristics, and show that a clumpy hour-glass shape is the most likely three-dimensional structure for NGC6369. Note that our knowledge of the nebular structure eliminates all uncertainties associated with classical distance determinations, and our method can be applied to any spatially resolved emission line nebula. We use the central star, nebular emission line, and optical+IR luminosities to show that NGC6369 is matter bound, as about 70% of the Lyman continuum flux escapes. Using evolutionary tracks from Bloecker(1995) we derive a central star mass of about 0.65 Msolar.

We have measured the proper motion (PM) and kinematics of the ansae in NGC 7009 using high dispersion echelle spectra and archive narrow band HST images. Assuming that the ansae are moving at equal and opposite velocities from the central star we obtain a system radial velocity of -53+-2km/s, the eastern ansa approaching and the western ansa receding at Vr=5.3+-1km/s with respect to this value. The PM of the eastern ansa is 28+-8mas/yr, which with our weighted distance to NGC 7009 of 0.86+-0.34kpc gives Vexp=112+-32km/s. The electron temperature and density in both ansae were determined to be Te=9000+-400K and Ne=2300+-400/cc. The dynamic age of the ansae is 925+-260yrs.and the implied PM of the central star is 1+-0.5mas/yr. This is in qualitative but not quantitative agreement with previous work.

In this work we investigate how a circumstellar disk affects the radiation emitted by an embedded star. We show correlations obtained from broad-band observations of bipolar nebulae indicating that an orientation effect is at play in these systems. The FIR radiation relative to total radiation increases with inclination while the NIR and BVR fractions decrease. This is an expected effect if we consider the system as being made up of a dense dusty disk being irradiated by a hot star. We calculate 2-D models to try and reproduce the observed behavior with different disk and star configurations.

We have measured the expansion velocities and proper motion of the ansae in NGC7009 using high dispersion echelle spectra and archive narrow band HST images. Assuming that the ansae are moving at equal and opposite velocities from the central star we obtain an average system radial velocity of -54 +-2 km/s, the eastern ansa approaching and the western ansa receding at Vr = 5.5 +-1 km/s relative to this value. Only the proper motion of the eastern ansa could be measured, leading to 2.8 +- 0.8 arcsec/century, or Vt=130 +- 40.d km/s, where d is the distance to the nebula in kpc. Additionally, the electron temperature and density for each ansa was measured using line intensity ratios. The results are Te = 9000 K and ne = 2000 /cm3 for both ansae within the errors.

We show that there are strong links between certain types of asymmetrical Planetary Nebulae (PNe) and symbiotic stars. Symbiotics are binaries and several have extended optical nebulae all of which are asymmetrical and ~-40% are bipolar. Bipolar PNe are likely to be formed by binaries and share many properties with symbiotic nebulae (SyNe). Some PNe show point symmetry which is naturally explained by precession in a binary system. M2-9 has both point and plane symmetry, and has been shown to have a binary central object. We show that inclination on the sky affects the observed properties of bipolar nebulae due to enhanced equatorial densities, and compare observations of a sample of BPNe with a simple model. Good agreement is obtained between model predicted and observed IR/optical flux ratios and apparent luminosities, which further confirms the binary hypothesis.

We use a 3D photoionization modeling tool to study the morpho-kinematic properties of the Planetary Nebula (PN) NGC 3132. We show that it is possible to reproduce the low resolution observations (spectra and images) with an ellipsoidal shell. However, high resolution observations, as those showing a density variation along the nebula and the [O III]5007 velocity profiles, definitively rule out this description. We show that a bipolar Diabolo shape with a 40 degres rotation of the symmetry axis relative to the line of sight successfully reproduces the observed images, as well as the high resolution observations.

We present a new modeling tool for planetary nebulae (PNe), based on 3D photoionization calculations. From models for two theoretical PNe, we show that the enhancement in the equatorial zone observed in several PNe is not necessarily due to a density gradient, as usually interpreted. We study the morpho-kinematic properties of the PN NGC 3132 and show that a bipolar Diabolo shape successfully reproduces the observed images, as well as the high resolution observations.