Source author record

M. Fabricius

M. Fabricius appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2025arXiv

Counterpart identification and classification for eRASS1 and characterisation of the AGN content

[abridged] Accurately accounting for the AGN phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the Western Galactic Hemisphere. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multiwavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a Machine Learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using Circlez. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; about 570,000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the surveys shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic Plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identify approximately 18,000 Gaia and 55,000 CatWISE2020 extragalactic sources. We release three high-quality counterpart catalogues, plus the training and validation sets, as a benchmark for the field. These datasets have many applications, but in particular empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universes most energetic engines.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

First Light for GRAVITY Wide: Large Separation Fringe Tracking for the Very Large Telescope Interferometer

GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth's turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth's atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning.

preprint2022arXiv

GRAVITY faint: reducing noise sources in GRAVITY$^+$ with a fast metrology attenuation system

With the upgrade from GRAVITY to GRAVITY$^+$ the instrument will evolve into an all-sky interferometer that can observe faint targets, such as high redshift AGN. Observing the faintest targets requires reducing the noise sources in GRAVITY as much as possible. The dominant noise source, especially in the blue part of the spectrum, is the backscattering of the metrology laser light onto the detector. To reduce this noise we introduce two new metrology modes. With a combination of small hardware changes and software adaptations, we can dim the metrology laser during the observation without losing the phase referencing. For single beam targets, we can even turn off the metrology laser for the maximum SNR on the detector. These changes lead to an SNR improvement of over a factor of two averaged over the whole spectrum and up to a factor of eight in the part of the spectrum currently dominated by laser noise.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2016arXiv

The angular momentum distribution and baryon content of star forming galaxies at z~1-3

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$α$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$λ$> ~ 0.037 and its dispersion ($σ_{log(λ)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $λ$ x (j_d/j_DM) and the nuclear stellar density $Σ_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

preprint2016arXiv

The SINFONI Black Hole Survey: The Black Hole Fundamental Plane revisited and the paths of (co-) evolution of supermassive black holes and bulges

We investigate the correlations between the black hole mass $M_{BH}$, the velocity dispersion $σ$, the bulge mass $M_{Bu}$, the bulge average spherical density $ρ_h$ and its spherical half mass radius $r_h$, constructing a database of 97 galaxies (31 core ellipticals, 17 power-law ellipticals, 30 classical bulges, 19 pseudo bulges) by joining 72 galaxies from the literature to 25 galaxies observed during our recent SINFONI black hole survey. For the first time we discuss the full error covariance matrix. We analyse the well known $M_{BH}-σ$ and $M_{BH}-M_{Bu}$ relations and establish the existence of statistically significant correlations between $M_{Bu}$ and $r_h$ and anti-correlations between $M_{Bu}$ and $ρ_h$. We establish five significant bivariate correlations ($M_{BH}-σ-ρ_h$, $M_{BH}-σ-r_h$, $M_{BH}-M_{Bu}-σ$, $M_{BH}-M_{Bu}-ρ_h$, $M_{BH}-M_{Bu}-r_h$) that predict $M_{BH}$ of 77 core and power-law ellipticals and classical bulges with measured and intrinsic scatter as small as $\approx 0.36$ dex and $\approx 0.33$ dex respectively, or 0.26 dex when the subsample of 45 galaxies defined by Kormendy and Ho (2013) is considered. In contrast, pseudo bulges have systematically lower $M_{BH}$, but approach the predictions of all the above relations at spherical densities $ρ_h\ge 10^{10} M_\odot/kpc^3$ or scale lengths $r_h\le 1$ kpc. These findings fit in a scenario of co-evolution of BH and classical-bulge masses, where core ellipticals are the product of dry mergers of power-law bulges and power-law Es and bulges the result of (early) gas-rich mergers and of disk galaxies. In contrast, the (secular) growth of BHs is decoupled from the growth of their pseudo bulge hosts, except when (gas) densities are high enough to trigger the feedback mechanism responsible for the existence of the correlations between $M_{BH}$ and galaxy structural parameters.

preprint2014arXiv

Evidence for Wide-Spread AGN Driven Outflows in the Most Massive z~1-2 Star Forming Galaxies

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

preprint2014arXiv

The Dynamical Fingerprint of Core Scouring in Massive Elliptical Galaxies

The most massive elliptical galaxies have low-density centers or cores that differ dramatically from the high-density centers of less massive ellipticals and bulges of disk galaxies. These cores have been interpreted as the result of mergers of supermassive black hole binaries, which depopulate galaxy centers by gravitationally slingshotting central stars toward large radii. Such binaries naturally form in mergers of luminous galaxies. Here, we analyze the population of central stellar orbits in 11 massive elliptical galaxies that we observed with the integral field spectrograph SINFONI at the European Southern Observatory Very Large Telescope. Our dynamical analysis is orbit-based and includes the effects of a central black hole, the mass distribution of the stars, and a dark matter halo. We show that the use of integral field kinematics and the inclusion of dark matter is important to conclude upon the distribution of stellar orbits in galaxy centers. Six of our galaxies are core galaxies. In these six galaxies, but not in the galaxies without cores, we detect a coherent lack of stars on radial orbits in the core region and a uniform excess of radial orbits outside of it: when scaled by the core radius, the radial profiles of the classical anisotropy parameter beta are nearly identical in core galaxies. Moreover, they match quantitatively the predictions of black hole binary simulations, providing the first convincing dynamical evidence for core scouring in the most massive elliptical galaxies.

preprint2014arXiv

The KMOS^3D Survey: design, first results, and the evolution of galaxy kinematics from 0.7<z<2.7

We present the KMOS^3D survey, a new integral field survey of over 600 galaxies at 0.7<z<2.7 using KMOS at the Very Large Telescope (VLT). The KMOS^3D survey utilizes synergies with multi-wavelength ground and space-based surveys to trace the evolution of spatially-resolved kinematics and star formation from a homogeneous sample over 5 Gyrs of cosmic history. Targets, drawn from a mass-selected parent sample from the 3D-HST survey, cover the star formation-stellar mass ($M_*$) and rest-frame $(U-V)-M_*$ planes uniformly. We describe the selection of targets, the observations, and the data reduction. In the first year of data we detect Halpha emission in 191 $M_*=3\times10^{9}-7\times10^{11}$ Msun galaxies at z=0.7-1.1 and z=1.9-2.7. In the current sample 83% of the resolved galaxies are rotation-dominated, determined from a continuous velocity gradient and $v_{rot}/σ>1$, implying that the star-forming 'main sequence' (MS) is primarily composed of rotating galaxies at both redshift regimes. When considering additional stricter criteria, the Halpha kinematic maps indicate at least ~70% of the resolved galaxies are disk-like systems. Our high-quality KMOS data confirm the elevated velocity dispersions reported in previous IFS studies at z>0.7. For rotation-dominated disks, the average intrinsic velocity dispersion decreases by a factor of two from 50 km/s at z~2.3 to 25 km/s at z~0.9 while the rotational velocities at the two redshifts are comparable. Combined with existing results spanning z~0-3, disk velocity dispersions follow an approximate (1+z) evolution that is consistent with the dependence of velocity dispersion on gas fractions predicted by marginally-stable disk theory.

preprint2013arXiv

Depleted Galaxy Cores and Dynamical Black Hole Masses

Shallow cores in bright, massive galaxies are commonly thought to be the result of scouring of stars by mergers of binary supermassive black holes. Past investigations have suggested correlations between the central black hole mass and the stellar light or mass deficit in the core, using proxy measurements of $M_{\rm BH}$ or stellar mass-to-light ratios ($Υ$). Drawing on a wealth of dynamical models which provide both $M_{\rm BH}$ and $Υ$, we identify cores in 23 galaxies, of which 20 have direct, reliable measurements of $M_{\rm BH}$ and dynamical stellar mass-to-light ratios ($Υ_{\star,{\rm dyn}}$). These cores are identified and measured using Core-Sérsic model fits to surface brightness profiles which extend out to large radii (typically more than the effective radius of the galaxy); for approximately one fourth of the galaxies, the best fit includes an outer (\sersic) envelope component. We find that the core radius is most strongly correlated with the black hole mass and that it correlates better with total galaxy luminosity than it does with velocity dispersion. The strong core-size-- $M_{\rm BH}$ correlation enables estimation of black hole masses (in core galaxies) with an accuracy comparable to the $M_{\rm BH}$--$σ$ relation (rms scatter of 0.30 dex in $\log M_{\rm BH}$), without the need for spectroscopy. The light and mass deficits correlate more strongly with galaxy velocity dispersion than they do with black hole mass. Stellar mass deficits span a range of 0.2--39 \mbh, with almost all (87%) being $< 10 \, M_{\rm BH}$; the median value is 2.2 $M_{\rm BH}$.

preprint2013arXiv

The Influence of Dark Matter Halos on Dynamical Estimates of Black Hole Mass: Ten New Measurements for High-sigma Early-Type Galaxies

Adaptive-Optics assisted SINFONI observations of the central regions of ten early-type galaxies are presented. Based primarily on the SINFONI kinematics, ten black hole masses occupying the high-mass regime of the M_BH-sigma relation are derived using three-integral Schwarzschild models. The effect of dark matter inclusion on the black hole mass is explored. The omission of a dark matter halo in the model results in a higher stellar mass-to-light ratio, especially when extensive kinematic data are used in the model. However, when the diameter of the sphere of influence -- computed using the black hole mass derived without a dark halo -- is at least 10 times the PSF FWHM during the observations, it is safe to exclude a dark matter component in the dynamical modeling, i.e. the change in black hole mass is negligible. When the spatial resolution is marginal, restricting the mass-to-light ratio to the right value returns the correct M_BH although dark halo is not present in the model. Compared to the M_BH-sigma and M_BH-L relations of McConnell et al. (2011), the ten black holes are all more massive than expected from the luminosities and seven black hole masses are higher than expected from the stellar velocity dispersions of the host bulges. Using new fitted relations which include the ten galaxies, we find that the space density of the most massive black holes (M_BH >~ 10^9 Msun) estimated from the M_BH-L relation is higher than the estimate based on the M_BH-sigma relation and the latter is higher than model predictions based on quasar counts, each by about an order of magnitude.

preprint2009arXiv

The old and heavy bulge of M31 I. Kinematics and stellar populations

We present new optical long-slit data along 6 position angles of the bulge region of M31. We derive accurate stellar and gas kinematics reaching 5 arcmin from the center, where the disk light contribution is always less than 30%, and out to 8 arcmin along the major axis, where the disk makes 55% of the total light. We show that the velocity dispersions of McElroy (1983) are severely underestimated (by up to 50 km/s) and previous dynamical models have underestimated the stellar mass of M31's bulge by a factor 2. Moreover, the light-weighted velocity dispersion of the galaxy grows to 166 km/s, thus reducing the discrepancy between the predicted and measured mass of the black hole at the center of M31. The kinematic position angle varies with distance, pointing to triaxiality. We detect gas counterrotation near the bulge minor axis. We measure eight emission-corrected Lick indices. They are approximately constant on circles. We derive the age, metallicity and alpha-element overabundance profiles. Except for the region in the inner arcsecs of the galaxy, the bulge of M31 is uniformly old (>12 Gyr, with many best-fit ages at the model grid limit of 15 Gyr), slightly alpha-elements overabundant ([alpha/Fe]~0.2) and at solar metallicity, in agreement with studies of the resolved stellar components. The predicted u-g, g-r and r-i Sloan color profiles match reasonably well the dust-corrected observations. The stellar populations have approximately radially constant mass-to-light ratios (M/L_R ~ 4-4.5 for a Kroupa IMF), in agreement with stellar dynamical estimates based on our new velocity dispersions. In the inner arcsecs the luminosity-weighted age drops to 4-8 Gyr, while the metallicity increases to above 3 times the solar value.