Open Access
Issue
A&A
Volume 692, December 2024
Article Number A124
Number of page(s) 12
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202452548
Published online 05 December 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model.

Open access funding provided by Max Planck Society.

1. Introduction

Massive early-type galaxies (ETGs) are commonly found at the centre of galaxy clusters and are the result of mostly dissipationless mergers. During these events, nuclear supermassive black hole (SMBH) binaries are formed. Gravitational slingshots eject stars on radial orbits from the center of the remnant galaxy, destroying the power-law surface brightness distributions found in lower-luminosity ellipticals (Faber et al. 1997). Gravitational wave recoil (Khonji et al. 2024) can enhance the scouring mechanism. Through this core-scouring mechanism, the surface brightness profile I(r) of most massive ETGs becomes almost constant within a break (or core) radius rb, and for r ≤ rb one finds I(r)∝rγ, with γ < 0.3 (Faber et al. 1997). The break radius rb is tightly correlated with the mass of the central black hole (BH); moreover, it is anti-correlated with the central surface brightness (Mehrgan et al. 2019). A broader correlation between rb and the luminosity or stellar mass of the galaxies is also established (Laine et al. 2003; Rusli et al. 2013a). Moreover, within rb, the distribution of orbits becomes tangentially anisotropic (Thomas et al. 2014), as only stars avoiding the center can survive the scouring (Milosavljević & Merritt 2001; Thomas et al. 2014; Rantala et al. 2018, 2019). Although alternative explanations for the formation of cores have been proposed, such as the “tidal deposition” discussed by Nasim et al. (2021) and the feedback by active galactic nuclei (see Teyssier et al. 2011; Martizzi et al. 2012; Choi et al. 2018), they fail to explain this tangential anisotropy signature.

Black holes with dynamically measured masses of greater than 1010M, that is, hypermassive black holes (HMBHs), are still rare. They cannot be found using the MBH − σ relation (Saglia et al. 2016): a dissipationless merger of equal-mass galaxies doubles the mass of the resulting BH, maintains or even reduces the velocity dispersion of the system (Lauer et al. 2007; Naab et al. 2009). A large fraction of brightest cluster galaxies (BCGs) in the local Universe have relatively low velocity dispersions: Kluge & Bender (2023) measure on average a velocity dispersion of 250 km s−1 for their large sample of BCGs, with only 10% of objects having σ > 300 km s−1. The MBH − σ relation translates σ = 250 km s−1 into BH masses of around only about 6 × 108M. Nevertheless, HMBHs are found in BCGs, the largest (4 × 1010M) known being in Holm 15A (Mehrgan et al. 2019). The most promising way to search for HMBHs is to select massive ETGs, in particular BCGs, with core radii of the order of or larger than 0.6 kpc (Holm 15A has a core radius of 4 kpc). The Euclid Wide and Deep Surveys will allow us to find these objects in large numbers and out to redshifts of around 1 thanks to their excellent spatial resolution, large area coverage, and depth. Here we report the detection of the 1 . 29 $ 1{{\overset{\prime\prime}{.}}}29 $ (or 0.45 kpc) core of NGC 1272, which was measured on the Euclid Early Release Observations (EROs) Visible Camera (VIS) (Euclid Collaboration 2024a) image of the Perseus cluster (Cuillandre et al. 2024a,b), and the dynamical determination of the mass of its BH.

The galaxy is the second brightest elliptical galaxy of Perseus. With a total magnitude in the V band corrected for Galactic absorption VT0 of 11.27 (de Vaucouleurs et al. 1991), we compute a luminosity of L = 1.3 × 1011L using 72 Mpc as the distance of the cluster (Kluge et al. 2024), with which 1″ translates to 0.35 kpc. The stellar mass of the galaxy is 9 × 1011M using our dynamically determined mass-to-light ratio of 7 M/L (see Sect. 3) and the effective radius quoted in de Vaucouleurs et al. (1991) is 57″ or 20 kpc. With these properties, NGC 1272 belongs to the class of cD galaxies, the most massive ellipticals. We expect such objects to be triaxial with a low V/σ parameter (the ratio between the mean stellar velocity V and the velocity dispersion σ of the galaxy), to be detected as radio sources, and to have extended X-ray emission (Bender et al. 1989). Consistent with this, Veale et al. (2017) present the stellar kinematics of NGC 1272 obtained with the VIRUS-P spectrograph and classify the galaxy kinematically as a slow rotator. Park et al. (2017) detect a faint radio source at its center and McBride & McCourt (2014) study the properties of the double jets emerging from the center of the galaxy (which are bent with a curvature radius of 2 kpc). Arakawa et al. (2019) detect and study the X-ray minicorona of the galaxy, measuring a temperature of 0.63 keV and a size of 1.2 kpc.

The structure of the paper is as follows. In Sect. 2 we describe the photometric and spectroscopic observations of NGC 1272. The dynamical modeling is presented in Sect. 3. We draw our conclusions in Sect. 4, where we discuss the prospect of exploiting the Euclid survey (Euclid Collaboration 2024b) to find large cores up to redshift 1 in order to probe the formation redshift of the most massive BHs in galaxies.

2. Observations

NGC 1272 was observed during the early days of the Euclid survey as part of the pointings covering the Perseus galaxy cluster (Cuillandre et al. 2024b), one of the objects selected for the EROs program. The VIS (Euclid Collaboration 2024a) and Near-Infrared Spectrometer and Photometer (NISP) (Euclid Collaboration 2024c) ERO images of the cluster were reduced as described in Cuillandre et al. (2024a). Based on this dataset, studies of the Perseus intracluster light and intracluster globular clusters are described in Kluge et al. (2024) and a study of its dwarf galaxy population is presented by Marleau et al. (2024).

In Sect. 2.1 we make quantitative use of the VIS image of NGC 1272, with pixel size and resolution of 0 . 1 $ 0{{\overset{\prime\prime}{.}}}1 $ and 0 . 17 $ 0{{\overset{\prime\prime}{.}}}17 $, respectively. We used the near-infrared images (with 0 . 3 $ 0{{\overset{\prime\prime}{.}}}3 $ pixels) to assess the absence of dust in the central regions of the galaxy. The complementary spectroscopic information is described in Sect. 2.2.

2.1. Photometry

Figure 1 shows a cutout of the Euclid VIS image of NGC 1272. The isophote shape analysis was performed following Bender & Moellenhoff (1987). Figure 2 shows the resulting surface brightness profile calibrated to the V band for compatibility with the results of Rusli et al. (2013a). We perform the calibration by integrating the profile in circular apertures, which we shift to reproduce the aperture photometry listed in Hyperleda1. Finally, we adopt the correction for Galactic absorption and cosmological surface brightness dimming adopted in de Vaucouleurs et al. (1991) by matching our aperture magnitude within 57″ –which is the half-luminosity radius– to VT0 + 2.5log102 –which is half the total luminosity of the galaxy. We measure the photometry out to 147″ from the center, down to 24.9 mag arcsec−2; this is five times the distance reached by our stellar kinematics. The galaxy is round, with ellipticities smaller than 0.15 and isophotes showing only small deviations from perfect ellipses. For radii of greater than about 45″, the center of the isophotes starts to drift toward the direction of NGC 1275 and the position angle twists by 70°. In the inner 1 . 2 $ 1{{\overset{\prime\prime}{.}}}2 $, the surface brightness increase toward the center slows down, pointing to the presence of a core.

thumbnail Fig. 1.

Cutout of the Euclid VIS image of NGC 1272.

The sizes of the cores of ETGs have been determined in the past using the Nuker law (Faber et al. 1997) and the core-Sérsic law (Graham et al. 2003). The advantages and disadvantages of the two approaches have been discussed at length in the literature and depend on how well the outer parts of a galaxy can be described by either law. In the following, we rely on both approaches as a way to estimate the systematic effects that affect our measurements.

We begin the process of deriving the size of the core of NGC 1272 by fitting a 5000 × 5000 pixel (250″ × 250″) image extracted from the VIS mosaic with the PSF-convolved core-Sérsic function provided by the Imfit code2 of Erwin (2015), using the image of a star extracted in the vicinity of the galaxy as the PSF. Fitting larger cutouts requires prohibitively large computing time without improving the determination of the core size. As implemented in Erwin (2015), the core-Sérsic function is

I CS ( r ) = I [ 1 + ( r b r ) α CS ] γ CS / α CS exp [ b n ( r α CS + r b α CS r e α CS ) 1 / n α CS ] , $$ \begin{aligned} I_{\rm CS}(r) = I^{\prime }\left[1+\left(\frac{r_{\rm b}}{r}\right)^{\alpha _{\rm CS}}\right]^{\gamma _{\rm CS}/\alpha _{\rm CS}} \exp \left[-b_n\left( \frac{r^{\alpha _{\rm CS}}+r_{\rm b}^{\alpha _{\rm CS}}}{r_{\rm e}^{\alpha _{\rm CS}}} \right)^{1/n\alpha _{\rm CS}}\right], \end{aligned} $$(1)

where

I = I CS , b 2 γ CS / α CS exp [ b n ( 2 1 / α CS r b r e ) 1 / n ] , $$ \begin{aligned} I^{\prime }=I_{\rm CS,b}2^{-\gamma _{\rm CS}/\alpha _{\rm CS}} \exp \left[b_n\left(2^{1/{\alpha _{\rm CS}}}\frac{r_{\rm b}}{r_{\rm e}}\right)^{1/n}\right], \end{aligned} $$(2)

and bn ∼ 2n − 1/3 + 4/405n. Similarly, the Sérsic function is

I S ( r ) = I e exp [ b n ( r r e ) 1 / n ] . $$ \begin{aligned} I_{\rm S}(r) = I_{\rm e}\exp \left[-b_n \left(\frac{r}{r_{\rm e}}\right)^{1/n}\right]. \end{aligned} $$(3)

Here, re is the half-luminosity radius, Ie the intensity at re, n the Sérsic index, rb the break radius, and ICS, b is the intensity at rb. Also, −γCS is the slope of the power-law inner profile, and αCS specifies the sharpness of the transition to the outer, Sérsic profile. In Tables 1 and 2, we provide the values of μV(re) and μV(rb) that calibrate the surface brightness profiles μS = −2.5 log10IS/Ie + μV(re) and μCS = −2.5 log10ICS/ICS, b + μV(rb) to the V band.

Table 1.

Parameters of the core-Sérsic best fits.

Table 2.

Parameters of the second Sérsic component.

The core-Sérsic model reproduces the surface brightness of the galaxy accurately, with residuals of less than 0.1 mag, even if it has a constant ellipticity and position angle. The resulting parameters of the fit are given in Table 1; in particular, the size of the core is perfectly resolved by the spatial resolution of the VIS image. According to Thomas et al. (2016), we expect this to match the size of the sphere of influence of the central black hole of the galaxy. The best fitting value of n (21.1) is unrealistically large, as is that of re, which is two orders of magnitudes larger than the size of the fitted image; this stems from the almost power-law behavior of the outer profile, which is typical of BCGs (Kluge & Bender 2023). Both parameters are to be considered as a convenient parametrization of the galaxy profile out to the limit of the image and increasing with the image size. More importantly, the values of rb and γCS show no significant dependency on this choice. The statistical errors listed in Table 1 (and further below in Tables 2 and 3) are minute, because there are a great many independent points in the image. We have rounded them up to the first or second digit. Fitting the surface brightness of Fig. 2 with a one-dimensional Sérsic profile without PSF convolution delivers similar results within the systematic errors estimated below.

thumbnail Fig. 2.

Photometry profiles of NGC 1272. Top Left: Surface brightness profile of NGC 1272 measured from the Euclid VIS image, calibrated to the V band, and corrected for Galactic absorption and cosmological dimming, as a function of the 1/4 power of the semi-major distance A on the sky in arcsec. Bottom Left: Difference ΔSB between the surface brightness profile of NGC 1272 and the surface brightness of the core-Sérsic+Sérsic model. Right: Ellipticity 1 − B/A, where B is the semi-minor axis length on the sky (top) and PA (bottom) as a function of the 1/4 power of A. The solid red lines show the core-Sérsic+Sérsic model. The dashed lines show its core radius. The blue line shows the ellipticity profile of the axisymmetric deprojection. The green line shows the PA profile of the triaxial deprojection.

Table 3.

Parameters of the Nuker best fit.

We further explore the systematic errors affecting the estimation of the core radius by fitting the same image with the Imfit implementation of the (PSF-convolved) Nuker-law:

I N ( r ) = I N , b 2 ( β N γ N ) / α N ( r b r ) γ N [ 1 + ( r r b ) α N ] ( γ N β N ) / α N . $$ \begin{aligned} I_{\rm N}(r) = I_{\rm N,b} 2^{(\beta _{\rm N}-\gamma _{\rm N})/\alpha _{\rm N}} \left(\frac{r_{\rm b}}{r}\right)^{\gamma _{\rm N}}\left[1+\left(\frac{r}{r_{\rm b}}\right)^{\alpha _{\rm N}}\right]^{(\gamma _{\rm N}-\beta _{\rm N})/\alpha _{\rm N}}. \end{aligned} $$(4)

Here, −γN is the asymptotic logarithmic slope inside rb, −βN is the asymptotic outer slope, and the αN parameter describes the sharpness of the break; IN, b is the intensity at rb and the surface brightness profile μN(r) =  − 2.5log10IN/IN, b + μV(rb) is calibrated to the V band through the value of μV(rb) given in Table 3. The Nuker fit delivers a core-size determination similar to that found using the core-Sérsic function (see Table 3) and, as noted above, describes the outer power-law behavior of the galaxy profile reasonably well.

Inspection of the NIR images (see Fig. 3) confirms that the central region of NGC 1272 is not strongly affected by dust: core sizes of between 1 . 25 $ 1{{\overset{\prime\prime}{.}}}25 $ and 1 . 29 $ 1{{\overset{\prime\prime}{.}}}29 $ are obtained when fitting these images. The slope of the surface brightness profile inside rb is between 0.1 (from the Nuker fit) and 0.2 (from the core-Sérsic fit), which is within the range expected for core ellipticals (Faber et al. 1997).

thumbnail Fig. 3.

The IE, YE, JE, and HE cutouts of the inner 6″ × 6″ of NGC 1272.

A more realistic estimate of the effective radius of the galaxy that better catches the varying ellipticity and Position Angle (PA) profiles (see Fig. 2) is obtained by fitting a two-component model: an inner core-Sérsic plus an outer Sérsic profile. The results are listed in Tables 1 and 2. The core radius is somewhat larger and the inner slope γCS of the profile somewhat shallower than above. Fig. 4 shows the fractional residuals between the image and the model, the absolute values of which are always smaller than 0.1. The surface brightness is reproduced with a root mean square (RMS) of 0.028 mag (see Fig. 2).

thumbnail Fig. 4.

Percentage residuals after subtraction of the core-Sérsic + Sérsic model.

de Rijcke et al. (2009) collected F555W and F814W ACS images of NGC 1272 with the Hubble Space Telescope (HST), obtaining a factor two better pixel size and resolution than our Euclid VIS images. We performed core-Sérsic fits to 179″ × 184″ images with Imfit, fixing the value of n to the result obtained from the VIS image. The results are given in Table 1, where we calibrate μb = −2.5log Ib to the V band as above. We measured r b = 1 . 27 $ r_{\mathrm{b}} = 1{{\overset{\prime\prime}{.}}}27 $ and 1 . 25 $ 1{{\overset{\prime\prime}{.}}}25 $ in the two bands, demonstrating that possible color gradients do not strongly affect the determination of rb.

The statistical errors reported in Tables 1, 2, and 3 are very small due to the large number of pixels fitted. More significant are the systematic errors that come from the different fitting functions used to measure rb. Averaging the five estimates of rb presented above, we obtain 1 . 29 $ 1{{\overset{\prime\prime}{.}}}29 $, or 0.45 kpc, with a RMS of 0 . 07 $ 0{{\overset{\prime\prime}{.}}}07 $, which we adopt as our measurement error.

We deproject the surface brightness profile using the axisymmetric deprojection code of Magorrian (1999), assuming that the galaxy is edge-on, as is usually done in such cases (Lipka & Thomas 2021). Other options are explored below, when triaxial deprojections are considered. The blue line in Fig. 2 shows that this deprojection reproduces the ellipticity profiles, but cannot reproduce the PA profile (assumed to be constant in axisymmetric deprojections). The intrinsic flattening profile q(r) –where q = c/a and a and c are the major and minor semi-axes of the galaxy– derived in this way is around 0.9 (see Fig. 5).

thumbnail Fig. 5.

Profiles of p(r) and q(r) (where p = b/a, q = c/a, and a, b, and c are the major, intermediate, and minor semi-axis of the galaxy) derived from the axisymmetric (dotted), triaxial (full line), spherical (dashed line), and prolate (dashed-dotted line) deprojections of the galaxy, as a function of distance r from the center. The red and blue shaded areas show the whole range of allowed deprojections with RMS ≤ 1.2 × RMSmin (de Nicola et al. 2020, 2022a,b).

We also explore the range of possible triaxial deprojections following de Nicola et al. (2020). The reconstructed p and q profiles (where p = b/a and b is the intermediate semi-axis of the galaxy) are shown in Fig. 5 and demonstrate that the galaxy is almost spherical and close to axisymmetric, with p ≈ 1 and q ≈ 0.9. For some viewing angles, the strong PA radial variation (see Fig. 5) forces a twist of the principal axis with radius, which explains why the p and q profiles can become larger than one at the same distance (de Nicola et al. 2020).

The deprojection with the lowest RMS in surface brightness is obtained at angles of (θ, ϕ, ψ) = (64 ° ,124 ° ,23 ° ), roughly 26° above the equatorial plane and about 34° away from the intermediate axis. However, it is clear that reconstructing the true orientation of the galaxy is almost impossible given its almost spherical geometry. The green line in Fig. 2 shows that this deprojection indeed reproduces the PA profile (in addition to the ellipticity profile). We further explore alternative deprojections with comparably good surface brightness RMS: a (mildly) prolate and an (almost) spherical deprojection (similar to the axisymmetric one, but matching the PA twist). Both are obtained assuming that the line of sight is along the major axis of the galaxy and starting the deprojection routine with constant profiles q(r) = 0.7 and q(r) = 0.95 in the prolate and spherical cases, respectively.

2.2. Spectroscopic observations and kinematics

We observed NGC 1272 spectroscopically with the Visible Integral-field Replicable Unit Spectrograph (VIRUS) at the Hobby-Eberly Telescope (HET) on 3 March 2022. The pointing of the telescope was optimized to observe NGC 1275; as a result, the integral field unit (IFU) covering NGC 1272 was slightly off-center and did not uniformly cover the galaxy. The seeing reported during the observations was FWHM = 2 . 36 $ \mathrm{FWHM} = 2{{\overset{\prime\prime}{.}}}36 $. The diameter of the single fiber is 1 . 5 $ 1{{\overset{\prime\prime}{.}}}5 $. Given the size of the core measured above (a diameter of 2 . 6 $ 2{{\overset{\prime\prime}{.}}}6 $), the spatial resolution of this data set is marginally sufficient to resolve the sphere of influence of the central BH of the galaxy. Rusli et al. (2013b) find that in such cases, an unbiased recovery of the BH mass is possible if the dark matter (DM) halo of the galaxy is taken into account in the dynamical modeling, as done here; see below. The data cover a wavelength interval ranging from 3470 Å to 5540 Å with a spectral resolution of 5.6 Å.

We used the Voronoi tessellation method of Cappellari & Copin (2003) to spatially bin the spectral data for a target average signal-to-noise ratio (S/N) of 40. With this target S/N, spectra were only binned together starting approximately 3″ from the center of the galaxy, thus maximizing the spatial resolution of our data within the core region. This resulted in a total of 110 spatial bins. We measured the stellar kinematics using WINGFIT (Thomas, in prep.), which delivers optimally smoothed non-parametric line-of-sight velocity distributions (LOSVDs) using the model optimization approach of Thomas & Lipka (2022). The stellar kinematic fits were performed using the MILES library (Sánchez-Blázquez et al. 2006) of stellar templates. Following the strategy laid out in Mehrgan et al. (2023), we performed a careful pre-selection of templates in order to minimize distortions of the LOSVDs due to template mismatch. To this end, we fitted the average spectrum of the central 2″ of the galaxy using all the templates of the MILES library with a Gauss-Hermite LOSVD that was fixed to be symmetric around a line-of-sight velocity of zero. We selected the set of 18 templates that in the best fit received a non-zero weight. Using this set, we fitted all bins of the galaxy with non-parametric LOSVDs without the symmetry costraint. Also following Mehrgan et al. (2023), we used no additive polynomials in the fit and only a minimal third-order multiplicative polynomial. Fits were performed in the wavelength interval between 4700 and 5400 Å. The resulting two-dimensional kinematic maps are shown in Fig. 6; the corresponding radial profiles can be seen in Fig. 7. We measured the stellar kinematics out to a maximum distance of 38″, or 0.66 times the effective radius quoted by de Vaucouleurs et al. (1991). The galaxy has small mean rotation (at most v = 20 km s−1), a relatively low velocity dispersion σ of around 250 km s−1, increasing to 270 km s−1 toward the center (admittedly with only one point within 1 arcsec of the center), an almost zero third-order Hermite parameter h3, and a zero fourth-order Hermite parameter h4, which decreases to about −0.05 toward the center. The data presented by Veale et al. (2017) match these findings, though the h4 obtained by these authors is always approximately zero. We averaged (v2 + σ2)0.5, with equal or luminosity weights, to get an estimate of σe, that is, the velocity dispersion within the half-luminosity radius (even if our stellar kinematics only reach out to two-thirds of re; see above). We find σe = 247 ± 3 km s−1, which we adopt in Sect. 4.

thumbnail Fig. 6.

Two-dimensional stellar kinematics of NGC 1272. The horizontal and vertical dashed lines show the major and minor axes of the galaxy, respectively. North is up and east is to the left.

thumbnail Fig. 7.

Radial stellar kinematics of NGC 1272. The red lines show the axisymmetric fit to the stellar kinematics of the galaxy (black data points with error bars) as a function of the distance R from the center of the galaxy on the sky.

For the subsequent dynamical analysis, we sampled and fit the non-parametric LOSVDs between ±1400 km s−1 with Nvel = 25 velocity bins, and no Hermite parameters. An example LOSVD measured near the center of the galaxy is shown in Fig. 8, where the red line connects the values produced by the best-fitting base model. Finally, to ensure that we are able to estimate the uncertainties of our dynamical models, we split our data into four quadrants (indicated by q1, q2, q3, and q4 in Fig. 6) along the minor and major axes of the galaxy for the axisymmetric dynamical models, and into two halves (quadrants q1/q4, northern, and q2/q3, southern) split by the major axis for the triaxial analysis. By modeling quadrants or halves independently, we can estimate the uncertainties of our best-fit modeling parameters from the scatter between them. However, the orientation and positioning of the IFU give us a much better coverage of the q3 and q4 quadrants, or the east side of the galaxy. Therefore, we expect the most reliable dynamical constraints to come from these regions.

thumbnail Fig. 8.

Line-of-sight velocity distribution measured at 4 arcsec from the center of NGC 1272 (filled circles with error bars). The red line connects the values provided by the base model.

3. Dynamical modeling

Given the results presented in Fig. 5 (NGC 1272 is almost spherical and axisymmetric, but triaxiality and a prolate shape cannot be excluded), we construct both axisymmetric and triaxial Schwarzschild models of the galaxy. The axisymmetric modeling is similar to that of Mehrgan et al. (2024), with the following modifications. We fit the four quadrants both independently and together, determining the mass of the central BH MBH and the stellar mass-to-light ratio Υ* (i.e., no radial variations of Υ* are considered because of the rather coarse and sparse sampling of our stellar kinematics). We use a spherical Zhao (1996) halo with α = 1 and β = 3, defined by ρ10; the DM density at 10 kpc; rs, the scale radius of the halo, which is allowed to vary up to the largest distance probed by our kinematics (Lipka et al. 2024); and γDM ≥ 0, the inner slope of the DM halo:

ρ DM ( r ) = k ( r / r s ) γ DM ( 1 + r / r s ) 3 γ DM , $$ \begin{aligned} \rho _{\rm DM}(r) = \frac{k}{\left(r/r_{\rm s}\right)^{\gamma _{\rm DM}} {\left(1+r/r_{\rm s}\right)^{3-\gamma _{\rm DM}}}}, \end{aligned} $$(5)

and k = ρ10(10 kpc/rs)γDM(1+10 kpc/rs)3 − γDM.

The triaxial modeling follows de Nicola et al. (2024) and uses the Schwarzschild code SMART (Neureiter et al. 2021) to determine MBH and Υ*, considering a DM halo that is triaxial, described by its shape parameters pDM and qDM plus ρ10 and γDM, fixing rs to a large value (158 kpc). We model the northern (quadrants q1 and q4 in Fig. 6) and southern (quadrants q2 and q3 in Fig. 6) halves of the galaxy separately to assess the systematic uncertainties.

In both the axisymmetric and triaxial cases, we maximize the quantity S ̂ = S α ̂ χ 2 $ \hat{S}=S-\hat{\alpha}\,\chi^2 $ to determine the orbital weights. Here χ2 is calculated from the model fit to the observed non-parametric LOSVDs, and S is the Boltzmann entropy (Thomas et al. 2004). The deprojected light distributions are used as a constraint and the parameter α ̂ $ \hat{\alpha} $ is the smoothing of the models, which is determined following the prescriptions of Lipka & Thomas (2021) and Thomas & Lipka (2022), which involve the determination of the effective degrees of freedom meff. The parameters MBH, Υ*, ρ10, and rs (in both the axisymmetric and triaxial cases), plus pDM and qDM in the triaxial case, are determined by minimizing the generalized Akaike information criterion AICp = χ2 + 2 meff over a grid of α ̂ $ \hat{\alpha} $ values.

The resulting axisymmetric best fits to the kinematics are shown in Fig. 7; the derived parameters for the different fit types are listed in Table 4. Our base result is the axisymmetric model of the entire stellar kinematic data set; it fits the kinematics very well (see red line in Fig. 7), delivering a reduced χ2 of χ2/(Ndata − meff) = 0.91. Figure 9 shows MBH, Υ*, V, ρ10, and γDM as a function of the quality of the fit measured by the AICp value. Every parameter is well constrained, with small statistical errors (so small that we do not quote them in Table 4). In particular, we detect a BH of 5 × 109M, a mass-to-light ratio Υ*, V of 7.1 M/L, in between the values derived from our stellar population analysis for the Kroupa or the Salpeter initial mass function (IMF, see below), a DM density at 10 kpc (approximately 30″) similar to the values reported for other massive elliptical galaxies (Mehrgan et al. 2024), and a cored DM density profile. We compute the radius rSOI of the sphere of influence of the BH as the distance from the center where the total mass (stellar plus DM without a BH) is equal to MBH (Thomas et al. 2016). Our rSOI matches the value of rb and is larger than half the FWHM of the seeing of the spectroscopic observations. This, together with the modeling of the DM halo of the galaxy, allows an unbiased estimate of the BH mass (Rusli et al. 2013b).

thumbnail Fig. 9.

Results of the axisymmetric modeling of NGC 1272 of the complete stellar kinematic dataset. As a function of the quality of the fits measured by the AICp we show: from left to right, and from top to bottom: the BH mass MBH; the dynamical V-band mass-to-light ratio Υ*, V; the DM density at 10 kpc ρ10; and the inner slope of the DM density profile γDM. The gray points show the individual models, the red dots show the best-fitting model, and the black lines show the lower envelope of the distributions of gray points.

Table 4.

Parameters of the axisymmetric and triaxial dynamical modeling.

We gauge our (systematic) errors by looking first at the axisymmetric modeling of the two quadrants covering larger parts of the galaxy, q3 and q4. Here the BH mass can be as low as 1.8 × 109M and Υ*, V as large as 8.9 M/L, with slightly larger DM densities. Further insights into our systematic errors are gained from the triaxial SMART modeling. All models fit the kinematic data well, with χ2/(Ndata − meff) between 0.6 and 1.0 The best-fitting triaxial model delivers MBH = (5.9 ± 1.7)×109M, averaging over the two halves of the galaxy; the smallest and largest values for the BH mass are obtained when fitting the southern half of the galaxy (where the kinematic coverage is relatively sparse) in the prolate and spherical cases, respectively. The dynamical stellar mass-to-light ratio Υ*, V ranges from 4 to 7.7 M/L. The density of the DM halo (log10ρ10/[M kpc−3] = 7.3 ± 0.3, averaging over all triaxial models and halves) agrees with the axisymmetric result. The slope γDM of the DM density profile is smaller than 1, but possibly steeper than the cored halo determined axisymmetrically. The errors quoted in Table 4 for our base result are the RMSs of the nine listed best-fitting models. Finally, the best-fitting shape of the DM halo is spherical (pDM = qDM = 1), with only the prolate case delivering pDM = 0.9.

Figure 10 summarizes the spherically averaged stellar and total mass profiles derived by the dynamical models we considered. The total mass distribution is robustly determined in the region probed by the measured kinematics, with only small deviations between the different models. Inside the sphere of influence (roughly the size of the core), the differences between the profiles reflect the observed scatter in the BH mass. The stellar mass profiles scale according to the derived Υ*, V values. The DM mass is equal to the stellar mass at approximately the outermost radius probed by our stellar kinematics. Using the best-fit value of Υ*, V = 7.1 M/L, we estimate the total stellar mass of the galaxy from the total luminosity L = 1.3 × 1011L quoted in Sect. 1 to be 9 × 1011M, which we use in Sect. 4.

thumbnail Fig. 10.

Spherically averaged mass profiles of NGC 1272. The solid and dashed lines show the total and stellar profiles, respectively. The triaxial, prolate, and spherical models are averaged over the two sides. The vertical dashed line shows the position of the core radius.

The anisotropy β profile (where β = 1 − σT2/σR2, and σT and σR are the spherical tangential and radial velocity dispersions, respectively) is not particularly well constrained, but displays the typical feature of core ellipticals (Thomas et al. 2014). Figure 11 shows that the β profile of the base model becomes tangentially anisotropic within the core radius (the result of core scouring) and more isotropic in the outer part, similarly to the spherical model. The triaxial model is overall mildly tangentially anisotropic, while the prolate model is radially anisotropic outside the core.

thumbnail Fig. 11.

Anisotropy profiles β (where β = 1 − σT2/σR2, and σT and σR are the spherical tangential and radial velocity dispersions, respectively) of the different models. The triaxial, prolate, and spherical models are averaged over the two sides. The distances to the center are in units of the core radius.

Measuring Lick indices and fitting them with the simple stellar population models of Thomas et al. (2003) and Maraston (2005), we find that the best fit has a simple stellar population as old as the Universe, a metallicity of slightly above solar, and is more than a factor of two overabundant in α-elements. The derived V-band mass-to-light ratio is 6 M/L with a Kroupa IMF and 8 M/L with a Salpeter IMF. This matches the dynamically determined Υ*, V, without unambiguously preferring one of the two options.

Finally, we estimate the mass that has been expelled from the core during its formation. We consider the core-Sérsic solution obtained with the second Sérsic component (see second line of Table 1) and consider the Sérsic function with n = 12.1, r e = 34 . 9 $ r_{\mathrm{e}} = 34{{\overset{\prime\prime}{.}}}9 $, and μV(re) = 23.57, which reproduces the core-Sérsic solution outside the core region. We integrate the luminosity difference, or the luminosity deficit Ldef, between the two functions out to 8″, finding Ldef = 2.7 × 109L. Using the dynamically determined Υ*, V = 7 M/L, this translates into a mass deficit of Mdef = 1.9 × 1010M, or 3.8 MBH, in the range found by Rusli et al. (2013a). According to the simulations of Gualandris & Merritt (2008), mass deficits of up to 5 × MBH can result from single dry mergers.

4. Conclusions

We present a measurement of the size (0.45 kpc) of the core of NGC 1272 based on the VIS image of the Perseus cluster taken as part of the Euclid ERO campaign. The dynamical modeling of the stellar kinematics collected with the VIRUS spectrograph at the HET allows us to measure the mass (5 ± 3)×109M of the BH at the center of the galaxy. While in line with expectations from the MBH − rb correlation of Thomas et al. (2016), the central surface brightness versus MBH correlation of Mehrgan et al. (2019), and the MBH − M* relation of Saglia et al. (2016), the BH mass of NGC 1272 is a factor of 8 larger than predicted by the MBH − σ relation of Saglia et al. (2016), or 1.8 times the 1σ error combined with the intrinsic scatter in the relation (see Fig. 12). This corroborates the conclusion that the velocity dispersion is not the best indicator of the BH mass for core galaxies with stellar masses of the order of or larger than 1012M: five out of the six galaxies with such a stellar mass in Fig. 12 have BH masses larger than predicted by the MBH − σ relation. Therefore, the most efficient and rapid method to search for galaxies harboring the most massive BHs is to look for passive objects with large cores and low central surface brightness. In the local Universe, a galaxy with a core size of 1 kpc contains a BH with a mass of 1010M.

thumbnail Fig. 12.

Position of NGC 1272 (in red) on the MBH − rb (left), the MBH − M* (middle), and the MBH − σ (right) relation. The black data points are from Rusli et al. (2013a), Saglia et al. (2016), Thomas et al. (2016), Mehrgan et al. (2019), Neureiter et al. (2023), and de Nicola et al. (2024). The blue data points are core ellipticals with a stellar mass of greater than 1012M. While the galaxy follows the MBH − rb given by Thomas et al. (2016) and the MBH − M* relation of Saglia et al. (2016) for the sample of CorePowerE, it deviates by a factor of 8.4 from the MBH − σ relation of Saglia et al. (2016), or by a factor of 1.8 of the 1σ error combined with the intrinsic scatter in the relation (shown by the dotted lines).

The Euclid VIS images in the IE band deliver a PSF with FWHM ≈ 0 . 17 $ 0{{\overset{\prime\prime}{.}}}17 $ and a pixel size of 0 . 1 $ 0{{\overset{\prime\prime}{.}}}1 $, with a depth of 24.5 mag in the Wide Survey (at 10σ for extended sources) and 2 mag deeper in the Deep Survey. Near-infrared YE, JE, and HE images provide photometry with 0 . 3 $ 0{{\overset{\prime\prime}{.}}}3 $ pixels. Combined with ground-based images, the surveys will deliver not only photometric redshifts for each detected source, but also physical parameters, such as stellar masses and sizes. At the end of the mission, the Wide Survey will cover about 14 000 deg2 of extragalactic sky, along with 50 deg2 at the Deep Survey. This unprecedented dataset will allow us to search for galaxies with cores that are larger than 2 kpc out to a redshift of 1 (where they will subtend an angle of 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $ on the sky) as a function of stellar mass. We plan to establish the redshift up to which the correlation between rb and total stellar mass exists (see Fig. 13) and study its possible evolution with a large statistical sample, indirectly probing the possible coevolution of BHs and galaxy properties at the highest BH mass end. For example, the blue cross in Fig. 13 shows the position of the BCG of the EDISCS cluster CL1216 at a redshift of 0.8 (Saglia et al. 2010). We measured the size of its core in the available HST images, deriving 1.5 kpc (or 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $) from a core-Sérsic fit, and 2.21 kpc (or 0 . 7 $ 0{{\overset{\prime\prime}{.}}}7 $) from a Nuker fit. Such a core will be measurable in the VIS mosaics of the Wide survey. With a stellar mass of log10M*/M = 11.82, the BCG appears to have a larger core than local core ellipticals of similar mass. Using the local rb − MBH relation, we estimate that an HMBH with a mass of greater than 1010M could be already in place at such a high redshift in this galaxy. Spectroscopic follow-up (possible at the Extremely Large Telescope) of selected galaxies with similarly large and bright cores will deliver the dynamical mass confirmation.

thumbnail Fig. 13.

Correlation between the core radius rb and stellar mass M*. The black datapoints are from Rusli et al. (2013a), Saglia et al. (2016), Thomas et al. (2016), Mehrgan et al. (2019), Neureiter et al. (2023), and de Nicola et al. (2024). NGC 1272 is shown in red. The blue cross shows the position of the BCG on the EDISCS cluster CL1216 (Saglia et al. 2010). The black solid line shows log10(M*/M) = 0.64log10(rb/kpc)+12.1; the dotted lines show the ±1σ scatter (0.26 dex) in the relation.


Acknowledgments

This work has made use of the Early Release Observations (ERO) data from the Euclid mission of the European Space Agency (ESA), 2024, https://doi.org/10.57780/esa-qmocze3. The Euclid Consortium acknowledges the European Space Agency and a number of agencies and institutes that have supported the development of Euclid, in particular the Agenzia Spaziale Italiana, the Austrian Forschungsförderungsgesellschaft funded through BMK, the Belgian Science Policy, the Canadian Euclid Consortium, the Deutsches Zentrum für Luft- und Raumfahrt, the DTU Space and the Niels Bohr Institute in Denmark, the French Centre National d’Etudes Spatiales, the Fundação para a Ciência e a Tecnologia, the Hungarian Academy of Sciences, the Ministerio de Ciencia, Innovación y Universidades, the National Aeronautics and Space Administration, the National Astronomical Observatory of Japan, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Research Council of Finland, the Romanian Space Agency, the State Secretariat for Education, Research, and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency. A complete and detailed list is available on the Euclid website (https://www.euclid-ec.org). RS, RB and MF acknowledge support by the Deutsches Zentrum für Luft- und Raumfahrt (DLR) grant 50 QE 1101. RS, RB, ML thank the Hobby Eberly Telescope (HET) project for allocating the observations and the technical support. The Hobby-Eberly Telescope is a joint project of the University of Texas at Austin, the Pennsylvania State University, Ludwig-Maximilians-Universität München, and Georg-August Universität Gottingen. The HET is named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly. The HET Collaboration acknowledges the support and resources from the Texas Advanced Computing Center. We thank the Resident Astronomers and Telescope Operators at the HET for the skillful execution of our observations with VIRUS. We would like to acknowledge that the HET is built on Indigenous land. Moreover, we would like to acknowledge and pay our respects to the Carrizo & Comecrudo, Coahuiltecan, Caddo, Tonkawa, Comanche, Lipan Apache, Alabama-Coushatta, Kickapoo, Tigua Pueblo, and all the American Indian and Indigenous Peoples and communities who have been or have become a part of these lands and territories in Texas, here on Turtle Island.

References

  1. Arakawa, N., Fabian, A. C., & Walker, S. A. 2019, MNRAS, 488, 894 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bender, R., & Moellenhoff, C. 1987, A&A, 177, 71 [NASA ADS] [Google Scholar]
  3. Bender, R., Surma, P., Doebereiner, S., Moellenhoff, C., & Madejsky, R. 1989, A&A, 217, 35 [Google Scholar]
  4. Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 [Google Scholar]
  5. Choi, E., Somerville, R. S., Ostriker, J. P., Naab, T., & Hirschmann, M. 2018, ApJ, 866, 91 [NASA ADS] [CrossRef] [Google Scholar]
  6. Cuillandre, J.-C., Bertin, E., Bolzonella, M., et al. 2024a, A&A, submitted [arXiv:2405.13496] [Google Scholar]
  7. Cuillandre, J.-C., Bolzonella, M., Boselli, A., et al. 2024b, A&A, submitted [arXiv:2405.13501] [Google Scholar]
  8. de Nicola, S., Saglia, R. P., Thomas, J., Dehnen, W., & Bender, R. 2020, MNRAS, 496, 3076 [CrossRef] [Google Scholar]
  9. de Nicola, S., Saglia, R. P., Thomas, J., et al. 2022a, ApJ, 933, 215 [NASA ADS] [CrossRef] [Google Scholar]
  10. de Nicola, S., Neureiter, B., Thomas, J., Saglia, R. P., & Bender, R. 2022b, MNRAS, 517, 3445 [NASA ADS] [CrossRef] [Google Scholar]
  11. de Nicola, S., Thomas, J., Saglia, R. P., et al. 2024, MNRAS, 530, 1035 [NASA ADS] [CrossRef] [Google Scholar]
  12. de Rijcke, S., Penny, S. J., Conselice, C. J., Valcke, S., & Held, E. V. 2009, MNRAS, 393, 798 [NASA ADS] [CrossRef] [Google Scholar]
  13. de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies (New York: Springer) [Google Scholar]
  14. Erwin, P. 2015, ApJ, 799, 226 [Google Scholar]
  15. Euclid Collaboration (Cropper, M. S., et al.) 2024a, A&A, in press, https://doi.org/10.1051/0004-6361/202450996 [Google Scholar]
  16. Euclid Collaboration (Mellier, Y., et al.) 2024b, A&A, in press, https://doi.org/10.1051/0004-6361/202450810 [Google Scholar]
  17. Euclid Collaboration (Jahnke, K., et al.) 2024c, A&A, in press, https://doi.org/10.1051/0004-6361/202450786 [Google Scholar]
  18. Faber, S. M., Tremaine, S., Ajhar, E. A., et al. 1997, AJ, 114, 1771 [NASA ADS] [CrossRef] [Google Scholar]
  19. Graham, A. W., Erwin, P., Trujillo, I., & Asensio Ramos, A. 2003, AJ, 125, 2951 [CrossRef] [Google Scholar]
  20. Gualandris, A., & Merritt, D. 2008, ApJ, 678, 780 [NASA ADS] [CrossRef] [Google Scholar]
  21. Khonji, N., Gualandris, A., Read, J. I., & Dehnen, W. 2024, ApJ, 974, 204 [NASA ADS] [CrossRef] [Google Scholar]
  22. Kluge, M., & Bender, R. 2023, ApJS, 267, 41 [NASA ADS] [CrossRef] [Google Scholar]
  23. Kluge, M., Hatch, N. A., Montes, M., et al. 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202450772 [Google Scholar]
  24. Laine, S., van der Marel, R. P., Lauer, T. R., et al. 2003, AJ, 125, 478 [Google Scholar]
  25. Lauer, T. R., Faber, S. M., Richstone, D., et al. 2007, ApJ, 662, 808 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lipka, M., & Thomas, J. 2021, MNRAS, 504, 4599 [NASA ADS] [CrossRef] [Google Scholar]
  27. Lipka, M., Thomas, J., Saglia, R., et al. 2024, ApJ, 976, 17 [NASA ADS] [CrossRef] [Google Scholar]
  28. Magorrian, J. 1999, MNRAS, 302, 530 [NASA ADS] [CrossRef] [Google Scholar]
  29. Maraston, C. 2005, MNRAS, 362, 799 [NASA ADS] [CrossRef] [Google Scholar]
  30. Marleau, F. R., Cuillandre, J.-C., Cantiello, M., et al. 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202450799 [Google Scholar]
  31. Martizzi, D., Teyssier, R., & Moore, B. 2012, MNRAS, 420, 2859 [CrossRef] [Google Scholar]
  32. McBride, J., & McCourt, M. 2014, MNRAS, 442, 838 [NASA ADS] [CrossRef] [Google Scholar]
  33. Mehrgan, K., Thomas, J., Saglia, R., et al. 2019, ApJ, 887, 195 [NASA ADS] [CrossRef] [Google Scholar]
  34. Mehrgan, K., Thomas, J., Saglia, R., Parikh, T., & Bender, R. 2023, ApJ, 948, 79 [NASA ADS] [CrossRef] [Google Scholar]
  35. Mehrgan, K., Thomas, J., Saglia, R., et al. 2024, ApJ, 961, 127 [NASA ADS] [CrossRef] [Google Scholar]
  36. Milosavljević, M., & Merritt, D. 2001, ApJ, 563, 34 [Google Scholar]
  37. Naab, T., Johansson, P. H., & Ostriker, J. P. 2009, ApJ, 699, L178 [Google Scholar]
  38. Nasim, I. T., Gualandris, A., Read, J. I., et al. 2021, MNRAS, 502, 4794 [NASA ADS] [CrossRef] [Google Scholar]
  39. Neureiter, B., Thomas, J., Saglia, R., et al. 2021, MNRAS, 500, 1437 [Google Scholar]
  40. Neureiter, B., Thomas, J., Rantala, A., et al. 2023, ApJ, 950, 15 [NASA ADS] [CrossRef] [Google Scholar]
  41. Park, S., Yang, J., Oonk, J. B. R., & Paragi, Z. 2017, MNRAS, 465, 3943 [NASA ADS] [CrossRef] [Google Scholar]
  42. Rantala, A., Johansson, P. H., Naab, T., Thomas, J., & Frigo, M. 2018, ApJ, 864, 113 [Google Scholar]
  43. Rantala, A., Johansson, P. H., Naab, T., Thomas, J., & Frigo, M. 2019, ApJ, 872, L17 [Google Scholar]
  44. Rusli, S. P., Erwin, P., Saglia, R. P., et al. 2013a, AJ, 146, 160 [CrossRef] [Google Scholar]
  45. Rusli, S. P., Thomas, J., Saglia, R. P., et al. 2013b, AJ, 146, 45 [NASA ADS] [CrossRef] [Google Scholar]
  46. Saglia, R. P., Sánchez-Blázquez, P., Bender, R., et al. 2010, A&A, 524, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Saglia, R. P., Opitsch, M., Erwin, P., et al. 2016, ApJ, 818, 47 [Google Scholar]
  48. Sánchez-Blázquez, P., Peletier, R. F., Jiménez-Vicente, J., et al. 2006, MNRAS, 371, 703 [Google Scholar]
  49. Teyssier, R., Moore, B., Martizzi, D., Dubois, Y., & Mayer, L. 2011, MNRAS, 414, 195 [NASA ADS] [CrossRef] [Google Scholar]
  50. Thomas, J., & Lipka, M. 2022, MNRAS, 514, 6203 [CrossRef] [Google Scholar]
  51. Thomas, D., Maraston, C., & Bender, R. 2003, MNRAS, 339, 897 [NASA ADS] [CrossRef] [Google Scholar]
  52. Thomas, J., Saglia, R. P., Bender, R., et al. 2004, MNRAS, 353, 391 [NASA ADS] [CrossRef] [Google Scholar]
  53. Thomas, J., Saglia, R. P., Bender, R., Erwin, P., & Fabricius, M. 2014, ApJ, 782, 39 [Google Scholar]
  54. Thomas, J., Ma, C.-P., McConnell, N. J., et al. 2016, Nature, 532, 340 [NASA ADS] [CrossRef] [Google Scholar]
  55. Veale, M., Ma, C.-P., Thomas, J., et al. 2017, MNRAS, 464, 356 [Google Scholar]
  56. Zhao, H. 1996, MNRAS, 278, 488 [Google Scholar]

All Tables

Table 1.

Parameters of the core-Sérsic best fits.

Table 2.

Parameters of the second Sérsic component.

Table 3.

Parameters of the Nuker best fit.

Table 4.

Parameters of the axisymmetric and triaxial dynamical modeling.

All Figures

thumbnail Fig. 1.

Cutout of the Euclid VIS image of NGC 1272.

In the text
thumbnail Fig. 2.

Photometry profiles of NGC 1272. Top Left: Surface brightness profile of NGC 1272 measured from the Euclid VIS image, calibrated to the V band, and corrected for Galactic absorption and cosmological dimming, as a function of the 1/4 power of the semi-major distance A on the sky in arcsec. Bottom Left: Difference ΔSB between the surface brightness profile of NGC 1272 and the surface brightness of the core-Sérsic+Sérsic model. Right: Ellipticity 1 − B/A, where B is the semi-minor axis length on the sky (top) and PA (bottom) as a function of the 1/4 power of A. The solid red lines show the core-Sérsic+Sérsic model. The dashed lines show its core radius. The blue line shows the ellipticity profile of the axisymmetric deprojection. The green line shows the PA profile of the triaxial deprojection.

In the text
thumbnail Fig. 3.

The IE, YE, JE, and HE cutouts of the inner 6″ × 6″ of NGC 1272.

In the text
thumbnail Fig. 4.

Percentage residuals after subtraction of the core-Sérsic + Sérsic model.

In the text
thumbnail Fig. 5.

Profiles of p(r) and q(r) (where p = b/a, q = c/a, and a, b, and c are the major, intermediate, and minor semi-axis of the galaxy) derived from the axisymmetric (dotted), triaxial (full line), spherical (dashed line), and prolate (dashed-dotted line) deprojections of the galaxy, as a function of distance r from the center. The red and blue shaded areas show the whole range of allowed deprojections with RMS ≤ 1.2 × RMSmin (de Nicola et al. 2020, 2022a,b).

In the text
thumbnail Fig. 6.

Two-dimensional stellar kinematics of NGC 1272. The horizontal and vertical dashed lines show the major and minor axes of the galaxy, respectively. North is up and east is to the left.

In the text
thumbnail Fig. 7.

Radial stellar kinematics of NGC 1272. The red lines show the axisymmetric fit to the stellar kinematics of the galaxy (black data points with error bars) as a function of the distance R from the center of the galaxy on the sky.

In the text
thumbnail Fig. 8.

Line-of-sight velocity distribution measured at 4 arcsec from the center of NGC 1272 (filled circles with error bars). The red line connects the values provided by the base model.

In the text
thumbnail Fig. 9.

Results of the axisymmetric modeling of NGC 1272 of the complete stellar kinematic dataset. As a function of the quality of the fits measured by the AICp we show: from left to right, and from top to bottom: the BH mass MBH; the dynamical V-band mass-to-light ratio Υ*, V; the DM density at 10 kpc ρ10; and the inner slope of the DM density profile γDM. The gray points show the individual models, the red dots show the best-fitting model, and the black lines show the lower envelope of the distributions of gray points.

In the text
thumbnail Fig. 10.

Spherically averaged mass profiles of NGC 1272. The solid and dashed lines show the total and stellar profiles, respectively. The triaxial, prolate, and spherical models are averaged over the two sides. The vertical dashed line shows the position of the core radius.

In the text
thumbnail Fig. 11.

Anisotropy profiles β (where β = 1 − σT2/σR2, and σT and σR are the spherical tangential and radial velocity dispersions, respectively) of the different models. The triaxial, prolate, and spherical models are averaged over the two sides. The distances to the center are in units of the core radius.

In the text
thumbnail Fig. 12.

Position of NGC 1272 (in red) on the MBH − rb (left), the MBH − M* (middle), and the MBH − σ (right) relation. The black data points are from Rusli et al. (2013a), Saglia et al. (2016), Thomas et al. (2016), Mehrgan et al. (2019), Neureiter et al. (2023), and de Nicola et al. (2024). The blue data points are core ellipticals with a stellar mass of greater than 1012M. While the galaxy follows the MBH − rb given by Thomas et al. (2016) and the MBH − M* relation of Saglia et al. (2016) for the sample of CorePowerE, it deviates by a factor of 8.4 from the MBH − σ relation of Saglia et al. (2016), or by a factor of 1.8 of the 1σ error combined with the intrinsic scatter in the relation (shown by the dotted lines).

In the text
thumbnail Fig. 13.

Correlation between the core radius rb and stellar mass M*. The black datapoints are from Rusli et al. (2013a), Saglia et al. (2016), Thomas et al. (2016), Mehrgan et al. (2019), Neureiter et al. (2023), and de Nicola et al. (2024). NGC 1272 is shown in red. The blue cross shows the position of the BCG on the EDISCS cluster CL1216 (Saglia et al. 2010). The black solid line shows log10(M*/M) = 0.64log10(rb/kpc)+12.1; the dotted lines show the ±1σ scatter (0.26 dex) in the relation.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.