Based on very deep photometry, Surface Brightness Fluctuations (SBF) have been traditionally used to determine galaxy distances. We have recently computed SBF spectra of stellar populations at moderately high resolution, which are fully based on empirical stars. We show that the SBF spectra provide an unprecedented potential for stellar population studies that, so far, have been tackled on the basis of the mean fluxes. We find that the SBFs are able to unveil metal-poor stellar components at the one percent level, which are not possible to disentangle with the standard analysis. As these metal-poor components correspond to the first stages of the chemical enrichment, the SBF analysis provides stringent constrains on the quenching epoch.

It is generally recognized that massive galaxies form through a combination of in-situ collapse and ex-situ accretion. The in-situ component forms early, where gas collapse and compaction leads to the formation of massive compact systems (blue and red "nuggets") seen at $z>1$. The subsequent accretion of satellites brings in ex-situ material, growing these nuggets in size and mass to appear as the massive early-type galaxies (ETGs) we see locally. Due to stochasticity in the accretion process, in a few rare cases a red nugget will evolve to the present day having undergone little ex-situ mass accretion. The resulting massive, compact and ancient objects have been termed "relic galaxies". Detailed stellar population and kinematic analyses are required to characterise these systems. However, an additional crucial aspect lies in determining the fraction of ex-situ mass they have accreted since their formation. Globular cluster systems can be used to constrain this fraction, since the oldest and most metal-poor globular clusters in massive galaxies are primarily an accreted, ex-situ population. Models for the formation of relic galaxies and their globular cluster systems suggest that, due to their early compaction and limited accretion of dark-matter dominated satellites, relic galaxies should have characteristically low dark-matter mass fractions compared to ETGs of the same stellar mass.

We examine the stellar population radial gradients (age, metallicity and [$\alpha/$Fe]) of $\sim$ 100 passive central galaxies up to $\sim 2 R_e$. The targeted groups have a halo mass range spanning from $11 < \log(M_{200}/M_{\odot}) < 15$, in the SAMI Galaxy Survey. The main goal of this work is to determine whether central galaxies have different stellar population properties when compared to similarly massive satellite galaxies. We find negative metallicity radial gradients, which become shallower with increasing stellar mass. The age and [$\alpha$/Fe] gradients are consistent with zero or slightly positive. [$\alpha$/Fe] gradients become more negative with increasing mass, while age gradients do not show any significant trend with mass. We do not observe a significant difference between the stellar population gradients of central and satellite galaxies, at fixed stellar mass. The mean metallicity gradients are $\overline{\Delta [Z/H]/\Delta \log(r/R_e)} = -0.25 \pm 0.03$ for central galaxies and $\overline{\Delta [Z/H]/\Delta \log(r/R_e)} = -0.30 \pm 0.01$ for satellites. The mean age and [$\alpha$/Fe] gradients are consistent between central and satellite galaxies, within the uncertainties, with a mean value of $\overline{\Delta \textrm{log (Age/Gyr)}/\Delta \log(r/R_e)} = 0.13 \pm 0.03$ for centrals and $\overline{\Delta \textrm{log (Age/Gyr)}/\Delta \log(r/R_e)} = 0.17 \pm 0.01$ for satellite and $\overline{\Delta [\alpha/Fe]/\Delta \log(r/R_e)} = 0.01 \pm 0.03$ for centrals and $\overline{\Delta [\alpha/Fe]/\Delta \log(r/R_e)} = 0.08 \pm 0.01$ for satellites. This evidence suggests that the central region of central passive galaxies form in a similar fashion to satellite passive galaxies, in agreement with a two-phase formation scenario.

We present a new set of Surface Brightness Fluctuation spectra computed with the E-MILES stellar population synthesis models. The model SBF spectra cover the range 1680-50000 Angstrom at moderately high resolution, all based on extensive empirical stellar libraries. The models span the metallicity range -2.3<[M/H<+0.26 for a suite of IMF types with varying slopes. These predictions can complement and aid fluctuation magnitudes studies, permitting a first order approximation by applying filter responses to the SBF spectra to obtain spectroscopic SBF magnitudes. We provide a recipe for obtaining the latter and discuss their uncertainties and limitations. We compare our spectroscopic SBF magnitudes to photometric data of a sample of early-type galaxies. We also show that the SBF spectra can be very useful for constraining relevant stellar population parameters. We find small (<5%) mass-fractions of extremely metal-poor components ([M/H]<-1) on the top of the dominant, old and metal-rich stellar population. These results put stringent constraints on the early stages of galaxy formation in massive elliptical galaxies. This is remarkable given the high degree of degeneracy of the standard spectral analysis to such metal-poor stellar populations in the visible and in the near-IR. The new SBF models show great potential for exploiting ongoing surveys, particularly those based on narrow-band filters.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

The 8.4m Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will start a ten-year survey of the southern hemisphere sky in 2023. LSST will revolutionise low surface brightness astronomy. It will transform our understanding of galaxy evolution, through the study of low surface brightness features around galaxies (faint shells, tidal tails, halos and stellar streams), discovery of low surface brightness galaxies and the first set of statistical measurements of the intracluster light over a significant range of cluster masses and redshifts.

The diffuse light in clusters of galaxies, or intracluster light, has attracted a lot of attention lately due to its potential in describing the assembly history of galaxy clusters and to explain the observed growth of the brightest cluster galaxy with time. The properties of this light (color, stellar populations, extent) give clues about its formation and, consequently, the processes that shape the cluster. Here, I will present a review on intracluster light, its history, properties and the particular observational problems and limitations associated with the study of this diffuse component in integrated light.

The ubiquitous presence of Galactic cirri in deep optical images represents a major obstacle to study the low surface brightness features of extragalactic sources. To address this issue, we have explored the optical properties of cirri using g, r, i and z bands in the Sloan Digital Sky Survey (SDSS) Stripe82 region. Using state-of-the-art $-$custom made$-$ image processing techniques, including the modeling and removal of the scattered light produced by the stars, we manage to isolate the optical diffuse emission by the cirri, allowing their photometric characterization. We find that their optical colors are driven by the dust column density: the cirri become redder as their 100 $\mu$m emission increases. This could explain the extended red emission previously found in dust clouds. In most cases, the optical colors of the Galactic cirri differ significantly from those of extragalactic sources. If future works confirm our findings, it would be possible the use of deep multi-band optical photometry (as the one produced by LSST) to identify the presence of cirri at a higher spatial resolution than those provided by far infrared observations. The combination of very deep data and multi-band photometry would make possible to build dust and extinction maps of unprecedented quality.

The claimed detection of a diffuse galaxy lacking dark matter represents a possible challenge to our understanding of the properties of these galaxies and galaxy formation in general. The galaxy, already identified in photographic plates taken in the summer of 1976 at the UK 48-in Schmidt telescope, presents normal distance-independent properties (e.g. colour, velocity dispersion of its globular clusters). However, distance-dependent quantities are at odds with those of other similar galaxies, namely the luminosity function and sizes of its globular clusters, mass-to-light ratio and dark matter content. Here we carry out a careful analysis of all extant data and show that they consistently indicate a much shorter distance (13 Mpc) than previously indicated (20 Mpc). With this revised distance, the galaxy appears to be a rather ordinary low surface brightness galaxy (R_e=1.4+-0.1 kpc; M*=6.0+-3.6x10^7 Msun) with plenty of room for dark matter (the fraction of dark matter inside the half mass radius is >75% and M_halo/M*>20) corresponding to a minimum halo mass >10^9 Msun. At 13 Mpc, the luminosity and structural properties of the globular clusters around the object are the same as those found in other galaxies.

High-resolution multi-wavelength photometry is crucial to explore the spatial distribution of star formation in galaxies and understand how these evolve. To this aim, in this paper we exploit the deep, multi-wavelength Hubble Space Telescope (HST) data available in the central parts of the GOODS fields and study the distribution of star formation activity and mass in galaxies located at different positions with respect to the Main Sequence (MS) of star-forming galaxies. Our sample consists of galaxies with stellar mass $\geq 10^{9.5} M_{\odot}$ in the redshift range 0.2 $ \leq z \leq 1.2$. Exploiting 10-band photometry from the UV to the near-infrared at HST resolution, we derive spatially resolved maps of galaxies properties, such as stellar mass and star formation rate and specific star formation rate, with a resolution of $\sim 0.16$ arcsec. We find that the star formation activity is centrally enhanced in galaxies above the MS and centrally suppressed below the MS, with quiescent galaxies (1 dex below the MS) characterised by the highest suppression. The sSFR in the outer region does not show systematic trends of enhancement or suppression above or below the MS. The distribution of mass in MS galaxies indicates that bulges are growing when galaxies are still on the MS relation. Galaxies below the MS are more bulge-dominated with respect to MS counterparts at fixed stellar mass, while galaxies in the upper envelope are more extended and have S\'ersic indexes that are always smaller than or comparable to their MS counterparts. The suppression of star formation activity in the central region of galaxies below the MS hints at \textit{inside-out} quenching, as star formation is still ongoing in the outer regions.

The bulk of stars in galaxy clusters are confined within their constituent galaxies. Those stars do not trace the extended distribution of dark matter well as they are located in the central regions of the cluster's dark matter sub-halos. A small fraction of stars is expected, however, to follow the global dark matter shape of the cluster. These are the stars whose extended spatial distribution results from the merging activity of galaxies and form the intracluster light (ICL). In this work, we compare the bi-dimensional distribution of dark matter in massive galaxy clusters (as traced by gravitational lensing models) with the distribution of the ICL. To do that, we use the superb data from the Hubble Frontier Fields Initiative. Using the Modified Hausdorff distance (MHD) as a way of quantifying the similarities between the mass and ICL distributions, we find an excellent agreement (MHD~25 kpc) between the two components. This result shows that the ICL exquisitely follows the global dark matter distribution, providing an accurate luminous tracer of dark matter. This finding opens up the possibility of exploring the distribution of dark matter in galaxy clusters in detail using only deep imaging observations.

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST).

Massive galaxies are thought to form in two phases: an initial, early collapse of gas and giant burst of central star formation, followed by the later accretion of material that builds up their stellar and dark matter haloes. The globular cluster systems of such galaxies are believed to form in a similar manner. The initial central burst forms metal-rich (red) clusters, while more metal-poor (blue) clusters are brought in by the later accretion of less massive satellites. This formation process is thought to lead the creation of the multimodal optical colour distributions seen in the globular cluster systems of massive galaxies. Here we report HST/ACS observations of the massive relic galaxy NGC 1277 and its globular clusters, a nearby unevolved example of a high redshift "red nugget". The g-z cluster colour distribution shows that the globular cluster system of the galaxy is unimodal and uniquely red. This is in strong contrast to normal galaxies of similar and larger stellar mass, whose cluster systems always exhibit (and are generally dominated by) blue clusters. We argue that the globular cluster system of NGC 1277 indicates that the galaxy has undergone little (if any) mass accretion after its initial collapse and use analytic merger trees to show that the total stellar mass accretion is likely less than ~ 10 %. These results confirm that NGC 1277 is a genuine relic galaxy and show that the blue, metal-poor globular clusters constitute an accreted population in present day massive galaxies.

A randomly chosen star in today's Universe is most likely to live in a galaxy with a stellar mass between that of the Milky Way and Andromeda. Yet it remains uncertain how the structural evolution of these bulge-disk systems proceeded. Most of the unobscured star formation we observe building Andromdeda progenitors at 0.7<z<1.5 occurs in disks, but >90% of their star formation is reprocessed by dust and remains unaccounted for. Here we map 500micron dust continuum emission in an Andromeda progenitor at z=1.25 to probe where it is growing through dust-obscured star formation. Combining resolved dust measurements from the NOEMA interferometer with Hubble Space Telescope Halpha maps and multicolor imaging (including new UV data from the HDUV survey), we find a bulge growing by dust-obscured star formation: while the unobscured star formation is centrally suppressed, the dust continuum is centrally concentrated, filling in the ring-like structures evident in the Halpha and UV emission. Reflecting this, the dust emission is more compact than the optical/UV tracers of star formation with r_e(dust)=3.4kpc, r_e(Halpha)/r_e(dust)=1.4, and r_e(UV)/r_e(dust)=1.8. Crucially, however, the bulge and disk of this galaxy are building simultaneously; although the dust emission is more compact than the rest-optical emission (r_e(optical)/r_e(dust)=1.4), it is somewhat less compact than the stellar mass (r_e(M_*)/r_e(dust)=0.9). Taking the 500micron emission as a tracer of star formation, the expected structural evolution of this galaxy can be accounted for by star formation: it will grow in size by Delta(r_e)/Delta(M_*)~0.3 and central surface density by Delta(Sigma_cen)/Delta(M_*)~0.9. Finally, our observations are consistent with a picture in which merging and disk instabilities drive gas to the center of galaxies, boosting global star formation rates above the main sequence and building bulges.

Multiwavelength deep observations are a key tool to understand the origin of the diffuse light in clusters of galaxies: the intra-cluster light (ICL). For this reason, we take advantage of the Hubble Frontier Fields survey to investigate the properties of the stellar populations of the ICL of its 6 massive intermediate redshift (0.3<z<0.6) clusters. We carry on this analysis down to a radial distance of ~120 kpc from the brightest cluster galaxy. We found that the average metallicity of the ICL is [Fe/H] ~-0.5, compatible with the value of the outskirts of the Milky Way. The mean stellar ages of the ICL are between 2 to 6 Gyr younger than the most massive galaxies of the clusters. Those results suggest that the ICL of these massive (> 10^15 Msol) clusters is formed by the stripping of MW-like objects that have been accreted at z<1, in agreement with current simulations. We do not find any significant increase in the fraction of light of the ICL with cosmic time, although the redshift range explored is narrow to derive any strong conclusion. When exploring the slope of the stellar mass density profile, we found that the ICL of the HFF clusters follows the shape of their underlying dark matter haloes, in agreement with the idea that the ICL is the result of the stripping of galaxies at recent times.

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]

Using the exquisite depth of the Hubble Ultra Deep Field (HUDF12 programme) dataset, we explore the ongoing assembly of the outermost regions of the most massive galaxies ($\rm M_{\rm stellar}\geq$ 5$\times$10$^{10}$ M$_{\odot}$) at $z \leq$ 1. The outskirts of massive objects, particularly Early-Types Galaxies (ETGs), are expected to suffer a dramatic transformation across cosmic time due to continuous accretion of small galaxies. HUDF imaging allows us to study this process at intermediate redshifts in 6 massive galaxies, exploring the individual surface brightness profiles out to $\sim$25 effective radii. We find that 5-20\% of the total stellar mass for the galaxies in our sample is contained within 10 $< R <$ 50 kpc. These values are in close agreement with numerical simulations, and higher than those reported for local late-type galaxies ($\lesssim$5\%). The fraction of stellar mass stored in the outer envelopes/haloes of Massive Early-Type Galaxies increases with decreasing redshift, being 28.7\% at $< z > =$ 0.1, 15.1\% at $< z > =$ 0.65 and 3.5\% at $< z > =$ 2. The fraction of mass in diffuse features linked with ongoing minor merger events is $>$ 1-2\%, very similar to predictions based on observed close pair counts. Therefore, the results for our small albeit meaningful sample suggest that the size and mass growth of the most massive galaxies have been solely driven by minor and major merging from $z =$ 1 to today.

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.

The ultra-deep multiwavelength HST Frontier Fields coverage of the Abell Cluster 2744 is used to derive the stellar population properties of its intra-cluster light (ICL). The restframe colors of the ICL of this intermediate redshift (z=0.3064) massive cluster are bluer (g-r=0.68+-0.04; i-J=0.56+-0.01) than those found in the stellar populations of its main galaxy members (g-r=0.83+-0.01; i-J=0.75+-0.01). Based on these colors, we derive the following mean metallicity Z=0.018+-0.007 for the ICL. The ICL age is 6+-3 Gyr younger than the average age of the most massive galaxies of the cluster. The fraction of stellar mass in the ICL component comprises at least 6% of the total stellar mass of the galaxy cluster. Our data is consistent with a scenario where the bulk of the ICL of Abell 2744 has been formed relatively recently (z<1). The stellar population properties of the ICL suggest that this diffuse component is mainly the result of the disruption of infalling galaxies with similar characteristics in mass (M*~ 3x10^10 Msolar) and metallicity than our own Milky Way. The amount of ICL mass in the central part of the cluster (<400 kpc) is equivalent to the disruption of 4-6 Milky Way-type galaxies.

We present a comprehensive multiwavelength photometric analysis of the innermost (3x3 square kpc) 110 globular clusters (GCs) of M87. Their spectral energy distributions (SEDs) were built taking advantage of new ground-based high resolution near-IR imaging aided by adaptive optics at the Very Large Telescope (VLT) combined with Hubble Space Telescope (HST) ultraviolet--optical archival data. These GC SEDs are among the best photometrically sampled extragalactic GC SEDs. To compare with our SEDs we constructed equally sampled SEDs of Milky Way GCs. Using both these Milky Way cluster templates and different stellar population models, ages of >10 Gyr and metallicities of [Fe/H] -0.6 dex are consistently inferred for the inner GCs of M87. In addition, the metallicity of these GCs is low (Dif([Fe/H]) 0.8 dex) compared to that of their host galaxy. These results agree with the idea that the GC formation in M87 ceased earlier than that of the bulk of the stars of the central part of the galaxy. The ages of the inner GCs of M87 support the idea that these central parts of the galaxy formed first. Our data do not support evidence of recent wet merging.

In order to probe the inside-out formation of the most massive galaxies in the Universe, we have explored the radial (0.1 < R < 8 kpc) variation of the spectral energy distribution (SED) of M87 from UV to IR. For this purpose, we have combined high resolution data in 16 different bands. Our analysis indicate that the age of the stellar population of M87 remains almost unchanged with radius. However, the metallicity ([Z/H]) profile presents three different zones: the innermost kpc shows a plateau with supersolar metallicity, followed by a decline in metallicity down to 5 kpc and another plateau afterwards. The size of the inner plateau is similar to the expected size (Re) of an object with the predicted mass of M87 at z=2. The global [Z/H] gradient is -0.26 +- 0.10, similar to those found in other nearby massive ellipticals. The observed change in the stellar population of M87 is consistent with a rapid formation of the central part (R<5 kpc) of this galaxy followed by the accretion of the outer regions through the infall of more metal-poor material.

The inner structure of AGNs is expected to change below a certain luminosity limit. The big blue bump, footprint of the accretion disk, is absent for the majority of low-luminosity AGNs (LLAGNs). Moreover, recent simulations suggest that the torus, a keystone in the Unified Model, vanishes for nuclei with L_bol < 10^42 erg/s. However, the study of LLAGN is a complex task due to the contribution of the host galaxy, which light swamps these faint nuclei. This is specially critical in the IR range, at the maximum of the torus emission, due to the contribution of the old stellar population and/or dust in the nuclear region. Adaptive optics imaging in the NIR (VLT/NaCo) together with diffraction limited imaging in the mid-IR (VLT/VISIR) permit us to isolate the nuclear emission for some of the nearest LLAGNs in the Southern Hemisphere. These data were extended to the optical/UV range (HST), radio (VLA, VLBI) and X-rays (Chandra, XMM-Newton, Integral), in order to build a genuine spectral energy distribution (SED) for each AGN with a consistent spatial resolution (< 0.5") across the whole spectral range. From the individual SEDs, we construct an average SED for LLAGNs sampled in all the wavebands mentioned before. Compared with previous multiwavelength studies of LLAGNs, this work covers the mid-IR and NIR ranges with high-spatial resolution data. The LLAGNs in the sample present a large diversity in terms of SED shapes. Some of them are very well described by a self-absorbed synchrotron (e.g. NGC 1052), while some other present a thermal-like bump at ~1 micron (NGC 4594). All of them are significantly different when compared with bright Seyferts and quasars, suggesting that the inner structure of AGNs (i.e. the torus and the accretion disk) suffers intrinsic changes at low luminosities.

High-spatial resolution near-infrared (NIR) images of the central 24 x 24 arcsec^2 (~ 2 x 2 kpc^2) of the elliptical galaxy NGC 1052 reveal a total of 25 compact sources randomly distributed in the region. Fifteen of them exhibit Halpha luminosities an order of magnitude above the estimate for an evolved population of extreme horizontal branch stars. Their Halpha equivalent widths and optical-to-NIR spectral energy distributions are consistent with them being young stellar clusters aged < 7 Myr. We consider this to be the first direct observation of spatially resolved star-forming regions in the central kiloparsecs of an elliptical galaxy. The sizes of these regions are ~< 11 pc and their median reddening is E(B - V) ~ 1 mag. According to previous works, NGC 1052 may have experienced a merger event about 1 Gyr ago. On the assumption that these clusters are spreaded with similar density over the whole galaxy, the fraction of galaxy mass (5 x 10^{-5}) and rate of star formation (0.01 Msun/yr) involved, suggest the merger event as the possible cause for the star formation we see today.

We present the results of monitoring the radio emission from the Type IIb supernova SN 2001gd between 2002 February 8 and 2002 October 28. Most of the data were obtained using the Very Large Array at the five wavelengths of $\lambda \lambda$1.3 cm (22.4 GHz), 2.0 cm (14.9 GHz), 3.6 cm (8.44 GHz), 6.2 cm (4.86 GHz), and 21 cm (1.4 GHz). Observations were also made with Giant Meterwave Radio Telescope at $\lambda$21 cm (1.4 GHz). The object was discovered optically well after maximum light, making any determination of the early radio evolution difficult. However, subsequent observations indicate that the radio emission has evolved regularly in both time and frequency and is well described by the SN shock/circumstellar medium interaction model.

Since 1997 the afterglow of gamma-ray bursting sources (GRBs) has occasionally been detected in the radio, as well in other wavelengths bands. In particular, the interesting and unusual gamma-ray burst GRB980425, thought to be related to the radio supernova SN1998bw, is a possible link between the two classes of objects. Analyzing the extensive radio emission data avaliable for SN1998bw, one can describe its time evolution within the well established framework available for the analysis of radio emission from supernovae. This then allows relatively detailed description of a number of physical properties of the object. The radio emission can best be explained as the interaction of a mildly relativistic shock with a dense preexplosion stellar wind-established circumstellar medium (CSM) that is highly structured both azimuthally, in clumps or filaments, and radially, with observed density enhancements. Because of its unusual characteristics for a Type Ib/c supernova, the relation of SN1998bw to GRB980425 is strengthened and suggests that at least some classes of GRBs originate in supernova (SN) massive star explosions. Thus, employing the formalism for describing the radio emission from supernovae (SNe) and following the link through SN1998bw/GRB980425, it is possible to model the gross properties of the radio and optical/infrared (OIR) emission from the half-dozen GRBs with extensive radio observations. From this we conclude that at least some members of the ``slow-soft'' class of GRBs can be attributed to the explosion of a massive star in a dense, highly structured CSM that was presumably established by the preexplosion stellar system.

The unusual supernova SN1998bw, which is thought to be related to the gamma-ray burster GRB980425, is a possible link between the two classes of objects. Analyzing the extensive radio emission data avaliable for SN1998bw, we are able to describe its time evolution within the well established framework available for the analysis of radio emission from supernovae. This then allows description of a number of physical properties of the object. The radio emission can best be explained as interaction of a mildly relativistic (Gamma about 1.6) shock with a dense pre-explosion stellar wind established circumstellar medium (CSM) which is highly structured both azimuthally, in clumps or filaments, and radially, with two observed density enhancements separated by about 3e17 cm. With assumptions as to pre-explosion stellar wind conditions, it is possible to estimate that the progenitor to SN 1998bw had a mass loss rate of about 3.5e-5 solar masses per yr with at least two approximately 30% increases in mass-loss rate; the most recent extending from about 1,600 - 4,700 yr before explosion and the oldest known having occurred, with possibly comparable length, about 12,000 yr before explosion. Because of its unusual characteristics for a Type Ib/c SN, the relation of SN1998bw to GRB980425 is strengthened with consequent improvement in our understanding of these poorly understood objects.

Study of radio supernovae (RSNe) over the past 20 years includes two dozen detected objects and more than 100 upper limits. From this work we are able to identify classes of radio properties, demonstrate conformance to and deviations from existing models, estimate the density and structure of the circumstellar material and, by inference, the evolution of the presupernova stellar wind, and reveal the last stages of stellar evolution before explosion. It is also possible to detect ionized hydrogen along the line of sight, to demonstrate binary properties of the stellar system, and to show clumpiness of the circumstellar material. More speculatively, it may be possible to provide distance estimates to radio supernovae. The interesting and unusual radio supernova SN 1998bw, which is thought to be related to the gamma-ray burst GRB 980425, is discussed in particular detail. Its radio properties are compared and contrasted with those of other known RSNe.

We present new radio observations of the supernova SN 1979C made with the VLA at 20, 6, 3.6, and 2 cm from 1991 July to 1998 October, which extend our previously published observations (Weiler et al.1986, 1991), beginning 8 days after optical maximum in 1979 April and continuing through 1990 December. We find that the radio emission from SN 1979C has stopped declining in flux density in the manner described by Weiler et al. (1992), and has apparently entered a new stage of evolution. The observed ``flattening,'' or possible brightening, of the radio light curves for SN 1979C is interpreted as due to the SN shock wave entering a denser region of material near the progenitor star and may be indicative of complex structure in the circumstellar medium established by the stellar wind from the red supergiant (RSG) progenitor.

Detailed radio observations of extragalactic supernovae are critical to obtaining valuable information about the nature and evolutionary phase of the progenitor star in the period of a few hundred to several tens-of-thousands of years before explosion. Additionally, radio observations of old supernovae (>20 years) provide important clues to the evolution of supernovae into supernova remnants, a gap of almost 300 years (SN ~1680=Cas A, to SN 1923A) in our current knowledge. Finally, new empirical relations indicate that it may be possible to use some types of radio supernovae as distance yardsticks, to give an independent measure of the distance scale of the Universe. However, the study of radio supernovae is limited by the sensitivity and resolution of current radio telescope arrays. Therefore, it is necessary to have more sensitive arrays, such as the Square Kilometer Array and the several other radio telescope upgrade proposals, to advance radio supernova studies and our understanding of supernovae, their progenitors, and the connection to supernova remnants.

The supernova SN 1978K has been noted for its lack of emission lines broader than a few thousand km/s since its discovery in 1990. Modeling of the radio spectrum of the peculiar SN 1978K indicates the existence of HII absorption along the line of sight. To determine the nature of this absorbing region, we have obtained a high-dispersion spectrum of SN 1978K at the wavelength range 6530--6610 \AA. The spectrum shows not only the moderately broad H-alpha emission of the supernova ejecta but also narrow nebular H-alpha and [N II] emission. The high [N II]6583/H-alpha ratio, 0.8-1.3, suggests that this radio absorbing region is a stellar ejecta nebula. The expansion velocity and emission measure of the nebula are consistent with those seen in ejecta nebulae of luminous blue variables. Previous low-dispersion spectra have detected a strong [N II]5755 line, indicating an electron density of 3-12x10^5 cm^{-3}. We argue that this stellar ejecta nebula is probably part of the pre-shock dense circumstellar envelope of SN 1978K. We further suggest that SN 1997ab may represent a young version of SN 1978K.

New observations of SN 1980K made with the VLA at 20 and 6 cm from 1994 April through 1996 October show that the supernova (SN) has undergone a significant change in its radio emission evolution, dropping by a factor of ~2 below the flux density S \propto t^{-0.73} power-law decline with time t observed earlier. However, although S at all observed frequencies has decreased significantly, its current spectral index of \alpha= -0.42\pm0.15 (S \propto \nu^{+\alpha}) is consistent with the previous spectral index of \alpha=-0.60_{-0.07}^{+0.04}. It is suggested that this decrease in emission may be due to the SN shock entering a new region of the circumstellar material which has a lower density than that expected for a constant speed (w), constant mass-loss rate (Mdot) wind from the progenitor. If such an interpretation is correct, the difference in wind and shock speeds appears to indicate a significant evolution in the mass-loss history of the SN progenitor ~10^4 years before explosion, with a change in circumstellar density (\propto Mdot/w) occurring over a time span of \lesssim 4 kyr. Such features could be explained in terms of a fast ``blue-loop'' evolutionary phase of a relatively massive pre-SN progenitor star. If so, we may, for the first time, provide a stringent constraint on the mass of the SN progenitor based solely on the SN's radio emission.

Long-term monitoring of the radio emission from supernovae with the Very Large Array (VLA) shows that the radio ``light curves'' evolve in a systematic fashion with a distinct peak flux density (and thus, in combination with a distance, a peak spectral luminosity) at each frequency and a well-defined time from explosion to that peak. Studying these two quantities at 6 cm wavelength, peak spectral luminosity (L_{6 cm peak}) and time after explosion date (t_0) to reach that peak (t_{6 cm peak} - t_0), we find that they appear related. In particular, based on two objects, Type Ib supernovae may be approximate radio ``standard candles'' with a 6 cm peak luminosity L_{6 cm peak} \approx 19.9 X 10^{26} erg s^{-1} Hz^{-1}; also based on two objects, Type Ic supernovae may be approximate radio ``standard candles'' with a 6 cm peak luminosity L_{6 cm peak} \approx 6.5 X 10^{26} erg s^{-1} Hz^{-1}; and, based on twelve objects, Type II supernovae appear to obey a relation L_{6 cm peak} \simeq 5.5 X 10^{23} (t_{6 cm peak} - t_0)^{1.4} erg s^{-1} Hz^{-1}, with time measured in days. If these relations are supported by further observations, they provide a means for determining distances to supernovae, and thus to their parent galaxies, from purely radio continuum observations.

Radio observations of SN 1986E have shown a clear detection of emission at 6 cm wavelength about 8 months after optical discovery. Combined with a number of new upper limits and a study of the possible models, it appears that SN 1986E was probably a fairly normal Type IIL supernova, somewhat similar to SN 1980K, with radio emission at roughly expected levels. This detection continues the correlation between radio detection and late time optical emission.

We extend the approximate radiative transfer analysis of Hershkowitz, Linder, and Wagoner (1986) to a more general class of supernova model atmospheres, using a simple fit to the effective continuum opacity produced by lines (Wagoner, Perez, and Vasu 1991). At the low densities considered, the populations of the excited states of hydrogen are governed mainly by photoionization and recombination, and scattering dominates absorptive opacity. We match the asymptotic expressions for the spectral energy density $\Jnu$ at the photosphere, whose location at each frequency is determined by a first-order calculation of the deviation of $\Jnu$ from the Planck function $\Bnu$. The emergent spectral luminosity then assumes the form $ L_\nu = 4\pi^2 \r*^2 \zeta^2 \Bnu(T_p) $, where $T_p(\nu)$ is the photospheric temperature, $\zeta$ is the dilution factor, and $\r*$ is a fiducial radius [ultimately taken to be the photospheric radius $r_p (\nu)$]. The atmosphere is characterized by an effective temperature $T_e$ ($\propto L^{1/4} r_*^{-1/2}$) and hydrogen density $n_H=n_* (r_*/r)^\alpha$ ; and less strongly by the heavy element abundance and velocity gradient. Our major result is the dependence of $\zeta$ on frequency $\nu$ and the parameters $T_p$, $r_p$, and $\alpha$.