In this proceeding, we present a short review of the fascinating nebulosities of symbiotic binary R Aquarii. The R Aquarii system, comprising the central binary and surrounding nebular material, has been the subject of near-continuous study since its discovery, with a few hundred papers listed in ADS. As such, it is impossible to provide here the comprehensive review that R Aquarii deserves, instead we chose to focus on the nebulosities -- covering both our own research and other relevant results from the literature.

We present an analysis of the evolution of circumstellar dust and molecules in the environment of the very late thermal pulse object V4334 Sgr (Sakurai's Object) over a $\sim20$-year period, drawing on ground-, airborne- and space-based infrared photometry and spectroscopy. The dust emission, which started in 1997, resembles a blackbody that cooled from $\sim1200$K in 1998 August to $\sim180$K in 2016 July. The dust mass, assuming amorphous carbon, was $\sim5\times10^{-10}$M$_\odot$ in 1998 August, and we estimate that the total dust mass was $\sim2\times10^{-5}$M$_\odot$ by $\sim2016$. The appearance of a near infrared excess in 2008 suggests a new episode of (or renewed) mass loss began then. We infer lower limits on the bolometric luminosity of the embedded star from that of the dust shell, which rose to $\sim16000$L$_\odot$ before declining to $\sim3000$L$_\odot$. There is evidence for weak 6-7$\mu$m absorption, which we attribute to hydrogenated amorphous carbon formed in material ejected by Sakurai's Object during a mass ejection phase that preceded the 1997 event. We detect small hydrocarbon and other molecules in the spectra, and trace the column densities in hydrogen cyanide (HCN) and acetylene (C$_2$H$_2$). We use the former to determine the $^{12}$C/$^{13}$C ratio to be $6.4\pm0.7$, 14 times smaller than the Solar System value.

In this PhD Thesis we study the occurrence of powerful stellar outbursts in the late phases of evolution of stars like the Sun. This is one of the most controversial topics in the field, as the extraordinary variety of shapes revealed by observations are not expected to be the result of the evolution of stars that are basically spherical all over their lives. The stellar ejecta that we study occur and develop on time-scales of years, even months. This allows us to study their evolution in real-time. A very rare opportunity in Astrophysics! To this aim, we have gathered over the years a unique collection of multi epoch images, from which we produce "movies" that reveal the expansion of the nebulae as projected in the plane of the sky. Images are complemented by spectra that measure the line-of-sight velocities, also required to build 3D models of the dynamical evolution of the outflows, a privileged information for theoretical modelling.

R Aquarii is a symbiotic binary surrounded by a large and complex nebula with a prominent curved jet. It is one of the closest known symbiotic systems, and therefore offers a unique opportunity to study the central regions of these systems and the formation and evolution of astrophysical jets. We studied the evolution of the central jet and outer nebula of R Aqr taking advantage of a long term monitoring campaign of optical imaging, as well as of high-resolution integral field spectroscopy. Narrow-band images acquired over a period of more than 21 years are compared in order to study the expansion and evolution of all components of the R Aqr nebula. The magnification method is used to derive the kinematic ages of the features that appear to expand radially. Integral field spectroscopy of the OIII 5007A emission is used to study the velocity structure of the central regions of the jet. New extended features, further out than the previously known hourglass nebula, are detected. The kinematic distance to R Aqr is calculated to be 178 pc using the expansion of the large hourglass nebula. This nebula of R Aqr is found to be roughly 650 years old, while the inner regions have ages ranging from 125 to 290 years. The outer nebula is found to be well described by a ballistic expansion, while for most components of the jet strong deviations from such behaviour are found. We find that the Northern jet is mostly red-shifted while its Southern part is blue-shifted, apparently at odds with findings from previous studies but almost certainly a consequence of the complex nature of the jet and variations in ionisation and illumination between observations.

The Galactic object MWC 137 was suggested to belong to the group of B[e] supergiants. However, with its large-scale optical bipolar ring nebula and the high velocity jet and knots, it is a rather atypical representative of this class. We performed multi-wavelength observations spreading from the optical to the radio regime. Based on optical imaging and long-slit spectroscopic data we found that the northern parts of the large-scale nebula are predominantly blue-shifted, while the southern regions appear mostly red-shifted. We developed a geometrical model consisting of two double-cones. While various observational features can be approximated with such a scenario, the observed velocity pattern is more complex. Using near-infrared integral-field unit spectroscopy we studied the hot molecular gas in the close vicinity of the star. The emission from the hot CO gas arises in a small-scale disk revolving around the star on Keplerian orbits. While the disk itself cannot be spatially resolved, its emission is reflected by dust arranged in arc-like structures and clumps surrounding MWC 137 on small scales. In the radio regime we mapped the cold molecular gas in the outskirts of the optical nebula. We found that large amounts of cool molecular gas and warm dust embrace the optical nebula in the east, south and west. No cold gas or dust were detected in the north and north-western regions. Despite the new insights on the nebula kinematics gained from our studies, the real formation scenario of the large-scale nebula remains an open issue.

The peculiar emission-line star MWC 137 with its extended optical nebula was recently classified as B[e] supergiant. To study the spatial distribution of its circumstellar molecular gas on small and large scales, we obtained near-infrared and radio observations using SINFONI and APEX, respectively. We find that the hot CO gas is arranged in moving clumpy ring and shell structures close to the star, while a cold CO envelope is encircling the borders of the optical nebula from the south to the west.

Context. To investigate the dynamical properties of protostellar jets. Aims. Determine the proper motion of protostellar jets around Class 0 and Class I sources in an active star forming region in Serpens. Methods. Multi-epoch deep images in the 2.122 $\mu$m line of molecular hydrogen, v=1-0 S(1), obtained with the near-infrared instrument NOTCam over a time-scale of 10 years, are used to determine proper motion of knots and jets. K-band spectroscopy of the brighter knots is used to supply radial velocities, estimate extinction, excitation temperature, and H$_2$ column densities towards these knots. Results. We measure the proper motion of 31 knots over different time scales (2, 4, 6, and 10 years). The typical tangential velocity is around 50 km/s for the 10 year base-line, but for shorter time-scales a maximum tangential velocity up to 300 km/s is found for a few knots. Based on morphology, velocity information and the locations of known protostars, we argue for the existence of at least three partly overlapping and deeply embedded flows, one Class 0 flow and two Class I flows. The multi-epoch proper motion results indicate time-variable velocities of the knots, for the first time directly measured for a Class 0 jet. We find in general higher velocities for the Class 0 jet than for the two Class I jets. While the bolometric luminosites of the three driving sources are about equal, the derived mass flow rate (dM/dt)$_{\rm out}$ is two orders of magnitude higher in the Class 0 flow than in the two Class I flows.

While analysing the images of the IPHAS H$\alpha$ survey, we noticed that the central star of the candidate planetary nebula IPHASXJ211420.0+434136 (also named Ou5) was clearly variable. This is generally considered as an indication of binarity. To confirm it, we performed a photometric monitoring of the central star, and obtained images and spectra of the nebula. The nebular spectrum confirms that IPHASXJ211420.0+434136 is a planetary nebula of moderately high excitation. It has a remarkable morphology with two nested pairs of bipolar lobes, and other unusual features. The light curve of the central star reveals that it is an eclipsing binary system with an orbital period of 8.74 hours. It also displays a strong irradiation effect with an amplitude of 1.5~mag. The presence of multiple bipolar outflows adds constraints to the formation of these nebulae, suggesting the occurrence of discrete ejection events during, or immediately before, the common-envelope phase. IPHASXJ211420.0+434136 also adds evidence to the hypothesis that a significant fraction of planetary nebulae with close binary central stars have a peculiar nebular chemistry and a relatively low nebular mass. This may point to low-mass, low-metallicity progenitors, with additional effects related to the binary evolution. We also suggest that these objects may be relevant to understand the abundance discrepancy problem in planetary nebulae.

Due to their large expansion speed, the apparent growth of nearby nova remnants such as GK Per can be easily resolved from ground-based optical imagery on a timescale of months. If the expansion in the plane of the sky is coupled with the Doppler shift velocities, an almost complete dynamical picture is drawn. We will discuss our latest results of such a study on GK Per.

We present a kinematical study of the optical ejecta of GK Per. It is based on proper motions measurements of 282 knots from ~20 images spanning 25 years. Doppler-shifts are also computed for 217 knots. The combination of proper motions and radial velocities allows a unique 3-D view of the ejecta to be obtained. The main results are: (1) the outflow is a thick shell in which knots expand with a significant range of velocities, mostly between 600 and 1000 km/s; (2) kinematical ages indicate that knots have suffered only a modest deceleration since their ejection a century ago; (3) no evidence for anisotropy in the expansion rate is found; (4) velocity vectors are generally aligned along the radial direction but a symmetric pattern of non-radial velocities is also observed at specific directions; (5) the total Halpha+[NII] flux has been linearly decreasing at a rate of 2.6 % per year in the last decade. The Eastern nebular side is fading at a slower rate than the Western one. Some of the knots displayed a rapid change of brightness during the 2004-2011 period. Over a longer timescale, a progressive circularization and homogenization of the nebula is taking place; (6) a kinematic distance of 400+-30 pc is determined. These results raise some problems to the previous interpretations of the evolution of GK Per. In particular, the idea of a strong interaction of the outflow with the surrounding medium in the Southwest quadrant is not supported by our data.

The nova outburst of V407 Cyg in 2010 Mar. 10 was the first observed for this star but its close resemblance to the well known symbiotic-like recurrent nova RS Oph suggests that it is also a member of this rare type of Galactic novae. The nova was the first detected at $\gamma$-ray energies and is the first known nova explosion for this system. The extensive multiwavelength coverage of this outburst makes it an ideal comparison with the few other outbursts known for similar systems. We extend our previous analysis of the Mira and the expanding shock from the explosion to detail the time development of the photoionized Mira wind, circumstellar medium, and shocked circumstellar environment to derive their physical parameters and how they relate to large scale structure of the environment, extending the previous coverage to more than 500 days after outburst. Absorption lines of Fe-peak ions formed in the Mira wind were visible as P Cyg profiles at low velocity before Day 69, around the time of the X-ray peak and we identified many absorption transitions without accompanying emission for metal lines. The H Balmer lines showed strong P Cyg absorption troughs that weakened during the 2010 observing period, through Day 128. We distinguish the components from the shock, the photoionized environment, and the chromosphere and inner Mira wind using spectra taken more than one year after outburst. The multiple shells and radiative excitation phenomenology are similar to those recently cited for GRBs and SNIa (severely truncated)

A proper determination of the abundance gradient in the Milky Way requires the observation of objects at large galactiocentric distances. With this aim, we are exploring the planetary nebula population towards the Galactic Anticentre. In this article, the discovery and physico-chemical study of a new planetary nebula towards the Anticentre direction, IPHASX J052531.19+281945.1 (PNG 178.1-04.0), is presented. The planetary nebula was discovered from the IPHAS survey. Long-slit follow-up spectroscopy was carried out to confirm its planetary nebula nature and to calculate its physical and chemical characteristics. The newly discovered planetary nebula turned out to be located at a very large galactocentric distance (D_GC=20.8+-3.8 kpc), larger than any previously known planetary nebula with measured abundances. Its relatively high oxygen abundance (12+log(O/H) = 8.36+-0.03) supports a flattening of the Galactic abundance gradient at large galactocentric distances rather than a linearly decreasing gradient.

On 2010 Mar 10, V407 Cyg was discovered in outburst, eventually reaching V< 8 and detected by Fermi. Using medium and high resolution ground-based optical spectra, visual and Swift UV photometry, and Swift X-ray spectrophotometry, we describe the behavior of the high-velocity profile evolution for this nova during its first three months. The peak of the X-ray emission occurred at about day 40 with a broad maximum and decline after day 50. The main changes in the optical spectrum began at around that time. The He II 4686A line first appeared between days 7 and 14 and initially displayed a broad, symmetric profile that is characteristic of all species before day 60. Low-excitation lines remained comparatively narrow, with v(rad,max) of order 200-400 km/s. They were systematically more symmetric than lines such as [Ca V], [Fe VII], [Fe X], and He II, all of which showed a sequence of profile changes going from symmetric to a blue wing similar to that of the low ionization species but with a red wing extended to as high as 600 km/s . The Na I D doublet developed a broad component with similar velocity width to the other low-ionization species. The O VI Raman features were not detected. We interpret these variations as aspherical expansion of the ejecta within the Mira wind. The blue side is from the shock penetrating into the wind while the red wing is from the low-density periphery. The maximum radial velocities obey power laws, v(rad,max) t^{-n} with n ~ 1/3 for red wing and ~0.8 for the blue. (truncated)

We report the discovery of peculiar features in the optical spectrum of 4C+22.25, a flat spectrum radio quasar at z=0.4183 observed in the SDSS and in a dedicated spectroscopic follow-up from the Nordic Optical Telescope. The Hbeta and Halpha lines show broad profiles (FWHM~12,000 km/s), faint fluxes and extreme offsets (Delta v=8,700+/-1,300 km/s) with respect to the narrow emission lines. These features show no significant variation in a time lag of ~3.1 yr (rest frame). We rule out possible interpretations based on the superposition of two sources or on recoiling black holes, and we discuss the virtues and limitations of a massive black hole binary scenario.

We present results from the first two quarters of a survey to search for pulsations in compact stellar objects with the Kepler spacecraft. The survey sample and the various methods applied in its compilation are described, and spectroscopic observations are presented to separate the objects into accurate classes. From the Kepler photometry we clearly identify nine compact pulsators, and a number of interesting binary stars. Of the pulsators, one shows the strong, rapid pulsations typical for a V361 Hya type sdB variable (sdBV), seven show long-period pulsations characteristic of V1093 Her type sdBVs, and one shows low-amplitude pulsations with both short and long periods. We derive effective temperatures and surface gravities for all the subdwarf B stars in the sample and demonstrate that below the boundary region where hybrid sdB pulsators are found, all our targets are pulsating. For the stars hotter than this boundary temperature a low fraction of strong pulsators (<10 per cent) is confirmed. Interestingly, the short-period pulsator also shows a low-amplitude mode in the long-period region, and several of the V1093 Her pulsators show low amplitude modes in the short-period region, indicating that hybrid behaviour may be common in these stars, also outside the boundary temperature region where hybrid pulsators have hitherto been found.

The dust ejecta of Main-Belt Comet P/2010 A2 (LINEAR) have been observed with several telescopes at the at the Observatorio del Roque de los Muchachos on La Palma, Spain. Application of an inverse dust tail Monte Carlo method to the images of the dust ejecta from the object indicates that a sustained, likely water-ice driven, activity over some eight months is the mechanism responsible for the formation of the observed tail. The total amount of dust released is estimated to be 5E7 kg, which represents about 0.3% of the nucleus mass. While the event could have been triggered by a collision, this cannot be decided from the currently available data.

We present time-resolved optical spectroscopy of V458 Vulpeculae (Nova Vul 2007 No. 1) spread over a period of 15 months starting 301 days after its discovery. Our data reveal radial velocity variations in the HeII {\lambda}5412 and HeII {\lambda}4686 emission lines. A period analysis of the radial velocity curves resulted in a period of 98.09647 \pm 0.00025 min (0.06812255 \pm 0.00000017 d) which we identify with the orbital period of the binary system. V458 Vul is therefore the planetary nebula central binary star with the shortest period known. We explore the possibility of the system being composed of a relatively massive white dwarf (M1 \gsim 1.0 Msun) accreting matter from a post-asymptotic giant branch star which produced the planetary nebula observed. In this scenario, the central binary system therefore underwent two common-envelope episodes. A combination of previous photoionisation modelling of the nebular spectra, post-asymptotic giant branch evolutionary tracks and the orbital period favour a mass of M2 \sim 0.6 Msun for the donor star. Therefore, the total mass of the system may exceed the Chandrasekhar mass, which makes V458 Vul a Type Ia supernova progenitor candidate.

Longtime monitoring of gravitational lens systems is often done using telescopes and recording equipment with modest resolution. Still, it would be interesting to get as much information as possible from the measured lightcurves. From high resolution images we know that the recorded quasar images are often blends and that the corresponding time series are not pure shifted replicas of the source variability. In this paper we will develop an algorithm to unscramble this kind of blended data. The proposed method is based on a simple idea. We use one of the photometric curves, which is supposedly a simple shifted replica of the source curve, to build different artificial combined curves. Then we compare these artificial curves with the blended curves. Proper solutions for a full set of time delays are then obtained by varying free input parameters and estimating statistical distances between the artificial and blended curves. We performed a check of feasibility and applicability of the new algorithm. For numerically generated data sets the time delay systems were recovered for a wide range of setups. Application of the new algorithm to the classical double quasar QSO 0957+561 A,B lightcurves shows a clear splitting of one of the images. This is an unexpected result and extremely interesting, especially in the context of the recent controversy about the exact time delay value for the system. The proposed method allows us to properly analyse the data from low resolution observations that have long time coverages. There are a number of gravitational lens monitoring programmes that can make use of the new algorithm.