This paper contains the results from the first 30 months of the NELIOTA project for impacts of NEOs/meteoroids on the lunar surface. Our analysis on the statistics concerning the efficiency of the campaign and the parameters of the projectiles and the impacts is presented. The parameters of the lunar impact flashes based on simultaneous observations in two wavelength bands are used to estimate the distributions of the masses, sizes and frequency of the impactors. These statistics can be used both in space engineering and science. The photometric fluxes of the flashes are measured using aperture photometry and their apparent magnitudes are calculated using standard stars. Assuming that the flashes follow a black body law of irradiation, the temperatures can be derived analytically, while the parameters of the projectiles are estimated using fair assumptions on their velocity and luminous efficiency of the impacts. 79 lunar impact flashes have been observed with the 1.2 m Kryoneri telescope in Greece. The masses of the meteoroids range between 0.7 g and 8 kg and their respective sizes between 1-20 cm depending on their assumed density, impact velocity, and luminous efficiency. We find a strong correlation between the observed magnitudes of the flashes and the masses of the meteoroids. Moreover, an empirical relation between the emitted energies of each band has been derived allowing the estimation of the physical parameters of the meteoroids that produce low energy impact flashes. The NELIOTA project has so far the highest detection rate and the faintest limiting magnitude for lunar impacts compared to other ongoing programs. Based on the impact frequency distribution on Moon, we estimate that sporadic meteoroids with typical masses less than 100 g and sizes less than 5 cm enter the mesosphere of the Earth with a rate ~108 meteoroids/hr and also impact Moon with a rate of ~8 meteoroids/hr.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Pismis 18 is a moderately populated, intermediate-age open cluster located within the solar circle at a Galactocentric distance of about 7 kpc. Few open clusters have been studied in detail in the inner disc region before the Gaia-ESO Survey. New data from the Gaia-ESO Survey allowed us to conduct an extended radial velocity membership study as well as spectroscopic metallicity and detailed chemical abundance measurements for this cluster. Gaia-ESO Survey data for 142 potential members, lying on the upper MS and on the red clump, yielded radial velocity measurements, which, together with proper motion measurements from the Gaia DR2, were used to determine the systemic velocity of the cluster and membership of individual stars. Photometry from Gaia DR2 was used to re-determine cluster parameters based on high confidence member stars only. Cluster abundance measurements of six radial-velocity member stars with UVES high-resolution spectroscopy are presented for 23 elements. According to the new estimates, based on high confidence members, Pismis 18 has an age of $700^{+40}_{-50}$ Myr, interstellar reddening of E(B-V) = $0.562^{+0.012}_{-0.026}$ mag and a de-reddened distance modulus of $DM_0 = 11.96^{+0.10}_{-0.24}$ mag. The median metallicity of the cluster (using the six UVES stars) is [Fe/H] = $+0.23 \pm 0.05$ dex, with [$\alpha$/Fe]= $0.07 \pm 0.13$ and a slight enhancement of s- and r- neutron-capture elements. With the present work, we fully characterized the open cluster Pismis 18, confirming its present location in the inner disc. We estimated a younger age than the previous literature values and gave, for the first time, its metallicity and its detailed abundances. Its [$\alpha$/Fe] and [s-process/Fe], both slightly super-solar, are in agreement with other inner-disc open clusters observed by the Gaia-ESO survey. [abridged]

NELIOTA (NEO Lunar Impacts and Optical TrAnsients) is an ESA-funded lunar monitoring project, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. A prime focus, high-speed, twin-camera system providing simultaneous observations in two photometric bands at a rate of 30 frames-per-second on the 1.2 m Kryoneri telescope of the National Observatory of Athens was commissioned for this purpose. A dedicated software processes the images and automatically detects candidate lunar impact flashes, which are then validated by an expert user. The four year observing campaign began in February 2017 and has so far detected more than 40 lunar impact events. The software routinely detects satellites, which typically appear as streaks or dots crossing the lunar disk. To avoid confusing these events with real flashes, we check different available catalogs with spacecraft orbital information and exclude spacecraft identifications.

We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrument was specially designed and installed at the 1.2 m Kryoneri telescope utilizing the fast-frame capabilities of scientific Complementary Metal-Oxide Semiconductor detectors (sCMOS). The system provides a wide field-of-view (17.0' $\times$ 14.4') and simultaneous observations in two photometric bands (R and I), reaching limiting magnitudes of 18.7 mag in 10 sec in both bands at a 2.5 signal-to-noise level. This makes it a unique instrument that can be used for the detection of NEO impacts on the Moon, as well as for any astronomy projects that demand high-cadence multicolor observations. The wide field-of-view ensures that a large portion of the Moon is observed, while the simultaneous, high-cadence, monitoring in two photometric bands makes possible, for the first time, the determination of the temperatures of the impacts on the Moon's surface and the validation of the impact flashes from a single site. Considering the varying background level on the Moon's surface we demonstrate that the NELIOTA system can detect NEO impact flashes at a 2.5 signal-to-noise level of ~12.4 mag in the I-band and R-band for observations made at low lunar phases ~0.1. We report 31 NEO impact flashes detected during the first year of the NELIOTA campaign. The faintest flash was at 11.24 mag in the R-band (about two magnitudes fainter than ever observed before) at lunar phase 0.32. Our observations suggest a detection rate of $1.96 \times 10^{-7}$ events $km^{-2} h^{-1}$.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We report the first scientific results from the NELIOTA (NEO Lunar Impacts and Optical TrAnsients) project, which has recently begun lunar monitoring observations with the 1.2-m Kryoneri telescope. NELIOTA aims to detect faint impact flashes produced by near-Earth meteoroids and asteroids and thereby help constrain the size-frequency distribution of near-Earth objects in the decimeter to meter range. The NELIOTA setup, consisting of two fast-frame cameras observing simultaneously in the $R$ and $I-$bands, enables - for the first time - direct analytical calculation of the flash temperatures. We present the first 10 flashes detected, for which we find temperatures in the range ~1,600-3,100 K, in agreement with theoretical values. Two of these flashes were detected on multiple frames in both filters and therefore yield the first measurements of the temperature drop for lunar flashes. In addition, we compute the impactor masses, which range between ~100 g and ~50 kg.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

A cluster finding method was developed and applied in four Local Group Galaxies (SMC, M31, M33 and NGC 6822). The aim is to study the young stellar population of these galaxies by identifying stellar structures in small and large scales. Also our aim is to assess the potential of using the observations of ESA's space mission Gaia for the study of nearby galaxies resolved in stars. The detection method used is a Hierarchical technique based on a modified friends of friends algorithm. The identified clusters are classified in five distinct categories according to their size. The data for our study were used from two ground based surveys, the Local Group Galaxy Survey and the Maggelanic Clouds Spectroscopic Survey. Relatively young main sequence stars were selected from the stellar catalogs and were used by the detection algorithm. Multiple young stellar structures were identified in all galaxies with size varying from very small scales of a few pc up to scales larger than 1 kpc. The same cluster finding method was used in six spiral galaxies observed with the Hubble Space Telescope in a previous study. The average size in each category of the identified structures in the Local Group galaxies presents values consistent with the identified structures in the relatively distant spiral galaxies. Most of the structures consist of stars within the observational limits of Gaia's instruments. It is expected that Gaia's observations will contribute significantly on the study of the young stellar population of nearby galaxies.

NELIOTA is a new ESA activity launched at the National Observatory of Athens in February 2015 aiming to determine the distribution and frequency of small near-earth objects (NEOs) via lunar monitoring. The project involves upgrading the 1.2m Kryoneri telescope at the National Observatory of Athens, procuring two fast-frame cameras, and developing a software system, which will control the telescope and the cameras, process the images and automatically detect NEO impacts. NELIOTA will provide a web-based user interface, where the impact events will be reported and made available to the scientific community and the general public. The objective of this 3.5 year activity is to design, develop and implement a highly automated lunar monitoring system, which will conduct an observing campaign for 2 years in search of NEO impact flashes on the Moon. The impact events will be verified, characterised and reported. The 1.2m telescope will be capable of detecting flashes much fainter than current, small-aperture, lunar monitoring telescopes. NELIOTA is therefore expected to characterise the frequency and distribution of NEOs weighing as little as a few grams.

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

Aims: To investigate a possible dependence between age and metallicity in the Magellanic Clouds (MCs) from a study of small open star clusters, using Str\"{o}mgren photometry. Our goal is to trace evidence of an age metallicity relation (AMR) and correlate it with the mutual interactions of the two MCs. Our aim is also to correlate the AMR with the spatial distribution of the clusters. In the Large Magellanic Cloud (LMC), the majority of the selected clusters are young (up to 1 Gyr) and our aim is to search for an AMR at this epoch which has not been much studied. Methods: We report on results for 15 LMC and 8 Small Magellanic Cloud (SMC) clusters, scattered all over the area of these galaxies, to cover a wide spatial distribution and metallicity range. The selected LMC clusters were observed with the 1.54m Danish Telescope in Chile, using the Danish Faint Object Spectrograph and Camera (DFOSC) with a single 2kx2k CCD, whereas the SMC clusters were observed with the ESO 3.6m Telescope also in Chile with the ESO Faint Object Spectrograph and Camera (EFOSC). The obtained frames were analysed with the conventional DAOPHOT and IRAF software. We used Str\"{o}mgren filters in order to achieve reliable metallicities from photometry. Isochrone fitting was used in order to determine the ages and metallicities. Results: The AMR for the LMC displays a metallicity gradient, with higher metallicities for the younger ages. The AMR for LMC-SMC star clusters shows a possible jump in metallicity and a considerable increase at about 6x10^8 yr. It is possible that this is connected to the latest LMC-SMC interaction. The AMR for the LMC also displays a metallicity gradient with distance from the centre. The metallicities in SMC are lower, as expected for a metal poor host galaxy.

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.

The photometric and spectroscopic observational campaign organized for the 2008/9 eclipse of EE Cep revealed features, which indicate that the eclipsing disk in the EE Cep system has a multi-ring structure. We suggest that the gaps in the disk can be related to the possible planet formation.

This paper is the second in a series, implementing a classification system for Gaia observations of unresolved galaxies. Our goals are to determine spectral classes and estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new extended library of synthetic galaxy spectra, (2) its comparison with various observations, and (3) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Using the PEGASE.2 code, based on galaxy evolution models that take account of metallicity evolution, extinction correction, and emission lines (with stellar spectra based on the BaSeL library), we improved our first library and extended it to cover the domain of most of the SDSS catalogue. We produce an extended library of 28885 synthetic galaxy spectra at zero redshift covering four general Hubble types of galaxies, over the wavelength range between 250 and 1050 nm at a sampling of 1 nm or less. The library is also produced for 4 random values of redshift in the range of 0-0.2. It is computed on a random grid of four key astrophysical parameters (infall timescale and 3 parameters defining the SFR) and, depending on the galaxy type, on two values of the age of the galaxy. The synthetic library was compared and found to be in good agreement with various observations. The first results from the SVM classifiers and parametrizers are promising, indicating that Hubble types can be reliably predicted and several parameters estimated with low bias and variance.

The star complexes (large scale star forming regions) of NGC 6822 were traced and mapped and their size distribution was compared with the size distribution of star complexes in the Magellanic Clouds (MCs). Furthermore, the spatial distributions of different age stellar populations were compared with each other. The star complexes of NGC 6822 were determined by using the isopleths, based on star counts, of the young stars of the galaxy, using a statistical cutoff limit in density. In order to map them and determine their geometrical properties, an ellipse was fitted to every distinct region satisfying this minimum limit. The Kolmogorov-Smirnov statistical test was used to study possible patterns in their size distribution. Isopleths were also used to study the stellar populations of NGC 6822. The star complexes of NGC 6822 were detected and a list of their positions and sizes was produced. Indications of hierarchical star formation, in terms of spatial distribution, time evolution and preferable sizes were found in NGC 6822 and the MCs. The spatial distribution of the various age stellar populations has indicated traces of an interaction in NGC 6822, dated before 350 +/- 50 Myr.

Aims:The Gaia astrometric survey mission will, as a consequence of its scanning law, obtain low resolution optical (330-1000 nm) spectrophotometry of several million unresolved galaxies brighter than V=22. We present the first steps in a project to design and implement a classification system for these data. The goal is both to determine morphological classes and to estimate intrinsic astrophysical parameters via synthetic templates. Here we describe (1) a new library of synthetic galaxy spectra, and (2) first results of classification and parametrization experiments using simulated Gaia spectrophotometry of this library. Methods:We have created a large grid of synthetic galaxy spectra using the PEGASE.2 code, which is based on galaxy evolution models that take into account metallicity evolution, extinction correction, emission lines (with stellar spectra based on the BaSeL library). Our classification and regression models are Support Vector Machines (SVMs), which are kernel-based nonlinear estimators. Results:We produce a basic library of about 4000 zero redshift galaxy spectra covering the main Hubble types over wavelength range 250 to 1050 nm at a sampling of 1 nm or less. It is computed on a regular grid of four key astrophysical parameters for each type and for intermediate random values of the same parameters. An extended library reproduces this at a series of redshifts. Initial results from the SVM classifiers and parametrizers are promising, indicating that Hubble types can be reliably predicted and several parameters estimated with low bias and variance. Comparing the colours of our synthetic library with Sloan Digital Sky Survey (SDSS) spectra we find good agreement over the full range of Hubble types and parameters.

We report the serendipitous discovery of a peculiar main sequence in archived Hubble Space Telescope WFPC2 observations of the young star cluster NGC 2011 in the Large Magellanic Cloud. The bright part of this main sequence exhibits a prominent double, fork-like feature, as if it consists of twin main sequences, one of them being redder. The color-magnitude diagram, constructed from the stars found in the only available WFPC2 field of the cluster, is used to distinguish the stars according to their membership to each of these sequences and to study their spatial distribution. We find that there are two well distinguished populations in the sense that the redder main sequence is dominated by stars that belong to the main body of the cluster, while the stars of the bluer main sequence belong to the surrounding region. Providing that NGC 2011 is a verified binary cluster, with the second companion unfortunately not observed, and taking into account the general region where this cluster is located, we discuss the possible scenarios from both star formation, and early dynamical evolution point-of-view that might explain this unique discovery.

UBVRI photometry observations of 10 symbiotic stars and related objects obtained in the period 2002-2003 are presented. Analysing differential light curves we found rapid light variations with timescales of tens of minutes and significant amplitudes in the well-known flickerers MWC 560, RS Oph, V407 Cyg and T CrB. MWC 560 and V407 Cyg demonstrate quasi periodic oscillations (QPO) with similar amplitudes and timescales. Flickering and unusual flare in V627 Cas as well as some indications of flickering presence in BX Mon are detected. The existence of 29 minutes oscillations in Z And with an amplitude about 0.02 mag in U band is confirmed. Only one symbiotic star, V471 Per, and both non symbiotic, CI Cam and V886 Her seem to be constant on flickering timescales. Nevertheless, small night to night changes in the brightness of V886 Her were observed as well.

We have observed S5 2007+777 and 3C371 in the B and I bands for 13 and 8 nights, respectively, during various observing runs in 2001, 2002 and 2004. The observations resulted in almost evenly sampled light curves, 6-9 hours long. We do not detect any flares within the observed light curves, but we do observe small amplitude, significant variations, in both bands, on time scales of hours and days. The average variability amplitude on time scales of minutes/hours is 2.5% and 1-1.5% in the case of S5 2007+777 and 3C371, respectively. The average amplitudes increase to 5-12% and 4-6%, respectively, on time scales of days. We find that the B and I band variations are highly correlated, on both short and long time scales. During the 2004 observations, which resulted in the longest light curves, we observe two well defined flux-decay and rising trends in the light curves of both objects. When the flux decays, we observe significant delays, with the B band flux decaying faster than the flux in the I band. As a result, we also observe significant, flux related spectral variations as well. The flux-spectral relation is rather complicated, with loop-like structures forming during the flux evolution. The presence of spectral variations imply that the observed variability is not caused by geometric effects. On the other hand, our results are fully consistent with the hypothesis that the observed variations are caused by perturbations which affect different regions in the jet of the sources.

Aims. This study focuses on very luminous Mbol<-6.0 mag AGB stars with J-Ks>1.5 mag and H-Ks>0.4 mag in the LMC, SMC, M31, and M33 from 2MASS data. Methods.The data were taken from the 2MASS All-Sky Point Source catalogue archive. We used Virtual Observatory tools and took advantage of its capabilities at various stages in the analysis. Results. It is well known that stars with the colors we selected correspond mainly to carbon stars. Although the most luminous AGBs detected here contain a large number of carbon stars,they are not included in existing catalogues produced from data in the optical domain, where they are not visible since they are dust-enshrouded. A comparison of the AGB stars detected with combined near and mid-infrared data from MSX and 2MASS in the LMC shows that 10% of the bright AGB stars are bright carbon stars never detected before and that the other 50% are OH/IR oxygen rich stars, whereas the 40% that remain were not cross-matched. Conclusions. The catalogues of the most luminous AGB stars compiled here are an important complement to existing data. In the LMC, these bright AGB stars are centrally located, whereas they are concentrated in an active star-formation ring in M31. In the SMC and M33, there are not enough of them to draw definite conclusions, although they tend to be centrally located. Their luminosity functions are similar for the four galaxies we studied.

There are 509 expanding neutral hydrogen shells catalogued in the Small Magellanic Cloud (SMC), all apparently very young, with dynamical ages of a few Myr. To examine their relationship with young stellar objects we cross-correlated the shell catalogue with various catalogues of OB associations, super giants, Cepheids, WR stars, supernova remnants, and star clusters. The incidence of chance line-ups was estimated via Monte-Carlo simulations, and found to be high. However, it is important that there are 1.5 times more shells that are {\em not} spatially correlated to an OB association, than shells that are. Moreover, 59 of the 509 shells lie mainly in low stellar density fields and have no young stellar objects associated with them, and therefore no obvious energy source. It is shown that, on the whole, the properties of these "empty" shells are very similar to the properties of the rest of the shells, once selection biases are taken into account. Theoretical implications of the results are discussed.

We present results from a deep mid-IR survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15um with the CAM instrument on board ISO. We found reliable optical/near-IR associations for 32 of the 35 sources detected in this field by Oliver et al. (2002, Paper I): eight of them were identified as stars, one is definitely an AGN, a second seems likely to be an AGN, too, while the remaining 22 appear to be normal spiral or starburst galaxies. Using model spectral energy distributions (SEDs) of similar galaxies, we compare methods for estimating the star formation rates (SFRs) in these objects, finding that an estimator based on integrated (3-1000um) IR luminosity reproduces the model SFRs best. Applying this estimator to model fits to the SEDs of our 22 spiral and starburst galaxies, we find that they are forming stars at rates of ~1-100 M_sol/yr, with a median value of ~40M_sol/yr, assuming an Einstein - de Sitter universe with a Hubble constant of 50 km/s/Mpc, and star formation taking place according to a Salpeter (1955) IMF across the mass range 0.1-100M_sol. We split the redshift range 0.0<z<0.6 into two equal-volume bins to compute raw estimates of the star formation rate density contributed by these sources, assuming the same cosmology and IMF as above and computing errors based on estimated uncertainties in the SFRs of individual galaxies. We compare these results with other estimates of the SFR density made with the same assumptions, showing them to be consistent with the results of Flores et al. (1999) from their ISO survey of the CFRS 1415+52 field. However, the relatively small volume of our survey means that our SFR density estimates suffer from a large sampling variance, implying that our results, by themselves, do not place tight constraints on the global mean SFR density.

We present results from a deep mid-infrared survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15 micron with the CAM instrument on board the Infrared Space Observatory (ISO). The final map in each band was constructed by the coaddition of four independent rasters, registered using bright sources securely detected in all rasters, with the absolute astrometry being defined by a radio source detected at both 7 and 15 micron. We sought detections of bright sources in a circular region of radius 2.5 arcmin at the centre of each map, in a manner that simulations indicated would produce highly reliable and complete source catalogues using simple selection criteria. Merging source lists in the two bands yielded a catalogue of 35 distinct sources, which we calibrated photometrically using photospheric models of late-type stars detected in our data. We present extragalactic source count results in both bands, and discuss the constraints they impose on models of galaxy evolution models, given the volume of space sampled by this galaxy population.