We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.9<z<4.6 and provide material to prepare future observations of the faint Universe. We use data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median M$_{UV}$ of -18 and down to M$_{UV}$ of -16), low stellar mass (median value of $10^{8.4}$ and down to $10^{7}M_{\odot}$) LAEs at redshift z>3. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.

Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time.

Here we describe a simple, efficient, and most importantly fully operational point-spread-function(PSF)-reconstruction approach for laser-assisted ground layer adaptive optics (GLAO) in the frame of the Multi Unit Spectroscopic Explorer (MUSE) Wide Field Mode. Based on clear astrophysical requirements derived by the MUSE team and using the functionality of the current ESO Adaptive Optics Facility we aim to develop an operational PSF-reconstruction (PSFR) algorithm and test it both in simulations and using on-sky data. The PSFR approach is based on a Fourier description of the GLAO correction to which the specific instrumental effects of MUSE Wide Field Mode (pixel size, internal aberrations, etc.) have been added. It was first thoroughly validated with full end-to-end simulations. Sensitivity to the main atmospheric and AO system parameters was analysed and the code was re-optimised to account for the sensitivity found. Finally, the optimised algorithm was tested and commissioned using more than one year of on-sky MUSE data. We demonstrate with an on-sky data analysis that our algorithm meets all the requirements imposed by the MUSE scientists, namely an accuracy better than a few percent on the critical PSF parameters including full width at half maximum and global PSF shape through the kurtosis parameter of a Moffat function. The PSFR algorithm is publicly available and is used routinely to assess the MUSE image quality for each observation. It can be included in any post-processing activity which requires knowledge of the PSF.

We present recent results obtained with the VLT/MUSE Integral Field Spectrograph fed by the 4LGSF and its laser tomography adaptive optics module GALACSI. While this so-called narrow-field mode of MUSE was not designed to perform directly imaging of exoplanets and outflows, we show that it can be a game changer to detect and characterize young exoplanets with a prominent emission lines (i.e H{\alpha}, tracer of accretion), at moderate contrasts. These performances are achieved thanks to the combo of a near-diffraction limited PSF and a medium resolution spectrograph and a cross-correlation approach in post-processing . We discuss this in the context of ground and space, infrared and visible wavelengths, preparing for missions like JWST and WFIRST in great synergy and as pathfinder for future ELT/GSMT (Extremely Large and/or Giant Segmented Mirror Telescopes) instruments.

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

We present spatially resolved maps of six individually-detected Lyman alpha haloes (LAHs) as well as a first statistical analysis of the Lyman alpha (Lya) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies using MUSE. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lya line profile. Using a three-dimensional two-component model for the Lya emission, we measure the full width at half maximum (FWHM), the peak velocity shift and the asymmetry of the Lya line in the core and in the halo of 19 galaxies. We find that the Lya line shape is statistically different in the halo compared to the core for ~40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lya line both at the interstellar and circum-galactic scales. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies show broader, more redshifted and less asymmetric Lya lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lya lines. The generally broad and red line shapes found in the halo component suggests that the Lya haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lya line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.

We aim to provide spectroscopic evidence regarding the nature of the putative star cluster in Eridanus 2 and to place constraints on the mass and abundance of massive astrophysical compact halo objects (MACHOs) as a constituent of dark matter. Methods. We present spectroscopic observations of the central square arcminute of Eridanus 2 from MUSE-Faint, a survey of ultra-faint dwarf galaxies with the Multi Unit Spectroscopic Explorer on the Very Large Telescope. We derive line-of-sight velocities for possible member stars of the putative cluster and for stars in the centre of Eridanus 2. We discuss the existence of the cluster and determine new constraints for MACHOs using the Fokker-Planck diffusion approximation. Results. Out of 182 extracted spectra, we identify 26 member stars of Eridanus 2, seven of which are possible cluster members. We find intrinsic mean line-of-sight velocities of $79.7^{+3.1}_{-3.8}\,\mathrm{km}\,\mathrm{s}^{-1}$ and $76.0^{+3.2}_{-3.7}\,\mathrm{km}\,\mathrm{s}^{-1}$ for the cluster and the bulk of Eridanus 2, respectively, and intrinsic velocity dispersions of ${<}7.6\,\mathrm{km}\,\mathrm{s}^{-1}$ (68-$\%$ upper limit) and $10.3^{+3.9}_{-3.2}\,\mathrm{km}\,\mathrm{s}^{-1}$, respectively. This indicates the cluster most likely exists as a distinct dynamical population hosted by Eridanus 2, without surplus of dark matter over the background distribution. Among the member stars in the bulk of Eridanus 2, we find possible carbon stars, alluding to the existence of an intermediate-age population. We derive constraints on the fraction of dark matter that can consist of MACHOs with a given mass between $1$-$10^5\,M_\mathrm{sun}$. For dark matter consisting purely of MACHOs, the mass of the MACHOs must be less than ${\sim}7.6\,M_\mathrm{sun}$ and ${\sim}44\,M_\mathrm{sun}$ at a $68$- and $95$-$\%$ confidence level, respectively. (Abridged)

One of the major science cases of the MUSE integral field spectrograph is the detection of Lyman-alpha emitters at high redshifts. The on-going and planned deep fields observations will allow for one large sample of these sources. An efficient tool to perform blind detection of faint emitters in MUSE datacubes is a prerequisite of such an endeavor. Several line detection algorithms exist but their performance during the deepest MUSE exposures is hard to quantify, in particular with respect to their actual false detection rate, or purity. {The aim of this work is to design and validate} an algorithm that efficiently detects faint spatial-spectral emission signatures, while allowing for a stable false detection rate over the data cube and providing in the same time an automated and reliable estimation of the purity. Results on simulated data cubes providing ground truth show that the method reaches its aims in terms of purity and completeness. When applied to the deep 30-hour exposure MUSE datacube in the Hubble Ultra Deep Field, the algorithms allows for the confirmed detection of 133 intermediate redshifts galaxies and 248 Lyman Alpha Emitters, including 86 sources with no HST counterpart. The algorithm fulfills its aims in terms of detection power and reliability. It is consequently implemented as a Python package whose code and documentation are available on GitHub and readthedocs.

NGC6240 is a well-studied nearby galaxy system in the process of merging. Based on optical, X-ray, and radio observations, it is thought to harbor two active nuclei We carried out a detailed optical 3D spectroscopic study to investigate the inner region of this system in connection with existing MERLIN and VLBA data. We observed NGC6240 with very high spatial resolution using the MUSE instrument in the Narrow-Field Mode with the four-laser GALACSI adaptive optics system on the ESO VLT under seeing conditions of 0.49 arcsec. Our 3D spectra cover the wavelength range from 4725 to 9350 AA at a spatial resolution of ~75 mas. We report the discovery of three nuclei in the final state of merging within a region of only 1 kpc in the NGC6240 system.Thanks to MUSE we are able to show that the formerly unresolved southern component actually consists of two distinct nuclei separated by only 198 pc. In combination with Gaia data we reach an absolute positional accuracy of only 30 mas that is essential to compare optical spectra with MERLIN and VLBA radio positions. The verification and detailed study of a system with three nuclei, two of which are active and each with a mass in excess of $9\times10^{7} M_{\odot}$, is of great importance for the understanding of hierarchical galaxy formation via merging processes since multiple mergers lead to a faster evolution of massive galaxies in comparison to binary mergers. So far it has been suggested that the formation of galactic nuclei with multiple supermassive black holes (SMBHs) is expected to be rare in the local universe.Triple massive black hole systems might be of fundamental importance for the coalescence of massive black hole binaries in less than a Hubble time leading to the loudest sources of gravitational waves in the millihertz regime.

We present results from our on-going MusE GAs FLOw and Wind (MEGAFLOW) survey, which consists of 22 quasar lines-of-sight, each observed with the integral field unit (IFU) MUSE and the UVES spectrograph at the ESO Very Large Telescopes (VLT). The goals of this survey are to study the properties of the circum-galactic medium around $z\sim1$ star-forming galaxies. The absorption-line selected survey consists of 79 strong \MgII\ absorbers (with rest-frame equivalent width (REW)$\gtrsim$0.3\AA) and, currently, 86 associated galaxies within 100 projected~kpc of the quasar with stellar masses ($M_\star$) from $10^9$ to $10^{11}$ \msun. We find that the cool halo gas traced by \MgII\ is not isotropically distributed around these galaxies, as we show the strong bi-modal distribution in the azimuthal angle of the apparent location of the quasar with respect to the galaxy major-axis. This supports a scenario in which outflows are bi-conical in nature and co-exist with a coplanar gaseous structure extending at least up to 60 to 80 kpc. Assuming that absorbers near the minor axis probe outflows, the current MEGAFLOW sample allowed us to select 26 galaxy-quasar pairs suitable for studying winds. From this sample, using a simple geometrical model, we find that the outflow velocity only exceeds the escape velocity when $M_{\star}\lesssim 4\times10^9$~\msun, implying the cool material is likely to fall back except in the smallest halos. Finally, we find that the mass loading factor $\eta$, the ratio between the ejected mass rate and the star formation rate (SFR), appears to be roughly constant with respect to the galaxy mass.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

We present an analysis of HI Lyman-alpha emission in deep VLT/MUSE observations of two highly magnified and extended galaxies at z=3.5 and 4.03, including a newly discovered, almost complete Einstein ring. While these Lyman-alpha haloes are intrinsically similar to the ones typically seen in other MUSE deep fields, the benefits of gravitational lensing allows us to construct exceptionally detailed maps of Lyman-alpha line properties at sub-kpc scales. By combining all multiple images, we are able to observe complex structures in the Lyman-alpha emission and uncover small (~ 120 km/s in Lyman-alpha peak shift), but significant at > 4 sigma, systematic variations in the shape of the Lyman-alpha line profile within each halo. Indeed, we observe a global trend for the line peak shift to become redder at large radii, together with a strong correlation between the peak wavelength and line width. This systematic intrahalo variation is markedly similar to the object-to-object variations obtained from the integrated properties of recent large samples. Regions of high surface brightness correspond to relatively small line shifts, which could indicate that Lyman-alpha emission escapes preferentially from regions where the line profile has been less severely affected by scattering of Lyman-alpha photons.

Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.

We present the first measurements of the Lyman-continuum photon production efficiency $\xi_{\textrm{ion,0}}$ at $z\sim4$-5 for galaxies fainter than 0.2 $L^*$ ($-$19 mag). $\xi_{\textrm{ion,0}}$ quantifies the production rate of ionizing photons with respect to the UV luminosity density assuming a fiducial escape fraction of zero. Extending previous measurements of $\xi_{\textrm{ion,0}}$ to the faint population is important, as ultra-faint galaxies are expected to contribute the bulk of the ionizing emissivity. We probe $\xi_{\textrm{ion,0}}$ to such faint magnitudes by taking advantage of 200-hour depth Spitzer/IRAC observations from the GREATS program and $\approx$300 3<$z$<6 galaxies with spectroscopic redshifts from the MUSE GTO Deep + Wide programs. Stacked IRAC [3.6]$-$[4.5] colors are derived and used to infer the H$\alpha$ rest-frame equivalent widths, which range from 403\r{A} to 2818\r{A}. The derived $\xi_{\textrm{ion,0}}$ is $\log_{10}(\xi_{\textrm{ion,0}} / \textrm{Hz erg}^{-1}) = 25.36 \pm 0.08$ over $-$20.5 < M$_{\textrm{UV}}$ < $-$17.5, similar to those derived for brighter galaxy samples at the same redshift and therefore suggesting that $\xi_{\textrm{ion}}$ shows no strong dependence on $M_{UV}$. The $\xi_{\textrm{ion,0}}$ values found in our sample imply that the Lyman-continuum escape fraction for $M_{\textrm{UV}} \approx -19$ star-forming galaxies cannot exceed $\approx$8-20\% in the reionization era.

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360<lambda<1330 nm with options for fiber-fed spectrographs at high resolution and a panoramic IFU at a separate focus. This would enable transformational progress via its ability to access a larger fraction of objects from Gaia, LSST, Euclid, and WFIRST than any currently funded or planned spectroscopic facility. An ESO-sponsored study (arXiv:1701.01976) discussed the scientific potential in ambitious new spectroscopic surveys in Galactic astronomy, extragalactic astronomy, and cosmology. The US community should establish links with European and other international communities to plan for such a powerful facility and maximize the potential of large aperture multi-object spectroscopy given the considerable investment in deep imaging surveys.

Newly forming proto-planets are expected to create cavities and substructures in young, gas-rich proto-planetary disks, but they are difficult to detect as they could be confused with disk features affected by advanced image-analysis techniques. Recently, a planet was discovered inside the gap of the transitional disk of the T-Tauri star PDS 70. Here we report on the detection of strong H-alpha emission from two distinct locations in the PDS 70 system, one corresponding to the previously discovered planet PDS 70 b, which confirms the earlier H$\alpha$ detection, and another located close to the outer-edge of the gap, coinciding with a previously identified bright dust spot in the disk and with a small opening in a ring of molecular emission. We identify this second H$\alpha$ peak as a second proto-planet in the PDS 70 system. The H$\alpha$ emission spectra of both proto-planets indicate ongoing accretion onto the proto-planets, which appear to be near a 2:1 mean motion resonance. Our observations show that adaptive-optics-assisted, medium-resolution, integral-field spectroscopy with MUSE targeting accretion signatures will be a powerful way to trace ongoing planet formation in transitional disks at different stages of their evolution. Finding more young planetary systems in mean motion resonance would give credibility to the Grand Tack hypothesis in which Jupiter and Saturn migrated in a resonance orbit during the early formation period of our Solar System.

Observations of Neptune, made in 2018 using the new Narrow Field Adaptive Optics mode of the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) from 0.48 - 0.93 micron, are analysed here to determine the latitudinal and vertical distribution of cloud opacity and methane abundance in Neptune's observable troposphere (0.1 - ~3 bar). Previous observations at these wavelengths in 2003 by HST/STIS (Karkoschka and Tomasko 2011, Icarus 205, 674-694) found that the mole fraction of methane above the cloud tops (at ~2 bar) varied from ~4% at equatorial latitudes to ~2% at southern polar latitudes, by comparing the observed reflectivity at wavelengths near 825 nm controlled primarily by either methane absorption or H2-H2/H2-He collision-induced absorption. We find a similar variation in cloud-top methane abundance in 2018, which suggests that this depletion of methane towards Neptune's pole is potentially a long-lived feature, indicative of long-term upwelling at mid-equatorial latitudes and subsidence near the poles. By analysing these MUSE observations along the central meridian with a retrieval model, we demonstrate that a broad boundary between the nominal and depleted methane abundances occurs at between 20 - 40S. We also find a small depletion of methane near the equator, perhaps indicating subsidence there, and a local enhancement near 60 - 70S, which we suggest may be associated with South Polar Features (SPFs) seen in Neptune's atmosphere at these latitudes. Finally, by the use of both a reflectivity analysis and a principal component analysis, we demonstrate that this depletion of methane towards the pole is apparent at all locations on Neptune's disc, and not just along its central meridian.

A nova is a cataclysmic event on the surface of a white dwarf in a binary system that increases the overall brightness by several orders of magnitude. Although binary systems with a white dwarf are expected to be overabundant in globular clusters (GCs) compared to the Galaxy, only two novae from Galactic globular clusters have been observed. We present the discovery of an emission nebula in the Galactic globular cluster M 22 (NGC 6656) in observations made with the integral-field spectrograph MUSE. We extract the spectrum of the nebula and use the radial velocity determined from the emission lines to confirm that the nebula is part of NGC 6656. Emission-line ratios are used to determine the electron temperature and density. It is estimated to have a mass of 1 to $17 \times 10^{-5}$ solar masses. This mass and the emission-line ratios indicate that the nebula is a nova remnant. Its position coincides with the reported location of a 'guest star', an ancient Chinese term for transients, observed in May 48 BCE. With this discovery, this nova may be one of the oldest confirmed extrasolar events recorded in human history.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

We present the MUSE-Wide survey, a blind, 3D spectroscopic survey in the CANDELS/GOODS-S and CANDELS/COSMOS regions. Each MUSE-Wide pointing has a depth of 1 hour and hence targets more extreme and more luminous objects over 10 times the area of the MUSE-Deep fields (Bacon et al. 2017). The legacy value of MUSE-Wide lies in providing "spectroscopy of everything" without photometric pre-selection. We describe the data reduction, post-processing and PSF characterization of the first 44 CANDELS/GOODS-S MUSE-Wide pointings released with this publication. Using a 3D matched filtering approach we detected 1,602 emission line sources, including 479 Lyman-$\alpha$ (Lya) emitting galaxies with redshifts $2.9 \lesssim z \lesssim 6.3$. We cross-match the emission line sources to existing photometric catalogs, finding almost complete agreement in redshifts and stellar masses for our low redshift (z < 1.5) emitters. At high redshift, we only find ~55% matches to photometric catalogs. We encounter a higher outlier rate and a systematic offset of $\Delta$z$\simeq$0.2 when comparing our MUSE redshifts with photometric redshifts. Cross-matching the emission line sources with X-ray catalogs from the Chandra Deep Field South, we find 127 matches, including 10 objects with no prior spectroscopic identification. Stacking X-ray images centered on our Lya emitters yielded no signal; the Lya population is not dominated by even low luminosity AGN. A total of 9,205 photometrically selected objects from the CANDELS survey lie in the MUSE-Wide footprint, which we provide optimally extracted 1D spectra of. We are able to determine the spectroscopic redshift of 98% of 772 photometrically selected galaxies brighter than 24th F775W magnitude. All the data in the first data release - datacubes, catalogs, extracted spectra, maps - are available on the website https://musewide.aip.de. [abridged]

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent.

With a projected size of about 450 kpc at z~2.3, the Slug Ly-alpha nebula is a rare laboratory to study, in emission, the properties of the intergalactic gas in the Cosmic Web. Since its discovery, the Slug has been the subject of several spectroscopic follow-ups to constrain the properties of the emitting gas. Here we report the results of a deep MUSE integral-field spectroscopic search for non-resonant, extended HeII1640 and metal emission. Extended HeII radiation is detected on scales of about 100 kpc, but only in some regions associated with the bright Ly-alpha emission and a continuum-detected source, implying large and abrupt variations in the line ratios across adjacent regions in projected space. The recent detection of associated H-alpha emission and similar abrupt variations in the Ly-alpha kinematics, strongly suggest that the HeII/Ly-alpha gradient is due to large variations in the physical distances between the associated quasar and these regions. This implies that the overall length of the emitting structure could extend to physical Mpc scales and be mostly oriented along our line of sight. At the same time, the relatively low HeII/Ly-alpha values suggest that the emitting gas has a broad density distribution that - if expressed in terms of a lognormal - implies dispersions as high as those expected in the interstellar medium of galaxies. These results strengthen the possibility that the density distribution of intergalactic gas at high-redshift is extremely clumpy and multiphase on scales below our current observational spatial resolution of a few physical kpc.

(Abridged) We investigate the Lyman $\alpha$ emitter luminosity function (LAE LF) within the redshift range $2.9 \leq z \leq 6$ from the first instalment of the blind integral field spectroscopic survey MUSE-Wide. This initial part of the survey probes a region of 22.2 arcmin$^2$ in the CANDELS/GOODS-S field. The dataset provided us with 237 LAEs from which we construct the LAE LF in the luminosity range $42.2 \leq \log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}]\leq 43.5$ within a volume of $2.3\times10^5$ Mpc$^3$. For the LF construction we utilise three different non-parametric estimators: The classical $1/V_\mathrm{max}$ method, the $C^{-}$ method, and an improved binned estimator for the differential LF. All three methods deliver consistent results, with the cumulative LAE LF being $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 43.5) \simeq 3\times 10^{-6}$ Mpc$^{-3}$ and $\Phi(\log L_\mathrm{Ly\alpha} [\mathrm{erg\,s}^{-1}] = 42.2) \simeq 2 \times 10^{-3}$ Mpc$^{-3}$ towards the bright- and faint-end of our survey, respectively. By employing a non-parametric statistical test, as well as by comparing the full sample to sub-samples in redshift bins, we find no supporting evidence for an evolving LAE LF over the probed redshift and luminosity range. We determine the best-fitting Schechter function parameters $\alpha = -1.84^{+0.42}_{-0.41}$ and $\log L^* [\mathrm{erg\,s}^{-1}] = 42.2^{+0.22}_{-0.16}$ with the corresponding normalisation $\log \phi^* [\mathrm{Mpc}^{-3}] = -2.71$. When correcting for completeness in the LAE LF determinations, we take into account that LAEs exhibit diffuse extended low surface-brightness haloes. We compare the resulting LF to one obtained where we apply a correction assuming compact point-like emission. We find that the standard correction underestimates the LAE LF at the faint end of our survey by a factor of 2.5.

Galaxies are surrounded by large reservoirs of gas, mostly hydrogen, fed by inflows from the intergalactic medium and by outflows due to galactic winds. Absorption-line measurements along the sightlines to bright and rare background quasars indicate that this circumgalactic medium pervades far beyond the extent of starlight in galaxies, but very little is known about the spatial distribution of this gas. A new window into circumgalactic environments was recently opened with the discovery of ubiquitous extended Lyman-alpha emission from hydrogen around high-redshift galaxies, facilitated by the extraordinary sensitivity of the MUSE instrument at the ESO Very Large Telescope. Due to the faintness of this emission, such measurements were previously limited to especially favourable systems or to massive statistical averaging. Here we demonstrate that low surface brightness Lyman-alpha emission surrounding faint galaxies at redshifts between 3 and 6 adds up to a projected sky coverage of nearly 100%. The corresponding rate of incidence (the mean number of Lyman-alpha emitters penetrated by any arbitrary line of sight) is well above unity and similar to the incidence rate of high column density absorbers frequently detected in the spectra of distant quasars. This similarity suggests that most circumgalactic atomic hydrogen at these redshifts has now been detected also in emission.

Using an ultra-deep blind survey with the MUSE integral field spectrograph on the ESO Very Large Telescope, we obtain spectroscopic redshifts to a depth never explored before: galaxies with observed magnitudes m > 30 - 32. Specifically, we detect objects via Lyman-alpha emission at 2.9 < z < 6.7 without individual continuum counterparts in areas covered by the deepest optical/near-infrared imaging taken by the Hubble Space Telescope, the Hubble Ultra Deep Field. In total, we find 102 such objects in 9 square arcminutes at these redshifts. Detailed stacking analyses confirm the Lyman-alpha emission as well as the 1216 Angstrom-breaks and faint UV continua (M_UV ~ -15). This makes them the faintest spectroscopically-confirmed objects at these redshifts, similar to the sources believed to reionize the universe. A simple model for the expected fraction of detected/undetected Lyman-alpha emitters as a function of luminosity is consistent with these objects being the high-equivalent width tail of the normal Lyman-alpha-emitter population at these redshifts.

Star-forming galaxies have been found to follow a relatively tight relation between stellar mass ($M_{*}$) and star formation rate (SFR), dubbed the `star formation sequence'. A turnover in the sequence has been observed, where galaxies with $M_{*} < 10^{10} {\rm M}_{\odot}$ follow a steeper relation than their higher mass counterparts, suggesting that the low-mass slope is (nearly) linear. In this paper, we characterise the properties of the low-mass end of the star formation sequence between $7 \leq \log M_{*}[{\rm M}_{\odot}] \leq 10.5$ at redshift $0.11 < z < 0.91$. We use the deepest MUSE observations of the Hubble Ultra Deep Field and the Hubble Deep Field South to construct a sample of 179 star-forming galaxies with high signal-to-noise emission lines. Dust-corrected SFRs are determined from H$\beta$ $\lambda 4861$ and H$\alpha$ $\lambda 6563$. We model the star formation sequence with a Gaussian distribution around a hyperplane between $\log M_{*}$, $\log {\rm SFR}$, and $\log (1+z)$, to simultaneously constrain the slope, redshift evolution, and intrinsic scatter. We find a sub-linear slope for the low-mass regime where $\log {\rm SFR}[{\rm M}_{\odot}/{\rm yr}] = 0.83^{+0.07}_{-0.06} \log M_{*}[{\rm M}_{\odot}] + 1.74^{+0.66}_{-0.68} \log (1+z)$, increasing with redshift. We recover an intrinsic scatter in the relation of $\sigma_{\rm intr} = 0.44^{+0.05}_{-0.04}$ dex, larger than typically found at higher masses. As both hydrodynamical simulations and (semi-)analytical models typically favour a steeper slope in the low-mass regime, our results provide new constraints on the feedback processes which operate preferentially in low-mass halos.

Aims. As a new approach to the study of resolved stellar populations in nearby galaxies, our goal is to demonstrate in NGC300 that integral field spectroscopy with high spatial resolution and excellent seeing conditions reaches an unprecedented depth in severely crowded fields. Methods. MUSE observations with seven pointings in NGC300 have resulted in datacubes that are analyzed in four ways: (1) PSF-fitting 3D spectroscopy with PampelMUSE yields deblended spectra of individually distinguishable stars. The technique also provides samples of planetary nebulae that are complete down to m5007=28. (2) pseudo-monochromatic images, created at the wavelengths of the most important emission lines and corrected for continuum light by using the P3D visualization tool, provide maps of HII regions, SNR, and the diffuse ISM at a high level of sensitivity, allowing for the discovery of planetary nebulae, WR stars etc. (3) The use of the P3D line-fitting tool yields emission line fluxes, surface brightness, and kinematic information for gaseous objects, corrected for absorption line profiles of the underlying stellar population. (4) Visual inspection of the datacubes is demonstrated to be effcient for data mining and the discovery of background galaxies and unusual objects. Results. We present a catalogue of luminous stars, rare stars such as WR and other emission line stars, carbon stars, symbiotic star candidates, planetary nebulae, HII regions, supernova remnants, giant shells, peculiar diffuse and filamentary emission line objects, and background galaxies, along with their spectra. Conclusions. The technique of crowded-field 3D spectroscopy is capable of deblending individual bright stars, the unresolved background of faint stars, gaseous nebulae, and the diffuse component of the interstellar medium, resulting in unprecedented legacy value for observations of nearby galaxies with MUSE.

The physical origin of the near-ultraviolet MgII emission remains an under-explored domain, contrary to more typical emission lines detected in the spectra of star-forming galaxies. We explore the nebular and physical properties for a sample of 381 galaxies between 0.70 < z < 2.34 drawn from the MUSE Hubble Ultra Deep Survey. The spectra of these galaxies show a wide variety of profiles of the MgII 2796,2803 resonant doublet, from absorption to emission. We present a study on the main drivers for the detection of MgII emission in galaxy spectra. By exploiting photoionization models we verified that the emission-line ratios observed in galaxies with MgII in emission are consistent with nebular emission from HII regions. From a simultaneous analysis of MUSE spectra and ancillary HST information via spectral energy distribution (SED) fitting, we find that galaxies with MgII in emission have lower stellar masses, smaller sizes, bluer spectral slopes and lower optical depth than those with absorption. This leads us to suggest that MgII emission is a potential tracer of physical conditions not merely related to those of the ionized gas. We show that these differences in MgII emission/absorption can be explained in terms of a higher dust and neutral gas content in the interstellar medium (ISM) of galaxies showing MgII in absorption, confirming the extreme sensitivity of MgII to the presence of the neutral ISM. We conclude with an analogy between the MgII doublet and the Ly-alpha line, due to their resonant nature. Further investigations with current and future facilities, including JWST, are promising as the detection of MgII emission and its potential connection with Ly-alpha could provide new insights on the ISM content in the early Universe.

The Lyman alpha (lya) line of Hydrogen is a prominent feature in the spectra of star-forming galaxies, usually redshifted by a few hundreds of km/s compared to the systemic redshift. This large offset hampers follow-up surveys, galaxy pair statistics and correlations with quasar absorption lines when only lya is available. We propose diagnostics that can be used to recover the systemic redshift directly from the properties of the lya line profile. We use spectroscopic observations of Lyman-Alpha Emitters (LAEs) for which a precise measurement of the systemic redshift is available. Our sample contains 13 sources detected between z~3 and z~6 as part of various Multi Unit Spectroscopic Explorer (MUSE) Guaranteed Time Observations (GTO). We also include a compilation of spectroscopic lya data from the literature spanning a wide redshift range (z~0-8). First, restricting our analysis to double-peaked lya spectra, we find a tight correlation between the velocity offset of the red peak with respect to the systemic redshift, Vpeak, and the separation of the peaks. Secondly, we find a correlation between Vpeak and the full width at half maximum of the lya line. Fitting formulas, to estimate systemic redshifts of galaxies with an accuracy of +-100 km/s when only the lya emission line is available, are given for the two methods.

The "Antennae Galaxy" (NGC 4038/39) is the closest major interacting galaxy system and therefore often taken as merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star-formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect HII regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of HII regions, showing fainter luminosities. By comparing HII region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each HII region against ionization-parameter sensitive emission line ratios. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that just from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking HII regions for the diffuse ionized gas in the Antennae.

We report the discovery of a 10^4 kpc^2 gaseous structure detected in [OII] in an over-dense region of the COSMOS-Gr30 galaxy group at z~0.725 thanks to deep MUSE Guaranteed Time Observations. We estimate the total amount of diffuse ionised gas to be of the order of (~5+-3)x10^10 Msun and explore its physical properties to understand its origin and the source(s) of the ionisation. The MUSE data allow the identification of a dozen of group members embedded in this structure from emission and absorption lines. We extracted spectra from small apertures defined for both the diffuse ionised gas and the galaxies. We investigated the kinematics and ionisation properties of the various galaxies and extended gas regions thanks to line diagnostics (R23, O32 and [OIII]/H\beta) available within the MUSE wavelength range. We compared these diagnostics to photo-ionisation models and shock models. The structure is divided in two kinematically distinct sub-structures. The most extended sub-structure of ionised gas is likely rotating around a massive galaxy and displays filamentary patterns linking some galaxies. The second sub-structure links another massive galaxy hosting an Active Galactic Nucleus to a low mass galaxy but also extends orthogonally to the AGN host disk over ~35 kpc. This extent is likely ionised by the AGN itself. The location of small diffuse regions in the R23 vs. O32 diagram is compatible with photo-ionisation. However, the location of three of these regions in this diagram (low O32, high R23) can also be explained by shocks, which is supported by their large velocity dispersions. One edge-on galaxy shares the same properties and may be a source of shocks. Whatever the hypothesis, the extended gas seems to be non primordial. We favour a scenario where the gas has been extracted from galaxies by tidal forces and AGN triggered by interactions between at least the two sub-structures.

We present the deepest study to date of the Lya luminosity function (LF) in a blank field using blind integral field spectroscopy from MUSE. We constructed a sample of 604 Lya emitters (LAEs) across the redshift range 2.91 < z < 6.64 using automatic detection software in the Hubble Ultra Deep Field. We calculate accurate total Lya fluxes capturing low surface brightness extended Lya emission now known to be a generic property of high-redshift star-forming galaxies. We simulated realistic extended LAEs to characterise the selection function of our samples, and performed flux-recovery experiments to test and correct for bias in our determination of total Lya fluxes. We find an accurate completeness correction accounting for extended emission reveals a very steep faint-end slope of the LF, alpha, down to luminosities of log10 L erg s^-1< 41.5, applying both the 1/Vmax and maximum likelihood estimators. Splitting the sample into three broad redshift bins, we see the faint-end slope increasing from -2.03+1.42-inf at z ~ 3.44 to -2.86+0.76-inf at z ~ 5.48, however no strong evolution is seen between the 68% confidence regions in L*-alpha parameter space. Using the Lya line flux as a proxy for star formation activity, and integrating the observed LFs, we find that LAEs' contribution to the cosmic SFRD rises with redshift until it is comparable to that from continuum-selected samples by z ~ 6. This implies that LAEs may contribute more to the star-formation activity of the early Universe than previously thought - any additional interglactic medium correction would act to further boost the Lya luminosities. Finally, assuming fiducial values for the escape of Lya and LyC radiation, and the clumpiness of the IGM, we integrated the maximum likelihood LF at 5.00 < z < 6.64 and find we require only a small extrapolation beyond the data (< 1 dex in L) for LAEs alone to maintain an ionised IGM at z ~ 6.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We provide, for the first time, robust observational constraints on the galaxy major merger fraction up to $z\approx 6$ using spectroscopic close pair counts. Deep Multi Unit Spectroscopic Explorer (MUSE) observations in the Hubble Ultra Deep Field (HUDF) and Hubble Deep Field South (HDF-S) are used to identify 113 secure close pairs of galaxies among a parent sample of 1801 galaxies spread over a large redshift range ($0.2<z<6$) and stellar masses ($10^7-10^{11} M_\odot$), thus probing about 12 Gyr of galaxy evolution. Stellar masses are estimated from spectral energy distribution (SED) fitting over the extensive UV-to-NIR HST photometry available in these deep Hubble fields, adding Spitzer IRAC bands to better constrain masses for high-redshift ($z\geqslant 3$) galaxies. These stellar masses are used to isolate a sample of 54 major close pairs with a galaxy mass ratio limit of 1:6. Among this sample, 23 pairs are identified at high redshift ($z\geqslant 3$) through their Ly$\alpha$ emission. The sample of major close pairs is divided into five redshift intervals in order to probe the evolution of the merger fraction with cosmic time. Our estimates are in very good agreement with previous close pair counts with a constant increase of the merger fraction up to $z\approx 3$ where it reaches a maximum of 20%. At higher redshift, we show that the fraction slowly decreases down to about 10% at $z\approx6$. The sample is further divided into two ranges of stellar masses using either a constant separation limit of $10^{9.5} M_\odot$ or the median value of stellar mass computed in each redshift bin. Overall, the major close pair fraction for low-mass and massive galaxies follows the same trend. These new, homogeneous, and robust estimates of the major merger fraction since $z\approx6$ are in good agreement with recent predictions of cosmological numerical simulations.

We report the detection of extended Lyman-alpha (Lya) haloes around 145 individual star-forming galaxies at redshifts 3<z<6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer at ESO-VLT. Our sample consists of continuum-faint (-15> M_{UV}> -22) Lya emitters (LAEs). Using a 2D, two-component decomposition of Lya emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lya haloes. We find that 80% of our objects having reliable Lya halo measurements show Lya emission that is significantly more extended than the UV continuum detected by HST (by a factor ~4 to >20). The median exponential scale length of the Lya haloes in our sample is ~4.5 kpc. By comparing the maximal detected extent of the Lya emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lya emission of our selected sample of LAEs probes a significant portion of the cold circum-galactic medium (CGM) of these galaxies (>50% in average). This result shows that there must be significant HI reservoirs in the CGM and reinforces the idea that Lya haloes are ubiquitous around high-redshift Lya emitting galaxies. Our characterization of the Lya haloes indicates that the majority of the Lya flux comes from the halo (~65%) and that their scale lengths seem to be linked to the UV properties of the galaxies. We do not observe a significant Lya halo size evolution with redshift. We also find that the Lya lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km/s. While the FWHM does not seem to be correlated to the Lya scale length, most compact Lya haloes and those that are not detected with high significance tend to have narrower Lya profiles. Finally, we investigate the origin of the extended Lya emission but we conclude that our data do not allow us to disentangle the possible processes.

Non-resonant FeII* 2365, 2396, 2612, 2626 emission can potentially trace galactic winds in emission and provide useful constraints to wind models. From the 3'x3' mosaic of the Hubble Ultra Deep Field (UDF) obtained with the VLT/MUSE integral field spectrograph, we identify a statistical sample of 40 FeII* emitters and 50 MgII 2796, 2803 emitters from a sample of 271 [OII] 3726, 3729 emitters with reliable redshifts from z = 0.85 - 1.5 down to 2E-18 (3 sigma) ergs/s/cm^2 (for [OII]), covering the stellar mass range 10^8 - 10^11 Msun. The FeII* and MgII emitters follow the galaxy main sequence, but with a clear dichotomy. Galaxies with masses below 10^9 Msun and star formation rates (SFRs) of <1 Msun/year have MgII emission without accompanying FeII* emission, whereas galaxies with masses above 10^10 Msun and SFRs >10 Msun/year have FeII* emission without accompanying MgII emission. Between these two regimes, galaxies have both MgII and FeII* emission, typically with MgII P-Cygni profiles. Indeed, the MgII profile shows a progression along the main sequence from pure emission to P-Cygni profiles to strong absorption, due to resonant trapping. Combining the deep MUSE data with HST ancillary information, we find that galaxies with pure MgII emission profiles have lower star formation rate surface densities than those with either MgII P-Cygni profiles or FeII* emission. These spectral signatures produced through continuum scattering and fluorescence, MgII P-Cygni profiles and FeII* emission, are better candidates for tracing galactic outflows than pure MgII emission, which may originate from HII regions. We compare the absorption and emission rest-frame equivalent widths for pairs of FeII transitions to predictions from outflow models and find that the observations consistently have less total re-emission than absorption, suggesting either dust extinction or non-isotropic outflow geometries.

We present spatially resolved stellar kinematic maps, for the first time, for a sample of 17 intermediate redshift galaxies (0.2 < z < 0.8). We used deep MUSE/VLT integral field spectroscopic observations in the Hubble Deep Field South (HDFS) and Hubble Ultra Deep Field (HUDF), resulting from ~30h integration time per field, each covering 1'x1' field of view, with ~0.65" spatial resolution. We selected all galaxies brighter than 25mag in the I band and for which the stellar continuum is detected over an area that is at least two times larger than the spatial resolution. The resulting sample contains mostly late-type disk, main-sequence star-forming galaxies with 10^8.5 - 10^10.5 Msun. Using a full-spectrum fitting technique, we derive two-dimensional maps of the stellar and gas kinematics, including the radial velocity V and velocity dispersion sigma. We find that most galaxies in the sample are consistent with having rotating stellar disks with roughly constant velocity dispersions and that the Vrms=sqrt{V^2+sigma^2} of the gas and stars, a scaling proxy for the galaxy gravitational potential, compare well to each other. These spatially resolved observations of intermediate redshift galaxies suggest that the regular stellar kinematics of disk galaxies that is observed in the local Universe was already in place 4 - 7 Gyr ago and that their gas kinematics traces the gravitational potential of the galaxy, thus is not dominated by shocks and turbulent motions. Finally, we build dynamical axisymmetric Jeans models constrained by the derived stellar kinematics for two specific galaxies and derive their dynamical masses. These are in good agreement (within 25%) with those derived from simple exponential disk models based on the gas kinematics. The obtained mass-to-light ratios hint towards dark matter dominated systems within a few effective radii.

This is the first of a series of papers presenting the results from our survey of 25 Galactic globular clusters with the MUSE integral-field spectrograph. In combination with our dedicated algorithm for source deblending, MUSE provides unique multiplex capabilities in crowded stellar fields and allows us to acquire samples of up to 20 000 stars within the half-light radius of each cluster. The present paper focuses on the analysis of the internal dynamics of 22 out of the 25 clusters, using about 500 000 spectra of 200 000 individual stars. Thanks to the large stellar samples per cluster, we are able to perform a detailed analysis of the central rotation and dispersion fields using both radial profiles and two-dimensional maps. The velocity dispersion profiles we derive show a good general agreement with existing radial velocity studies but typically reach closer to the cluster centres. By comparison with proper motion data we derive or update the dynamical distance estimates to 14 clusters. Compared to previous dynamical distance estimates for 47 Tuc, our value is in much better agreement with other methods. We further find significant (>3sigma) rotation in the majority (13/22) of our clusters. Our analysis seems to confirm earlier findings of a link between rotation and the ellipticities of globular clusters. In addition, we find a correlation between the strengths of internal rotation and the relaxation times of the clusters, suggesting that the central rotation fields are relics of the cluster formation that are gradually dissipated via two-body relaxation.

The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer.

We tested the performance of photometric redshifts for galaxies in the Hubble Ultra Deep field down to 30th magnitude. We compared photometric redshift estimates from three spectral fitting codes from the literature (EAZY, BPZ and BEAGLE) to high quality redshifts for 1227 galaxies from the MUSE integral field spectrograph. All these codes can return photometric redshifts with bias |Dzn|=|z-z_phot|/(1+z)<0.05 down to F775W=30 and spectroscopic incompleteness is unlikely to strongly modify this statement. We have, however, identified clear systematic biases in the determination of photometric redshifts: in the 0.4<z<1.5 range, photometric redshifts are systematically biased low by as much as Dzn=-0.04 in the median, and at z>3 they are systematically biased high by up to Dzn = 0.05, an offset that can in part be explained by adjusting the amount of intergalactic absorption applied. In agreement with previous studies we find little difference in the performance of the different codes, but in contrast to those we find that adding extensive ground-based and IRAC photometry actually can worsen photo-z performance for faint galaxies. We find an outlier fraction, defined through |Dzn|>0.15, of 8% for BPZ and 10% for EAZY and BEAGLE, and show explicitly that this is a strong function of magnitude. While this outlier fraction is high relative to numbers presented in the literature for brighter galaxies, they are very comparable to literature results when the depth of the data is taken into account. Finally, we demonstrate that while a redshift might be of high confidence, the association of a spectrum to the photometric object can be very uncertain and lead to a contamination of a few percent in spectroscopic training samples that do not show up as catastrophic outliers, a problem that must be tackled in order to have sufficiently accurate photometric redshifts for future cosmological surveys.

We have conducted a two-layered spectroscopic survey (1'x1' ultra deep and 3'x3' deep regions) in the Hubble Ultra Deep Field (HUDF) with the Multi Unit Spectroscopic Explorer (MUSE). The combination of a large field of view, high sensitivity, and wide wavelength coverage provides an order of magnitude improvement in spectroscopically confirmed redshifts in the HUDF; i.e., 1206 secure spectroscopic redshifts for HST continuum selected objects, which corresponds to 15% of the total (7904). The redshift distribution extends well beyond z>3 and to HST/F775W magnitudes as faint as ~30 mag (AB, 1-sigma). In addition, 132 secure redshifts were obtained for sources with no HST counterparts that were discovered in the MUSE data cubes by a blind search for emission-line features. In total, we present 1338 high quality redshifts, which is a factor of eight increase compared with the previously known spectroscopic redshifts in the same field. We assessed redshifts mainly with the spectral features [OII] at z<1.5 (473 objects) and Lya at 2.9<z<6.7 (692 objects). With respect to F775W magnitude, a 50% completeness is reached at 26.5 mag for ultra deep and 25.5 mag for deep fields, and the completeness remains >~20% up to 28-29 mag and ~27 mag, respectively. We used the determined redshifts to test continuum color selection (dropout) diagrams of high-z galaxies. The selection condition for F336W dropouts successfully captures ~80% of the targeted z~2.7 galaxies. However, for higher redshift selections (F435W, F606W, and F775W dropouts), the success rates decrease to ~20-40%. We empirically redefine the selection boundaries to make an attempt to improve them to ~60%. The revised boundaries allow bluer colors that capture Lya emitters with high Lya equivalent widths falling in the broadbands used for the color-color selection. Along with this paper, we release the redshift and line flux catalog.

MUSE (Multi Unit Spectroscopic Explorer) is an integral-field spectrograph mounted on the Very Large Telescope (VLT) in Chile and made available to the European community since October 2014. The Centre de Recherche Astrophysique de Lyon has developed a dedicated software to help MUSE users analyze the reduced data. In this paper we introduce MPDAF, the MUSE Python Data Analysis Framework, based on several well-known Python libraries (Numpy, Scipy, Matplotlib, Astropy) which offers new tools to manipulate MUSE-specific data. We present different examples showing how this Python package may be useful for MUSE data analysis.

We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90\% of the entire HUDF region with a 10-hour deep exposure time, plus a deeper 31-hour exposure in a single 1.15 arcmin2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 arcsec at 7000 A) and accurate astrometry (0.07 arcsec rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to mAB=28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 sigma emission line detection limit for a point source is 1.5 and 3.1 10-19 erg.s-1.cm-2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts (mAB>31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 hours of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

Recent theoretical models suggest that the early phase of galaxy formation could involve an epoch when galaxies are gas-rich but inefficient at forming stars: a "dark galaxy" phase. Here, we report the results of our MUSE (Multi Unit Spectroscopic Explorer) survey for dark galaxies fluorescently illuminated by quasars at $z>3$. Compared to previous studies which are based on deep narrow-band (NB) imaging, our integral field survey provides a nearly uniform sensitivity coverage over a large volume in redshift space around the quasars as well as full spectral information at each location. Thanks to these unique features, we are able to build control samples at large redshift distances from the quasars using the same data taken under the same conditions. By comparing the rest-frame equivalent width (EW$_{0}$) distributions of the Ly$\alpha$ sources detected in proximity to the quasars and in control samples, we detect a clear correlation between the locations of high EW$_{0}$ objects and the quasars. This correlation is not seen in other properties such as Ly$\alpha$ luminosities or volume overdensities, suggesting the possible fluorescent nature of at least some of these objects. Among these, we find 6 sources without continuum counterparts and EW$_{0}$ limits larger than $240\,\mathrm{\AA}$ that are the best candidates for dark galaxies in our survey at $z>3.5$. The volume densities and properties, including inferred gas masses and star formation efficiencies, of these dark galaxy candidates are similar to previously detected candidates at $z\approx2.4$ in NB surveys. Moreover, if the most distant of these are fluorescently illuminated by the quasar, our results also provide a lower limit of $t=60$ Myr on the quasar lifetime.

We present an analysis of MUSE observations obtained on the massive Frontier Fields cluster Abell 2744. This new dataset covers the entire multiply-imaged region around the cluster core. We measure spectroscopic redshifts for HST-selected continuum sources together with line emitters blindly detected in the datacube. The combined catalog consists of 514 spectroscopic redshifts (with 414 new identifications), including 156 cluster members and 326 magnified background sources. We use this redshift information to perform a strong-lensing analysis of all multiple images previously found in the deep Frontier Field images, and add three new MUSE-detected multiply-imaged systems with no obvious HST counterpart. The combined strong lensing constraints include a total of 60 systems producing 188 images altogether, out of which 29 systems and 83 images are spectroscopically confirmed, making Abell 2744 one of the most well-constrained clusters to date. A parametric mass model including two cluster-scale components in the core and several group-scale substructures at larger radii accurately reproduces all the spectroscopic multiple systems, reaching an rms of 0.67" in the image plane. Overall, the large number of spectroscopic redshifts gives us a robust model and we estimate the systematics on the mass density and magnification within the cluster core to be typically ~9%.

We present a first instalment of the MUSE-Wide survey, covering an area of 22.2 arcmin$^2$ (corresponding to $\sim$20% of the final survey) in the CANDELS/Deep area of the Chandra Deep Field South. We use the MUSE integral field spectrograph at the ESO VLT to conduct a full-area spectroscopic mapping at a depth of 1h exposure time per 1 arcmin$^2$ pointing. We searched for compact emission line objects using our newly developed LSDCat software based on a 3-D matched filtering approach, followed by interactive classification and redshift measurement of the sources. Our catalogue contains 831 distinct emission line galaxies with redshifts ranging from 0.04 to 6. Roughly one third (237) of the emission line sources are Lyman $\alpha$ emitting galaxies with $3 < z < 6$, only four of which had previously measured spectroscopic redshifts. At lower redshifts 351 galaxies are detected primarily by their [OII] emission line ($0.3 \lesssim z \lesssim 1.5$), 189 by their [OIII] line ($0.21 \lesssim z \lesssim 0.85$), and 46 by their H$\alpha$ line ($0.04 \lesssim z \lesssim 0.42$). Comparing our spectroscopic redshifts to photometric redshift estimates from the literature, we find excellent agreement for $z<1.5$ with a median $\Delta z$ of only $\sim 4 \times 10^{-4}$ and an outlier rate of 6%, however a significant systematic offset of $\Delta z = 0.26$ and an outlier rate of 23% for Ly$\alpha$ emitters at $z>3$. Together with the catalogue we also release 1D PSF-weighted extracted spectra and small 3D datacubes centred on each of the 831 sources.

We present a clustering analysis of a sample of 238 Ly{$\alpha$}-emitters at redshift 3<z<6 from the MUSE-Wide survey. This survey mosaics extragalactic legacy fields with 1h MUSE pointings to detect statistically relevant samples of emission line galaxies. We analysed the first year observations from MUSE-Wide making use of the clustering signal in the line-of-sight direction. This method relies on comparing pair-counts at close redshifts for a fixed transverse distance and thus exploits the full potential of the redshift range covered by our sample. A clear clustering signal with a correlation length of r0 = 2.9(+1.0/-1.1) Mpc (comoving) is detected. Whilst this result is based on only about a quarter of the full survey size, it already shows the immense potential of MUSE for efficiently observing and studying the clustering of Ly{$\alpha$}-emitters.

Cosmological simulations suggest that most of the matter in the Universe is distributed along filaments connecting galaxies. Illuminated by the cosmic UV background (UVB), these structures are expected to glow in fluorescent Lyman alpha emission with a Surface Brightness (SB) that is well below current observational limits for individual detections. Here, we perform a stacking analysis of the deepest MUSE/VLT data using three-dimensional regions (subcubes) with orientations determined by the position of neighbouring Lyman alpha galaxies (LAEs) at 3<z<4. Our method should increase the probability of detecting filamentary Lyman alpha emission, provided that these structures are Lyman Limit Systems (LLSs). By stacking 390 oriented subcubes we reach a 2 sigma sensitivity level of SB ~ 0.44e-20 erg/s/cm^2/arcsec^2 in an aperture of 1 arcsec^2 x 6.25 Angstrom, which is three times below the expected fluorescent Lyman alpha signal from the Haardt-Madau 2012 (HM12) UVB at z~3.5. No detectable emission is found on intergalactic scales, implying that at least two thirds of our subcubes do not contain oriented LLSs for a HM12 UVB. On the other hand, significant emission is detected in the circum-galactic medium (CGM) of galaxies in the direction of the neighbours. The signal is stronger for galaxies with a larger number of neighbours and appears to be independent of any other galaxy properties such as luminosity, redshift and neighbour distance. We estimate that preferentially oriented satellite galaxies cannot contribute significantly to this signal, suggesting instead that gas densities in the CGM are typically larger in the direction of neighbouring galaxies on cosmological scales.

Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using deep observations (27 hours) of the Hubble Deep Field South from the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M* = 8 x 10^9 Msun, star formation rate SFR = 77 Msun/yr, and star formation rate surface brightness 1.6 Msun/kpc^2 within the [OII] half-light radius R_1/2,[OII] = 2.76 +- 0.17 kpc. From a component of the strong resonant MgII and FeII absorptions at -350 km/s, we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant FeII* emission, at lambda 2626, 2612, 2396, and 2365, at 1.2-2.4-1.5-2.7 x 10^-18 egs s-1 cm-2 respectively. These flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant FeII* emission lines, we spatially map the FeII* emission from an individual galaxy for the first time. The FeII* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius R_1/2,FeII* = 4.1 +- 0.4 kpc is 50% larger than the stellar continuum (R_1/2,* = 2.34 +- 0.17 kpc) or the [OII] nebular line. Moreover, the FeII* emission shows a blue wing extending up to -400 km/s, which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p-v diagram along that axis. These features are consistent with a bi-conical outflow.

We present a MUSE and KMOS dynamical study 405 star-forming galaxies at redshift z=0.28-1.65 (median redshift z=0.84). Our sample are representative of star-forming, main-sequence galaxies, with star-formation rates of SFR=0.1-30Mo/yr and stellar masses M=10^8-10^11Mo. For 49+/-4% of our sample, the dynamics suggest rotational support, 24+/-3% are unresolved systems and 5+/-2% appear to be early-stage major mergers with components on 8-30kpc scales. The remaining 22+/-5% appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination corrected rotational velocities and show these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j*=J/M*\propto M^(2/3) but with a redshift evolution that scales as j*\propto M^{2/3}(1+z)^(-1). We identify a correlation between specific angular momentum and disk stability such that galaxies with the highest specific angular momentum, log(j*/M^(2/3))>2.5, are the most stable, with Toomre Q=1.10+/-0.18, compared to Q=0.53+/-0.22 for galaxies with log(j*/M^(2/3))<2.5. At a fixed mass, the HST morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas disks: at z~1, massive galaxies that have disks with low specific angular momentum, appear to be globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular have evolved in to stable disks with spiral structures.

The Multi Unit Spectroscopic Explorer (MUSE) is an integral-field spectrograph operating in the visible wavelength range, and installed at the Very Large Telescope (VLT). The official MUSE pipeline is available from ESO. However, for the data reduction of the Deep Fields program (Bacon et al., in prep.), we have built a more sophisticated reduction pipeline, with additional reduction tasks, to extend the official pipeline and produce cubes with fewer instrumental residuals.

The physical properties of galactic winds are one of the keys to understand galaxy formation and evolution. These properties can be constrained thanks to background quasar lines of sight (LOS) passing near star-forming galaxies (SFGs). We present the first results of the MusE GAs FLOw and Wind (MEGAFLOW) survey obtained of 2 quasar fields which have 8 MgII absorbers of which 3 have rest-equivalent width greater than 0.8 \AA. With the new Multi Unit Spectroscopic Explorer (MUSE) spectrograph on the Very Large Telescope (VLT), we detect 6 (75$\%$) MgII host galaxy candidates withing a radius of 30 arcsec from the quasar LOS. Out of these 6 galaxy--quasar pairs, from geometrical arguments, one is likely probing galactic outflows, two are classified as "ambiguous", two are likely probing extended gaseous disks and one pair seems to be a merger. We focus on the wind$-$pair and constrain the outflow using a high resolution quasar spectra from Ultraviolet and Visual Echelle Spectrograph (UVES). Assuming the metal absorption to be due to gas flowing out of the detected galaxy through a cone along the minor axis, we find outflow velocities of the order of $\approx$ 150 km/s (i.e. smaller than the escape velocity) with a loading factor, $\eta =\dot M_{\rm out}/$SFR, of $\approx$ 0.7. We see evidence for an open conical flow, with a low-density inner core. In the future, MUSE will provide us with about 80 multiple galaxy$-$quasar pairs in two dozen fields.

(Abridged) We make use of the deepest VLT/MUSE observations performed so far on the Hubble Deep Field South (HDFS) to characterize the low-mass (< $10^{10}$M$_\odot$) galaxy population at intermediate redshift. We identify a sample of 28 spatially-resolved emission-line galaxies in the deep (27h integration time) MUSE data cube, spread over a redshift interval of 0.2 < z < 1.4. The public HST images and multi-band photometry over the HDFS are used to constrain the stellar mass and star formation rate (SFR) of the galaxies and to perform a morphological analysis. We derive the resolved ionized gas properties of these galaxies from the MUSE data and model the disk (both in 2D and with GalPaK$^{\rm 3D}$) to retrieve their intrinsic gas kinematics. We build a sample of resolved emission-line galaxies of much lower stellar mass and SFR (by $\sim$1-2 orders of magnitude) than previous 3D spectroscopic surveys. Most of the spatially-resolved MUSE-HDFS galaxies have gas kinematics consistent with disk-like rotation, but about 20% have velocity dispersions larger than the rotation velocities, and 30% are part of a close pair and/or show clear signs of recent gravitational interactions. In the high-mass regime, the MUSE-HDFS galaxies follow the Tully-Fisher relation defined from previous surveys in a similar redshift range. This scaling relation extends also to lower masses/velocities but with a higher dispersion. The MUSE-HDFS galaxies follow the scaling relations defined in the local universe between the specific angular momentum and the stellar mass. However, we find that intermediate-redshift star-forming galaxies fill a continuum transition from the spiral to elliptical local scaling relations, according to the dynamical state of the gas. This indicates that some galaxies may lose their angular momentum and become dispersion-dominated prior to becoming passive.

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

We present the first estimate of the Ly{\alpha} luminosity function using blind spectroscopy from the Multi Unit Spectroscopic Explorer, MUSE, in the Hubble Deep Field South. Using automatic source-detection software, we assemble a homogeneously-detected sample of 59 Ly{\alpha} emitters covering a flux range of -18.0 < log10 (F) < -16.3 (erg s^-1 cm^-2), corresponding to luminosities of 41.4 < log10 (L) < 42.8 (erg s^-1). As recent studies have shown, Ly{\alpha} fluxes can be underestimated by a factor of two or more via traditional methods, and so we undertake a careful assessment of each object's Ly{\alpha} flux using a curve-of-growth analysis to account for extended emission. We describe our self-consistent method for determining the completeness of the sample, and present an estimate of the global Ly{\alpha} luminosity function between redshifts 2.91 < z < 6.64 using the 1/Vmax estimator. We find the luminosity function is higher than many number densities reported in the literature by a factor of 2 - 3, although our result is consistent at the 1{\sigma} level with most of these studies. Our observed luminosity function is also in good agreement with predictions from semi-analytic models, and shows no evidence for strong evolution between the high- and low-redshift halves of the data. We demonstrate that one's approach to Ly{\alpha} flux estimation does alter the observed luminosity function, and caution that accurate flux assessments will be crucial in measurements of the faint end slope. This is a pilot study for the Ly{\alpha} luminosity function in the MUSE deep-fields, to be built on with data from the Hubble Ultra Deep Field which will increase the size of our sample by almost a factor of 10.

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

Direct Ly $\alpha$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $\alpha$ nebulae around 17 bright radio-quiet quasars at $3<z<4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant Ly $\alpha$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx-1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively "cold" (T$\sim$10$^4$K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

We use a background quasar to detect the presence of circum-galactic gas around a $z=0.91$ low-mass star forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the VLT show that the host galaxy has a dust-corrected star-formation rate (SFR) of 4.7$\pm$0.2 Msun/yr, with no companion down to 0.22 Msun/yr (5 $\sigma$) within 240 kpc (30"). Using a high-resolution spectrum (UVES) of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle $\alpha$ of only $15^\circ$), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold flow disk" extending at least 12 kpc ($3\times R_{1/2}$). We estimate the mass accretion rate $\dot M_{\rm in}$ to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the HI column density of $\log N_{\rm HI} \simeq 20.4$ obtained from a {\it HST}/COS NUV spectrum. From a detailed analysis of the low-ionization lines (e.g. ZnII, CrII, TiII, MnII, SiII), the accreting material appears to be enriched to about 0.4 $Z_\odot$ (albeit with large uncertainties: $\log Z/Z_\odot=-0.4~\pm~0.4$), which is comparable to the galaxy metallicity ($12+\log \rm O/H=8.7\pm0.2$), implying a large recycling fraction from past outflows. Blue-shifted MgII and FeII absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The MgII and FeII doublet ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent FeII* emission.

We present the results obtained with MUSE on the core of the lensing cluster A1689. Integral-field observations with MUSE provide a unique view of the central region, allowing us to conduct a complete census on both cluster galaxies and lensed background sources, identified based on their spectral features without preselection. We investigate the multiple-image configuration for all known sources in the field. Previous to our survey, 28 different lensed galaxies displaying 46 multiple images were known in the MUSE field of view, most of them based on photometric redshifts and lensing considerations. Among them, we spectroscopically confirm 12 images based on their emission-lines, corresponding to 7 different lensed galaxies between z = 0.95 and 5.0. In addition, 14 new galaxies have been spectroscopically identified in this area, with redshifts ranging between 0.8 and 6.2. All background sources within the MUSE field of view correspond to multiple-imaged systems lensed by A1689. 17 sources in total are found at z > 3 based on their Lyman-alpha emission, with Lyman-alpha luminosities ranging between 40.5 < log(Ly{\alpha}) < 42.5 after correction for magnification. This sample is particularly sensitive to the slope of the LF toward the faintest-end. The density of sources obtained in this survey is consistent with a steep value of {\alpha} < -1.5, although this result still needs further investigation. These results illustrate the efficiency of MUSE in the characterization of lensing clusters on one hand, and the study of faint and distant populations of galaxies on the other hand. In particular, our current survey of lensing clusters should provide a unique census of sources responsible for the reionization in a representative volume at z ~ 4-7.

We introduce Zurich Atmosphere Purge (ZAP), an approach to sky subtraction based on principal component analysis (PCA) that we have developed for the Multi Unit Spectrographic Explorer (MUSE) integral field spectrograph. ZAP employs filtering and data segmentation to enhance the inherent capabilities of PCA for sky subtraction. Extensive testing shows that ZAP reduces sky emission residuals while robustly preserving the flux and line shapes of astronomical sources. The method works in a variety of observational situations from sparse fields with a low density of sources to filled fields in which the target source fills the field of view. With the inclusion of both of these situations the method is generally applicable to many different science cases and should also be useful for other instrumentation. ZAP is available for download at http://muse-vlt.eu/science/tools.

Aims. We demonstrate the high multiplex advantage of crowded field 3D spectroscopy using the new integral field spectrograph MUSE by means of a spectroscopic analysis of more than 12,000 individual stars in the globular cluster NGC 6397. Methods. The stars are deblended with a PSF fitting technique, using a photometric reference catalogue from HST as prior, including relative positions and brightnesses. This catalogue is also used for a first analysis of the extracted spectra, followed by an automatic in-depth analysis using a full-spectrum fitting method based on a large grid of PHOENIX spectra. Results. With 18,932 spectra from 12,307 stars in NGC 6397 we have analysed the largest sample so far available for a single globular cluster. We derived a mean radial velocity of vrad=17.84+-0.07 km/s and a mean metallicity of [Fe/H]=-2.120+-0.002, with the latter seemingly varying with temperature for stars on the RGB. We determine T_eff and [Fe/H] from the spectra, and log g from HST photometry. This is the first very comprehensive HRD for a globular cluster based on the analysis of several thousands of stellar spectra, ranging from the main sequence to the tip of the RGB. Furthermore, two interesting objects were identified with one being a post-AGB star and the other a possible millisecond-pulsar companion.

We present a detailed analysis of the kinematics of the Galactic globular cluster NGC 6397 based on more than ~18,000 spectra obtained with the novel integral field spectrograph MUSE. While NGC 6397 is often considered a core collapse cluster, our analysis suggests a flattening of the surface brightness profile at the smallest radii. Although it is among the nearest globular clusters, the low velocity dispersion of NGC 6397 of <5km/s imposes heavy demands on the quality of the kinematical data. We show that despite its limited spectral resolution, MUSE reaches an accuracy of 1km/s in the analysis of stellar spectra. We find slight evidence for a rotational component in the cluster and the velocity dispersion profile that we obtain shows a mild central cusp. To investigate the nature of this feature, we calculate spherical Jeans models and compare these models to our kinematical data. This comparison shows that if a constant mass-to-light ratio is assumed, the addition of an intermediate-mass black hole with a mass of 600M_sun brings the model predictions into agreement with our data, and therefore could be at the origin of the velocity dispersion profile. We further investigate cases with varying mass-to-light ratios and find that a compact dark stellar component can also explain our observations. However, such a component would closely resemble the black hole from the constant mass-to-light ratio models as this component must be confined to the central ~5arcsec of the cluster and must have a similar mass. Independent constraints on the distribution of stellar remnants in the cluster or kinematic measurements at the highest possible spatial resolution should be able to distinguish the two alternatives.

Since 2014, MUSE, the Multi-Unit Spectroscopic Explorer, is in operation at the ESO-VLT. It combines a superb spatial sampling with a large wavelength coverage. By design, MUSE is an integral-field instrument, but its field-of-view and large multiplex make it a powerful tool for multi-object spectroscopy too. Every data-cube consists of 90,000 image-sliced spectra and 3700 monochromatic images. In autumn 2014, the observing programs with MUSE have commenced, with targets ranging from distant galaxies in the Hubble Deep Field to local stellar populations, star formation regions and globular clusters. This paper provides a brief summary of the key features of the MUSE instrument and its complex data reduction software. Some selected examples are given, how multi-object spectroscopy for hundreds of continuum and emission-line objects can be obtained in wide, deep and crowded fields with MUSE, without the classical need for any target pre-selection.

We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.

Spatially resolved studies of high redshift galaxies, an essential insight into galaxy formation processes, have been mostly limited to stacking or unusually bright objects. We present here the study of a typical (L$^{*}$, M$_\star$ = 6 $\times 10^9$ $M_\odot$) young lensed galaxy at $z=3.5$, observed with MUSE, for which we obtain 2D resolved spatial information of Ly$\alpha$ and, for the first time, of CIII] emission. The exceptional signal-to-noise of the data reveals UV emission and absorption lines rarely seen at these redshifts, allowing us to derive important physical properties (T$_e\sim$15600 K, n$_e\sim$300 cm$^{-3}$, covering fraction f$_c\sim0.4$) using multiple diagnostics. Inferred stellar and gas-phase metallicities point towards a low metallicity object (Z$_{\mathrm{stellar}}$ = $\sim$ 0.07 Z$_\odot$ and Z$_{\mathrm{ISM}}$ $<$ 0.16 Z$_\odot$). The Ly$\alpha$ emission extends over $\sim$10 kpc across the galaxy and presents a very uniform spectral profile, showing only a small velocity shift which is unrelated to the intrinsic kinematics of the nebular emission. The Ly$\alpha$ extension is $\sim$4 times larger than the continuum emission, and makes this object comparable to low-mass LAEs at low redshift, and more compact than the Lyman-break galaxies and Ly$\alpha$ emitters usually studied at high redshift. We model the Ly$\alpha$ line and surface brightness profile using a radiative transfer code in an expanding gas shell, finding that this model provides a good description of both observables.

We present a pilot study of the z=2.923 radio galaxy MRC0943-242, where we for the first time combine information from ALMA and MUSE data cubes. Even with modest integration times, we disentangle an AGN and a starburst dominated set of components. These data reveal a highly complex morphology, as the AGN, starburst, and molecular gas components show up as widely separated sources in dust continuum, optical continuum and CO line emission observations. CO(1-0) and CO(8-7) line emission suggest that there is a molecular gas reservoir offset from both the dust and the optical continuum that is located ~90kpc from the AGN. The UV line emission has a complex structure in emission and absorption. The line emission is mostly due to i) a large scale ionisation cone energised by the AGN, ii) a Ly-alpha emitting bridge of gas between the radio galaxy and a heavily star-forming set of components. Strangely, the ionisation cone has no Ly-alpha emission. We find this is due to an optically thick layer of neutral gas with unity covering fraction spread out over a region of at least ~100kpc from the AGN. Other, less thick absorption components are associated with Ly-alpha emitting gas within a few tens of kpc from the radio galaxy and are connected by a bridge of emission. We speculate that this linear structure of dust, Ly-alpha and CO emission, and the redshifted absorption seen in the circum-nuclear region may represent an accretion flow feeding gas into this massive AGN host galaxy.

[abridged] The polar disk galaxy NGC4650A was observed during the commissioning of the MUSE at the ESO VLT to obtain the first 2D map of the velocity and velocity dispersion for both stars and gas. The new MUSE data allow the analysis of the structure and kinematics towards the central regions of NGC4650A, where the two components co-exist. These regions were unexplored by the previous long-slit literature data available for this galaxy. The extended view of NGC~4650A given by the MUSE data is a galaxy made of two perpendicular disks that remain distinct and drive the kinematics right into the very centre of this object. In order to match this observed structure for NGC4650A, we constructed a multicomponent mass model made by the combined projection of two disks. By comparing the observations with the 2D kinematics derived from the model, we found that the modelled mass distribution in these two disks can, on average, account for the complex kinematics revealed by the MUSE data, also in the central regions of the galaxy where the two components coexist. This result is a strong constraint on the dynamics and formation history of this galaxy; it further supports the idea that polar disk galaxies like NGC~4650A were formed through the accretion of material that has different angular momentum.

We present a new integral-field spectroscopic dataset of the central part of the Orion Nebula (M 42), observed with the MUSE instrument at the ESO VLT. We reduced the data with the public MUSE pipeline. The output products are two FITS cubes with a spatial size of ~5.9'x4.9' (corresponding to ~0.76 pc x 0.63 pc) and a contiguous wavelength coverage of 4595...9366 Angstrom, spatially sampled at 0.2". We provide two versions with a sampling of 1.25 Angstrom and 0.85 Angstrom in dispersion direction. Together with variance cubes these files have a size of 75 and 110 GiB on disk. They represent one of the largest integral field mosaics to date in terms of information content. We make them available for use in the community. To validate this dataset, we compare world coordinates, reconstructed magnitudes, velocities, and absolute and relative emission line fluxes to the literature and find excellent agreement. We derive a two-dimensional map of extinction and present de-reddened flux maps of several individual emission lines and of diagnostic line ratios. We estimate physical properties of the Orion Nebula, using the emission line ratios [N II] and [S III] (for the electron temperature $T_e$) and [S II] and [Cl III] (for the electron density $N_e$), and show two-dimensional images of the velocity measured from several bright emission lines.

MUSE, a giant integral field spectrograph, is about to become the newest facility instrument at the VLT. It will see first light in February 2014. Here, we summarize the properties of the instrument as built and outline functionality of the data reduction system, that transforms the raw data that gets recorded separately in 24 IFUs by 4k CCDs, into a fully calibrated, scientifically usable data cube. We then describe recent work regarding geometrical calibration of the instrument and testing of the processing pipeline, before concluding with results of the Preliminary Acceptance in Europe and an outlook to the on-sky commissioning.

MUSE observations of NGC5813 reveal a complex structure in the velocity dispersion map, previously hinted by SAURON observations. The structure is reminiscent of velocity dispersion maps of galaxies comprising two counter-rotating discs, and may explain the existence of the kinematically distinct core (KDC). Further evidence for two counter-rotating components comes from the analysis of the higher moments of the stellar line-of-sight velocity distributions and fitting MUSE spectra with two separate Gaussian line-of-sight velocity distributions. The emission-line kinematics show evidence of being linked to the present cooling flows and the buoyant cavities seen in X-rays. We detect ionised gas in a nuclear disc-like structure, oriented like the KDC, which is, however, not directly related to the KDC. We build an axisymmetric Schwarzschild dynamical model, which shows that the MUSE kinematics can be reproduced well with two counter-rotating orbit families, characterised by relatively low angular momentum components, but clearly separated in integral phase space and with radially varying contributions. The model indicates that the counter-rotating components in NGC5813 are not thin discs, but dynamically hot structures. Our findings give further evidence that KDCs in massive galaxies should not necessarily be considered as structurally or dynamically decoupled regions, but as the outcomes of the mixing of different orbital families, where the balance in the distribution of mass of the orbital families is crucial. We discuss the formation of the KDC in NGC5813 within the framework of gas accretion, binary mergers and formation of turbulent thick discs from cold streams at high redshift.

We present deep MUSE integral-field unit (IFU) spectroscopic observations of the giant (~150 x 80 kpc) Ly-alpha halo around the z=4.1 radio galaxy TNJ J1338-1942. This 9-hr observation maps the two-dimensional kinematics of the Ly-alpha emission across the halo. We identify two HI absorbers which are seen against the Ly-alpha emission, both of which cover the full 150 x 80 kpc extent of the halo and so have covering fractions ~1. The stronger and more blue-shifted absorber (dv~1200 km/s) has dynamics that mirror that of the underlying halo emission and we suggest that this high column material (n(HI) ~ 10^19.4 /cm^2), which is also seen in CIV absorption, represents an out-flowing shell that has been driven by the AGN (or star formation) within the galaxy. The weaker (n(HI)~10^14 /cm^2) and less blue shifted (dv~500 km/s) absorber most likely represents material in the cavity between the out-flowing shell and the Ly-alpha halo. We estimate that the mass in the shell must be of order 10^10 Msol -- a significant fraction of the ISM from a galaxy at z=4. The large scales of these coherent structures illustrate the potentially powerful influence of AGN feedback on the distribution and energetics of material in their surroundings. Indeed, the discovery of high-velocity (~1000 km/s), group-halo-scale (i.e. >150 kpc) and mass-loaded winds in the vicinity of the central radio source are broadly in agreement with the requirements of models that invoke AGN-driven outflows to regulate star formation and black-hole growth in massive galaxies at early times.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

We present new observations of the lensing cluster SMACSJ2031.8-4036 obtained with the MUSE integral field spectrograph as part of its commissioning on the Very Large Telescope. By providing medium-resolution spectroscopy over the full 4750-9350 Angstroms domain and a 1x1 arcmin2 field of view, MUSE is ideally suited for identifying lensed galaxies in the cluster core, in particular multiple-imaged systems. We perform a redshift analysis of all sources in the datacube and identify a total of 12 systems ranging from $z=1.46$ to $z=6.4$, with all images of each system confirmed by a spectroscopic redshift. This allows us to accurately constrain the cluster mass profile in this region. We foresee that future MUSE observations of cluster cores should help us discover very faint Lyman-alpha emitters thanks to the strong magnification and the high sensitivity of this instrument.

The next generation of wide-field spectroscopic redshift surveys will map the large-scale galaxy distribution in the redshift range 0.7< z<2 to measure baryonic acoustic oscillations (BAO). The primary optical signature used in this redshift range comes from the [OII] emission line doublet, which provides a unique redshift identification that can minimize confusion with other single emission lines. To derive the required spectrograph resolution for these redshift surveys, we simulate observations of the [OII] (3727,3729) doublet for various instrument resolutions, and line velocities. We foresee two strategies about the choice of the resolution for future spectrographs for BAO surveys. For bright [OII] emitter surveys ([OII] flux ~30.10^{-17} erg /cm2/s like SDSS-IV/eBOSS), a resolution of R~3300 allows the separation of 90 percent of the doublets. The impact of the sky lines on the completeness in redshift is less than 6 percent. For faint [OII] emitter surveys ([OII] flux ~10.10^{-17} erg /cm2/s like DESi), the detection improves continuously with resolution, so we recommend the highest possible resolution, the limit being given by the number of pixels (4k by 4k) on the detector and the number of spectroscopic channels (2 or 3).

Following a first study of the central regions of M32 that illustrated the power of integral-field spectroscopy (IFS) in detecting and measuring the [O III]{\lambda}5007 emission of PNe against a strong stellar background, we turn to the very nuclear PN population of M31, within 80 pc of its centre. We show that PNe can also be found in the presence of emission from diffuse gas and further illustrate the excellent sensitivity of IFS in detecting extragalactic PNe through a comparison with narrowband images obtained with the Hubble Space Telescope. We find the nuclear PNe population of M31 is only marginally consistent with the generally adopted form of the PNe luminosity function (PNLF). In particular, this is due to a lack of PNe with absolute magnitude M5007 brighter than -3, which would only result from a rather unfortunate draw from such a model PNLF. We suggest that the observed lack of bright PNe in the nuclear regions of M31 is due to a horizontal-branch population that is more tilted toward less massive and hotter He-burning stars, so that its progeny consists mostly of UV-bright stars that fail to climb back up the asymptotic giant branch (AGB) and only of few, if any, bright PNe powered by central post-AGB stars. These results are also consistent with recent reports on a dearth of bright post-AGB stars towards the nucleus of M31, and lend further support to the idea that the metallicity of a stellar population has an impact on the way the horizontal branch is populated and to the loose anticorrelation between the strength of the UV-upturn and the specific number of PNe that is observed in early-type galaxies. Finally, our investigation also serves to stress the importance of considering the same spatial scales when comparing the PNe population of galaxies with the properties of their stellar populations.

MUSE, the Multi Unit Spectroscopic Explorer, is a 2nd generation integral-field spectrograph under final assembly to see first light at the Very Large Telescope in 2013. By capturing ~ 90000 optical spectra in a single exposure, MUSE represents a challenge for data reduction and analysis. We summarise here the main features of the Data Reduction System, as well as some of the tools under development by the MUSE consortium and the DAHLIA team to handle the large MUSE datacubes (about 4x?10^8 pixels) to recover the original astrophysical signal.

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population.

The unexpected rising flux of early-type galaxies at decreasing ultraviolet (UV) wavelengths is a long-standing mystery. One important observational constraint is the correlation between UV-optical colours and Mg2 line strengths found by Burstein et al. (1988). The simplest interpretation of this phenomenon is that the UV strength is related to the Mg line strength. Under this assumption, we expect galaxies with larger Mg gradients to have larger UV colour gradients. By combining UV imaging from GALEX, optical imaging from MDM and SAURON integral-field spectroscopy, we investigate the spatially-resolved relationships between UV colours and stellar population properties of 34 early-type galaxies from the SAURON survey sample. We find that galaxies with old stellar populations show tight correlations between the FUV colours (FUV-V and FUV-NUV) and the Mgb index, H{\beta} index and metallicity [Z/H]. The equivalent correlations for the Fe5015 index, {\alpha}-enhancement [{\alpha}/Fe] and age are present but weaker. We have also derived logarithmic internal radial colour, measured line strength and derived stellar population gradients for each galaxy and again found a strong dependence of the FUV-V and FUV-NUV colour gradients on both the Mg b line strength and the metallicity gradients for galaxies with old stellar populations. In particular, global gradients of Mg b and [Z/H] with respect to the UV colour across galaxies are consistent with their local gradients within galaxies, suggesting that the global correlations also hold locally. From a simple model based on multi-band colour fits of UV upturn and UV-weak galaxies, we have identified a plausible range of parameters that reproduces the observed radial colour profiles. In these models, the centers of elliptical galaxies, where the UV flux is strong, are enhanced in metals by roughly 60% compared to UV-weak regions.

We investigate the [3.6] - [4.5] Spitzer-IRAC colour behaviour of the early-type galaxies of the SAURON survey, a representative sample of 48 nearby ellipticals and lenticulars. We investigate how this colour, which is unaffected by dust extinction, can be used to constrain the stellar populations in these galaxies. We find a tight relation between the [3.6]-[4.5] colour and effective velocity dispersion, a good mass-indicator in early-type galaxies. Contrary to other colours in the optical and near-infrared, we find that the colours become bluer for larger galaxies. The relations are tighter when using the colour inside r_e, rather than the much smaller r_e/8 aperture, due to the presence of young populations in the central regions. We also obtain strong correlations between the [3.6]-[4.5] colour and 3 strong absorption lines (H beta, Mg b and Fe 5015). Comparing our data with the models of Marigo et al., which show that more metal rich galaxies are bluer, we can explain our results in a way consistent with results from the optical, by stating that larger galaxies are more metal rich. The blueing is caused by a strong CO absorption band, whose line strength increases strongly with decreasing temperature and which covers a considerable fraction of the 4.5 micron filter. In galaxies that contain a compact radio source, the [3.6]-[4.5] colour is generally slightly redder than in the other galaxies, indicating small amounts of either hot dust, non-thermal emission, or young stars near the center. We find that the large majority of the galaxies show redder colours with increasing radius. abbreviated...

We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical $\Lambda$CDM scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via $\lambda_R$. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the \atl survey, assuming that fast rotating ETGs have at least 10% of their total stellar mass in a disc component. In agreement with \atl observations we find that slow rotators are predominantly galaxies with $ M_* > 10^{10.5}$ M$_{\odot}$ contributing $\sim 20%$ to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between $\sim 50% -90%$ of their stellar mass from satellites and their most massive progenitors have on average up to 3 major mergers during their evolution. Fast rotators in contrast, accrete less than 50% and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stops, galaxies grow via infall from satellites. Frequent minor mergers thereby, destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering $\lambda_R$ into the slow rotator regime. Abridged...

[Abridged] We present ground-based MDM V-band and Spitzer/IRAC 3.6um-band photometric observations of the 72 representative galaxies of the SAURON Survey. In combination with the SAURON stellar velocity dispersion measured within an effective radius (se), this allows us to explore the location of our galaxies in the main scaling relations. We investigate the dependence of these relations on our recent kinematical classification of early-type galaxies (i.e. Slow/Fast Rotators) and the stellar populations. Slow Rotator and Fast Rotator E/S0 galaxies do not populate distinct locations in the scaling relations, although Slow Rotators display a smaller intrinsic scatter. Surprisingly, extremely young objects do not display the bluest (V-[3.6]) colours in our sample, as is usually the case in optical colours. This can be understood in the context of the large contribution of TP-AGB stars to the infrared, even for young populations, resulting in a very tight (V-[3.6]) - se relation that in turn allows us to define a strong correlation between metallicity and velocity dispersion. Many Sa galaxies appear to follow the Fundamental Plane defined by E/S0 galaxies. Galaxies that appear offset from the relations correspond mostly to objects with extremely young populations, with signs of on-going, extended star formation. We correct for this effect in the Fundamental Plane, by replacing luminosity with stellar mass using an estimate of the stellar mass-to-light ratio, so that all galaxies are part of a tight, single relation. The new estimated coefficients are consistent in both photometric bands and suggest that differences in stellar populations account for about half of the observed tilt with respect to the virial prediction. After these corrections, the Slow Rotator family shows almost no intrinsic scatter around the best-fit Fundamental Plane.

Extragalactic Planetary Nebulae (PNe) are not only useful as distance signposts or as tracers of the dark-matter content of their host galaxies, but constitute also good indicators of the main properties of their parent stellar populations. Yet, so far, the properties of PNe in the optical regions of galaxies where stellar population gradients can be more extreme have remained largely unexplored, mainly because the detection of PNe with narrow-band imaging or slit-less spectroscopy is considerably hampered by a strong stellar background. Integral-field spectroscopy (IFS) can overcome this limitation, and here we present a study of the PN population in the nearby compact elliptical M32. Using SAURON data taken with just two 10-minutes-long pointings we have doubled the number of known PNe within the effective radius of M32, detecting PNe five times fainter than previously found in narrow-band images that collected nearly the same number of photons. Furthermore, by carefully accounting for the incompleteness of our survey we could conclude, despite having only 15 sources, that the central PNe population of M32 is consistent with the generally adopted shape for the PNe Luminosity Function and its typical normalization observed in early-type galaxies. Finally, owing to the proximity of M32 and to UV images taken with HST, we could identify the most likely candidates for the central star of a subset of our detected PNe and conclude that these stars are affected by substantial amounts of circumstellar dust extinction, a finding that could reconcile the intriguing discrepancy previously reported in M32 between model predictions and observations for the later stages of stellar evolution. Considering the modest time investment on a 4m-class telescope that delivered these results, this work illustrates the potential of future IFS studies for the central PNe population of early-type galaxies.

Using far (FUV) and near (NUV) ultraviolet photometry from guest investigator programmes on the Galaxy Evolution Explorer (GALEX) satellite, optical photometry from the MDM Observatory and optical integral-field spectroscopy from SAURON, we explore the UV-linestrength relations of the 48 nearby early-type galaxies in the SAURON sample. Identical apertures are used for all quantities, avoiding aperture mismatch. We show that galaxies with purely old stellar populations show well-defined correlations of the integrated FUV-V and FUV-NUV colours with the integrated Mgb and Hbeta absorption linestrength indices, strongest for FUV-NUV. Correlations with the NUV-V colour, Fe5015 index and stellar velocity dispersion are much weaker. These correlations put stringent constraints on the origin of the UV-upturn phenomenon in early-type galaxies, and highlight its dependence on age and metallicity. In particular, despite recent debate, we recover the negative correlation between FUV-V colour and Mg linestrength originally publicised by Burstein et al. (1988), which we refer to as the "Burstein relation", suggesting a positive dependence of the UV-upturn on metallicity. We argue that the scatter in the correlations is real, and present mild evidence that a strong UV excess is preferentially present in slow-rotating galaxies. We also demonstrate that most outliers in the correlations are galaxies with current or recent star formation, some at very low levels. We believe that this sensitivity to weak star formation, afforded by the deep and varied data available for the SAURON sample, explains why our results are occasionally at odds with other recent but shallower surveys. This is supported by the analysis of a large, carefully-crafted sample of more distant early-type galaxies from the Sloan Digital Sky Survey (SDSS), more easily comparable with current and future large surveys.

The Atlas3D project is a multi-wavelength survey combined with a theoretical modeling effort. The observations span from the radio to the millimeter and optical, and provide multi-colour imaging, two-dimensional kinematics of the atomic (HI), molecular (CO) and ionized gas (Hbeta, [OIII] and [NI]), together with the kinematics and population of the stars (Hbeta, Fe5015 and Mgb), for a carefully selected, volume-limited (1.16*10^5 Mpc^3) sample of 260 early-type (elliptical E and lenticular S0) galaxies (ETGs). The models include semi-analytic, N-body binary mergers and cosmological simulations of galaxy formation. Here we present the science goals for the project and introduce the galaxy sample and the selection criteria. The sample consists of nearby (D<42 Mpc) morphologically-selected ETGs extracted from a parent sample of 871 galaxies (8% E, 22% S0 and 70% spirals) brighter than M_K<-21.5 mag (stellar mass M_Star>6*10^9 M_Sun). We analyze possible selection biases and we conclude that the parent sample is essentially complete and statistically representative of the nearby galaxy population. We present the size-luminosity relation for the spirals and ETGs and show that the ETGs in the Atlas3D sample define a tight red sequence in a colour-magnitude diagram, with few objects in the transition from the blue cloud. We describe the strategy of the SAURON integral-field observations and the extraction of the stellar kinematics with the pPXF method. We give an overview of the characteristics of the other main datasets already available for our sample and of the ongoing modelling projects.

(Abridged) We present a stellar population analysis of the absorption line strength maps for 48 early-type galaxies from the SAURON sample. Using the line strength index maps of Hbeta, Fe5015, and Mgb, measured in the Lick/IDS system and spatially binned to a constant signal-to-noise, together with predictions from up-to-date stellar population models, we estimate the simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio [alpha/Fe] over a two-dimensional field extending up to approximately one effective radius. We find a large range of SSP-equivalent ages in our sample, of which ~40% of the galaxies show signs of a contribution from a young stellar population. The most extreme cases of post-starburst galaxies, with SSP-equivalent ages of <=3 Gyr observed over the full field-of-view, and sometimes even showing signs of residual star-formation, are restricted to low mass systems(sigma_e <= 100 k/ms or ~2x10^10 M_sol). Spatially restricted cases of young stellar populations in circumnuclear regions can almost exclusively be linked to the presence of star-formation in a thin, dusty disk/ring, also seen in the near-UV or mid-IR. The flattened components with disk-like kinematics previously identified in all fast rotators (Krajnovi\'c et al.) are shown to be connected to regions of distinct stellar populations. These range from the young, still star-forming circumnuclear disks and rings with increased metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures with extended disk-like kinematics, which are observed to have an increased metallicity and mildly depressed [alpha/Fe] ratio compared to the main body of the galaxy. The slow rotators generally show no stellar population signatures over and above the well known metallicity gradients and are largely consistent with old (>=10 Gyr) stellar populations.

We investigate the accuracy of mass determinations M_BH of supermassive black holes in galaxies using dynamical models of the stellar kinematics. We compare 10 of our M_BH measurements, using integral-field OASIS kinematics, to published values. For a sample of 25 galaxies we confront our new M_BH derived using two modeling methods on the same OASIS data.

We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/IRAC imaging in order to investigate the presence of trace star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by PAH emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disk galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disk or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, co-rotating inner stellar disks ubiquitous in this population.

Following our study on the incidence, morphology and kinematics of the ionised gas in early-type galaxies we now address the question of what is powering the observed nebular emission. To constrain the likely sources of gas excitation, we resort to a variety of ancillary data, draw from complementary information on the gas kinematics, stellar populations and galactic potential from the SAURON data, and use the SAURON-specific diagnostic diagram juxtaposing the [OIII]/Hb and [NI]/Hb line ratios. We find a tight correlation between the stellar surface brightness and the flux of the Hb recombination line across our sample, which points to a diffuse and old stellar source as the main contributor of ionising photons in early-type galaxies, with post-asymptotic giant branch (pAGB) stars being still the best candidate based on ionising-balance arguments. Other ionising sources such as a central AGN, OB-stars, shocks or the interaction between the hot and warm phases of the interestellar medium are found to play only a limited or localised role in powering the diffuse nebular emission observed in our sample galaxies. These results lead us to investigate the relative importance of stellar and AGN photoionisation in explaining the ionised-gas emission observed in early-type galaxies by the Sloan Digital Sky Survey (SDSS). By simulating how our sample galaxies would appear if placed at further distance and targeted by the SDSS, we conclude that only in very few, if any, of the SDSS early-type galaxies that display modest values for the equivalent width of the [OIII] line (less than ~2.4\AA) and LINER-like [OIII]/Hb values, the nebular emission is truly powered by an AGN.

We present GALEX far (FUV) and near (NUV) ultraviolet imaging of 34 nearby early-type galaxies from the SAURON representative sample of 48 E/S0 galaxies, all of which have ground-based optical imaging from the MDM Observatory. The surface brightness profiles of nine galaxies show regions with blue UV-optical colours suggesting recent star formation. Five of these show blue integrated UV-optical colours that set them aside in the NUV integrated colour-magnitude relation. They also have other properties confirming they have had recent star formation, in particular H_beta absorption higher than expected for a quiescent population and a higher CO detection rate. NUV-blue galaxies are generally drawn from the lower stellar velocity dispersion and thus lower dynamical mass part of the sample. We have also constructed the first UV Fundamental Planes and show that NUV blue galaxies bias the slopes and increase the scatters. If they are eliminated the fits get closer to expectations from the virial theorem. Although our analysis is based on a limited sample, it seems that a dominant fraction of the tilt and scatter of the UV Fundamental Planes is due to the presence of young stars in preferentially low-mass early-type galaxies. Interestingly, the UV-optical radial colour profiles reveal a variety of behaviours, with many galaxies showing signs of recent star formation, a central UV-upturn phenomenon, smooth but large-scale age and metallicity gradients, and in many cases a combination of these. In addition, FUV-NUV and FUV-V colours even bluer than those normally associated with UV-upturn galaxies are observed at the centre of some quiescent galaxies. Four out of the five UV-upturn galaxies are slow rotators. These objects should thus pose interesting challenges to stellar evolutionary models of the UV-upturn.

We present the results of an investigation of the local escape velocity (Vesc) - line strength index relationship for 48 early type galaxies from the SAURON sample, the first such study based on a large sample of galaxies with both detailed integral field observations and extensive dynamical modelling. Values of Vesc are computed using Multi Gaussian Expansion (MGE) photometric fitting and axisymmetric, anisotropic Jeans' dynamical modelling simultaneously on HST and ground-based images. We determine line strengths and escape velocities at multiple radii within each galaxy, allowing an investigation of the correlation within individual galaxies as well as amongst galaxies. We find a tight correlation between Vesc and the line-strength indices. For Mgb we find that this correlation exists not only between different galaxies but also inside individual galaxies - it is both a local and global correlation. The relation within individual galaxies has the same slope and offset as the global relation to a good level of agreement, though there is significant intrinsic scatter in the local gradients. We transform our line strength index measurements to the single stellar population (SSP) equivalent age (t), metallicity ([Z/H]) and enhancement ([$\alpha$/Fe]) and carry out a principal component analysis of our SSP and Vesc data. We find that in this four-dimensional parameter space the galaxies in our sample are to a good approximation confined to a plane, given by $\log \mathrm({V}_{\mathrm{esc}}/500\mathrm{km/s}) = 0.85 \mathrm{[Z/H]} + 0.43 \log (\mathrm{t}/\mathrm{Gyrs})$ - 0.20. It is surprising that it is a combination of age and metallicity that is conserved; this may indicate a 'conspiracy' between age and metallicity or a weakness in the SSP models.

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight velocity distribution and absorption line strengths out to four effective radii (Re) in the early-type galaxies NGC 3379 and NGC 821. With our newly developed observing technique we can now probe these faint regions in galaxies that were previously not accessible with traditional long-slit spectroscopy. We make optimal use of the large field-of-view and high throughput of the spectrograph: by adding the signal of all ~1400 lenslets into one spectrum, we obtain sufficient signal-to-noise in a few hours of observing time to reliably measure the absorption line kinematics and line strengths out to large radius. We find that the line strength gradients previously observed within 1 Re remain constant out to at least 4 Re, which puts constraints on the merger histories of these galaxies. The stellar halo populations are old and metal-poor. By constructing orbit-based Schwarzschild dynamical models we find that dark matter is necessary to explain the observed kinematics in NGC 3379 and NGC 821, with 30 - 50 per cent of the total matter being dark within 4 Re. The radial anisotropy in our best-fit halo models is less than in our models without halo, due to differences in orbital structure. The halo also has an effect on the Mgb - Vesc relation: its slope is steeper when a dark matter halo is added to the model.

[Abridged] We analysed two-dimensional maps of 48 early-type galaxies obtained with the SAURON and OASIS integral-field spectrographs using kinemetry, a generalisation of surface photometry to the higher order moments of the line-of-sight velocity distribution (LOSVD). In the SAURON sample, we find that 31% of early-type galaxies are single component systems. 91% of the multi-components systems have two kinematic subcomponents, the rest having three. In addition, 29% of galaxies have kinematically decoupled components, nuclear components with significant kinematic twists. We find that the velocity maps of fast rotators closely resemble those of inclined disks, except in the transition regions between kinematic subcomponents. In terms of E/S0 classification, this means that 74% of Es and 92% of S0s have components with disk-like kinematics. For the majority of fast rotators, the kinematic axial ratios are equal to or less than their photometric axial ratios, contrary to what is predicted with isotropic Jeans models viewed at different inclinations. The position angles of fast rotators are constant, while they vary abruptly in slow rotators. Velocity dispersion maps of face-on galaxies have shapes similar to the distribution of light. We constructed local (bin-by-bin) h3 - V/sigma and h4 - V/sigma diagrams from SAURON observations. We confirm the classical anti-correlation of h3 and V\sigma, but we also find that h3 is almost zero in some objects or even weakly correlated with V/sigma. The distribution of h4 for fast and slow rotators is mildly positive on average. The difference between slow and fast rotators is traceable throughout all moments of the LOSVD, with evidence for different intrinsic shapes and orbital contents and, hence, likely different evolutionary paths.

We discuss some recent integral field spectroscopy using the SAURON instrument of a sample consisting of 24 early-type spirals, part of the SAURON Survey, and 18 late-type spirals. Using 2-dimensional maps of their stellar radial velocity, velocity dispersion, and absorption line strength, it is now much easier to understand the nature of nearby galactic bulges. We discuss a few highlights of this work, and point out some new ideas about the formation of galactic bulges.

Using the unique dataset obtained within the course of the SAURON project, a radically new view of the structure, dynamics and stellar populations of early-type galaxies has emerged. We show that galaxies come in two broad flavours (slow and fast rotators), depending on whether or not they exhibit clear large-scale rotation, as indicated via a robust measure of the specific angular momentum of baryons. This property is also linked with other physical characteristics of early-type galaxies, such as: the presence of dynamically decoupled cores, orbital structure and anisotropy, stellar populations and dark matter content. I here report on the observed link between this baryonic angular momentum and a mass sequence, and how this uniquely relates to the building of the red sequence via dissipative/dissipationless mergers and secular evolution.

Supermassive black holes are a key element in our understanding of how galaxies form. Most of the progress in this very active field of research is based on just ~30 determinations of black hole mass, accumulated over the past decade. We illustrate how integral-field spectroscopy, and in particular our OASIS modeling effort, can help improve the current situation.

Motivated by recent progress in the study of early-type galaxies owing to technological advances, the launch of new space telescopes and large ground-based surveys, we attempt a short review of our current understanding of the recent star-formation activity in such intriguing galactic systems.

In this proceeding we make use of the two-dimensional stellar kinematics of a representative sample of E and S0 galaxies obtained with the SAURON integral-field spectrograph to reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We measure the level of rotation via a new parameter LambdaR and use it as a basis for a new kinematic classification that separates early-type galaxies into slow and fast rotators. With the aid of broad-band imaging we will reinforce this finding by comparing our kinematic results to the photometric properties of these two classes.

In this proceeding we look at the relationship between the photometric nuclear properties of early-type galaxies from Hubble Space Telescope imaging and their overall kinematics as observed with the SAURON integral-field spectrograph. We compare the inner slope of their photometric profiles and the Slow/Fast rotator classes, defined by the amplitude of a newly defined LambdaR parameter, to show that slow rotators tend to be more massive systems and display shallower inner profiles and fast rotators steper ones. It is important to remark, however, that there is not a one-to-one relationship between the two photometric and kinematic groups.

2D stellar kinematics of 48 representative E and S0 galaxies obtained with the SAURON IFS reveal that early-type galaxies appear in two broad flavours, depending on whether they exhibit clear large-scale rotation or not. We define a new parameter LambdaR which involves luminosity weighted averages over the full 2D kinematic field, as a proxy to quantify the observed projected stellar angular momentum per unit mass. We use it as a basis for a new kinematic classification: early-type galaxies are separated into slow and fast rotators (SRs, FRs), depending on whether they have LambdaR values within their effective radius Re below or above 0.1, respectively. SRs and FRs are shown to be physically distinct classes of galaxies, a result which cannot simply be the consequence of a biased viewing angle. FRs tend to be relatively low luminosity galaxies. SRs tend to be brighter and more massive galaxies, but are still spread over a wide range of absolute magnitude. 3 slow rotators of our sample, among the most massive ones, are consistent with zero rotation. Remarkably, all other SRs contain a large kpc-scale KDC. All FRs show well aligned photometric and kinemetric axes, and small velocity twists, in contrast with most SRs which exhibit significant misalignments and velocity twists. In a companion paper (Paper X), we also show that FRs and SRs are distinct classes in terms of their orbital distribution. We suggest that gas is a key ingredient in the formation and evolution of FRs, and that the slowest rotators are the extreme evolutionary end point reached deep in gravitational potential wells where dissipationless mergers had a major role in the evolution, and for which most of the baryonic angular momentum was expelled outwards. (abridged)

We present absorption line-strength maps for a sample of 18 Sb-Sd galaxies observed using the integral-field spectrograph SAURON. The SAURON spectral range allows the measurement of the Lick/IDS indices Hbeta, Fe5015 and Mgb, which can be used to estimate the stellar population parameters. We present here the two-dimensional line-strength maps for each galaxy. From the maps, we learn that late-type spiral galaxies tend to have high Hbeta and low Fe5015 and Mgb values, and that the Hbeta index has often a positive gradient over the field, while the metal indices peak in the central region. We investigate the relations between the central line-strength indices and their correlations with morphological type and central velocity dispersion, and compare the observed behaviour with that for ellipticals, lenticulars and early-type spirals from the SAURON survey. We find that our galaxies lie below the Mg - sigma relation determined for elliptical galaxies and that the indices show a clear trend with morphological type. From the line-strength maps we calculate age, metallicity and abundance ratio maps and discuss the results from a one-SSP approach and from a two-SSP approach. Late-type galaxies are generally younger and more metal poor than ellipticals and have abundance ratios closer to solar values. We also explore a continuous star formation scenario, and try to recover the star formation history using the evolutionary models of Bruzual & Charlot (2003), assuming constant or exponentially declining star formation rate (SFR). We find a correlation between the e-folding time-scale tau of the starburst and the central velocity dispersion: more massive galaxies tend to have shorter tau, suggesting that the star formation happened long ago and has now basically ended, while for smaller objects with larger values of tau it is still active.

We present absorption line strength maps of a sample of 24 representative early-type spiral galaxies, mostly of type Sa, obtained as part of the SAURON survey of nearby galaxies using our custom-built integral-field spectrograph. Using high-quality spectra, spatially binned to a constant signal-to-noise, we measure several key age, metallicity and abundance ratio sensitive indices from the Lick/IDS system over a contiguous two-dimensional field including bulge and inner disc. We present maps of H beta, Fe 5015 and Mg b, for each galaxy The absorption line maps show that many galaxies contain some younger populations (<= 1 Gyr), distributed in small or large inner discs, or in circumnuclear star forming rings. In many cases these young stars are formed in circumnuclear mini-starbursts, which are dominating the light in the centres of some of the early-type spirals. These mini-starburst cause a considerable scatter in index-index diagrams such as Mg b- H beta and Mg b -Fe 5015, more than is measured for early-type galaxies. We find that the central regions of Sa galaxies display a wide range in ages, even within the galaxies. 50% of the sample show velocity dispersion drops in their centres. All of the galaxies of our sample lie on or below the Mg b- $\sigma$ relation for elliptical galaxies in the Coma cluster, and above the H beta absorption line - $\sigma$ relation for elliptical galaxies. If those relations are considered to be relations for the oldest local galaxies we see that our sample of spirals has a considerable scatter in age, with the largest scatter at the lowest $\sigma$. This is in disagreement with highly inclined samples, in which generally only old stellar populations are found in the central regions. All this can be understood if ... (see paper for rest of abstract)

We analyse the orbital distribution of E and S0 galaxies using SAURON integral-field stellar kinematics. We construct the anisotropy (V/sigma,epsilon) diagram, for the 48 E/S0 galaxies from the SAURON survey. For a subsample of 24 galaxies consistent with axisymmetry, we use three-integral axisymmetric Schwarzschild dynamical models to recover the detailed orbital distribution and we find good agreement with the anisotropy derived from the (V/sigma,epsilon) diagram. In the companion Paper IX we show that the early-type galaxies can be subdivided into two classes of systems with or without a significant amount of specific stellar angular momentum. Here we show that the two classes have different distributions on the (V/sigma,epsilon) diagram. The slow rotators are more common among the most massive systems and are generally classified as E from photometry alone. Those in our sample tend to be fairly round (epsilon<0.3), but can have significant kinematical misalignments, indicating that as a class they are moderately triaxial, and span a range of anisotropies (delta<0.3). The fast rotators are generally fainter and are classified either E or S0. They can appear quite flattened (epsilon<0.7), do not show significant kinematical misalignments (unless barred or interacting), indicating they are nearly axisymmetric, and span an even larger range of anisotropies (delta<0.5). These results are confirmed when we extend our analysis to 18 additional E/S0 galaxies observed with SAURON. All these results support the idea that fast rotators are nearly oblate and contain disk-like components. The role of gas must have been important for their formation. The slow rotators are weakly triaxial. Current collisionless merger models seem unable to explain their detailed observed properties. (Abridged)

We discuss SAURON absorption line strength maps of a sample of 24 early-type spirals, mostly Sa. From the Lick indices H$\beta$, Mg b and Fe 5015 we derive SSP-ages and metallicities. By comparing the scaling relations of Mg b and H beta and central velocity dispersion with the same relation for the edge-on sample of Falcon-Barroso et al. (2002) we derive a picture in which the central regions of Sa galaxies contain at least 2 components: one (or more) thin, disc-like component, often containing recent star formation, and another, elliptical-like component, consisting of old stars and rotating more slowly, dominating the light above the plane. If one defines a bulge to be the component responsible for the light in excess of the outer exponential disc, then many Sa-bulges are dominated by a thin, disc-like component containing recent star formation.

By detecting ionised-gas emission in 75% of the cases, the SAURON integral-field spectroscopic survey has further demonstrated that early-type galaxies often display nebular emission. Furthermore, the SAURON data have shown that such emission comes with an intriguing variety of morphologies, kinematic behaviours and line ratios. Perhaps most puzzling was the finding that round and slowly rotating objects generally display uncorrelated stellar and gaseous angular momenta, consistent with an external origin for the gas, whereas flatter and fast rotating galaxies host preferentially co-rotating gas and stars, suggesting internal production of gas. Alternatively, a bias against the internal production of ionised gas and against the acquisition of retrograde material may be present in these two kinds of objects, respectively. In light of the different content of hot gas in these systems, with slowly rotating objects being the only systems capable of hosting massive X-ray halos, we suggest that a varying importance of evaporation of warm gas in the hot interstellar medium can contribute to explain the difference in the relative behaviour of gas and stars in these two kinds of objects. Namely, whereas in X-ray bright and slowly rotating galaxies stellar-loss material would quickly evaporate in the hot medium, in X-ray faint and fast rotating objects such material would be allowed to lose angular momentum and settle in a disk, which could also obstruct the subsequent acquisition of retrograde gas. Evidence for a connection between warm and hot gas phases, presumably driven by heat conduction, is presented for four slowly rotating galaxies with CHANDRA observations.

We present high spatial resolution integral-field spectroscopy of 28 elliptical (E) and lenticular (S0) galaxies from the SAURON representative survey obtained with the OASIS spectrograph during its operation at the CFHT. These seeing-limited observations explore the central 8"x10" (typically one kiloparsec diameter) regions of these galaxies using a spatial sampling four times higher than SAURON (0.27" vs. 0.94" spatial elements), resulting in almost a factor of two improvement in the median PSF. These data allow accurate study of the central regions to complement the large-scale view provided by SAURON. Here we present the stellar and gas kinematics, stellar absorption-line strengths and nebular emission-line strengths for this sample. We also characterise the stellar velocity maps using the "kinemetry" technique, and derive maps of the luminosity-weighted stellar age, metallicity and abundance ratio via stellar population models. We give a brief review of the structures found in our maps, linking also to larger-scale structures measured with SAURON. We present two previously unreported kinematically-decoupled components (KDCs) in the centres of NGC3032 and NGC4382. We compare the intrinsic size and luminosity-weighted stellar age of all the visible KDCs in the full SAURON sample, and find two types of components: kiloparsec-scale KDCs, which are older than 8 Gyr, and are found in galaxies with little net rotation; and compact KDCs, which have intrinsic diameters of less than a few hundred parsec, show a range of stellar ages from 0.5 - 15 Gyr (with 5/6 younger than 5 Gyr), are found exclusively in fast-rotating galaxies, and are close to counter-rotating around the same axis as their host. (abridged)

The Multi Unit Spectroscopic Explorer (MUSE) is a second-generation VLT panoramic integral-field spectrograph under preliminary design study. MUSE has a field of 1x1 arcmin**2 sampled at 0.2x0.2 arcsec**2 and is assisted by the VLT ground layer adaptive optics ESO facility using four laser guide stars. The simultaneous spectral range is 465-930 nm, at a resolution of R~3000. MUSE couples the discovery potential of a large imaging device to the measuring capabilities of a high-quality spectrograph, while taking advantage of the increased spatial resolution provided by adaptive optics. This makes MUSE a unique and tremendously powerful instrument for discovering and characterizing objects that lie beyond the reach of even the deepest imaging surveys. MUSE has also a high spatial resolution mode with 7.5x7.5 arcsec**2 field of view sampled at 25 milli-arcsec. In this mode MUSE should be able to obtain diffraction limited data-cubes in the 600-930 nm wavelength range. Although the MUSE design has been optimized for the study of galaxy formation and evolution, it has a wide range of possible applications; e.g. monitoring of outer planets atmosphere, environment of young stellar objects, super massive black holes and active nuclei in nearby galaxies or massive spectroscopic surveys of stellar fields in the Milky Way and nearby galaxies.

(Abridged) We present observations of the stellar and gas kinematics for a representative sample of 24 Sa galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. Our maps typically cover the bulge dominated region. We find a significant fraction of kinematically decoupled components (12/24), many of them displaying central velocity dispersion minima. They are mostly aligned and co-rotating with the main body of the galaxies, and are usually associated with dust discs and rings detected in unsharp-masked images. Almost all the galaxies in the sample (22/24) contain significant amounts of ionised gas which, in general, is accompanied by the presence of dust. The kinematics of the ionised gas is consistent with circular rotation in a disc co-rotating with respect to the stars. The distribution of mean misalignments between the stellar and gaseous angular momenta in the sample suggest that the gas has an internal origin. The [OIII]/Hbeta ratio is usually very low, indicative of current star formation, and shows various morphologies (ring-like structures, alignments with dust lanes or amorphous shapes). The star formation rates in the sample are comparable with that of normal disc galaxies. Low gas velocity dispersion values appear to be linked to regions of intense star formation activity. We interpret this result as stars being formed from dynamically cold gas in those regions. In the case of NGC5953, the data suggest that we are witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC5954.

We present results from our ongoing effort to understand the nature and evolution of nearby galaxies using the SAURON integral-field spectrograph. In this proceeding we focus on the study of the particular case formed by the interacting galaxies NGC5953 and NGC5954. We present stellar and gas kinematics of the central regions of NGC5953. We use a simple procedure to determine the age of the stellar populations in the central regions and argue that we may be witnessing the formation of a kinematically decoupled component from cold gas being acquired during the ongoing interaction with NGC5954.

We present results from a series of follow-up observations of a sub-sample of the representative SAURON survey elliptical (E) and lenticular (S0) galaxies using the OASIS integral-field spectrograph. These observations focus on the central 10" x 10", with roughly double the spatial resolution of the SAURON observations. This increased spatial resolution reveals a number of interesting and previously unresolved features in the measured stellar kinematics and absorption-line strengths. We find that galaxies exhibiting the youngest global stellar populations (as measured with SAURON) often contain a distinctly young central region (on scales of a few hundred parsec or less) compared to the rest of the galaxy. Moreover, these compact, young components are found to be mostly counter-rotating with respect to the rest of the galaxy. Given that there is no well-established reason for such young components to `prefer' counter- over co-rotation, this finding raises the following questions: How common are these small KDCs as a function of age? Why are there more young than old compact KDCs? Where are the equivalent co-rotating components? We explore these questions using simple simulated velocity fields and stellar population models, and find that the fading of the young component as it evolves, coupled with the fact that counter-rotating components are more easily detected in the velocity field, may help explain the observed trends.

We present absorption line strength maps of 48 representative elliptical and lenticular galaxies obtained as part of the SAURON survey of nearby galaxies. Using high-quality spectra, spatially binned to a constant signal-to-noise, we measure four key age, metallicity and abundance ratio sensitive indices from the Lick/IDS system over a two-dimensional field extending up to approximately one effective radius. We modify the classical Fe5270 index to define a new index, Fe5270S, which maximizes the usable spatial coverage of SAURON. Maps of Hbeta, Fe5015, Mgb and Fe5270S are presented for each galaxy. We use the maps to compute average line strengths integrated over circular apertures of one-eighth effective radius, and compare the resulting relations of index versus velocity dispersion with previous long-slit work. The metal line strength maps show generally negative gradients with increasing radius roughly consistent with the morphology of the light profiles. Remarkable deviations from this general trend exist, particularly the Mgb isoindex contours appear to be flatter than the isophotes of the surface brightness for about 40% of our galaxies without significant dust features. Generally these galaxies exhibit significant rotation. We infer from this that the fast-rotating component features a higher metallicity and/or an increased Mg/Fe ratio as compared to the galaxy as a whole. The Hbeta maps are typically flat or show a mild positive outwards radial gradient, while a few galaxies show strong central peaks and/or elevated overall Hbeta-strength likely connected to recent star-formation activity. For the most prominent post-starburst galaxies even the metal line strength maps show a reversed gradient. (abridged)

We present selected results from integral-field spectroscopy of 48 early-type galaxies observed as part of the SAURON survey. Maps of the Hbeta, Fe5015, Mgb and Fe5270 indices in the Lick/IDS system were derived for each of the survey galaxies. The metal line strength maps show generally negative gradients with increasing radius roughly consistent with the morphology of the light profiles. Remarkable deviations from this general trend exist, particularly the Mgb isoindex contours appear to be flatter than the isophotes of the surface brightness for about 40% of our galaxies without significant dust features. Generally these galaxies exhibit significant rotation. We infer from this that the fast-rotating component features a higher metallicity and/or an increased Mg/Fe ratio as compared to the galaxy as a whole. We also use the line strengths maps to compute average values integrated over circular apertures of one effective radius, and derive luminosity weighted ages and metallicities. The lenticular galaxies show a wide range in age and metallicity estimates, while elliptical galaxies tend to occupy regions of older stellar populations.

We investigate with unprecedented accuracy the correlations between the dynamical mass-to-light ratio M/L and other global observables of E and S0 galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics. We find a tight correlation of the form (M/L)=(3.80+/-0.14)*(sigma/200 km/s)^(0.84+/-0.07) between the dynamical M/L (in the I-band) and the luminosity-weighted second moment (sigma) of the line-of-sight velocity-distribution within Re. The observed rms scatter in M/L for our sample is 18%, while the inferred intrinsic scatter is ~13%. The (M/L)-sigma relation can be included in the remarkable series of tight correlations between sigma and other galaxy global observables. The comparison of the observed correlations with the predictions of the Fundamental Plane (FP), and with simple virial estimates, shows that the `tilt' of the FP of early-type galaxies, is due to a real M/L variation, while structural and orbital non-homology have a negligible effect. The virial mass is a reliable estimator of the mass in the central regions of galaxies. The best-fitting virial relation has the form (M/L)_vir=(5.0+/-0.1)*Re*sigma^2/(L*G). The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates a median dark matter fraction in early-type galaxies of ~30% of the total mass inside one Re. (Abridged)

We present results from our ongoing effort to understand the morphological and kinematical properties of early-type galaxies using the integral-field spectrograph SAURON. We discuss the relation between the stellar and gas morphology and kinematics in our sub-sample of 24 representative Sa spiral bulges. We focus on the frequency of kinematically decoupled components and on the presence of star formation in circumnuclear rings.

We investigate the well-known correlations between the dynamical mass-to-light ratio M/L and other global observables of elliptical (E) and lenticular (S0) galaxies. We construct two-integral Jeans and three-integral Schwarzschild dynamical models for a sample of 25 E/S0 galaxies with SAURON integral-field stellar kinematics to about one effective (half-light) radius R_e. The comparison of the dynamical M/L with the (M/L)_pop inferred from the analysis of the stellar population, indicates that dark matter in early-type galaxies contributes ~30% of the total mass inside one R_e, in agreement with previous studies, with significant variations from galaxy to galaxy. Our results suggest a variation in M/L at constant (M/L)_pop, which seems to be linked to the galaxy dynamics. We speculate that fast rotating galaxies have lower dark matter fractions than the slow rotating and generally more massive ones.

We present the emission-line fluxes and kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph SAURON. Hb, [OIII], and [NI] emission lines were measured using a new procedure that simultaneously fits both the stellar spectrum and the emission lines. Using this technique we can detect emission lines down to an equivalent width of 0.1A set by the current limitations in describing galaxy spectra with synthetic and real stellar templates, rather than by the quality of our spectra. Emission is detected in 75% of our sample galaxies, and comes in a variety of resolved spatial distributions and kinematic behaviours. The ionised-gas kinematics is rarely consistent with simple coplanar circular motions. However, the gas almost never displays completely irregular kinematics, generally showing coherent motions with smooth variations in angular momentum. In the majority of the cases the gas kinematics is decoupled from the stellar kinematics, and in half of the objects this decoupling implies a recent acquisition of gaseous material. Over the entire sample however, the distribution of the mean misalignment values between stellar and gaseous angular momenta is inconsistent with a purely external origin. Consistent with previous studies, the presence of dust features is always accompanied by gas emission while the converse is not always true. A considerable range of values for the [OIII]/Hb ratio is found both across the sample and within single galaxies. Despite the limitations of this ratio as an emission-line diagnostic, this finding suggests either that a variety of mechanisms is responsible for the gas excitation in E and S0 galaxies or that the metallicity of the interstellar material is quite heterogeneous.

We use integral-field observation of the stellar kinematics obtained with SAURON in combination with Schwarzschild dynamical models to revisit our understanding of the classic (V/sigma)-epsilon anisotropy diagram of early-type galaxies.

We summarise the results and achievements of integral-field spectroscopy of early-type galaxies, observed as part of a survey using both the SAURON and OASIS spectrographs. From the perspective of integral-field spectroscopy, these otherwise smooth and featureless objects show a wealth of structure, both in their stellar kinematics and populations. We focus on the stellar content, and examine properties on both kiloparsec scales with SAURON, and scales of 100's of parsecs with OASIS. These complementary studies reveal two types of kinematically distinct components (KDCs), differing primarily in their intrinsic sizes. In previous studies, KDCs and their host galaxies have generally been found to be unremarkable in other aspects. We show that large KDCs, typical of the well-studied cases, indeed show little or no age differences with their host galaxy. The KDCs detected with the higher spatial-resolution of OASIS are intrinsically smaller and include, in contrast, a significant fraction of young stars. We speculate on the relationship between KDCs and their host galaxies, and the implications for young populations in early-type galaxies.

We present results of our ongoing study of the morphology and kinematics of the ionised gas in 48 representative nearby elliptical and lenticular galaxies using the SAURON integral-field spectrograph on the 4.2m William Herschel Telescope. Making use of a recently developed technique, emission is detected in 75% of the galaxies. The ionised-gas distributions display varied morphologies, ranging from regular gas disks to filamentary structures. Additionally, the emission-line kinematic maps show, in general, regular motions with smooth variations in kinematic position angle. In most of the galaxies, the ionised-gas kinematics is decoupled from the stellar counterpart, but only some of them present signatures of recent accretion of gaseous material. The presence of dust is very common in our sample and is usually accompanied by gas emission. Our analysis of the [OIII]/Hbeta emission-line ratios, both across the whole sample as well as within the individual galaxies, suggests that there is no unique mechanism triggering the ionisation of the gas.

High-velocity galactic outflows, driven by intense bursts of star formation and black hole accretion, are invoked by current theories of galaxy formation to terminate star formation in the most massive galaxies and to deposit heavy elements in the intergalactic medium. From existing observational evidence on high-redshift galaxies, it is unclear whether such outflows are localized to regions of intense star formation just a few kiloparsecs in extent, or whether they instead have a significant impact on the entire galaxy and its surroundings. Here we present two-dimensional spectroscopy of a star-forming galaxy at redshift z=3.09 (seen 11.5 Gyr ago, when the Universe was 20 per cent of its current age): its spatially extended Ly-alpha emission appears to be absorbed by HI in a foreground screen covering the entire galaxy, with a lateral extent of at least 100 kpc and remarkable velocity coherence. It was plausibly ejected from the galaxy during a starburst several 1E8 yr earlier and has subsequently swept up gas from the surrounding intergalactic medium and cooled. This demonstrates the galaxy-wide impact of high-redshift superwinds.

We present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

Supermassive black holes in the centres of giant elliptical galaxies are thought to induce chaos and eliminate triaxiality in their hosts. We address whether this process operates in real systems, by modeling the stellar kinematics of the old elliptical NGC 4365. This galaxy has a mean stellar population age > 12 Gyr and is known for its kinematically decoupled core and skew rotation at larger radii. We fit the two-dimensional mean velocity field obtained by the SAURON integral-field spectrograph, and the isophotal ellipticity and position-angle profiles, using the velocity field fitting approach. We find NGC 4365 to be strongly triaxial (triaxiality parameter T approximately 0.45) and somewhat flatter than it appears. Axisymmetry or near axisymmetry (T < 0.1) is ruled out at > 95% confidence. In the most probable orientation the long axis points roughly toward the observer, extending to the southwest in projection. Strong triaxiality has persisted for hundreds of dynamical times, ruling out black holes above 3 x 10^9 solar masses. The M-sigma relation predicts 4 x 10^8 solar masses, which would not preclude long-lived triaxiality and is consistent with the observations. NGC 4365 lends support to a picture in which supermassive black holes, while omnipresent in luminous giant elliptical galaxies, are not massive enough to alter their global structure through chaos.

We present the stellar kinematics of 48 representative elliptical and lenticular galaxies obtained with our custom-built integral-field spectrograph SAURON operating on the William Herschel Telescope. The data were homogeneously processed through a dedicated reduction and analysis pipeline. All resulting SAURON datacubes were spatially binned to a constant minimum signal-to-noise. We have measured the stellar kinematics with an optimized (penalized pixel-fitting) routine which fits the spectra in pixel space, via the use of optimal templates, and prevents the presence of emission lines to affect the measurements. We have thus generated maps of the mean stellar velocity, the velocity dispersion, and the Gauss-Hermite moments h3 and h4 of the line-of-sight velocity distributions. The maps extend to approximately one effective radius. Many objects display kinematic twists, kinematically decoupled components, central stellar disks, and other peculiarities, the nature of which will be discussed in future papers of this series.

We have used the SAURON panoramic integral field spectrograph to study the structure of the Ly-alpha emission-line halo, LAB1, surrounding the sub-millimeter galaxy SMM J221726+0013. This emission-line halo was discovered during a narrow-band imaging survey of the z=3.1 large-scale structure in the SSA22 region. Our observations trace the emission halo out to almost 100 kpc from the sub-millimeter source and identify two distinct Ly-alpha `mini-haloes' around the nearby Lyman-break galaxies. The main emission region has a broad line profile, with variations in the line profile seeming chaotic and lacking evidence for a coherent velocity structure. The data also suggests that Ly-alpha emission is suppressed around the sub-mm source. Interpretation of the line structure needs care because Ly-alpha may be resonantly scattered, leading to complex radiative transfer effects, and we suggest that the suppression in this region arises because of such effects. We compare the structure of the central emission-line halo with local counter-parts, and find that the emission line halo around NGC 1275 in the Perseus cluster may be a good local analogue, although the high redshift halo is factor of ~100 more luminous and appears to have higher velocity broadening. Around the Lyman-break galaxy C15, the emission line is narrower, and a clear shear in the emission wavelength is seen. A plausible explanation for the line profile is that the emission gas is expelled from C15 in a bipolar outflow, similar to that seen in M82.

We present SAURON integral-field observations of the S0 galaxy NGC7332. Existing broad-band ground-based and HST photometry reveals a double disk structure and a boxy bulge interpreted as a bar viewed close to edge-on. The SAURON two-dimensional stellar kinematic maps confirm the existence of the bar and inner disk but also uncover the presence of a cold counter-rotating stellar component within the central 250 pc. The Hbeta and [OIII] emission line maps show that the ionised gas has a complex morphology and kinematics, including both a component counter-rotating with respect to the stars and a fainter co-rotating one. Analysis of the absorption line-strength maps show that NGC7332 is young everywhere. The presence of a large-scale bar can explain most of those properties, but the fact that we see a significant amount of unsettled gas, together with a few peculiar features in the maps, suggest that NGC7332 is still evolving. Interactions as well as bar-driven processes must thus have played an important role in the formation and evolution of NGC7332, and presumably of S0 galaxies in general.

The Nearby Supernova Factory (SNfactory) is an ambitious project to find and study in detail approximately 300 nearby Type Ia supernovae (SNe~Ia) at redshifts 0.03<z<0.08. This program will provide an exceptional data set of well-studied SNe in the nearby smooth Hubble flow that can be used as calibration for the current and future programs designed to use SNe to measure the cosmological parameters. The first key ingredient for this program is a reliable supply of Hubble-flow SNe systematically discovered in unprecedented numbers using the same techniques as those used in distant SNe searches. In 2002, 35 SNe were found using our test-bed pipeline for automated SN search and discovery. The pipeline uses images from the asteroid search conducted by the Near Earth Asteroid Tracking group at JPL. Improvements in our subtraction techniques and analysis have allowed us to increase our effective SN discovery rate to ~12 SNe/month in 2003.

We present a summary of high-spatial resolution follow-up observations of the elliptical (E) and lenticular (S0) galaxies in the SAURON survey using the OASIS integral field spectrograph. The OASIS observations explore the central 8x10" regions of these galaxies using a spatial sampling four times higher than SAURON, often revealing previously undiscovered features. Around 75% (31/48) of the SAURON E/S0s with central velocity dispersion >= 120 km/s were observed with OASIS, covering well the original SAURON representative sample. We present here an overview of this follow-up survey, and some preliminary results on individual objects, including a previously unreported counter-rotating core in NGC 4382; the decoupled stellar and gas velocity fields of NGC 2768; and the strong age gradient towards the centre of NGC 3489.

Kinematic and population studies show that bulges are generally rotationally flattened systems similar to low-luminosity ellipticals. However, observations with state-of-the-art integral field spectrographs, such as SAURON, indicate that the situation is much more complex, and allow us to investigate phenomena such as triaxiality, kinematic decoupling and population substructure, and to study their connection to current formation and evolution scenarios for bulges of early-type galaxies. We present the examples of two S0 bulges from galaxies in our sample of nearby galaxies: one that shows all the properties expected from classical bulges (NGC5866), and another case that presents kinematic features appropriate for barred disk galaxies (NGC7332).

We describe our program for the dynamical modeling of early-type galaxies observed with the panoramic integral-field spectrograph SAURON. We are using Schwarzschild's numerical orbit superposition method to reproduce in detail all kinematical and photometric observables, and recover the intrinsic orbital structure of the galaxies. Since catastrophes are the most prominent features in the orbital observables, two-dimensional kinematical coverage is essential to constrain the dynamical models.

We have measured the two-dimensional kinematics and line-strength distributions of 72 representative nearby early-type galaxies, out to approximately one effective radius, with our panoramic integral-field spectrograph SAURON. The resulting maps reveal a rich variety in kinematical structures and linestrength distributions, indicating that early-type galaxies are more complex systems than often assumed. We are building detailed dynamical models for these galaxies, to derive their intrinsic shape and dynamical structure, and to determine the mass of the supermassive central black hole. Here we focus on two examples, the compact elliptical M32 and the E3 galaxy NGC4365. These objects represent two extreme cases: M32 has very regular kinematics which can be represented accurately by an axisymmetric model in which all stars rotate around the short axis, while NGC4365 is a triaxial galaxy with a prominent kinematically decoupled core, with an inner core that rotates about an axis that is nearly perpendicular to the rotation axis of the main body of the galaxy. Our dynamical models for these objects demonstrate that two-dimensional observations are essential for deriving the intrinsic orbital structure and dark matter content of galaxies.

We present results from a new and unique integral-field spectrograph, SAURON. It has a large field of view and high throughput and is primarily built for the study of stellar & gaseous kinematics and stellar populations in galaxies. Its aim is to carry out a systematic survey of the velocity fields, velocity dispersions, and line-strength distributions of nearby ellipticals, lenticular galaxies and spiral bulges. Its wide field is especially useful for the study of complicated velocity structures. Together with other spectroscopic data, images, and dynamical modelling, SAURON will help to constrain the intrinsic shapes, mass-to-light ratios, and stellar populations of early-type galaxies and spiral bulges.

We present dynamical models of the nearby compact elliptical galaxy M32, using high quality kinematical measurements, obtained with the integral-field spectrograph SAURON mounted on the William Herschel Telescope on La Palma. We also include STIS data obtained by Joseph et al. We find a best-fit black hole mass of (2.5 +/- 0.5) million solar masses and a stellar I-band mass-to-light ratio of (1.85 +/- 0.15) in solar units. For the first time, we are also able to constrain the inclination along which M32 is observed to (70 +/- 5) degrees. Combined with an averaged observed flattening of 0.73, this corresponds to an intrinsic flattening of approximately 0.68 +/- 0.03. These tight constraints are mainly caused by the use of integral-field data. We show this quantitatively by comparing with models that are constrained by multiple slits only. We show the phase-space distribution and intrinsic velocity structure of the best-fit model and investigate the effect of regularisation on the orbit distribution.

We present the SAURON project, which is aimed at studying the morphology, two-dimensional kinematics and stellar populations of a representative sample of elliptical galaxies and spiral bulges. SAURON, a dedicated integral-field spectrograph that is optimized for wide-field observations and has high throughput, was built in Lyon and is now operated at the WHT 4.2m telescope. At present, we have observed approximately two thirds of the seventy-two sample galaxies with SAURON. A comparison with published long-slit measurements demonstrates that the SAURON-data is of equal or better quality, and provides full two-dimensional coverage. The velocity and velocity dispersion fields exhibit a large variety of morphologies: from simple rotating systems to cylindrical, disky and triaxial velocity fields, bars and decoupled cores. Most of these kinematical signatures do not have counterparts in the light distribution. While some galaxies are consistent with axisymmetry, most are more complex systems than assumed previously. This suggests that the kinematical properties of nearby E/S0 galaxies do not agree with the often assumed simplistic two-family model, in which the giant non-rotating triaxial ellipticals are opposed to the fast-rotating axisymmetric faint ellipticals and S0s.

The panoramic integral-field spectrograph SAURON is currently being used to map the stellar kinematics, gaseous kinematics, and stellar populations of a large number of early-type galaxies and bulges. Here, we describe SAURON observations of cold stellar disks embedded in spheroids (NGC3384, NGC4459, NGC4526), we illustrate the kinematics and ionization state of large-scale gaseous disks (NGC4278, NGC7742), and we show preliminary comparisons of SAURON data with barred galaxy N-body simulations (NGC3623).

The SAURON project will deliver two-dimensional spectroscopic data of a sample of nearby early-type galaxies with unprecedented quality. In this paper, we focus on the mapping of their stellar populations using the SAURON data, and present some preliminary results on a few prototypical cases.

Early results are reported from the SAURON survey of the kinematics and stellar populations of a representative sample of nearby E, S0 and Sa galaxies. The survey is aimed at determining the intrinsic shape of the galaxies, their orbital structure, the mass-to-light ratio as a function of radius, the age and metallicity of the stellar populations, and the frequency of kinematically decoupled cores and nuclear black holes. The construction of the representative sample is described, and its properties are illustrated. A comparison with long-slit spectroscopic data establishes that the SAURON measurements are comparable to, or better than, the highest-quality determinations. Comparisons are presented for NGC 3384 and NGC 4365 where stellar velocities and velocity dispersions are determined to a precision of 6 km/s, and the h3 and h4 parameters of the line-of-sight velocity distribution to a precision of better than 0.02. Extraction of accurate gas emission-line intensities, velocities and line widths from the datacubes is illustrated for NGC 5813. Comparisons with published line-strengths for NGC 3384 and NGC 5813 reveal uncertainties of < 0.1 A on the measurements of the Hbeta, Mgb and Fe5270 indices. Integral-field mapping uniquely connects measurements of the kinematics and stellar populations to the galaxy morphology. The maps presented here illustrate the rich stellar kinematics, gaseous kinematics, and line-strength distributions of early-type galaxies. The results include the discovery of a thin, edge-on, disk in NGC 3623, confirm the axisymmetric shape of the central region of M32, illustrate the LINER nucleus and surrounding counter-rotating star-forming ring in NGC 7742, and suggest a uniform stellar population in the decoupled core galaxy NGC 5813.

A summary of the SAURON project and its current status is presented. SAURON is a panoramic integral-field spectrograph designed to study the stellar kinematics, gaseous kinematics, and stellar populations of spheroids. Here, the sample of galaxies and its properties are described. The instrument is detailed and its capabilities illustrated through observational examples. These includes results on the structure of central stellar disks, the kinematics and ionization state of gaseous disks, and the stellar populations of galaxies with decoupled cores.

We describe a new technique, kinemetry, to quantify kinematic maps of early-type galaxies in an efficient way. We present the first applications to velocity fields obtained with the integral-field spectrograph SAURON.

We describe MIFS, a second generation integral-field spectrograph for the VLT, operating in the visible wavelength range. It combines a 1'x1' field of view with the improved spatial resolution provided by multi-conjugate adaptive optics and covers a large simultaneous spectral range (0.6-1.0 microns). A separate mode exploits the highest spatial resolution provided by adaptive optics. With this unique combination of capabilities, MIFS has a wide domain of application and a large discovery potential. The MIFS low-spatial resolution mode (sampled at 0.2") combined with the initial MCAO capabilities planned for the VLT will provide ultra-deep fields with a limiting magnitude for spectroscopy of R~28. MIFS will improve the present day detection limit of Lyman alpha emitters by a factor of 100, and will detect low-mass star-forming galaxies to z~7. The MIFS high-spatial resolution mode (3"x3" field sampled at 0.01") is optimized for the next step in (MC)AO. It will probe, e.g., the relationship between supermassive central black holes and their host galaxy and the physics of winds from accretion disks in young stellar objects at unprecedented spatial resolution. MIFS will extend Europe's lead in integral-field spectroscopy. It capitalizes on new developments in adaptive optics, and is a key step towards instrumentation for OWL.

Early results are reported from the SAURON survey of the kinematics and stellar populations of nearby cluster and field E, S0, and Sa galaxies. We present maps of the Hbeta and [OIII] emission-line distribution and kinematics for NGC 5813 and NGC 7742.

A new integral-field spectrograph, SAURON, is described. It is based on the Tiger principle, and uses a lenslet array. SAURON has a large field of view and high throughput, and allows simultaneous sky subtraction. Its design is optimized for studies of the stellar kinematics, gas kinematics, and line-strength distributions of nearby early-type galaxies. The instrument design and specifications are described, as well as the extensive analysis software which was developed to obtain fully calibrated spectra, and the associated kinematic and line-strength measurements. A companion paper reports on the first results obtained with SAURON on the William Herschel Telescope.

We present observations with the adaptive optics assisted integral field spectrograph OASIS of the M 31 double nucleus at a spatial resolution better than 0.5 arcsec FWHM. These data are used to derive the two-dimensional stellar kinematics within the central 2 arcsec. Archival WFPC2/HST images are revisited to perform a photometric decomposition of the nuclear region. We also present STIS/HST kinematics obtained from the archive. The luminosity distribution of the central region is well separated into the respective contributions of the bulge, the nucleus including P1 and P2, and the so-called UV peak. We then show that the axis joining P1 and P2, the two local surface brightness maxima, does not coincide with the kinematic major-axis, which is also the major-axis of the nuclear isophotes (excluding P1). We also confirm that the velocity dispersion peak is offset by ~ 0.2 arcsec from the UV peak, assumed to mark the location of the supermassive black hole. The newly reduced STIS/HST velocity and dispersion profiles are then compared to OASIS and other published kinematics. We find significant offsets with previously published data. Simple parametric models are then built to successfully reconcile all the available kinematics. We finally interpret the observations using new N-body simulations. The nearly keplerian nuclear disk of M31 is subject to a natural m=1 mode, with a very slow pattern speed (3 km/s/pc for M_BH = 7 10^7~\Msun), that can be maintained during more than a thousand dynamical times. The resulting morphology and kinematics of the mode can reproduce the M~31 nuclear-disk photometry and mean stellar velocity, including the observed asymmetries. It requires a central mass concentration and a cold disk system representing between 20 and 40% of its mass. Abridged..

We report the first wide-field mapping of the kinematics and stellar populations in the E3 galaxy NGC 4365. The velocity maps extend previous long-slit work. They show two independent kinematic subsystems: the central 300 x 700 pc rotates about the projected minor axis, and the main body of the galaxy, 3 x 4 kpc, rotates almost at right angles to this. The line-strength maps show that the metallicity of the stellar population decreases from a central value greater than solar, to one-half solar at a radius of 2 kpc. The decoupled core and main body of the galaxy have the same luminosity-weighted age, of ~14 Gyr, and the same elevated magnesium-to-iron ratio. The two kinematically distinct components have thus shared a common star formation history. We infer that the galaxy underwent a sequence of mergers associated with dissipative star formation that ended >12 Gyr ago. The misalignment between the photometric and kinematic axes of the main body is unambiguous evidence of triaxiality. The similarity of the stellar populations in the two components suggests that the observed kinematic structure has not changed substantially in 12 Gyr.

We briefly describe the SAURON project, aimed at determining the intrinsic shape and internal dynamics of spheroids. We focus here on the ability of SAURON to identify gaseous and stellar disks and to measure their morphology and kinematics. We illustrate some of our results with complete maps of NGC3377, NGC3623, and NGC4365.

In this paper, The central black hole masses and the Doppler factors are derived for PKS 0528+134, PKS 0537-441, 3C279, PKS 1406-074, PKS 1622-297, Q1633+382, Mkn 501, and BL Lacertae. The masses obtained are in the range of (1 -7)$\times 10^{7}M_{\odot}$ and compared with that obtained with the Klein-Nishina cross section considered (Dermer & Gehrels 1995). If we considered only the Thomson cross section, the masses are in the range of 2.6$\times 10^{6}M_{\odot}$ - 2$\times 10^{11}M_{\odot}$. The masses obtained from our method are less sensitive to the flux than those obtained from Dermer & Gehrels (1995) method. The masses obtained from two flares (1991 and 1996 flares) of 3C279 are almost the same. For 3C279 and BL Lacertae, viewing angle, $\theta$, and Lorentz factor, $\Gamma$, are estimated from the derived Doppler factor and the measured superluminal velocity. For 3C279, $\theta = 10^{\circ}.9-15^{\circ}.6$, $\Gamma$ = 2.4-14.4 for $\delta$ = 3.37; $\theta = 8^{\circ}.45-9^{\circ}.7$, $\Gamma$ = 2.95-11.20 for $\delta$ = 4.89; For BL Lacertae, $\theta = 25^{\circ}-29^{\circ}.4$, $\Gamma$ = 2.0-4.0.

We present new dynamical models of the S0 galaxy N3115, making use of the available published photometry and kinematics as well as of two-dimensional TIGER spectrography. We first examined the kinematics in the central 40 arcsec in the light of two integral f(E,J) models. Jeans equations were used to constrain the mass to light ratio, and the central dark mass whose existence was suggested by previous studies. The even part of the distribution function was then retrieved via the Hunter & Qian formalism. We thus confirmed that the velocity and dispersion profiles in the central region could be well fit with a two-integral model, given the presence of a central dark mass of ~10^9 Msun. However, no two integral model could fit the h_3 profile around a radius of 25 arcsec where the outer disc dominates the surface brightness distribution. Three integral analytical models were therefore built using a Quadratic Programming technique. These models showed that three integral components do indeed provide a reasonable fit to the kinematics, including the higher Gauss-Hermite moments. Again, models without a central dark mass failed to reproduce the observed kinematics in the central arcseconds. This clearly supports the presence of a nuclear black hole of at least 6.5 10^8 Msun in the centre of NGC 3115. These models were finally used to estimate the importance of the dark matter in the outer part of NGC 3115, suggested by the flat stellar rotation curve observed by Capaccioli et al. (1993).

Sub-arcsecond imagery (HRCAM, 0".35 - 0".57 FWHM) and two-dimensional spectrography (TIGER, 0".9 FWHM) of the central nucleus of M31 have been obtained at CFHT. The photometric data clearly show the double-peaked nucleus, in excellent agreement with a recent HST image by Lauer et al. 1993. We built deconvolved surface brightness models, using the multi-Gaussian expansion method. We then perform a detailed morphological analysis of the three central photometric components (bulge, nucleus and bright secondary peak) and derive various spatial luminosity models (oblate and triaxial). Stellar velocity and velocity dispersion fields were derived from the TIGER data: the former displays an extremely rapid rotation around the true center of the galaxy, while the latter exhibits a peaked structure offset in the opposite direction of the brightest light peak. Neglecting these offsets,both extended versions of the virial theorem and detailed hydrodynamical models confirm the classical strong central mass concentration, perhaps a supermassive black hole of about 7.10**7 Msun, as well as a large excess of circular motions. These offsets, however, suggest that the nucleus presently undergoes a strong stellar oscillation, with non-stationary, non-axisymmetric dynamics. This too reopens the case for a strong central mass concentration.