Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics.

We use single-epoch spectroscopy of three gravitationally lensed quasars, HE0435-1223, WFI2033-4723, and HE2149-2745, to study their inner structure (BLR and continuum source). We detect microlensing-induced magnification in the wings of the broad emission lines of two of the systems (HE0435-1223 and WFI2033-4723). In the case of WFI2033-4723, microlensing affects two "bumps" in the spectra which are almost symmetrically arranged on the blue (coincident with an AlIII emission line) and red wings of CIII]. These match the typical double-peaked profile that follows from disk kinematics. The presence of microlensing in the wings of the emission lines indicates the existence of two different regions in the BLR: a relatively small one with kinematics possibly related to an accretion disk, and another one that is substantially more extended and insensitive to microlensing. There is good agreement between the estimated size of the region affected by microlensing in the emission lines, $r_s=10^{+15}_{-7} \sqrt{M/M_{\odot}}$ light-days (red wing of CIV in HE0435-1223) and $r_s=11^{+28}_{-7} \sqrt{M/M_{\odot}}$ light-days (CIII] bumps in WFI2033-4723) with the sizes inferred from the continuum emission, $r_s=13^{+5}_{-4} \sqrt{M/M_{\odot}}$ light-days (HE0435-1223) and $r_s=10^{+3}_{-2} \sqrt{M/M_{\odot}}$ light-days (WFI2033-4723). For HE2149-2745 we measure an accretion disk size $r_s=8^{+11}_{-5} \sqrt{M/M_{\odot}}$ light-days. The estimates of $p$, the exponent of the size vs. wavelength ($r_s\propto\lambda^p$), are $1.2\pm0.6$, $0.8\pm0.2$, and $0.4\pm0.3$ for HE0435-1223, WFI2033-4723, and HE2149-2745, respectively. In conclusion, the continuum microlensing amplitude in the three quasars and chromaticity in WFI2033-4723 and HE2149-2745 are below expectations for the thin disk model. The disks are larger and their temperature gradients are flatter than predicted by this model.

We present new mid-IR observations of the quadruply lensed quasar Q2237+0305 taken with CanariCam on the Gran Telescopio Canarias. Mid-IR emission by hot dust, unlike the optical and near-IR emission from the accretion disk, is unaffected by the interstellar medium (extinction/scattering) or stellar microlensing. We compare these "true" ratios to the (stellar) microlensed flux ratios observed in the optical/near-IR to constrain the structure of the quasar accretion disk. We find a half-light radius of $R_{1/2}=3.4_{-2.1}^{+5.3}\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days at $\lambda_{rest}=1736$ {\AA}, and an exponent for the temperature profile $R \propto \lambda^{p}$ of $p=0.79\pm0.55$, where $p=4/3$ for a standard thin-disk model. If we assume that the differences in the mid-IR flux ratios measured over the years are due to microlensing variability, we find a lower limit for the size of the mid-IR-emitting region of $R_{1/2} \gtrsim 200\,\sqrt{\langle M \rangle/0.3\,\rm{M_{\odot}}}$ light-days. We also test for the presence of substructure/satellites by comparing the observed mid-IR flux ratios with those predicted from smooth lens models. We can explain the differences if the surface density fraction in satellites near the lensed images is $\alpha = 0.033_{-0.019}^{+0.046}$ for a singular isothermal ellipsoid plus external shear mass model or $\alpha = 0.013_{-0.008}^{+0.019}$ for a mass model combining ellipsoidal NFW and de Vaucouleurs profiles in an external shear.

We propose to use the flux variability of lensed quasar images induced by gravitational microlensing to measure the transverse peculiar velocity of lens galaxies over a wide range of redshift. Microlensing variability is caused by the motions of the observer, the lens galaxy (including the motion of the stars within the galaxy), and the source; hence, its frequency is directly related to the galaxy's transverse peculiar velocity. The idea is to count time-event rates (e.g., peak or caustic crossing rates) in the observed microlensing light curves of lensed quasars that can be compared with model predictions for different values of the transverse peculiar velocity. To compensate for the large time-scale of microlensing variability we propose to count and model the number of events in an ensemble of gravitational lenses. We develop the methodology to achieve this goal and apply it to an ensemble of 17 lensed quasar systems. In spite of the shortcomings of the available data, we have obtained tentative estimates of the peculiar velocity dispersion of lens galaxies at $z\sim 0.5$, $\sigma_{\rm pec}(0.53\pm0.18)\simeq(638\pm213)\sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$. Scaling at zero redshift we derive, $\sigma_{\rm pec}(0)\simeq(491\pm164) \sqrt{\langle m \rangle/0.3 M_\odot} \, \rm km\, s^{-1}$, consistent with peculiar motions of nearby galaxies and with recent $N$-body nonlinear reconstructions of the Local Universe based on $\Lambda$CDM. We analyze the different sources of uncertainty of the method and find that for the present ensemble of 17 lensed systems the error is dominated by Poissonian noise, but that for larger ensembles the impact of the uncertainty on the average stellar mass may be significant.

We present estimates for the size and the logarithmic slope of the disk temperature profile of the lensed quasar Q2237+0305, independent of the component velocities. These estimates are based on six epochs of multi-wavelength narrowband images from the Nordic Optical Telescope. For each pair of lensed images and each photometric band, we determine the microlensing amplitude and chromaticity using pre-existing mid-IR photometry to define the baseline for no microlensing magnification. A statistical comparison of the combined microlensing data (6 epochs $\times$ 5 narrow bands $\times$ 6 image pairs) with simulations based on microlensing magnification maps gives Bayesian estimates for the half-light radius of $R_{1/2}=8.5^{+7.5}_{-4.0}\sqrt{ \langle M \rangle/0.3\, M_\odot}$ light-days, and $p=0.95\pm0.33$ for the exponent of the logarithmic temperature profile $T\propto R^{ -1/p}$. This size estimate is in good agreement with most recent studies. Other works based on the study of single microlensing events predict smaller sizes, but could be statistically biased by focusing on high-magnification events.

We study three high magnification microlensing events, generally recognized as probable caustic crossings, in the optical light curves of the multiply imaged quasar Q 2237+0305. We model the light curve of each event as the convolution of a standard thin disk luminosity profile with a straight fold caustic. We also allow for a linear gradient that can account for an additional varying background effect of microlensing. This model not only matches noticeably well the global shape of each of the three independent microlensing events but also gives remarkably similar estimates for the disk size parameter. The measured average half-light radius, $R_{1/2}=(3.0\pm 1.5)\sqrt{M/0.3M\odot}$ light-days, agrees with previous estimates. In the three events, the core of the magnification profile exhibits "fine structure" related to the innermost region of the accretion disk (located at a radial distance of $2.7\pm 1.4$ Schwarzschild radii according to our measurement). Relativistic beaming at the internal rim of the accretion disk can explain the shape and size of the fine structure, although alternative explanations are also possible. This is the first direct measurement of the size of a structure, likely the innermost stable circular orbit, at $\sim 3$ Schwarzschild radii in a quasar accretion disk. The monitoring of thousands of lensed quasars with future telescopes will allow the study of the event horizon environment of black holes in hundreds of quasars in a wide range of redshifts $(0.5<z<5)$.

We use the statistics of caustic crossings induced by microlensing in the lens system Q~2237+0305 to study the lens galaxy peculiar velocity. We calculate the caustic crossing rates for a comprehensive family of stellar mass functions and find a dependence of the average number of caustic crossings with the effective transverse velocity and the average mass, $\langle n \rangle \propto {v_{eff} / \sqrt{\langle m \rangle}}$, equivalent to the theoretical prediction for the case of microlenses with identical masses. We explore the possibilities of the method to measure $v_{eff}$ using the $\sim$12 years of OGLE monitoring of the four images of Q 2237+0305. To determine a lower limit for $v_{eff}$ we count, conservatively, a single caustic crossing for each one of the 4 high magnification events identified in the literature (plus one additional proposed by us) obtaining $v_{eff} \gtrsim 240\sqrt{\langle m \rangle/0.17M_\odot}\rm\,km\, s^{-1}$ at 68\% of confidence. From this value and the average $FWHM$ of 4 high magnification events we obtain a lower limit of $r_s \gtrsim 1.4 \sqrt{\langle m \rangle/0.17M_\odot}$ light-days for the radius of the source ($r_s=FWHM/2.35$). Tentative identification of 3 additional caustic crossing events leads to estimates of $v_{eff}\simeq (493\pm 246)\sqrt{\langle m \rangle/0.17M_\odot}\rm\,km\, s^{-1}$ for the effective transverse velocity and of $r_s \simeq (2.7\pm 1.3)\sqrt{\langle m \rangle/0.17M_\odot}$ light-days for the source size. The estimated transverse peculiar velocity of the galaxy is $v_t \simeq(429\pm246)\sqrt{\langle m \rangle/0.17M_\odot}\rm\,km\, s^{-1}$.

We measure the differential microlensing of the UV Fe II and Fe III emission line blends between 14 quasar image pairs in 13 gravitational lenses. We find that the UV iron emission is strongly microlensed in 4 cases with amplitudes comparable to that of the continuum. Statistically modeling the magnifications we infer a typical size of r ~ 4*sqrt(M/Msun) light-days for the Fe line emitting regions which is comparable to the size of the region generating the UV continuum (3 ~ 7 light-days). This may indicate that a significant part of the UV Fe II and Fe III emission originates in the quasar accretion disk.

We measure the differential microlensing of the broad emission lines between 18 quasar image pairs in 16 gravitational lenses. We find that high ionization lines such as CIV are more strongly microlensed than low ionization lines, indicating that the high ionization line emission regions are more compact. If we statistically model the distribution of microlensing magnifications, we obtain estimates for the broad line region radius of 24 (-15/+22) and 55 (-35/+150) light-days (90% confidence) for the high and low ionization lines, respectively. When the sample is divided attending to quasar luminosity, we find that the line emission regions of more luminous quasars are larger, with a slope consistent with the expected scaling from photoionization models. Our estimates also agree well with the results from local reveberation mapping studies.

Using microlensing measurements from a sample of 27 image-pairs of 19 lensed quasars we determine a maximum likelihood estimate for the accretion disk size of an {{\em}average} quasar of $r_s=4.0^{+2.4}_{-3.1} $ light days at rest frame $<\lambda>=1736$\AA\ for microlenses with a mean mass of $<M>=0.3M_\odot$. This value, in good agreement with previous results from smaller samples, is roughly a factor of 5 greater than the predictions of the standard thin disk model. The individual size estimates for the 19 quasars in our sample are also in excellent agreement with the results of the joint maximum likelihood analysis.

We analyse new integral field spectroscopy of the inner region (central 2.5 kpc) of the spiral galaxy NGC 4321 to study the peculiar kinematics of this region. Fourier analysis of the velocity residuals obtained by subtracting an axisymmetric rotation model from the $\rm H\alpha$ velocity field, indicates that the distortions are {\em global} features generated by an $m=2$ perturbation of the gravitational potential which can be explained by the nuclear bar. This bar has been previously observed in the near-infrared but not in the optical continuum dominated by star formation. We detect the optical counterpart of this bar in the 2D distribution of the old stellar population (inferred from the equivalent width map of the stellar absorption lines). We apply the Tremaine--Weinberg method to the stellar velocity field to calculate the pattern speed of the inner bar, obtaining a value of $\Omega_b$=160$\pm70\rm km s^{-1} kpc^{-1} $. This value is considerably la rger than the one obtained when a simple bar model is considered. However the uncertainties in the pattern speed determination prevent us to give support to alternative scenarios.

We present a homogeneous set of stellar atmospheric parameters Teff, log g, [Fe/H] for MILES, a new spectral stellar library covering the range 3525 - 7500 angstrom at 2.3 angstrom (FWHM) spectral resolution. The library consists of 985 stars spanning a large range in atmospheric parameters, from super metal-rich, cool stars to hot, metal-poor stars. The spectral resolution, spectral type coverage and number of stars represent a substantial improvement over previous libraries used in population synthesis models. The atmospheric parameters that we present here are the result of a previous, extensive compilation from the literature. In order to construct a homogeneous dataset of atmospheric parameters we have taken the sample of stars of Soubiran, Katz & Cayrel, which has very well determined fundamental parameters, as the standard reference system for our field stars, and have calibrated and bootstrapped the data from other papers against it. The atmospheric parameters for our cluster stars have also been revised and updated according to recent metallicity scales, colour-temperature relations and improved set of isochrones.

A new stellar library developed for stellar population synthesis modeling is presented. The library consist of 985 stars spanning a large range in atmospheric parameters. The spectra were obtained at the 2.5m INT telescope and cover the range 3525-7500A at 2.3A (FWHM) spectral resolution. The spectral resolution, spectral type coverage, flux calibration accuracy and number of stars represent a substantial improvement over previous libraries used in population synthesis models.

We present integral field spectroscopy of NGC 5033, a low luminosity Seyfert galaxy. The observations were made with INTEGRAL, a fiber based system operating at the WHT. The intensity map of the H$\beta$ emission line represents a spiral or ring-like pattern of HII regions. On the contrary, the [OIII] intensity map morphology is markedly anisotropic. The strong morphological differences imply that the [OIII] emitters represent highly ionized gas illuminated by the central source. The [OIII] map morphology is compatible with a biconical structure of ionization induced by strong extinction in the galaxy disc that also obscures half of the spheroidal stellar bulge. We identify the spectrum corresponding to the Seyfert 1 nucleus from the presence of H$\beta$ broad emission lines. This spectrum is located in a region where strong extinction is expected but exhibits the bluest spectral energy distribution. The Seyfert 1 nucleus seems to be offcenter with respect to the stellar rotation center. This result has been also found in other Seyfert galaxies and interpreted in terms of a past merger. The offcentering could indicate the presence of nonsymmetric departures in the gravitational potential which could be fueling the active nucleus. The kinematics of the [OIII] emitters show important deviations at a kpc scale with respect to the stellar velocity field and show features related to the asymmetrical morphology of the high ionization region.

We propose a hypothesis of the truncation of stellar discs based on the magnetic model of the rotation curve of spiral galaxies. Once the disc had formed and acquired its present structure, approximately, three balanced forces were acting on the initial gas: gravity and magnetic forces, inwards, and centrifugal force. When stars are formed from this gas, the magnetic force is suddenly suppressed. Gravitation alone cannot retain the newly-formed stars and at birth places beyond a certain galactocentric radius they escape to intergalactic space. This radius is the so-called ``truncation radius'', which is predicted to be at about 4-5 disc radial scale lengths, in promising agreement with observations.

We present a statistical study of the presence of galaxy warps in the Hubble deep fields. Among a complete sample of 45 edge-on galaxies above a diameter of 1.''3, we find 5 galaxies to be certainly warped and 6 galaxies as good candidates. In addition, 4 galaxies reveal a characteristic U-warp. Compared to statistical studies of local warps, and taking into account the strong bias against observing the outer parts of galaxies at high redshift, these numbers point towards a very high frequency of warps at z \sim 1: almost all galaxy discs might be warped. Furthermore, the amplitude of warps are stronger than for local warps. This is easily interpreted in terms of higher galaxy interactions and matter accretion in the past. This result supports these two mechanisms as the best candidates for the origin of early warps. The mean observed axis ratio of our sample of edge-on galaxies is significantly larger in the high-z sample than is found for samples of local spiral galaxies. This might be due to disk thickening due to more frequent galaxy interactions.

We present spectroscopic observations of the stellar motions in the disk of the superthin edge-on spiral galaxy IC 5249 and re-analyse synthesis observations of the HI. We find that the HI rotation curve rises initially to about 90-100 km/s, but contrary to the conclusion of Abe et al. (1999) flattens well before the edge of the optical disk. Over most part of the optical disk we have been able to establish that the (tangential) stellar velocity dispersion is 25-30 km/s. From earlier surface photometry we adopt a value for the radial scalelength of the disk of 7 +/- 1 kpc, a vertical scaleheight of 0.65 +/- 0.05 kpc and a disk truncation radius of 17 +/- 1 kpc. The very thin appearance of IC 5249 on the sky is the result of a combination of a low (face-on) surface brightness, a long scalelength and a a sharp truncation at only about 2.5 scalelengths. From various arguments we derive the stellar velocity dispersions at a position one radial scalelength as sigma_R about 35 km/s, sigma_{theta} about 30 km/s and sigma_z about 20 km/s. This is comparable to the values for the disk of our Galaxy in the solar neighborhood.

We present a first study of truncation of the stellar disks of spiral galaxies in the near infrared. Observations of NGC4013, NGC4217, NGC6504 and NGC5981 were made with the CAIN NIR camera on the CST in Tenerife. This wavelength range provides the best description of the phenomenon, not only because extinction effects are minimized, but also because the distribution of the old stellar population is directly obtained. The four galaxies are edge-on and an inversion method was developed to obtain the deprojected profiles. We did not assume any model of the different galactic components. The ``truncation curve'', i.e. T(R)= mu(R)- mu_{D}(R), where mu is the actual surface brightness in mag/arcsec^2 and mu_{D} the exponential disk surface brightness, has been obtained with unprecedented precision. It is suggested that T(R) is proportional to (R_{t}-R)^{-1}, where R_{t} is the truncation radius, i.e. the radius beyond which no star is observed.

We have observed the nearly face-on spiral galaxy NGC 5668 with the TAURUS II Fabry-Perot interferometer at the William Herschel Telescope using the $H\alpha$ line to study the kinematics of the ionized gas. From the extracted data cube we construct intensity, velocity and velocity dispersion maps. We calculate the rotation curve in the innermost 2 arcmin and we use the residual velocity field to look for regions with important vertical motions. By comparing the geometry of these regions in the residual velocity field with the geometry in the intensity and velocity dispersion maps we are able to select some regions which are very likely to be shells or chimneys in the disk. The geometry and size of these regions are very similar to the shells or chimneys detected in other galaxies by different means. Moreover, it is worth noting than this galaxy has been reported to have a population of neutral hydrogen high velocity clouds (Schulman et al. 1996) which, according to these observations, could have been originated by chimneys similar to those reported in this paper.

The nearly face-on spiral galaxy NGC 3938 has been observed in the $H_\alpha$ line with the TAURUS II Fabry-Perot interferometer at the William Herschel Telescope in order to study the kinematics of the ionized gas. We are able to construct intensity, velocity and velocity dispersion maps for this galaxy. The rotation curve of the galaxy is calculated up to 4.5 radial scale lengths from the galactic centre. The residual velocity field shows very small values with no systematic pattern. The mean velocity dispersion is approximately constant with radius at about 11 km/s as previously reported for the neutral and molecular gas. We have also studied the relation between intensity and velocity dispersion for the ionized gas. We have found that this distribution is compatible with a turbulent gas relaxing to a Kolmogorov type turbulence as the stationary regime. The average dispersion varies with intensity as $\sigma \propto I^{1/8}$ although it becomes much steeper at higher intensities, for which the dispersion is kept almost constant at a value of about 19 km/s.

Assuming the magnetic hypothesis for the formation of warps, we deduce a very simple formula for the warp curve in an idealized scenario. According to this formula the warp rises as the third power of radius in the innermost warped region, reaches the maximum slope at intermediate radii and has an asymptotic slope coincident with the direction of the extragalactic magnetic field. In most cases, however, the galaxy's limited size prevents the observation of the full curve. Even though the model is highly simplified, it basically reproduces real warp curves, in particular the 21 cm warp curve of NGC 5907. If the magnetic model were considered to be correct, the fitting of warp curves could allow rough estimations of the strength and direction of the magnetic field. We also propose a magnetic field distribution for the outermost part of the galaxy, matching the boundary conditions of being azimuthal inside, and constant at infinity. We use this magnetic field distribution to show that the assumptions made to obtain the warp curve with our simple model cannot introduce important errors.

Observations of the Local Arm have been carried out in the near-infrared with the 1.5 m ``Carlos S\'anchez'' telescope in Tenerife. A model of the disc with adjustable parameters fitted to reproduce the DIRBE-COBE survey, has been subtracted from the observational data in order to obtain a clean map of the Local Arm, uncontaminated by other components of the disc. The arm is more than 70 pc over the plane and is wider than 200 pc. At a latitude of about $80^\circ$ the deviation of the arm from the the galactic plane is so large that was only partially observed with our observational window of $\triangle b = \pm 6^\circ$. The elevation over the plane vanishes and the width decreases as the Arm comes closer to the Sun.

The truncation of stellar discs is not abrupt but characterized by a continuous distancing from the exponential profile. There exists a truncation curve, $t(r)$, ending at a truncation radius, $r_t$. We present here a theoretical model in which it is assumed that the magnetic hypothesis explaining the flat rotation curve also explains the truncation. Once stars are born, the centripetal magnetic force previously acting on the progenitor gas cloud is suddenly interrupted, and stars must move to larger orbits or escape. The agreement between theoretical and observed truncation curves is very satisfactory. Parameters defining the disc gas rotation curve should therefore be related to those defining the truncation. It is predicted that rotation curves that quickly reach the asymptotic value $\theta_0 = \theta (r=\infty)$ would have small truncation radii. On the contrary, $r_t$ and $\theta_0$ itself, would be uncorrelated quantities.

We propose new parameters to describe the geometry of a warped disc when viewed edge-on, a global warp parameter $w$, and a family of three parameters $A$, $B$ and $C$ that describe, independently, the shape of the warp. These parameters are useful to detect some key effects in the understanding of warps. We have also developed software (WIG) which is able to calculate these parameters for a given image of a galaxy in any wavelength and which is described here with some examples, namely ESO 235-53 in the optical, NGC 4013 in 21 cm and NGC 4565 in 1.2 mm continuum.

The radial scale length of the exponential component of the disc of the Milky Way has been determined in the near infrared. We have used the TMGS (Two Micron Galactic Survey) database which contains positions and K-magnitudes of about 700.000 stars distributed in serveral regions along the galactic plane. From those we have selected areas more than $5^\circ$ off the plane to minimize the effect of extinction and the contributions of the young disc; and with longitudes ranging from $30^\circ$ to $70^\circ$, to avoid contaminations from the central bulge, the bar and the molecular ring, in the inner end, and from the local arm, the warp and the truncation of the disc, in the outer end. We observed stars with magnitude $ 9 \leq m_K \leq 10 $. The use of the NIR K-band also reduces the effect of extinction. In the observation region, m_K=9.5 mag stars are K2-K5 III stars, with an absolute magnitude that is nearly constant, which also greatly simplifies the problem. We have obtained the value $2.1 \pm 0.3 kpc$ for the radial scale length, which is a typical value when compared with external galaxies of similar types.

We consider that no mean magnetic field exists during this epoch, but that there is a mean magnetic energy associated with large-scale magnetic inhomogeneities. We study the evolution of these inhomogeneities and their influence on the large scale density structure, by introducing linear perturbations in Maxwell equations, the conservation of momentum-energy equation, and in Einstein field equations. The primordial magnetic field structure is time independent in the linear approximation, only being diluted by the general expansion, so that $\vec{B}R^2$ is conserved in comoving coordinates. Magnetic fields have a strong influence on the formation of large-scale structure. Firstly, relatively low fields are able to generate density structures even if they were inexistent at earlier times. Second, magnetic fields act anisotropically more recently, modifying the evolution of individual density clouds. Magnetic flux tubes have a tendency to concentrate photons in filamentary patterns.