We present new CRIRES spectroscopic observations of BrGamma in the nuclear region of the Circinus galaxy, obtained with the aim of measuring the black hole (BH) mass with the spectroastrometric technique. The Circinus galaxy is an ideal benchmark for the spectroastrometric technique given its proximity and secure BH measurement obtained with the observation of its nuclear H2O maser disk. The kinematical data have been analyzed both with the classical method based on the analysis of the rotation curves and with the new method developed by us and based on spectroastrometry. The classical method indicates that the gas disk rotates in the gravitational potential of an extended stellar mass distribution and a spatially unresolved mass of (1.7 +- 0.2) 10^7 Msun, concentrated within r < 7 pc. The new method is capable of probing gas rotation at scales which are a factor ~3.5 smaller than those probed by the rotation curve analysis. The dynamical mass spatially unresolved with the spectroastrometric method is a factor ~2 smaller, 7.9 (+1.4 -1.1) 10^6 Msun indicating that spectroastrometry has been able to spatially resolve the nuclear mass distribution down to 2 pc scales. This unresolved mass is still a factor ~4.5 larger than the BH mass measurement obtained with the H2O maser emission indicating that it has not been possible to resolve the sphere of influence of the BH. Based on literature data, this spatially unresolved dynamical mass distribution is likely dominated by molecular gas and it has been tentatively identified with the circum-nuclear torus which prevents a direct view of the central BH in Circinus. This mass distribution, with a size of ~2pc, is similar in shape to that of the star cluster of the Milky Way suggesting that a molecular torus, forming stars at a high rate, might be the earlier evolutionary stage of the nuclear star clusters which are common in late type spirals.

We measure the black hole mass in the nearby active galaxy Centaurus A (NGC 5128) using a new method based on spectroastrometry of a rotating gas disk. The spectroastrometric approach consists in measuring the photocenter position of emission lines for different velocity channels. In a previous paper we focused on the basic methodology and the advantages of the spectroastrometric approach with a detailed set of simulations demonstrating the possibilities for black hole mass measurements going below the conventional spatial resolution. In this paper we apply the spectroastrometric method to multiple longslit and integral field near infrared spectroscopic observations of Centaurus A. We find that the application of the spectroastrometric method provides results perfectly consistent with the more complex classical method based on rotation curves: the measured BH mass is nearly independent of the observational setup and spatial resolution and the spectroastrometric method allows the gas dynamics to be probed down to spatial scales of ~0.02", i.e. 1/10 of the spatial resolution and ~1/50 of BH sphere of influence radius. The best estimate for the BH mass based on kinematics of the ionized gas is then log(MBH (sin i)^2/M\odot)=7.5 \pm 0.1 which corresponds to MBH = 9.6(+2.5-1.8) \times 10^7 M\odot for an assumed disk inclination of i = 35deg. The complementarity of this method with the classic rotation curve method will allow us to put constraints on the disk inclination which cannot be otherwise derived from spectroastrometry. With the application to Centaurus A, we have shown that spectroastrometry opens up the possibility of probing spatial scales smaller than the spatial resolution, extending the measured MBH range to new domains which are currently not accessible: smaller BHs in the local universe and similar BHs in more distant galaxies.

Galaxy dynamical masses are important physical quantities to constrain galaxy evolutionary models, especially at high redshifts. However, at z~2 the limited signal to noise ratio and spatial resolution of the data usually do not allow spatially resolved kinematical modeling and very often only virial masses can be estimated from line widths. But even such estimates require a good knowledge of galaxy size, which may be smaller than the spatial resolution. Spectroastrometry is a technique which combines spatial and spectral resolution to probe spatial scales significantly smaller than the spatial resolution of the observations. Here we apply it to the case of high-z galaxies and present a method based on spectroastrometry to estimate dynamical masses of high z galaxies, which overcomes the problem of size determination with poor spatial resolution. We construct and calibrate a "spectroastrometric" virial mass estimator, modifying the "classical" virial mass formula. We apply our method to the [O III] or H{\alpha} emission line detected in z~2-3 galaxies from AMAZE, LSD and SINS samples and we compare the spectroastrometric estimator with dynamical mass values resulting from full spatially resolved kinematical modeling. The spectroastrometric estimator is found to be a good approximation of dynamical masses, presenting a linear relation with a residual dispersion of only 0.15 dex. This is a big improvement compared to the "classical" virial mass estimator which has a non linear relation and much larger dispersion (0.47 dex) compared to dynamical masses. By applying our calibrated estimator to 16 galaxies from the AMAZE and LSD samples, we obtain masses in the ~10^7-10^10 M\odot range extending the mass range attainable with dynamical modeling.

We present a SINFONI integral field kinematical study of 33 galaxies at z~3 from the AMAZE and LSD projects which are aimed at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively large sample of galaxies at z~3. 11 galaxies show ordered rotational motions (~30% of the sample), in these cases we estimate dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses in the range 2 \times 10^9 M\odot - 2 \times 10^11 M\odot with a mean value of ~ 2 \times 10^10 M\odot. By comparing observed gas velocity dispersion with that expected from models, we find that most rotating objects are dynamically "hot", with intrinsic velocity dispersions of the order of ~90 km s-1. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion for the rotating objects is 1.6, much lower than observed in local galaxies value (~10) and slightly lower than the z~2 value (2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our measurements indicate that z~3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet "in place" at these early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking into account the velocity dispersion with the use of the S_0.5 indicator, suggesting that turbulent motions might have an important dynamical role.

It has recently been suggested that galaxies in the early Universe can grow through the accretion of cold gas, and that this may have been the main driver of star formation and stellar mass growth. Because the cold gas is essentially primordial, it has a very low abundance of elements heavier than helium (metallicity). As it is funneled to the centre of a galaxy, it will lead the central gas having an overall lower metallicity than gas further from the centre, because the gas further out has been enriched by supernovae and stellar winds, and not diluted by the primordial gas. Here we report chemical abundances across three rotationally-supported star-forming galaxies at z~3, only 2 Gyr after the Big Bang. We find an 'inverse' gradient, with the central, star forming regions having a lower metallicity than less active ones, opposite to what is seen in local galaxies. We conclude that the central gas has been diluted by the accretion of primordial gas, as predicted by 'cold flow' models.

We show that the mass-metallicity relation observed in the local universe is due to a more general relation between stellar mass M*, gas-phase metallicity and SFR. Local galaxies define a tight surface in this 3D space, the Fundamental Metallicity Relation (FMR), with a small residual dispersion of ~0.05 dex in metallicity, i.e, ~12%. At low stellar mass, metallicity decreases sharply with increasing SFR, while at high stellar mass, metallicity does not depend on SFR. High redshift galaxies, up to z~2.5 are found to follow the same FMR defined by local SDSS galaxies, with no indication of evolution. The evolution of the mass-metallicity relation observed up to z=2.5 is due to the fact that galaxies with progressively higher SFRs, and therefore lower metallicities, are selected at increasing redshifts, sampling different parts of the same FMR. By introducing the new quantity mu_alpha=log(M*)-alpha log(SFR), with alpha=0.32, we define a projection of the FMR that minimizes the metallicity scatter of local galaxies. The same quantity also cancels out any redshift evolution up to z~2.5, i.e, all galaxies have the same range of values of mu_0.32. At z>2.5, evolution of about 0.6 dex off the FMR is observed, with high-redshift galaxies showing lower metallicities. The existence of the FMR can be explained by the interplay of infall of pristine gas and outflow of enriched material. The former effect is responsible for the dependence of metallicity with SFR and is the dominant effect at high-redshift, while the latter introduces the dependence on stellar mass and dominates at low redshift. The combination of these two effects, together with the Schmidt-Kennicutt law, explains the shape of the FMR and the role of mu_0.32. The small metallicity scatter around the FMR supports the smooth infall scenario of gas accretion in the local universe.

This is the first in a series of papers in which we study the application of spectroastrometry in the context of gas kinematical studies aimed at measuring the mass of supermassive black holes. The spectroastrometrical method consists in measuring the photocenter of light emission in different wavelength or velocity channels. In particular we explore the potential of spectroastrometry of gas emission lines in galaxy nuclei to constrain the kinematics of rotating gas disks and to measure the mass of putative supermassive black holes. By means of detailed simulations and test cases, we show that the fundamental advantage of spectroastrometry is that it can provide information on the gravitational potential of a galaxy on scales significantly smaller (~ 1/10) than the limit imposed by the spatial resolution of the observations. We then describe a simple method to infer detailed kinematical informations from spectroastrometry in longslit spectra and to measure the mass of nuclear mass concentrations. Such method can be applied straightforwardly to integral field spectra, which do not have the complexities due to a partial spatial covering of the source in the case of longslit spectra.

We present the first results of a project, LSD, aimed at obtaining spatially-resolved, near-infrared spectroscopy of a complete sample of Lyman-Break Galaxies at z~3. Deep observations with adaptive optics resulted in the detection of the main optical lines, such as [OII], Hbeta and [OIII], which are used to study sizes, SFRs, morphologies, gas-phase metallicities, gas fractions and effective yields. Optical, near-IR and Spitzer/IRAC photometry is used to measure stellar mass. We obtain that morphologies are usually complex, with the presence of several peaks of emissions and companions that are not detected in broad-band images. Typical metallicities are 10-50% solar, with a strong evolution of the mass-metallicity relation from lower redshifts. Stellar masses, gas fraction, and evolutionary stages vary significantly among the galaxies, with less massive galaxies showing larger fractions of gas. In contrast with observations in the local universe, effective yields decrease with stellar mass and reach solar values at the low-mass end of the sample. This effect can be reproduced by gas infall with rates of the order of the SFRs. Outflows are present but are not needed to explain the mass-metallicity relation. We conclude that a large fraction of these galaxies are actively creating stars after major episodes of gas infall or merging.