The intermediate-mass pre-main sequence Herbig Ae/Be stars are key to understanding the differences in formation mechanisms between low- and high-mass stars. The study of the general properties of these objects is hampered by the fact that few and mostly serendipitously discovered sources are known. Our goal is to identify new Herbig Ae/Be candidates to create a homogeneous and well defined catalogue of these objects. We have applied machine learning techniques to 4,150,983 sources with data from Gaia DR2, 2MASS, WISE, and IPHAS or VPHAS+. Several observables were chosen to identify new Herbig Ae/Be candidates based on our current knowledge of this class, which is characterised by infrared excesses, photometric variabilities, and H$\alpha$ emission lines. Classical techniques are not efficient for identifying new Herbig Ae/Be stars mainly because of their similarity with classical Be stars, with which they share many characteristics. By focusing on disentangling these two types of objects, our algorithm has also identified new classical Be stars. We have obtained a large catalogue of 8470 new pre-main sequence candidates and another catalogue of 693 new classical Be candidates with a completeness of $78.8\pm1.4\%$ and $85.5\pm1.2\%$, respectively. Of the catalogue of pre-main sequence candidates, at least 1361 sources are potentially new Herbig Ae/Be candidates according to their position in the Hertzsprung-Russell diagram. In this study we present the methodology used, evaluate the quality of the catalogues, and perform an analysis of their flaws and biases. For this assessment, we make use of observables that have not been accounted for by the algorithm and hence are selection-independent, such as coordinates and parallax based distances. The catalogue of new Herbig Ae/Be stars that we present here increases the number of known objects of the class by an order of magnitude.

We use Gaia Data Release 2 (DR2) to place 252 Herbig Ae/Be stars in the HR diagram and investigate their characteristics and properties. For all known Herbig Ae/Be stars with parallaxes in Gaia DR2, we collected their atmospheric parameters and photometric and extinction values from the literature. To these data we added near- and mid-infrared photometry, collected H$\alpha$ equivalent widths and line profiles, and their binarity status. In addition, we developed a photometric variability indicator from Gaia's DR2 information. We provide masses, ages, luminosities, distances, variabilities and infrared excesses homogeneously derived for the most complete sample of Herbig Ae/Be stars to date. We find that high mass stars have a much smaller infrared excess and have much lower optical variabilities compared to lower mass stars, with the break at around 7M$_{\odot}$. H$\alpha$ emission is generally correlated with infrared excess, with the correlation being stronger for infrared emission at wavelengths tracing the hot dust closest to the star. The variability indicator as developed by us shows that $\sim$25% of all Herbig Ae/Be stars are strongly variable. We observe that the strongly variable objects display doubly peaked H$\alpha$ line profiles, indicating an edge-on disk. The fraction of strongly variable Herbig Ae stars is close to that found for A-type UX Ori stars. It had been suggested that this variability is in most cases due to asymmetric dusty disk structures seen edge-on. The observation here is in strong support of this hypothesis. Finally, the difference in dust properties occurs at 7M$_{\odot}$, while various properties traced at UV/optical wavelengths differ at 3M$_{\odot}$. The latter has been linked to different accretion mechanisms at work whereas the differing infrared properties and variabilities are related to different (dust-)disk dispersal mechanisms.

Microvariability (intra-night variability) is a low amplitude flux change at short time scales (i.e. hours). It has been detected in unobscured type 1 AGNs and blazars. However in type 2 AGNs, the detection is hampered by the low contrast between the presumably variable nucleus and the host galaxy. In this paper, we present a search for microvariability in a sample of four type 2 quasars as an astrostatistical problem. We are exploring the use of a newly introduced enhanced F-test, proposed by de Diego 2014. The presented results show that out of our four observed targets, we were able to apply this statistical method to three of them. Evidence of microvariations is clear in the case of quasar J0802+2552 in all used filters (g',r' and i') during both observing nights, the microvariations are present in one of the nights of observations of J1258+5239 in one filter (i'), while for the J1316+4452, there is evidence for microvariability within our detection levels during one night and two filters (r' and i'). We demonstrate the feasibility of the enhanced F-test to detect microvariability in obscured type 2 quasars. At the end of this paper, we discuss possible causes of microvariability. One of the options is the misclassification of the targets. A likely scenario for explanation of the phenomenon involves optically thin gaps in a clumpy obscuring medium, in accordance with the present view of the circumnuclear region. There is a possible interesting connection between the merging state of the targets and detection of microvariability.

Fast brightness variations are a unique tool to probe the innermost regions of active galactic nuclei (AGN). These variations are called microvariability or intra-night variability, and this phenomenon has been monitored in samples of blazars and unobscured AGNs. Detecting optical microvariations in targets hidden by the obscuring torus is a challenging task because the region responsible for the variations is hidden from our sight. However, there have been reports of fast variations in obscured Seyfert galaxies in X-rays, which rises the question whether microvariations can also be detected in obscured AGNs in the optical regime. Because the expected variations are very small and can easily be lost within the noise, the analysis requires a statistical approach. We report the use of a one-way analysis of variance, ANOVA, with which we searched for microvariability. ANOVA was successfully employed in previous studies of unobscured AGNs. As a result, we found microvariable events during three observing blocks: in two we observed the same object (Mrk 477), and in another, J0759+5050. The results on Mrk 477 confirm previous findings. However, since Mrk 477 is quite a peculiar target with hidden broad-line regions, we cannot rule out the possibility that we have serendipitously chosen a target prone to variations.

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

We take advantage of the exceptional photometric coverage provided by the combination of GALEX data in the UV and the ALHAMBRA survey in the optical and near-IR to analyze the physical properties of a sample of 1225 GALEX-selected Lyman break galaxies (LBGs) at $0.8 \lesssim z \lesssim 1.2$ located in the COSMOS field. This is the largest sample of LBGs studied at that redshift range so far. According to a spectral energy distribution (SED) fitting with synthetic stellar population templates, we find that LBGs at $z \sim 1$ are mostly young galaxies with a median age of 341 Myr and have intermediate dust attenuation, $\ < E_s (B-V) \ > \sim 0.20$. Due to their selection criterion, LBGs at $z \sim 1$ are UV-bright galaxies and have high dust-corrected total SFR, with a median value of 16.9 $M_\odot {\rm yr}^{-1}$. Their median stellar mass is $\log{\left(M_*/M_\odot \right)} = 9.74$. We obtain that the dust-corrected total SFR of LBGs increases with stellar mass and the specific SFR is lower for more massive galaxies. Only 2% of the galaxies selected through the Lyman break criterion have an AGN nature. LBGs at $z \sim 1$ are mostly located over the blue cloud of the color-magnitude diagram of galaxies at their redshift, with only the oldest and/or the dustiest deviating towards the green valley and red sequence. Morphologically, 69% of LBGs are disk-like galaxies, with the fraction of interacting, compact, or irregular systems being much lower, below 12%. LBGs have a median effective radius of 2.5 kpc and bigger galaxies have higher total SFR and stellar mass. Comparing to their high-redshift analogues, we find evidence that LBGs at lower redshifts are bigger, redder in the UV continuum, and have a major presence of older stellar populations in their SEDs. However, we do not find significant difference in the distributions of stellar mass or dust attenuation.

Lyman break galaxies (LBGs) represent one of the kinds of star-forming galaxies that are found in the high-redshift universe. The detection of LBGs in the FIR domain can provide very important clues on their dust attenuation and total SFR, allowing a more detailed study than those performed so far. In this work we explore the FIR emission of a sample of 16 LBGs at z ~ 3 in the GOODS-North and GOODS-South fields that are individually detected in PACS-100um or PACS-160um. These detections demonstrate the possibility of measuring the dust emission of LBGs at high redshift. We find that PACS-detected LBGs at z ~ 3 are highly obscured galaxies which belong to the Ultra luminous IR galaxies or Hyper luminous IR galaxies class. Their total SFR cannot be recovered with the dust attenuation factors obtained from their UV continuum slope or their SED-derived dust attenuation employing Bruzual & Charlot (2003) templates. Both methods underestimate the results for most of the galaxies. Comparing with a sample of PACS-detected LBGs at z ~ 1 we find evidences that the FIR emission of LBGs might have changed with redshift in the sense that the dustiest LBGs found at z ~ 3 have more prominent FIR emission, are dustier for a given UV slope, and have higher SFR for a given stellar mass than the dustiest LBGs found at z ~ 1.

The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.

Aimed at understanding the evolution of galaxies in clusters, the GLACE survey is mapping a set of optical lines ([OII]3727, [OIII]5007, Hbeta and Halpha/[NII] when possible) in several galaxy clusters at redshift around 0.40, 0.63 and 0.86, using the Tuneable Filters (TF) of the OSIRIS instrument (Cepa et al. 2005) at the 10.4m GTC telescope. This study will address key questions about the physical processes acting upon the infalling galaxies during the course of hierarchical growth of clusters. GLACE is already ongoing: we present some preliminary results on our observations of the galaxy cluster Cl0024+1654 at z = 0.395; on the other hand, GLACE@0.86 has been approved as ESO/GTC large project to be started in 2011.