This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We present new ATCA multi-wavelength radio measurements (range 2.1-21.2 GHz) of the early-type magnetic star rho Oph A, performed in March 2019 during 3 different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from rho Oph A at 2.1 GHz. We collected high-resolution optical spectra of rho Oph A acquired by several instruments over a time span of about ten years. We also report new magnetic field measurements of rho Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the HeI line (lambda=5015 Angstrom), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of rho Oph A's magnetosphere. Past XMM measurements showed periodic X-ray pulses from rho Oph A. We correlate the X-ray light curve with the magnetic field geometry of rho Oph A. The already published XMM data have been re-analyzed showing that the X-ray spectra of rho Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of rho Oph A as having an auroral origin.

We present a catalogue of a large sample of extended radio sources in the SCORPIO field, observed and resolved by the Australia Telescope Compact Array. SCORPIO, a pathfinder project for addressing the early operations of the Australia SKA Pathfinder, is a survey of ~5 square degrees between 1.4 and 3.1 GHz, centered at l=343.5{\deg}, b=0.75{\deg} and with an angular resolution of about 10 arcsec. It is aimed at understanding the scientific and technical challenges to be faced by future Galactic surveys. With a mean sensitivity around 100 $\mu$Jy/beam and the possibility to recover angular scales at least up to 4 arcmin, we extracted 99 extended sources, 35 of them detected for the first time. Among the 64 known sources 55 had at least a tentative classification in literature. Studying the radio morphology and comparing the radio emission with infrared we propose as candidates 6 new H II regions, 2 new planetary nebulae, 2 new luminous blue variable or Wolf--Rayet stars and 3 new supernova remnants. This study provides an overview of the potentiality of future radio surveys in terms of Galactic source extraction and characterization and a discussion on the difficulty to reduce and analyze interferometric data on the Galactic plane.

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

A new era in radioastronomy will begin with the upcoming large-scale surveys planned at the Australian Square Kilometre Array Pathfinder (ASKAP). ASKAP started its Early Science program in October 2017 and several target fields were observed during the array commissioning phase. The SCORPIO field was the first observed in the Galactic Plane in Band 1 (792-1032 MHz) using 15 commissioned antennas. The achieved sensitivity and large field of view already allow to discover new sources and survey thousands of existing ones with improved precision with respect to previous surveys. Data analysis is currently ongoing to deliver the first source catalogue. Given the increased scale of the data, source extraction and characterization, even in this Early Science phase, have to be carried out in a mostly automated way. This process presents significant challenges due to the presence of extended objects and diffuse emission close to the Galactic Plane. In this context we have extended and optimized a novel source finding tool, named CAESAR , to allow extraction of both compact and extended sources from radio maps. A number of developments have been done driven by the analysis of the SCORPIO map and in view of the future ASKAP Galactic Plane survey. The main goals are the improvement of algorithm performances and scalability as well as of software maintainability and usability within the radio community. In this paper we present the current status of CAESAR and report a first systematic characterization of its performance for both compact and extended sources using simulated maps. Future prospects are discussed in light of the obtained results.

The luminous blue variable (LBV) RMC143 is located in the outskirts of the 30~Doradus complex, a region rich with interstellar material and hot luminous stars. We report the $3\sigma$ sub-millimetre detection of its circumstellar nebula with ALMA. The observed morphology in the sub-millimetre is different than previously observed with HST and ATCA in the optical and centimetre wavelength regimes. The spectral energy distribution (SED) of RMC143 suggests that two emission mechanisms contribute to the sub-mm emission: optically thin bremsstrahlung and dust. Both the extinction map and the SED are consistent with a dusty massive nebula with a dust mass of $0.055\pm0.018~M_{\odot}$ (assuming $\kappa_{850}=1.7\rm\,cm^{2}\,g^{-1}$). To date, RMC143 has the most dusty LBV nebula observed in the Magellanic Clouds. We have also re-examined the LBV classification of RMC143 based on VLT/X-shooter spectra obtained in 2015/16 and a review of the publication record. The radiative transfer code CMFGEN is used to derive its fundamental stellar parameters. We find an effective temperature of $\sim 8500$~K, luminosity of log$(L/L_{\odot}) = 5.32$, and a relatively high mass-loss rate of $1.0 \times 10^{-5}~M_{\odot}$~yr$^{-1}$. The luminosity is much lower than previously thought, which implies that the current stellar mass of $\sim8~M_{\odot}$ is comparable to its nebular mass of $\sim 5.5~M_{\odot}$ (from an assumed gas-to-dust ratio of 100), suggesting that the star has lost a large fraction of its initial mass in past LBV eruptions or binary interactions. While the star may have been hotter in the past, it is currently not hot enough to ionize its circumstellar nebula. We propose that the nebula is ionized externally by the hot stars in the 30~Doradus star-forming region.

CU Vir has been the first main sequence star that showed regular radio pulses that persist for decades, resembling the radio lighthouse of pulsars and interpreted as auroral radio emission similar to that found in planets. The star belongs to a rare group of magnetic chemically peculiar stars with variable rotational period. We study the ultraviolet (UV) spectrum of CU Vir obtained using STIS spectrograph onboard the Hubble Space Telescope (HST) to search for the source of radio emission and to test the model of the rotational period evolution. We used our own far-UV and visual photometric observations supplemented with the archival data to improve the parameters of the quasisinusoidal long-term variations of the rotational period. We predict the flux variations of CU Vir from surface abundance maps and compare these variations with UV flux distribution. We searched for wind, auroral, and interstellar lines in the spectra. The UV and visual light curves display the same long-term period variations supporting their common origin. New updated abundance maps provide better agreement with the observed flux distribution. The upper limit of the wind mass-loss rate is about $10^{-12}\,M_\odot\,\rm{yr}^{-1}$. We do not find any auroral lines. We find rotationally modulated variability of interstellar lines, which is most likely of instrumental origin. Our analysis supports the flux redistribution from far-UV to near-UV and visual domains originating in surface abundance spots as the main cause of the flux variability in chemically peculiar stars. Therefore, UV and optical variations are related and the structures leading to these variations are rigidly confined to the stellar surface. The radio emission of CU Vir is most likely powered by a very weak presumably purely metallic wind, which leaves no imprint in spectra.

We report the detection of the auroral radio emission from the early-type magnetic star HD142301. New VLA observations of HD142301 detected highly polarized amplified emission occurring at fixed stellar orientations. The coherent emission mechanism responsible for the stellar auroral radio emission amplifies the radiation within a narrow beam, making the star where this phenomenon occurs similar to a radio lighthouse. The elementary emission process responsible for the auroral radiation mainly amplifies one of the two magneto-ionic modes of the electromagnetic wave. This explains why the auroral pulses are highly circularly polarized. The auroral radio emission of HD142301 is characterized by a reversal of the sense of polarization as the star rotates. The effective magnetic field curve of HD142301 is also available making it possible to correlate the transition from the left to the right-hand circular polarization sense (and vice-versa) of the auroral pulses with the known orientation of the stellar magnetic field. The results presented in this letter have implications for the estimation of the dominant magneto-ionic mode amplified within the HD142301 magnetosphere.

At radio wavelengths, solar-type stars emit thermal free-free and gyroresonance, gyrosynchrotron, and impulsive coherent emission. Thermal free-free emission originates at layers where the optical depth is close to unit, while high brightness temperature, variable emission, can be due to flares via gyrosynchrotron emission. We observed the alpha Cen system with the Australian Telescope Compact Array at 2 GHz for three days and 17 GHz for one day. Both stars have been detected at 17 GHz, while only an upper limit has been obtained at low frequency despite the longer integration time. The brightness temperatures are consistent with the temperature of the upper chromosphere of the Sun. Inverting the formulae of the free-free emission, the average electron density of the plasma has been inferred. The same procedure was applied to the data in the millimetre recently acquired with ALMA. A comparison with the atmospheric solar models reveals a higher level of activity in alpha Cen B rather than in alpha Cen A, even if still at quiescent level. The non detection at 2 GHz allows us to put a lower limit in the filling factor of active regions. The claimed detection of an Earth size planet in close orbit to alpha Cen B, although doubtful, opens the opportunity to check the existence of Star-Planet Magnetic Interaction (SPMI). We constructed dynamic spectra in the 1.3 - 2.9 GHz of the 2 stars to search for time-variable coherent emission but obtained a null result.

We present X-ray observations of CU Vir performed with XMM-Newton and Chandra. With Lx = 3 x 10^28 erg/s the source is moderately X-ray bright, but its spectrum is extremely hard compared to other Ap stars. Spectral modelling requires multi-component models with predominant hot plasma at temperatures of about Tx = 25 MK or, alternatively, a nonthermal spectral component. The Chandra observations was performed six years later than the one by XMM-Newton, yet the source has similar X-ray flux and spectrum, suggesting a steady and persistent X-ray emission. To explain its full X-ray properties, a generating mechanism beyond standard explanations like the presence of a low-mass companion or magnetically confined wind-shocks is required. Magnetospheric activity might be present or, as proposed for fast rotating strongly magnetic Bp stars, the X-ray emission of CU Vir is predominantly auroral in nature.

We present new radio/millimeter measurements of the hot magnetic star HR5907 obtained with the VLA and ALMA interferometers. We find that HR5907 is the most radio luminous early type star in the cm-mm band among those presently known. Its multi-wavelength radio light curves are strongly variable with an amplitude that increases with radio frequency. The radio emission can be explained by the populations of the non-thermal electrons accelerated in the current sheets on the outer border of the magnetosphere of this fast rotating magnetic star. We classify HR5907 as another member of the growing class of strongly magnetic fast rotating hot stars where the gyro-synchrotron emission mechanism efficiently operates in their magnetospheres. The new radio observations of HR5907 are combined with archival X-ray data to study the physical condition of its magnetosphere. The X-ray spectra of HR5907 show tentative evidence for the presence of non-thermal spectral component. We suggest that non-thermal X-rays originate a stellar X-ray aurora due to streams of non-thermal electrons impacting on the stellar surface. Taking advantage of the relation between the spectral indices of the X-ray power-law spectrum and the non-thermal electron energy distributions, we perform 3-D modeling of the radio emission for HR5907. The wavelength-dependent radio light-curves probe magnetospheric layers at different heights above the stellar surface. A detailed comparison between simulated and observed radio light-curves leads us to conclude that the stellar magnetic field of HR5907 is likely non-dipolar, providing further indirect evidence of the complex magnetic field topology of HR5907.

In this paper, we present the first extended catalogue of far-infrared fluxes of Galactic bubbles. Fluxes were estimated for 1814 bubbles, defined here as the `golden sample', and were selected from the Milky Way Project First Data Release (Simpson et al.) The golden sample was comprised of bubbles identified within the Wide-field Infrared Survey Explorer (WISE) dataset (using 12- and 22-$\mu$m images) and Herschel data (using 70-, 160-, 250-, 350- and 500-$\mu$m wavelength images). Flux estimation was achieved initially via classical aperture photometry and then by an alternative image analysis algorithm that used active contours. The accuracy of the two methods was tested by comparing the estimated fluxes for a sample of bubbles, made up of 126 H II regions and 43 planetary nebulae, which were identified by Anderson et al. The results of this paper demonstrate that a good agreement between the two was found. This is by far the largest and most homogeneous catalogue of infrared fluxes measured for Galactic bubbles and it is a step towards the fully automated analysis of astronomical datasets.

In the next few years the classification of radio sources observed by the large surveys will be a challenging problem, and spectral index is a powerful tool for addressing it. Here we present an algorithm to estimate the spectral index of sources from multiwavelength radio images. We have applied our algorithm to SCORPIO (Umana et al. 2015), a Galactic Plane survey centred around 2.1 GHz carried out with ATCA, and found we can measure reliable spectral indices only for sources stronger than 40 times the rms noise. Above a threshold of 1 mJy, the source density in SCORPIO is 20 percent greater than in a typical extra-galactic field, like ATLAS (Norris et al. 2006), because of the presence of Galactic sources. Among this excess population, 16 sources per square degree have a spectral index of about zero, suggesting optically thin thermal emission such as Hii regions and planetary nebulae, while 12 per square degree present a rising spectrum, suggesting optically thick thermal emission such as stars and UCHii regions.

Direct measurements of the stellar magnetic fields are based on the splitting of spectral lines into polarized Zeeman components. With few exceptions, Zeeman signatures are hidden in data noise and a number of methods have been developed to measure the average, over the visible stellar disk, of longitudinal components of the magnetic field. As to faint stars, at present observable only with low resolution spectropolarimetry, a method is based on the regression of the Stokes V signal against the first derivative of Stokes I. Here we present an extension of this method to obtain a direct measurement of the transverse component of stellar magnetic fields by the regression of high resolution Stokes Q and U as a function of the second derivative of Stokes I. We also show that it is possible to determine the orientation in the sky of the rotation axis of a star on the basis of the periodic variability of the transverse component due to its rotation. The method is applied to data, obtained with the Catania Astrophysical Observatory Spectropolarimeter, along the rotational period of the well known magnetic star \b{eta} CrB.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

We present a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.

We present ALMA and ATCA observations of the luminous blue variable \rmc. The radio maps show for the first time the core of the nebula and evidence that the nebula is strongly asymmetric with a Z-pattern shape. Hints of this morphology are also visible in the archival \emph{HST} $\rm H\alpha$ image, which overall resembles the radio emission. The emission mechanism in the outer nebula is optically thin free-free in the radio. At high frequencies, a component of point-source emission appears at the position of the star, up to the ALMA frequencies. The rising flux density distribution ($S_{\nu}\sim \nu^{0.78\pm0.05}$) of this object suggests thermal emission from the ionized stellar wind and indicates a departure from spherical symmetry with $n_{e}(r)\propto r^{-2}$. We examine different scenarios to explain this excess of thermal emission from the wind and show that this can arise from a bipolar outflow, supporting the suggestion by other authors that the stellar wind of \rmc is aspherical. We fit the data with two collimated ionized wind models and we find that the mass-loss rate can be a factor of two or more smaller than in the spherical case. We also fit the photometry obtained by IR space telescopes and deduce that the mid- to far-IR emission must arise from extended, cool ($\sim80\,\rm K$) dust within the outer ionized nebula. Finally we discuss two possible scenarios for the nebular morphology: the canonical single star expanding shell geometry, and a precessing jet model assuming presence of a companion star.

In this paper we simulate the cyclic circularly-polarised pulses of the ultra-cool dwarf TVLM513-46546, observed with the VLA at 4.88 and 8.44 GHz on May 2006, by using a 3D model of the auroral radio emission from the stellar magnetosphere. During this epoch, the radio light curves are characterised by two pulses left-hand polarised at 4.88 GHz, and one doubly-peaked (of opposite polarisations) pulse at 8.44 GHz. To take into account the possible deviation from the dipolar symmetry of the stellar magnetic field topology, the model described in this paper is also able to simulate the auroral radio emission from a magnetosphere shaped like an offset-dipole. To reproduce the timing and pattern of the observed pulses, we explored the space of parameters controlling the auroral beaming pattern and the geometry of the magnetosphere. Through the analysis of the TVLM513-46546 auroral radio emission, we derive some indications on the magnetospheric field topology that is able to simultaneously reproduce the timing and patterns of the auroral pulses measured at 4.88 and 8.44 GHz. Each set of model solutions simulates two auroral pulses (singly or doubly peaked) per period. To explain the presence of only one 8.44 GHz pulse per period, we analyse the case of auroral radio emission limited only to a magnetospheric sector activated by an external body, like the case of the interaction of Jupiter with its moons.

We present new JVLA multi-frequency measurements of a set of stars in transition from the post-AGB to the Planetary Nebula phase monitored in the radio range over several years. Clear variability is found for five sources. Their light curves show increasing and decreasing patterns. New radio observations at high angular resolution are also presented for two sources. Among these is IRAS 18062+2410, whose radio structure is compared to near-infrared images available in the literature. With these new maps, we can estimate inner and outer radii of 0.03$"$ and 0.08$"$ for the ionised shell, an ionised mass of $3.2\times10^{-4}$ M$_\odot$, and a density at the inner radius of $7.7\times 10^{-5}$ cm$^{-3}$, obtained by modelling the radio shell with the new morphological constraints. The combination of multi-frequency data and, where available, spectral-index maps leads to the detection of spectral indices not due to thermal emission, contrary to what one would expect in planetary nebulae. Our results allow us to hypothesise the existence of a link between radio variability and non-thermal emission mechanisms in the nebulae. This link seems to hold for IRAS 22568+6141 and may generally hold for those nebulae where the radio flux decreases over time.

In this paper we investigate the multiwavelengths properties of the magnetic early B-type star HR7355. We present its radio light curves at several frequencies, taken with the Jansky Very Large Array, and X-ray spectra, taken with the XMM X-ray telescope. Modeling of the radio light curves for the Stokes I and V provides a quantitative analysis of the HR7355 magnetosphere. A comparison between HR7355 and a similar analysis for the Ap star CUVir, allows us to study how the different physical parameters of the two stars affect the structure of the respective magnetospheres where the non-thermal electrons originate. Our model includes a cold thermal plasma component that accumulates at high magnetic latitudes that influences the radio regime, but does not give rise to X-ray emission. Instead, the thermal X-ray emission arises from shocks generated by wind stream collisions close to the magnetic equatorial plane. The analysis of the X-ray spectrum of HR7355 also suggests the presence of a non-thermal radiation. Comparison between the spectral index of the power-law X-ray energy distribution with the non-thermal electron energy distribution indicates that the non-thermal X-ray component could be the auroral signature of the non-thermal electrons that impact the stellar surface, the same non-thermal electrons that are responsible for the observed radio emission. On the basis of our analysis, we suggest a novel model that simultaneously explains the X-ray and the radio features of HR7355 and is likely relevant for magnetospheres of other magnetic early type stars.

We present new observations of the nebula around the Magellanic candidate Luminous Blue Variable S61. These comprise high-resolution data acquired with the Australia Telescope Compact Array (ATCA), the Atacama Large Millimetre/Submillimetre Array (ALMA), and VISIR at the Very Large Telescope (VLT). The nebula was detected only in the radio, up to 17 GHz. The 17 GHz ATCA map, with 0.8 arcsec resolution, allowed a morphological comparison with the H$\alpha$ Hubble Space Telescope image. The radio nebula resembles a spherical shell, as in the optical. The spectral index map indicates that the radio emission is due to free-free transitions in the ionised, optically thin gas, but there are hints of inhomogeneities. We present our new public code RHOCUBE to model 3D density distributions, and determine via Bayesian inference the nebula's geometric parameters. We applied the code to model the electron density distribution in the S61 nebula. We found that different distributions fit the data, but all of them converge to the same ionised mass, ~0.1 $\rm M\odot$, which is an order of magnitude smaller than previous estimates. We show how the nebula models can be used to derive the mass-loss history with high-temporal resolution. The nebula was probably formed through stellar winds, rather than eruptions. From the ALMA and VISIR non-detections, plus the derived extinction map, we deduce that the infrared emission observed by space telescopes must arise from extended, diffuse dust within the ionised region.

We present a multi-wavelength study of the Galactic Luminous Blue Variable HR Carinae, based on new high resolution mid-infrared (IR) and radio images obtained with the Very Large Telescope (VLT) and the Australia Telescope Compact Array (ATCA), which have been complemented by far-infrared Herschel-PACS observations and ATCA archive data. The Herschel images reveal the large-scale distribution of the dusty emitting nebula, which extends mainly to the North-East direction, up to 70 arcsec from the central star, and is oriented along the direction of the space motion of the star. In the mid-infrared images, the brightness distribution is characterized by two arc-shaped structures, tracing an inner envelope surrounding the central star more closely. At radio wavelengths, the ionized gas emission lies on the opposite side of the cold dust with respect to the position of the star, as if the ionized front was confined by the surrounding medium in the North-South direction. The comparison with previous data indicates significant changes in the radio nebula morphology and in the mass-loss rate from the central star, which has increased from 6.1$\times10^{-6}$ M$_{\odot}$yr$^{-1}$ in 1994-1995 to $1.17\times10^{-5}$ M$_{\odot}$yr$^{-1}$ in 2014. We investigate possible scenarios that could have generated the complex circumstellar environment revealed by our multi-wavelength data.

We present radio observations of 18 MIPSGAL bubbles performed at 5 GHz (6 cm) with the Karl G. Jansky Very Large Array in configuration B and BnA. The observations were aimed at understanding what kind of information high-resolution and high-sensitivity radio maps can supply on the circumstellar envelopes of different kinds of evolved stars and what their comparison with infrared images with similar resolution can tell us. We found that the 18 bubbles can be grouped into five categories according to their radio morphology. The three bubbles presenting a central point source in the radio images all correspond to luminous blue variable star candidates. Eleven bubbles show an elliptical shape and the total lack of a central object in the radio, and are likely associated with planetary nebulae. Under this assumption we derive their distance, their ionized mass and their distribution on the Galactic plane. We discuss the possibility that the MIPSGAL bubbles catalogue (428 objects) may contain a large fraction of all Galactic planetary nebulae located at a distance between 1.4 kpc and 6.9 kpc and lying in the MIPSGAL field of view. Among the remaining bubbles we identify also a H II region and a proto-planetary nebula candidate.

Automated source extraction and parameterization represents a crucial challenge for the next-generation radio interferometer surveys, such as those performed with the Square Kilometre Array (SKA) and its precursors. In this paper we present a new algorithm, dubbed CAESAR (Compact And Extended Source Automated Recognition), to detect and parametrize extended sources in radio interferometric maps. It is based on a pre-filtering stage, allowing image denoising, compact source suppression and enhancement of diffuse emission, followed by an adaptive superpixel clustering stage for final source segmentation. A parameterization stage provides source flux information and a wide range of morphology estimators for post-processing analysis. We developed CAESAR in a modular software library, including also different methods for local background estimation and image filtering, along with alternative algorithms for both compact and diffuse source extraction. The method was applied to real radio continuum data collected at the Australian Telescope Compact Array (ATCA) within the SCORPIO project, a pathfinder of the ASKAP-EMU survey. The source reconstruction capabilities were studied over different test fields in the presence of compact sources, imaging artefacts and diffuse emission from the Galactic plane and compared with existing algorithms. When compared to a human-driven analysis, the designed algorithm was found capable of detecting known target sources and regions of diffuse emission, outperforming alternative approaches over the considered fields.

The electron cyclotron maser is the coherent emission process that gives rise to the radio lighthouse effect observed in the hot magnetic chemically peculiar star CU Virginis. It has also been proposed to explain the highly circularly polarized radio pulses observed on some ultra cool dwarfs, with spectral type earlier than M7. Such kind of coherent events resemble the auroral radio emission from the magnetized planets of the solar system. In this paper, we present a tridimensional model able to simulate the timing and profile of the pulses emitted by those stars characterized by a dipolar magnetic field by following the hypothesis of the laminar source model, used to explain the beaming of the terrestrial auroral kilometric radiation. This model proves to be a powerful tool to understand the auroral radio-emission phenomenon, allowing us to derive some general conclusions about the effects of the model's free parameters on the features of the coherent pulses, and to learn more about the detectability of such kind of pulsed radio emission.

We conducted radio observations searching for OH 18-cm maser emission from a sample of 169 unclassified MIPSGAL compact Galactic bubbles. These sources are thought to be the circumstellar envelopes of different kinds of evolved stars. Our observations were aimed at shedding light on the nature of MIPSGAL bubbles, since their characterisation is a fundamental aid for the development of accurate physical models of stellar and Galaxy evolution. The maser emission is observatively linked to the last stages of the life of low- and intermediate-mass stars, which may constitute a significant fraction of the MIPSGAL bubbles. In particular OH masers are usually observed towards post-AGB stars. Our observations were performed with the Green Bank Telescope and, for each source, produced spectra around the four OH 18-cm transitions. The observations were compared with archive interferometer data in order to exclude possible contamination from nearby sources. The main result is that the OH maser emission is not a common feature among the MIPSGAL bubbles, with only one certain detection. We conclude that among the MIPSGAL bubbles the post-AGB stars could be very rare.

The detection of an excess of emission at microwave frequencies with respect to the predicted free-free emission has been reported for several Galactic HII regions. Here, we investigate the case of RCW 49, for which the Cosmic Background Imager tentatively (~ 3 sigma) detected Anomalous Microwave Emission at 31 GHz on angular scales of 7'. Using the Australia Telescope Compact Array, we carried out a multi-frequency (5 GHz, 19 GHz and 34 GHz) continuum study of the region, complemented by observations of the H109$\alpha$ radio recombination line. The analysis shows that: 1) the spatial correlation between the microwave and IR emission persists on angular scales from 3.4' to 0.4", although the degree of the correlation slightly decreases at higher frequencies and on smaller angular scales, 2) the spectral indices between 1.4 and 5 GHz are globally in agreement with optically thin free-free emission, however, ~ 30% of these are positive and much greater than -0.1, consistently with a stellar wind scenario, 3) no major evidence for inverted free-free radiation is found, indicating that this is likely not the cause of the Anomalous Emission in RCW 49. Although our results cannot rule out the spinning dust hypothesis to explain the tentative detection of Anomalous Microwave emission in RCW 49, they emphasize the complexity of astronomical sources very well known and studied such as HII regions, and suggest that, at least in these objects, the reported excess of emission might be ascribed to alternative mechanisms such as stellar winds and shocks.

Radio emission has been detected in a broad variety of stellar objects from all stages of stellar evolution. However, most of our knowledge originates from targeted observations of small samples, which are strongly biased to sources which are peculiar at other wavelengths. In order to tackle this problem we have conducted a deep 1.4 GHz survey by using the Australian Telescope Compact Array (ATCA), following the same observing setup as that used for the Australia Telescope Large Area Survey (ATLAS) project, this time choosing a region more appropriate for stellar work. In this paper, the SCORPIO project is presented as well as results from the pilot experiment. The achieved rms is about 30 /uJy and the angular resolution ~10 arcsec. About six hundred of point-like sources have been extracted just from the pilot field. A very small percentage of them are classified in SIMBAD or the NASA/IPAC Extragalactic Database (NED). About 80 % of the extracted sources are reported in one of the inspected catalogues and 50 % of them appears to belong to a reddened stellar/galactic population. The evaluation of extragalactic contaminants is very difficult without further investigations. Interesting results have been obtained for extended radio sources that fall in the SCORPIO field. Many bubble-like structures have been found, some of which are classified at other wavelengths. However, for all of these sources, our project has provided us with images of unprecedented sensitivity and angular resolution at 2.1 GHz.

The SKA will be a state of the art radiotelescope optimized for both large area surveys as well as for deep pointed observations. In this paper we analyze the impact that the SKA will have on Galactic studies, starting from the immense legacy value of the all-sky survey proposed by the continuum SWG but also presenting some areas of Galactic Science that particularly benefit from SKA observations both surveys and pointed. The planned all-sky survey will be characterized by unique spatial resolution, sensitivity and survey speed, providing us with a wide-field atlas of the Galactic continuum emission. Synergies with existing, current and planned radio Galactic Plane surveys will be discussed. SKA will give the opportunity to create a sensitive catalog of discrete Galactic radio sources, most of them representing the interaction of stars at various stages of their evolution with the environment: complete census of all stage of HII regions evolution; complete census of late stages of stellar evolution such as PNe and SNRs; detection of stellar winds, thermal jets, Symbiotic systems, Chemically Peculiar and dMe stars, active binary systems in both flaring and quiescent states. Coherent emission events like Cyclotron Maser in the magnetospheres of different classes of stars can be detected. Pointed, deep observations will allow new insights into the physics of the coronae and plasma processes in active stellar systems and single stars, enabling the detection of flaring activity in larger stellar population for a better comprehension of the mechanism of energy release in the atmospheres of stars with different masses and age.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

Late stages of stellar evolution are characterized by copious mass-loss events whose signature is the formation of circumstellar envelopes (CSE). Planck multi-frequency measurements have provided relevant information on a sample of Galactic planetary nebulae (PNe) in the important and relatively unexplored observational band between 30 and 857GHz. Planck enables the assembly of comprehensive PNe spectral energy distributions (SEDs) from radio {\bf to} far-infrared frequencies. Modelling of the derived SEDs provides us with information on physical properties of CSEs and the mass content of both main components: ionised gas, traced by the free-free emission at cm--mm waves; and thermal dust, traced by the millimetre and far-IR emission. In particular, the amount of ionised gas and dust has been derived here. Such quantities have also been estimated for the very young PN CRL618, where the strong variability observed in its radio and millimetre emission has previously prevented the construction of its SED. A morphological study of the Helix Nebula has also been performed. Planck maps reveal, for the first time, the spatial distribution of the dust inside the envelope, allowing us to identify different components, the most interesting of which is a very extended component (up to 1pc) that may be related to a region where the slow expanding envelope is interacting with the surrounding interstellar medium.

We present radio observations of two supernova remnants and four planetary nebulae with the Very Large Array and the Green Bank Telescope. These objects are part of a larger sample of radio sources, discussed in a previous paper, counterpart of the MIPSGAL 24-micron compact bubbles. For the two supernova remnants we combined the interferometric observations with single-dish data to obtain both a high resolution and a good sensitivity to extended structures. We discuss in detail the entire combination procedure adopted and the reliability of the resulting maps. For one supernova remnant we pose a more stringent upper limit for the flux density of its undetected pulsar, and we also show prominent spectral index spatial variations, probably due to inhomogeneities in the magnetic field and in its ejecta or to an interaction between the supernova shock and molecular clouds. We eventually use the 5-GHz maps of the four planetary nebulae to estimate their distance and their ionized mass.

We present a multiwavelength analysis of the nebula around the candidate luminous blue variable G79.29+0.46. The study is based on our radio observations performed at the Expanded Very Large Array and at the Green Bank Telescope and on archival infrared datasets, including recent images obtained by the Herschel Space Observatory. We confirm that the radio central object is characterized by a stellar wind and derive a current mass-loss rate of about 1.4x10-6 Msun yr-1. We find the presence of a dusty compact envelope close to the star, with a temperature between 40 and 1200 K. We estimate for the outer ejecta an ionised gas mass of 1.51 Msun and a warm (60--85 K) dust mass of 0.02 Msun. Diagnostics of the far-infrared spectra indicate the presence of a photo-dissociation region around the ionised gas. Finally, we model the nebula with the photo-ionization code CLOUDY, using as input parameters those estimated from our analysis. We find for the central star a luminosity of 10^5.4 Lsun and an effective temperature of 20.4 kK.

MWC 930 is a star just ~2{\deg} above the Galactic plane whose nature is not clear and that has not been studied in detail so far. While a post-Asymptotic Giant Branch (AGB) classification was proposed in the past, studies of its optical spectrum and photometry pointed toward strong variability, therefore the object was reclassified as a Luminous Blue Variable (LBV) candidate. LBVs typically undergo phases of strong mass loss in the form of eruptions that can create shells of ejecta around the star. Our goal is to search for the presence of such a circumstellar nebula in MWC 930 and investigate its properties. To do so, we make use of space-based infrared data from our Spitzer campaign performed with the InfraRed Array Camera (IRAC) and the InfraRed Spectrograph (IRS) as well as data from optical and infrared (IR) surveys. In our Spitzer images, we clearly detect an extended shell around MWC 930 at wavelengths longer than 5 um. The mid-infrared spectrum is dominated by the central star and mostly shows forbidden lines of [FeII], with an underlying continuum that decreases with wavelength up to ~15 um and then inverts its slope, displaying a second peak around 60 um, evidence for cold dust grains formed in a past eruption. By modeling the SED, we identify two central components, besides the star and the outer shell. These extra sources of radiation are interpreted as material close to the central star, maybe due to a recent ejection. Features of C-bearing molecules or grains are not detected.

We report the radio observations of a sub-sample of the 428 galactic compact bubbles discovered at 24 $\mu$m with the MIPSGAL survey. Pervasive through the entire Galactic plane, these objects are thought to be different kinds of evolved stars. The very large majority of the bubbles (~ 70%) are however not yet classified. We conducted radio observations with the EVLA at 6 cm and 20 cm in order to obtain the spectral index of 55 bubbles. We found that at least 70 per cent of the 31 bubbles for which we were effectively able to compute the spectral index (or its lower limit) are likely to be thermal emitters. We were also able to resolve some bubbles, obtaining that the size of the radio nebula is usually similar to the IR size, although our low resolution (with respect to IR images) did not allow further morphological studies. Comparisons between radio flux densities and IR archive data from Spitzer and IRAS suggest that at least 3 unclassified bubbles can be treated as planetary nebula candidates.

The nebulae associated to four Luminous Blue Variables (LBVs) in the Large Magellanic Cloud (LMC) have been observed at 5.5 and 9 GHz using the Australia Telescope Compact Array, and radio emission has been detected for first time in these sources, R127, R143, S61 and S119. The radio maps of the nebulae have an angular resolution of \sim 1.5" and a sensitivity of 1.5-3.0\times10-2 mJy beam-1, and show a very similar morphology to that observed in H{\alpha}. This similarity permit us to assume that the H{\alpha} emission is not affected by strong intrinsic extinction due to dust within the nebulae. We estimate the masses of the ionized gas in the LBVs nebulae and their values are consistent with those measured in Galactic LBVs.

The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and \gamma-ray flux apparently fades in time, likely because of an increasing optical to \gamma-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor \delta. Under the hypothesis that in the period 2008-2011 all the \gamma-ray and optical variability on a one-week timescale were due to changes in \delta, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle \theta\ only, then \theta\ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

In this paper we present new and archive radio measurements obtained with the Very Large Array of the magnetic chemically peculiar (MCP) star \sigma Ori E. The radio data have been obtained at different frequencies and are well distributed along the rotational phases. We analyze in detail the radio emission from \sigma Ori E with the aim to search evidence of circularly polarized radio pulses. Up to now, among the MCP stars only CU Virginis shows 100% polarized time-stable radio pulses, explained as highly directive electron cyclotron maser emission, visible from Earth at particular rotational phases, like a pulsar. Our analysis shows that there is no hint of coherent emission at frequencies below 15 GHz. We conclude that the presence of a quadrupolar component of the magnetic field, dominant within few stellar radii from the star, where the maser emission should be generated, inhibits the onset of the cyclotron maser instability in \sigma Ori E.

CU Vir and V901 Ori belong among these few magnetic chemically peculiar stars whose rotation periods vary on timescales of decades. We aim to study the stability of the periods in CU Vir and V901 Ori using all accessible observational data containing phase information. We collected all available relevant archived observations supplemented with our new measurements of these stars and analysed the period variations of the stars using a novel method that allows for the combination of data of diverse sorts. We found that the shapes of their phase curves were constant, while the periods were changing. Both stars exhibit alternating intervals of rotational braking and acceleration. The rotation period of CU Vir was gradually shortening until the year 1968, when it reached its local minimum of 0.52067198 d. The period then started increasing, reaching its local maximum of 0.5207163 d in the year 2005. Since that time the rotation has begun to accelerate again. We also found much smaller period changes in CU Vir on a timescale of several years. The rotation period of V901 Ori was increasing for the past quarter-century, reaching a maximum of 1.538771 d in the year 2003, when the rotation period began to decrease. A theoretically unexpected alternating variability of rotation periods in these stars would remove the spin-down time paradox and brings a new insight into structure and evolution of magnetic upper-main-sequence stars.

CU Virginis is a rapidly rotating Magnetic Chemically Peculiar star with at present unique characteristics as radio emitter. The most intriguing one is the presence of intense, 100% circularly polarized radiation ascribed to Cyclotron Maser. Each time the star rotates, this highly beamed emission points two times toward the Earth, like a pulsar. We observed CU Vir in April 2010 with the EVLA in two bands centered at 1450 and 1850 MHz. We covered nearly the whole rotational period, confirming the presence of the two pulses at a flux density up to 20 mJy. Dynamical spectra, obtained with unprecedented spectral and temporal sensitivity, allow us to clearly see the different time delays as a function of the frequency. We interpret this behaviour as a propagation effect of the radiation inside the stellar magnetosphere. The emerging scenario suggests interesting similarities with the auroral radio emission from planets, in particular with the Auroral Kilometric Radiation (AKR) from Earth, which originates at few terrestrial radii above the magnetic poles and was only recently discovered to be highly beamed. We conclude that the magnetospheres of CU Vir, Earth and other planets, maybe also exoplanets, could have similar geometrical and physical characteristics in the regions where the cyclotron maser is generated. In addition, the pulses are perfect "markers" of the rotation period. This has given us for the first time the possibility to measure with extraordinary accuracy the spin down of a star on or near the main sequence.

We have observed the radio nebula surrounding the Galactic LBV candidate G79.29+0.46 with the EVLA at 6 cm. These new radio observations allow a morphological comparison between the radio emission, which traces the ionized gas component, and the mid-IR emission, a tracer of the dust component. The IRAC (8 \mu m) and MIPS (24 \mu m and 70 \mu m) images have been reprocessed and compared with the EVLA map. We confirm the presence of a second shell at 24 \mu m and also provide evidence for its detection at 70 \mu m. The differences between the spatial morphology of the radio and mid-IR maps indicate the existence of two dust populations, the cooler one emitting mostly at longer wavelengths. Analysis of the two dusty, nested shells have provided us with an estimate of the characteristic timescales for shell ejection, providing important constraints for stellar evolutionary models. Finer details of the ionized gas distribution can be appreciated thanks to the improved quality of the new 6 cm image, most notably the highly structured texture of the nebula. Evidence of interaction between the nebula and the surrounding interstellar medium can be seen in the radio map, including brighter features that delineate regions where the shell structure is locally modified. In particular, the brighter filaments in the south-west region appear to frame the shocked southwestern clump reported from CO observations.

Among its great findings, the IRAS mission showed the existence of an unidentified mid-IR feature around 21 um. Since its discovery, this feature has been detected in all C-rich proto-PNe of intermediate spectral type (A-G) and - weakly - in a few PNe and AGB stars, but the nature of its carriers remains unknown. In this paper, we show the detection of this feature in the spectra of three new stars transiting from the AGB to the PN stage obtained with the Spitzer Space Telescope. Following a recent suggestion, we try to model the SEDs of our targets with amorphous carbon and FeO, which might be responsible for the unidentified feature. The fit thus obtained is not completely satisfactory, since the shape of the feature is not well matched. In the attempt to relate the unidentified feature to other dust features, we retrieved mid-IR spectra of all the 21-um sources currently known from ISO and Spitzer on-line archives and noticed a correlation between the flux emitted in the 21-um feature and that emitted at 7 and 11 um (PAH bands and HAC broad emission). Such a correlation may point to a common nature of the carriers.

The Simultaneous Medicina-{\it Planck} Experiment (SiMPlE) is aimed at observing a selected sample of 263 extragalactic and Galactic sources with the Medicina 32-m single dish radio telescope in the same epoch as the Planck satellite observations. The data acquired with a frequency coverage down to 5 GHz, also combined with Planck at frequencies above 30 GHz, will constitute a useful reference catalogue of bright sources over the whole Northern hemisphere. Furthermore, source observations performed in different epochs and comparison with other catalogues allow the investigation of source variabilities on different timescales. In this work, we describe the sample selection, the on-going data acquisition campaign, the data reduction procedures, the developed tools, and the comparison with other data-sets. We present the data at 5 and 8.3 GHz for the SiMPlE Northern Sample consisting of 79 sources with $\delta \geq 45^\circ$ selected in our catalogue and observed during the first 6 months of the project. A first analysis of their spectral behaviour and long-term variability is also presented.

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak.

Photometric variability of chemically peculiar (CP) stars of the upper main sequence is closely connected to their local stellar magnetic field and their rotational period. Long term investigations, as presented here, help us to identify possible stellar cycles (as in the Sun). Furthermore, these data provide a basis for detailed surface mapping techniques. Photoelectric Stroemgren uvby time series for 27 CP stars within the boundaries of open clusters are presented. In addition, Hipparcos photometric data (from 1989 to 1993) are used for our analysis. Our observations cover a time period of about six years (1986 to 1992) with typically fifteen measurements for each objects. These observations help us to determine the rotational periods of these objects. A standard reduction procedure was applied to the data. When possible, we merged our data sets with already published ones to obtain a more significant result. A detailed time series analysis was performed, involving five different methods to minimize spurious detections. We established, for the first time, variability for fourteen CP stars. For additional two stars, a merging of already published data sets, resulted in more precise periods, whereas for six objects, the published periods could be confirmed. Last, but not least, no significant variations were found for five stars. Apart from six stars, all targets seem to be members of their host open clusters.

We present first results of our study of a sample of Galactic LBV, aimed to contribute to a better understanding of the LBV phenomenon, by recovering the mass-loss history of the central object from the analysis of its associated nebula. Mass-loss properties have been derived by a synergistic use of different techniques, at different wavelengths, to obtain high-resolution, multi-wavelength maps, tracing the different emitting components coexisting in the stellar ejecta: the ionized/neutral gas and the dust. Evidence for asymmetric mass-loss and observational evidence of possible mutual interaction between gas and dust components have been observed by the comparison of mid-IR (Spitzer/IRAC, VLT/VISIR) and radio (VLA) images of the nebulae, while important information on the gas and dust composition have been derived from Spitzer/IRS spectra.

We present some preliminary results from our program of intensive near-infrared photometric monitoring of a sample of confirmed and candidate Luminous Blue Variables (LBVs) conducted from 2008 to 2010. Clear long-term variability has been observed for Wray 17-96 and V481 Sct, with overall brightness variation greater than 1 mag in the J band. Other sources, such as LBV 1806-20 showed detectable variability with amplitudes of few tenths of a magnitude with time-scale of about 60 days.

Searching for variability, we have observed a sample of hot post-AGB stars and young Planetary Nebulae candidates with the Very Large Array at 4.8, 8.4, and 22.4 GHz. The sources had been previously detected in the radio continuum, which is a proof that the central stars have started ionising their circumstellar envelopes and an increase in radio flux with time can be expected as a result of the progression of the ionisation front. Such a behaviour has been found in IRAS 18062+2410, whose radio modelling has allowed us to determine that its ionised mass has increased from 10^{-4} to 3.3 10^{-4} M_sun in 8 years and its envelope has become optically thin at lower frequencies. Different temporal behaviours have been found for three other sources. IRAS 17423-1755 has shown a possibly periodic pattern and an inversion of its radio spectral index, as expected from a varying stellar wind. We estimate that the radio flux arises from a very compact region around the central star (10^{15} cm) with an electron density of 2 10^6 cm^{-3}. IRAS 22568+6141 and 17516-2525 have decreased their radio flux densities of about 10% per year over 4 years. While a linear increase of the flux density with time points out to the progression of the ionisation front in the envelope, decreases as well as quasi-periodic patterns may indicate the presence of unstable stellar winds/jets or thick dusty envelopes absorbing ionising photons.

The transition from Asymptotic Giant Branch star to Planetary Nebula is short-lived and mysterious. Though it lasts only a few thousand years, it is thought to be the time when the asymmetries observed in subsequent phases arise. During this epoch, the star is shrouded in thick clouds of dust and molecular gas; infrared observations are needed to reveal these objects at their most pivotal moment. I present preliminary results of a Spitzer study of targets spanning the range from post-AGB stars to Planetary Nebulae with the goal of determining the genesis of asymmetry in these objects.

The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 at radio, near-IR, and optical frequencies. During this period, high-energy observations were performed by XMM-Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to interpret the source broad-band variability. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, as well as moderate flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the line of sight, producing an increase of the Doppler beaming factor. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability.

High spatial and sensitivity images of the Luminous Blue Variable IRAS 18576+0341 were obtained using the mid infrared imager VISIR at the Very Large Telescope and the Very Large Array interferometer. The resulting mid-infrared continuum maps show a similar clumpy and approximately circular symmetric nebula, which contrasts sharply with the asymmetry that characterizes the ionized component of the envelope, as evidenced from the radio and [Ne II] line images obtained with comparable spatial resolution. In particular, there is excellent overall agreement between the 12.8 micron map and the radio images, consistent with free-free emission from circumstellar ionized material surrounding a central stellar wind. The color temperature and optical depth maps obtained from mid-infrared images show only slight fluctuations, suggesting quite uniform dust characteristics over the dust shell. We explore various possibilities to understand the cause of the different morphology of the dusty and gaseous component of the circumstellar envelope which are compatible with the observations.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

Ultra compact HII (UCHII) regions are indicators of high-mass star formation sites and are distributed mainly in the Galactic plane. We intend to investigate the possible contribution of the forthcoming ESA Planck mission to the science of UCHII regions by evaluating the possibility of detecting UCHIIs that are bright in the radio regime. We performed new 7 mm observations of a sample of UCHII regions. For each source in our sample, the free-free radio spectrum has been modeled. Along with far-IR measurements, our spectra allow us to estimate the flux densities of the sources in the millimeter and sub-millimeter bands. The possibility of Planck detecting the selected sources can be assessed by comparing the estimated flux densities to the expected sensitivity in each Planck channel. We conclude that, in the case of the present sample, located close to the Galactic center, Planck will have a very low detection rate. In contrast, assuming that our sample is representative of the whole UCHII-region population, we derive a very high probability of detecting this kind of source with Planck if located instead close to the anticenter.

We report the detection by the AGILE satellite of a rapid gamma-ray flare from the powerful gamma-ray quasar PKS 1510-089, during a pointing centered on the Galactic Center region from 1 March to 30 March 2008. This source has been continuosly monitored in the radio-to-optical bands by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT). Moreover, the gamma-ray flaring episode triggered three ToO observations by the Swift satellite in three consecutive days, starting from 20 March 2008. In the period 1-16 March 2008, AGILE detected gamma-ray emission from PKS 1510-089 at a significance level of 6.2-sigma with an average flux over the entire period of (84 +/- 17) x 10^{-8} photons cm^{-2} s^{-1} for photon energies above 100 MeV. After a predefined satellite re-pointing, between 17 and 21 March 2008, AGILE detected the source at a significance level of 7.3-sigma, with an average flux (E > 100 MeV) of (134 +/- 29) x 10^{-8} photons cm^{-2} s^{-1} and a peak level of (281 +/- 68) x 10^{-8} photons cm^{-2} s^{-1} with daily integration. During the observing period January-April 2008, the source also showed an intense and variable optical activity, with several flaring episodes and a significant increase of the flux was observed at millimetric frequencies. Moreover, in the X-ray band the Swift/XRT observations seem to show an harder-when-brighter behaviour of the source spectrum. The spectral energy distribution of mid-March 2008 is modelled with a homogeneous one-zone synchrotron self Compton emission plus contributions from inverse Compton scattering of external photons from both the accretion disc and the broad line region. Indeed, some features in the optical-UV spectrum seem to indicate the presence of Seyfert-like components, such as the little blue bump and the big blue bump.

We have studied the growth of the shell-like radio structure of supernova SN1993J in M81 from September 1993 through October 2003 with very-long-baseline interferometry (VLBI) observations at the wavelengths of 3.6, 6, and 18cm. For this purpose, we have developed a method to accurately determine the outer radius (R) of any circularly symmetric compact radio structure like SN1993J. The source structure of SN1993J remains circularly symmetric (with deviations from circularity under 2%) over almost 4000 days. We characterize the decelerated expansion of SN 1993J through approximately day 1500 after explosion with an expansion parameter $m= 0.845\pm0.005$ ($R \propto t^{m}$). However, from that day onwards the expansion is different when observed at 6 and 18cm. Indeed, at 18cm, the expansion can be well characterized by the same $m$ as before day 1500, while at 6cm the expansion appears more decelerated, and is characterized by another expansion parameter, $m_{6}= 0.788\pm0.015$. Therefore, since about day 1500 on, the radio source size has been progressively smaller at 6cm than at 18cm. These findings are in stark contrast to previous reports by other authors on the details of the expansion. In our interpretation the supernova expands with a single expansion parameter, $m= 0.845\pm0.005$, and the 6cm results beyond day 1500 are due to physical effects, perhaps also coupled to instrumental limitations. Two physical effects may be involved: (a) a changing opacity of the ejecta to the 6cm radiation, and (b) a radial decrease of the magnetic field in the emitting region. (Long abstract cut. Please, read full abstract in manuscript).

We have observed a small sample of hot post-AGB stars with the InfraRed Array Camera (IRAC) and the InfraRed Spectrograph (IRS) on-board the Spitzer Space Telescope. The stars were selected from the literature on the basis of their far-Infrared excess (i.e., post-AGB candidates) and B spectral type (i.e., close to the ionization of the envelope). The combination of our IRAC observations with 2MASS and IRAS catalog data, along with previous radio observations in the cm range (where available) allowed us to model the SEDs of our targets and find that in almost all of them at least two shells of dust at different temperatures must be present, the hot dust component ranging up to 1000 K. In several targets grains larger than 1 micron are needed to match the far-IR data points. In particular, in IRAS 17423-1755 grains up to 100 micron must be introduced to match the emission in the mm range. We obtained IRS spectra to identify the chemistry of the envelopes and found that more than 1/3 of the sources in our sample have mixed chemistry, showing both mid-IR bands attributed to Polycyclic Aromatic Hydrocarbons (PAH) and silicate features. The analysis of the PAH features indicates that these molecules are located in the outflows, far away from the central stars. We consider the larger than expected percentage of mixed-chemistry targets as a selection bias towards stars with a disk or torus around them. Our results strengthen the current picture of mixed chemistry being due to the spatial segregation of different dust populations in the envelopes.

We present mid-IR observations of the Galactic Luminous Blue Variable (LBV) HR Car and its associated nebula carried out with the Spitzer Space Telescope using both IRAC and IRS, as part of a GTO program aimed to study stellar ejecta from evolved stars. Our observations reveal a rich mid-IR spectrum of the inner nebula showing both solid state and atomic gas signatures. Strong low-excitation atomic fine structure lines such as $ 26.0 \mu$m [\ion{Fe}{2}] and $ 34.8 \mu$m [\ion{Si}{2}], indicate, for the first time, the presence of a PDR in this object class. While the physics and chemistry of the low-excitation gas appears to be dominated by photodissociation, a possible contribution due to shocks can be inferred from the evidence of gas phase Fe abundance enhancement. The presence of amorphous silicates, inferred from the observed characteristic broad feature at $10 \mu$m located in the inner nebula, suggests that dust has formed during the LBV outburst. This is in contrast with the detection of crystalline dust in other probably more evolved Galactic LBVs, which is similar to the crystalline dust observed in red supergiants. This has been considered to be evidence of dust production during evolutionary phases prior to the outburst.

We have selected a small sample of post-AGB stars in transition towards the planetary nebula and present new Very Large Array multi-frequency high-angular resolution radio observations of them. The multi-frequency data are used to create and model the targets' radio continuum spectra, proving that these stars started their evolution as very young planetary nebulae. In the optically thin range, the slopes are compatible with the expected spectral index (-0.1). Two targets (IRAS 18062+2410 and 17423-1755) seem to be optically thick even at high frequency, as observed in a handful of other post-AGB stars in the literature, while a third one (IRAS 20462+3416) shows a possible contribution from cold dust. In IRAS 18062+2410, where we have three observations spanning a period of four years, we detect an increase in its flux density, similar to that observed in CRL 618. High-angular resolution imaging shows bipolar structures that may be due to circumstellar tori, although a different hypothesis (i.e., jets) could also explain the observations. Further observations and monitoring of these sources will enable us to test the current evolutionary models of planetary nebulae.

We have performed new laboratory experiments which gave us the possibility to obtain an estimate of the amount of carbon chain oxides (namely C3O2, C2O, and C3O) formed after irradiation (with 200 keV protons) of pure CO ice, at 16 K. The analysis of laboratory data indicates that in dense molecular clouds, when high CO depletion occurs, an amount of carbon chain oxides as high as 2-3x10^-3 with respect to gas phase carbon monoxide can be formed after ion irradiation of icy grain mantles. Then we have searched for gas phase C2O and C3O towards ten low-mass young stellar objects. Among these we have detected the C3O line at 38486.891 MHz towards the low-mass protostar Elias 18. On the basis of the laboratory results we suggest that in dense molecular clouds gas phase carbon chain oxides are formed in the solid phase after cosmic ion irradiation of CO-rich icy mantles and released to the gas phase after desorption of icy mantles. We expect that the Atacama Large Millimeter Array (ALMA), thanks to its high sensitivity and resolution, will increase the number of carbon chain oxides detected in dense molecular clouds.

IC 4406 is a large (about 100'' x 30'') southern bipolar planetary nebula, composed of two elongated lobes extending from a bright central region, where there is evidence for the presence of a large torus of gas and dust. We show new observations of this source performed with IRAC (Spitzer Space Telescope) and the Australia Telescope Compact Array. The radio maps show that the flux from the ionized gas is concentrated in the bright central region and originates in a clumpy structure previously observed in H_alpha, while in the infrared images filaments and clumps can be seen in the extended nebular envelope, the central region showing toroidal emission. Modeling of the infrared emission leads to the conclusion that several dust components are present in the nebula.

The fast rotating star CU Virginis is a magnetic chemically peculiar star with an oblique dipolar magnetic field. The continuum radio emission has been interpreted as gyrosyncrotron emission arising from a thin magnetospheric layer. Previous radio observations at 1.4 GHz showed that a 100% circular polarized and highly directive emission component overlaps to the continuum emission two times per rotation, when the magnetic axis lies in the plane of the sky. This sort of radio lighthouse has been proposed to be due to cyclotron maser emission generated above the magnetic pole and propagating perpendicularly to the magnetic axis. Observations carried out with the Australia Telescope Compact Array at 1.4 and 2.5 GHz one year after this discovery show that this radio emission is still present, meaning that the phenomenon responsible for this process is steady on a timescale of years. The emitted radiation spans at least 1 GHz, being observed from 1.4 to 2.5 GHz. On the light of recent results on the physics of the magnetosphere of this star, the possibility of plasma radiation is ruled out. The characteristics of this radio lighthouse provides us a good marker of the rotation period, since the peaks are visible at particular rotational phases. After one year, they show a delay of about 15 minutes. This is interpreted as a new abrupt spinning down of the star. Among several possibilities, a quick emptying of the equatorial magnetic belt after reaching the maximum density can account for the magnitude of the breaking. The study of the coherent emission in stars like CU Vir, as well as in pre main sequence stars, can give important insight into the angular momentum evolution in young stars. This is a promising field of investigation that high sensitivity radio interferometers such as SKA can exploit.

We present new millimetre 43 GHz observations of a sample of radio-bright Planetary Nebulae. Such observations were carried out to have a good determination of the high-frequency radio spectra of the sample in order to evaluate, together with far-IR measurements (IRAS), the fluxes emitted by the selected source in the millimetre and sub-millimetre band. This spectral range, even very important to constraint the physics of circumstellar environment, is still far to be completely exploited. To estimate the millimetre and sub-millimetre fluxes, we extrapolated and summed together the ionized gas (free-free radio emission) and dust (thermal emission) contributions at this frequency range. By comparison of the derived flux densities to the foreseen sensitivity we investigate the possible detection of such source for all the channels of the forthcoming ESA's PLANCK mission. We conclude that almost 80% of our sample will be detected by PLANCK, with the higher detection rate in the higher frequency channels, where there is a good combination of brighter intrinsic flux from the sources and reduced extended Galactic foregrounds contamination despite a worst instrumental sensitivity. From the new 43 GHz, combined with single-dish 5 GHz observations from the literature, we derive radio spectral indexes, which are consistent with optically thin free-free nebula. This result indicates that the high frequency radio spectrum of our sample sources is dominated by thermal free-free and other emission, if present, are negligible.

The radio emission from the youngest known Planetary nebula, SAO244567, has been mapped at 1384, 2368, 4800, 8640, 16832 and 18752 MHz by using the Australian Telescope Compact Array (ATCA). These observations constitute the first detailed radio study of this very interesting object, as they allow us to obtain the overall radio morphology of the source and to compute, for the first time, the radio spectrum up to millimetre range. Radio emission is consistent with free-free from a wind-like shell, which is also the region where most of the [OIII] comes from as revealed by HST images. Physical parameters of the radio nebula and of the central star were derived, all consistent with SAO 244567 being a very young Planetary Nebula still embedded in the dusty remnant of the AGB phase. The optically thin radio flux density appear to decrease when compared to data from the literature. Even very appealing, the variability of the radio emission, probably related to the evolution of the central object, needs further investigations.

IC 4406 is a large (about 100'' x 30'') southern bipolar planetary nebula, composed by two elongated lobes, extending from a bright central region, where there is evidence for the presence of a large torus of gas and dust. In this poster we show new observations of this source performed with IRAC (Spitzer Space Telescope) and the Australia Telescope Compact Array. Although the possibility for faint extended emission to be missing in the radio maps cannot be ruled out, flux from the ionized gas appears to be concentrated in the bright central region. Comparing ATCA to IRAC images, it seems that, like in other planetary nebulae, ionized and neutral components spatially co-exist in IC 4406.

Observing objects in transition from pre- to young Planetary Nebula (PN), when the central star radiation starts to excite the envelope, can help us to understand the evolution of the circumstellar ejecta and their shaping mechanism/s. In our project we have selected a sample of hot post-AGB stars as Transition Phase candidates. Radio observations have led to detect free-free radiation from an ionized shell in about half of our targets, providing us with two sub-samples of ionized and non ionized Transition Objects. We are now using IRAC and IRS on the Spitzer Space Telescope to determine if extended emission is present (IRAC) and to study our targets' chemistry (IRS). In particular, by comparing spectra from the two sub-samples, the IRS observations will enable us to check how the presence of an ionization front effects the circumstellar envelope. The IRAC measurements, combined with previous ones in the literature, will give us information on the extent and physical conditions of the dust components.

Context: Being dominated by non-thermal emission from aligned relativistic jets, blazars allow us to elucidate the physics of extragalactic jets, and, ltimately, how the energy is extracted from the central black hole in radio-loud active galactic nuclei. Aims: Crucial information is provided by broad-band spectral energy distributions (SEDs), their trends with luminosity and correlated multi-frequency variability. With this study we plan to obtain a database of contemporaneous radio-to-optical spectra of a sample of blazars, which are and will be observed by current and future high-energy satellites. Methods: Since December 2004 we are performing a monthly multi-frequency radio monitoring of a sample of 35 blazars at the antennas in Medicina and Noto. Contemporaneous near-IR and optical observations for all our observing epochs are organised. Results: Until June 2006 about 4000 radio measurements and 5500 near-IR and optical measurements were obtained. Most of the sources show significant variability in all observing bands. Here we present the multi-frequency data acquired during the first eighteen months of the project, and construct the SEDs for the best-sampled sources.

We performed a millimetric survey of a sample of 24 post-AGB stars aimed to search for emission from circumstellar matter, in order to investigate the physical properties of the outer parts of the envelopes. The observations were conducted using the 37-channel Max-Planck Millimeter Bolometer array at the 30-meter IRAM telescope. The continuum emission toward the detected sources was used to quantify the mass of the emitting dust. We combined our observations with data available in literature to construct the spectral energy distribution (SED) of the sources. When the observational data cover a spectral range wide enough, some properties of circumstellar envelopes were derived by comparison with spectra computed using a radiative transfer code. Of the 24 objects in our sample, we detected millimetric continuum emission toward 11 sources. Two other sources were detected at flux level close to 3$\sigma$. The derived circumstellar dust masses range between 0.4 and $24 \times 10^{-4}$ M$_{\odot}$, but these results are affected from the uncertain about the source distances. The parameters derived from the SED fits are consistent with the values characteristic for these kind of object. As confirmed from the flux density extrapolated in the first light channels of the Atacama Large Millimetric Array, such sources would be good targets for future high resolution mapping with the ALMA facility.

In order to fully understand the physical processes in the magnetospheres of the Magnetic Chemically Peculiar stars, we performed multi-frequency radio observations of CU Virginis. The radio emission of this kind of stars arises from the interaction between energetic electrons and magnetic field. The radio data were acquired with the VLA and cover the whole rotational period of CU Virginis. For each observed frequency the radio light curves of the total flux density and fraction of circular polarization were fitted using a three-dimensional MCP magnetospheric model simulating the stellar radio emission as a function of the magnetospheric physical parameters. The observations show a clear correlation between the radio emission and the orientation of the magnetosphere of this oblique rotator. Radio emission is explained as the result of the acceleration of the wind particles in the current sheets just beyond the Alfv\'en radius, that eventually return toward the star following the magnetic field and emitting radiation by gyrosyncrotron mechanism. The accelerated electrons have a hard energetic spectrum ($N(E)\propto E^{-2}$) and the acceleration process has an efficiency of about $10^{-3}$. The Alfv\'en radius we determined is in the range of $12-17 R_\ast$ and, for a dipolar field of 3000 Gauss at the magnetic pole of the star, we determine a mass loss from the star of about $10^{-12}$ M$_{\sun}$ yr$^{-1}$. In the inner magnetosphere a detectable X-ray emission is expected.

In the last few years we have started a systematic study aimed at detecting very young Planetary Nebulae (YPNs). The project, based on both radio (multi-frequency and high resolution) and optical observations, has already provided interesting results, including the identification of new YPNs and the observation, for the first time in PN in this very early phase of their formation, of radio bipolar morphologies. In this paper we present radio and optical results on one of our targets, namely IRAS 19590-1249 (LS IV -12 111). The analysis of spectroscopic observations, obtained with the Anglo Australian Telescope (AAT) in different epochs, has allowed us to conclude that the optical spectrum of this object has evolved quite rapidly, showing, very recently, strong forbidden nebular emission lines consistent with a young PN. This result confirms that IRAS 19590-1249, has turned into a PN within the last decades and makes this source a perfect target for studying the early evolution and structure of Planetary Nebulae.

Despite numerous efforts, the transition from Asymptotic Giant Branch (AGB) stars to Planetary Nebulae (PN) is a poorly understood phase of stellar evolution. We have therefore carried out interferometric (VLA) radio observations of a sample of hot post-AGB stars, selected on the basis of their optical and infrared properties. Ten sources out of the 16 observed were detected. This indicates that most of our targets are surrounded by a nebula where the ionization has already started. This definitively determines the evolutionary status of the selected sources and provides us with a unique sample of very young Planetary Nebulae (yPNe). A comparison with another sample of yPNe confirms our working hypothesis that our targets are indeed very young, probably just in the transition toward PN. Finding transition objects is extremely important as they can provide unique clues for a better understanding of this important phase of stellar evolution.

In this paper we present a three-dimensional numerical model for the radio emission of Magnetic Chemically Peculiar stars, on the hypothesis that energetic electrons emit by the gyrosynchrotron mechanism. For this class of radio stars, characterized by a mainly dipolar magnetic field whose axis is tilted with respect to the rotational axis, the geometry of the magnetosphere and its deformation due to the stellar rotation are determined. The radio emitting region is determined by the physical conditions of the magnetosphere and of the stellar wind. Free-free absorption by the thermal plasma trapped in the inner magnetosphere is also considered. Several free parameters are involved in the model, such as the size of the emitting region, the energy spectrum and the number density of the emitting electrons, and the characteristics of the plasma in the inner magnetosphere. By solving the equation of radiative transfer, along a path parallel to the line of sight, the radio brightness distribution and the total flux density as a function of stellar rotation are computed. As the model is applied to simulate the observed 5 GHz lightcurves of HD37479 and HD37017, several possible magnetosphere configurations are found. After simulations at other frequencies, in spite of the large number of parameters involved in the modeling, two solutions in the case of HD37479 and only one solution in the case of HD37017 match the observed spectral indices. The results of our simulations agree with the magnetically confined wind-shock model in a rotating magnetosphere. The X-ray emission from the inner magnetosphere is also computed, and found to be consistent with the observations.

The radio emission from the youngest known Planetary nebula, SAO244567, has been mapped at 20,13, 6, 3.6 and 1.2 cm by using the Australian Telescope Compact Array (ATCA). These observations constitute the first detailed radio study of this very interesting object, as they allow us to obtain, for the first time, the radio morphology of the source and to compute the radio spectrum up to 18752 MHz.

Despite of numerous efforts, the stage from Asymptotic Giant Branch (AGB) to Planetary Nebulae (PN) is a poorly understood phase of stellar evolution. We have therefore carried out interferometric (VLA) observations of a sample of hot post-AGB stars, selected on the basis of their optical and infrared properties. Ten sources, out of the 16 been observed, were detected. This indicates that most of our targets are surrounded by a nebula where the ionization has already started. This definitively determines the evolutionary status of the selected sources and provides us with a unique sample of very young Planetary Nebulae (yPNs).

SN1993J is to date the radio supernova whose evolution has been monitored in greatest detail and the one which holds best promise for a comprehensive theoretical-observational analysis. The shell-like radio structure of SN1993J has expanded in general accord with models of shock excited emission, showing almost circular symmetry for over 8 years, except for a bright feature at the south-eastern region of the shell that has been observed at every epoch. The spectrum of SN1993J has flattened from alpha =-1 to alpha =-0.67 (S_(\nu) propto nu**(alpha)). The decelerated expansion can be modeled well with a single slope but apparently better with two slopes. There are also intriguing hints of structure in the expansion curve. The results by the two VLBI groups carrying out this research show general agreement, but also some differences. A comparison of the optical and VLBI results about the details of the deceleration show some discrepancies.

In this paper we present new observational evidence that supports the presence of an extra source of continuum emission in the binary system beta Lyrae. New VLA and IRAM observations, together with published data from the literature and ISO archive data, allow us to build the Spectral Energy Distribution of the binary between 5x10^9 Hz and 5x10^15Hz. The radio-millimeter part of the spectrum is consistent with free-free emission from a symbiotic-like wind associated with the primary component and ionized by the radiation field of the hidden companion. Furthermore, we also consider the possibility that the observed radio flux originates from collimated radio structures associated with the mass gaining component and its disk (Conical thermal jets). An extrapolation of this emission to the far-IR part of the spectrum indicates that in both cases the free-free contribution at these frequencies cannot explain the observations and that the observed infrared excess flux is due principally to the secondary component and its associated disk.

In this paper we present the results of a 4-year (1996 - 1999) radio flux density monitoring program for a sample of X- and $\gamma$-ray loud blazars. Our program started in January 1996 and was carried out on monthly basis at the frequencies of 5 GHz and 8.4 GHz with the 32-m antennas located in Medicina (Bologna, Italy) and Noto (Siracusa, Italy). 22 GHz data collected in Medicina from January 1996 to June 1997 will also be presented. The sample of selected sources comprises most radio loud blazars with $\delta \ge -10^{\circ}$ characterised by emission in the X- and $\gamma$-ray regimes, and target sources for the BeppoSAX X-ray mission. All sources in the sample, except J1653+397 (MKN 501), are variable during the four years of our monitoring program. We classified the type of variability in each source by means of a structure function analysis. We also computed th$\alpha_{5}^{8.4}$ for all epochs with nearly simultaneous observations ate spectral index and found that $\alpha_{5}^{8.4}$ starts flattening at the very beginning o f a radio flare, or flux density increase.

We present multifrequency VLA and VLBA observations at 8.4 GHz of the RS CVn system AR Lac, that were performed in autumn 1997 simultaneously with X-ray observations obtained from Rodono` et al. (1999). Our VLBA data indicate a resolved source with dimension close to the system separation, while the study of the flux density curve evidences a small amplitude outside of the eclipse variability. The derived five-frequencies spectra, combined with the size information from VLBA data, are compared with gyrosynchrotron emission from a two component structured source. A comparison with the results of the X-ray observations allow us to exclude the possibility that thermal gyrosynchrotron is responsible for the radio emission, but it is compatible with the hypothesis of co-spatial X-ray and radio emitting sources.

We obtained EVN observations of the X-ray transient source GRS1915+105 in its radio-loud state. In 2 images obtained before a major flare, the source is resolved, showing two opposite asymmetric jets. This indicates that steady jets exist also before a major flare, with an estimated velocity of 0.2-0.6c.

Radio observations of the magnetic chemically peculiar star CU Vir, carried out with the VLA in three different days, show that the radio emission at 20 cm is characterized by a strong enhancement at particular rotational phases. This radio emission is found to be right hand polarized with a degree of polarization close to 100 %. As common for this class of stars, the magnetic axis of CU Vir is oblique with respect to the rotational axis. By comparing the 20 cm radio light curve with the effective magnetic field available from the literature, a coincidence of the main peaks of the radio emission with the magnetic nulls has been found. This happens when the magnetic axis lies in the plane of the sky. We suggest that the high degree of polarization, together with the high directivity of the radiation, can be explained in terms of coherent radio emission. The data have been interpreted on the basis of the Electron Cyclotron Maser Emission from electrons accelerated in current sheets out of the Alfven radius toward the stellar surface and eventually reflected outward by magnetic mirroring.

In this paper we present high spatial resolution radio images of the puzzling binary system beta Lyrae obtained with MERLIN at 5 GHz. We find a nebula surrounding the binary with a brightness temperature of 11000+-700K approximately 40AU across. This definitively confirms the thermal origin of the radio emission, which is consistent with emission from the wind of the B6-8II component (mass loss of order of 10^-7 Msun per year), ionized by the radiation field of the hotter companion. This nebula, surrounding the binary, is the proof that beta Layrae evolved in a non-conservative way, i. e. not all the mass lost by the primary is accretted by the secondary, and present measurements indicate that almost 0.015Msun had been lost from the system since the onset of the Roche lobe overflow phase. Moreover, the nebula is aligned with the jet-like structures inferred from recent optical measurements, indicating a possible connection among them.

A rarity among supernova, SN 1993J in M81 can be studied with high spatial resolution. Its radio power and distance permit VLBI observations to monitor the expansion of its angular structure. This radio structure was previously revealed to be shell-like and to be undergoing a self-similar expansion at a constant rate. From VLBI observations at the wavelengths of 3.6 and 6 cm in the period 6 to 42 months after explosion, we have discovered that the expansion is decelerating. Our measurement of this deceleration yields estimates of the density profiles of the supernova ejecta and circumstellar material in standard supernova explosion models.