Researcher profile

Chung-Pei Ma

Chung-Pei Ma contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

preprint2022arXiv

Gas inflows in the polar ring of NGC 4111: the birth of an AGN

We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ Åvelocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$γ$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring.

preprint2022arXiv

The MASSIVE Survey -- XVII. A Triaxial Orbit-based Determination of the Black Hole Mass and Intrinsic Shape of Elliptical Galaxy NGC 2693

We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph (GMOS) with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central $\sim 150$ pc out to $\sim2.5$ effective radii. Observations from Hubble Space Telescope (HST) WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6-D galaxy model parameter space to determine NGC 2693's SMBH mass ($M_\text{BH}$), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We find $M_\text{BH} = \left(1.7\pm 0.4\right)\times 10^{9}\ M_\odot$ and a triaxial intrinsic shape with axis ratios $p=b/a=0.902 \pm 0.009$ and $q=c/a=0.721^{+0.011}_{-0.010}$, triaxiality parameter $T = 0.39 \pm 0.04$. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 prefer $M_\text{BH} = \left(2.4\pm 0.6\right)\times 10^{9}\ M_\odot$ and $M_\text{BH} = \left(2.9\pm 0.3\right)\times 10^{9}\ M_\odot$, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data.

preprint2022arXiv

Toward determining the number of observable supermassive black hole shadows

We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ${\sim}0.1$ $μ$as resolution and ${\sim}1$ $μ$Jy sensitivity could access ${>}10^6$ SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ${\lesssim}2$ $μ$as resolution and ${\lesssim}0.5$ mJy sensitivity, while doing the same for second-order photon rings requires ${\lesssim}0.1$ $μ$as resolution and ${\lesssim}5$ $μ$Jy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), while a next-generation EHT observing at 345 GHz should have access to ${\sim}$3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings.

preprint2021arXiv

The MASSIVE SURVEY XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-Type Galaxies

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work we focus on the most massive galaxies in the nearby universe $\log(M_{\star}/M_{\odot})>11.2$. We obtain high quality Magellan/LDSS-3 long slit spectroscopy with a wide wavelength coverage of $0.4μ{\rm m}-1.01μ{\rm m}$ for 41 early-type galaxies (ETGs) in the MASSIVE survey, and derive high S/N spectra within an aperture of $R_{\rm e}/8$. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, $α_{\rm IMF}=(M/L)/(M/L)_{\rm MW}$, ranges from 1.12 to 3.05, with an average of $\langle α_{\rm IMF} \rangle=1.84$, suggesting that on average, the IMF is more bottom-heavy than Salpeter. Comparing the estimated stellar mass with the dynamical mass, we find that most galaxies have stellar masses smaller than their dynamical masses within the $1σ$ uncertainty. We complement our sample with lower-mass galaxies from the literature, and confirm that $\log(α_{\rm IMF})$ is positively correlated with $\log(σ)$, $\log(M_{\star})$, and $\log(M_{\rm dyn})$. The IMF in the centers of more massive ETGs is more bottom-heavy. In addition, we find that $\log(α_{\rm IMF})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density $Σ_{\rm Kroupa}$ might be responsible for the variation of $α_{\rm IMF}$. We conclude that $σ$, [Mg/Fe] and [Z/H] are the primary drivers of the global stellar IMF variation.

preprint2021arXiv

Triaxial Orbit-based Dynamical Modeling of Galaxies with Supermassive Black Holes and an Application to Massive Elliptical Galaxy NGC 1453

Most stellar-dynamical determinations of the masses of nearby supermassive black holes (SMBHs) have been obtained with the orbit superposition technique under the assumption of axisymmetry. However, few galaxies -- in particular massive early-type galaxies -- obey exact axisymmetry. Here we present a revised orbit superposition code and a new approach for dynamically determining the intrinsic shapes and mass parameters of triaxial galaxies based on spatially-resolved stellar kinematic data. The triaxial TriOS code described here corrects an error in the original van den Bosch et al. (2008) code that gives rise to incorrect projections for most orbits in triaxial models and can significantly impact parameter search results. The revised code also contains significant improvements in orbit sampling, mass constraints, and run time. Furthermore, we introduce two new parameter searching strategies -- a new set of triaxial shape parameters and a novel grid-free sampling technique -- that together lead to a remarkable gain in efficiency in locating the best-fit model. We apply the updated code and search method to NGC 1453, a fast-rotating massive elliptical galaxy. A full 6D parameter search finds $p=b/a=0.933^{+0.014}_{-0.015}$ and $q=c/a=0.779\pm0.012$ for the intrinsic axis ratios and $T=0.33\pm0.06$ for the triaxiality parameter. Despite the deviations from axisymmetry, the best-fit SMBH mass, stellar mass-to-light ratio, and dark matter enclosed mass for NGC 1453 are consistent with the axisymmetric results. More comparisons between axisymmetric and triaxial modeling are needed before drawing general conclusions.

preprint2020arXiv

The MASSIVE Survey -- XV. A Stellar Dynamical Mass Measurement of the Supermassive Black Hole in Massive Elliptical Galaxy NGC 1453

We present a new stellar dynamical mass measurement ($M_\mathrm{BH}$) of the supermassive black hole (SMBH) in NGC 1453, a fast-rotating massive elliptical galaxy in the MASSIVE survey. We measure stellar kinematics in 135 spatial bins in the central 1.5 kpc by 2 kpc region of the galaxy using high signal-to-noise ($S/N \sim 130$) spectra from the Gemini-North GMOS integral-field spectrograph (IFS). Combining with wide-field IFS kinematics out to $\sim 3$ effective radii and stellar light distributions from Hubble Space Telescope (HST) WFC3 images, we perform Schwarzschild orbit-based mass modeling in the axisymmetric limit to constrain the mass components in NGC 1453. The best-fit black hole mass is $M_\mathrm{BH} =(2.9 \pm 0.4) 10^9 M_\odot$; the mass models without a central black hole are excluded at the $8.7σ$ level. The NGC 1453 black hole lies within the intrinsic scatter of the SMBH and galaxy scaling relations, unlike three other galaxies hosting $\gtrsim 10^{10} M_\odot$ SMBHs in the MASSIVE sample. The high-$S/N$ GMOS spectra enable us to determine 8 moments of the Gauss-Hermite expansion of the line-of-sight velocity distributions (LOSVDs), which are used as constraints in the orbit modeling. The stellar orbits in the mass models are further constrained to produce negligible $h_9$ through $h_{12}$ to minimize spurious behavior in the LOSVDs. We show that truncating the series at $h_4$, as was often done in prior work, leads to a much weaker constraint on the inferred $M_\mathrm{BH}$ for NGC 1453. Furthermore, we discuss precautions and modifications that are needed to achieve axisymmetry in triaxial orbit codes that use the Schwarzschild method to sample the start space of stellar orbits in triaxial gravitational potentials.

preprint2020arXiv

The MASSIVE Survey XIV -- Stellar Velocity Profiles and Kinematic Misalignments from 200 pc to 20 kpc in Massive Early-type Galaxies

We use high spatial resolution stellar velocity maps from the Gemini GMOS integral-field spectrograph (IFS) and wide-field velocity maps from the McDonald Mitchell IFS to study the stellar velocity profiles and kinematic misalignments from $\sim 200$ pc to $\sim 20$ kpc in 20 early-type galaxies with stellar mass $M_* > 10^{11.7} M_\odot$ in the MASSIVE survey. While 80% of the galaxies have low spins ($λ< 0.1$) and low rotational velocities ($< 50$ km/s) in both the central region and the main body, we find a diverse range of velocity features and misalignment angles. For the 18 galaxies with measurable central kinematic axes, 10 have well aligned kinematic axis and photometric major axis, and the other 8 galaxies have misalignment angles that are distributed quite evenly from $15^\circ$ to the maximal value of $90^\circ$. There is a strong correlation between central kinematic misalignment and galaxy spin, where all 4 galaxies with significant spins have well aligned kinematic and photometric axes, but only 43% of the low-spin galaxies are well aligned. The central and main-body kinematic axes within a galaxy are not always aligned. When the two kinematic axes are aligned ($\sim 60$% of the cases), they are either also aligned with the photometric major axis or orthogonal to it. We find 13 galaxies to also exhibit noticeable local kinematic twists, and one galaxy to have a counter-rotating core. A diverse assembly history consisting of multiple gas-poor mergers of a variety of progenitor mass ratios and orbits is likely to be needed to account for the predominance of low spins and the wide range of central and main-body velocity features reported here for local massive ETGs.

preprint2019arXiv

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

preprint2012arXiv

Galaxies in X-ray Groups. II. A Weak Lensing Study of Halo Centering

Locating the centers of dark matter halos is critical for understanding the mass profiles of halos as well as the formation and evolution of the massive galaxies that they host. The task is observationally challenging because we cannot observe halos directly, and tracers such as bright galaxies or X-ray emission from hot plasma are imperfect. In this paper we quantify the consequences of miscentering on the weak lensing signal from a sample of 129 X-ray selected galaxy groups in the COSMOS field with redshifts 0<z<1 and halo masses in the range 10^13 - 10^14 M_sun. By measuring the stacked lensing signal around eight different candidate centers (such as the brightest member galaxy, the mean position of all member galaxies, or the X-ray centroid), we determine which candidates best trace the center of mass in halos. In this sample of groups, we find that massive galaxies near the X-ray centroids trace the center of mass to <~75 kpc, while the X-ray position and centroids based on the mean position of member galaxies have larger offsets primarily due to the statistical uncertainties in their positions (typically ~50-150 kpc). Approximately 30% of groups in our sample have ambiguous centers with multiple bright or massive galaxies, and these groups show disturbed mass profiles that are not well fit by standard models, suggesting that they are merging systems. We find that halo mass estimates from stacked weak lensing can be biased low by 5-30% if inaccurate centers are used and the issue of miscentering is not addressed.

preprint2011arXiv

Two ten-billion-solar-mass black holes at the centres of giant elliptical galaxies

Observational work conducted over the last few decades indicates that all massive galaxies have supermassive black holes at their centres. Although the luminosities and brightness fluctuations of quasars in the early Universe suggest that some are powered by black holes with masses greater than 10 billion solar masses, the remnants of these objects have not been found in the nearby Universe. The giant elliptical galaxy Messier 87 hosts the hitherto most massive known black hole, which has a mass of 6.3 billion solar masses. Here we report that NGC 3842, the brightest galaxy in a cluster at a distance from Earth of 98 megaparsecs, has a central black hole with a mass of 9.7 billion solar masses, and that a black hole of comparable or greater mass is present in NGC 4889, the brightest galaxy in the Coma cluster (at a distance of 103 megaparsecs). These two black holes are significantly more massive than predicted by linearly extrapolating the widely-used correlations between black hole mass and the stellar velocity dispersion or bulge luminosity of the host galaxy. Although these correlations remain useful for predicting black hole masses in less massive elliptical galaxies, our measurements suggest that different evolutionary processes influence the growth of the largest galaxies and their black holes.

preprint2010arXiv

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo &#39;destruction&#39; rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite &#39;over-quenching&#39; in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

preprint2010arXiv

The Merger Rates and Mass Assembly Histories of Dark Matter Haloes in the Two Millennium Simulations

We construct merger trees of dark matter haloes and quantify their merger rates and mass growth rates using the joint dataset from the Millennium and Millennium-II simulations. The finer resolution of the Millennium-II Simulation has allowed us to extend our earlier analysis of halo merger statistics to an unprecedentedly wide range of descendant halo mass (10^10 < M0 < 10^15 Msun), progenitor mass ratio (10^-5 < xi < 1), and redshift (0 < z < 15). We update our earlier fitting form for the mean merger rate per halo as a function of M_0, xi, and z. The overall behavior of this quantity is unchanged: the rate per unit redshift is nearly independent of z out to z~15; the dependence on halo mass is weak (M0^0.13); and it is nearly a power law in the progenitor mass ratio (xi^-2). We also present a simple and accurate fitting formula for the mean mass growth rate of haloes as a function of mass and redshift. This mean rate is 46 Msun/yr for 10^12 Msun haloes at z=0, and it increases with mass as M^{1.1} and with redshift as (1+z)^2.5 (for z > 1). When the fit for the mean mass growth rate is integrated over a halo&#39;s history, we find excellent match to the mean mass assembly histories of the simulated haloes. By combining merger rates and mass assembly histories, we present results for the number of mergers over a halo&#39;s history and the statistics of the redshift of the last major merger.

preprint2008arXiv

Conditional Mass Functions and Merger Rates of Dark Matter Halos in the Ellipsoidal Collapse Model

Analytic models based on spherical and ellipsoidal gravitational collapse have been used to derive the mass functions of dark matter halos and their progenitors (the conditional mass function). The ellipsoidal model generally provides a better match to simulation results, but there has been no simple analytic expression in this model for the conditional mass function that is accurate for small time steps, a limit that is important for generating halo merger trees and computing halo merger rates. We remedy the situation by deriving accurate analytic formulae for the first-crossing distribution, the conditional mass function, and the halo merger rate in the ellipsoidal collapse model in the limit of small look-back times. We show that our formulae provide a closer match to the Millennium simulation results than those in the spherical collapse model and the ellipsoidal model of Sheth & Tormen (2002).

preprint2008arXiv

How to Grow a Healthy Merger Tree

We investigate seven Monte Carlo algorithms -- four old and three new -- for constructing merger histories of dark matter halos using the extended Press-Schechter (EPS) formalism based on both the spherical and ellipsoidal collapse models. We compare, side-by-side, the algorithms&#39; abilities at reproducing the analytic EPS conditional (or progenitor) mass function over a broad range of mass and redshift (z=0 to 15). Among the four old algorithms (Lacey & Cole 1993, Kauffmann & White 1993, Somerville & Kolatt 1999, Cole et al 2000), we find that only KW93 produces a progenitor mass function that is consistent with the EPS prediction for all look-back redshifts. The origins of the discrepancies in the other three algorithms are discussed. Our three new algorithms are designed to generate the correct progenitor mass function at each timestep. We show that this is a necessary and sufficient condition for consistency with EPS at any look-back time. We illustrate the differences among the three new algorithms and KW93 by investigating two other conditional statistics: the mass function of the i_{th} most massive progenitors and the mass function for descendants with N_p progenitors.

preprint2005arXiv

Impact of Dark Matter Substructure on the Matter and Weak Lensing Power Spectra

We explore the effect of substructure in dark matter halos on the power spectrum and bispectrum of matter fluctuations and weak lensing shear. By experimenting with substructure in a cosmological N = 512^3 simulation, we find that when a larger fraction of the host halo mass is in subhalos, the resulting power spectrum has less power at 1 < k < 100 h Mpc^{-1} and more power at k > 100 h Mpc^{-1}. We explain this effect using an analytic halo model including subhalos, which shows that the 1 < k < 100 h Mpc^{-1} regime depends sensitively on the radial distribution of subhalo centers while the interior structure of subhalos is important at k > 100 h Mpc^{-1}. The corresponding effect due to substructures on the weak lensing power spectrum is up to about 11% at angular scale l < 10^4. Predicting the nonlinear power spectrum to a few percent accuracy for future surveys would therefore require large cosmological simulations that also have exquisite numerical resolution to model accurately the survivals of dark matter subhalos in the tidal fields of their hosts.