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\sigma$ 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.

We use observations from the GEMINI-N/GMOS integral-field spectrograph (IFS) to obtain spatially resolved stellar kinematics of the central $\sim 1$ kpc of 20 early-type galaxies (ETGs) with stellar masses greater than $10^{11.7} M_\odot$ in the MASSIVE survey. Together with observations from the wide-field Mitchell IFS at McDonald Observatory in our earlier work, we obtain unprecedentedly detailed kinematic maps of local massive ETGs, covering a scale of $\sim 0.1-30$ kpc. The high ($\sim 120$) signal-to-noise of the GMOS spectra enable us to obtain two-dimensional maps of the line-of-sight velocity, velocity dispersion $\sigma$, as well as the skewness $h_3$ and kurtosis $h_4$ of the stellar velocity distributions. All but one galaxy in the sample have $\sigma(R)$ profiles that increase towards the center, whereas the slope of $\sigma(R)$ at one effective radius ($R_e$) can be of either sign. The $h_4$ is generally positive, with 14 of the 20 galaxies having positive $h_4$ within the GMOS aperture and 18 having positive $h_4$ within $1 R_e$. The positive $h_4$ and rising $\sigma(R)$ towards small radii are indicative of a central black hole and velocity anisotropy. We demonstrate the constraining power of the data on the mass distributions in ETGs by applying Jeans anisotropic modeling (JAM) to NGC~1453, the most regular fast rotator in the sample. Despite the limitations of JAM, we obtain a clear $\chi^2$ minimum in black hole mass, stellar mass-to-light ratio, velocity anisotropy parameters, and the circular velocity of the dark matter halo.

Short-duration GRBs with extended emission form a subclass of short GRBs, comprising about 15% of the short-duration sample. Afterglow detections of short GRBs are also rare (about 30%) due to their smaller luminosity. We present a multi-band data set of the short burst with extended emission GRB 150424A, comprising of GROND observations, complemented with data from Swift/UVOT, Swift/XRT, HST, Keck/LRIS and data points from the literature. The GRB 150424A afterglow shows an extended plateau phase, lasting about 8hrs. The analysis of this unique GRB afterglow might shed light on the understanding of afterglow plateau emission, the nature of which is still under debate. We present a phenomenological analysis by applying fireball closure relations, and interpret the findings in the context of the fireball model. We discuss the plausibility of a magnetar as a central engine, being responsible for additional and prolonged energy injection into the fireball. We find convincing evidence for energy injection into the afterglow of GRB 150424A. We find that a magnetar spin down as source for a prolonged energy injection requires that at least 4% of the spin-down energy is converted to radiation.

We analyse the environmental properties of 370 local early-type galaxies (ETGs) in the MASSIVE and ATLAS3D surveys, two complementary volume-limited integral-field spectroscopic (IFS) galaxy surveys spanning absolute $K$-band magnitude $-21.5 > M_K > -26.6$, or stellar mass $8 \times 10^{9} < M_* < 2 \times 10^{12} M_\odot$. We find these galaxies to reside in a diverse range of environments measured by four methods: group membership (whether a galaxy is a brightest group/cluster galaxy, satellite, or isolated), halo mass, large-scale mass density (measured over a few Mpc), and local mass density (measured within the $N$th neighbour). The spatially resolved IFS stellar kinematics provide robust measurements of the spin parameter $\lambda_e$ and enable us to examine the relationship among $\lambda_e$, $M_*$, and galaxy environment. We find a strong correlation between $\lambda_e$ and $M_*$, where the average $\lambda_e$ decreases from $\sim 0.4$ to below 0.1 with increasing mass, and the fraction of slow rotators $f_{\rm slow}$ increases from $\sim 10$% to 90%. We show for the first time that at fixed $M_*$, there are almost no trends between galaxy spin and environment; the apparent kinematic morphology-density relation for ETGs is therefore primarily driven by $M_*$ and is accounted for by the joint correlations between $M_*$ and spin, and between $M_*$ and environment. A possible exception is that the increased $f_{\rm slow}$ at high local density is slightly more than expected based only on these joint correlations. Our results suggest that the physical processes responsible for building up the present-day stellar masses of massive galaxies are also very efficient at reducing their spin, in any environment.

We present the first systematic investigation of the existence, spatial distribution, and kinematics of warm ionized gas as traced by the [O II] 3727AA emission line in 74 of the most massive galaxies in the local Universe. All of our galaxies have deep integral field spectroscopy from the volume- and magnitude-limited MASSIVE survey of early-type galaxies with stellar mass log(M_*/M_sun) > 11.5 (M_K < -25.3 mag) and distance D < 108 Mpc. Of the 74 galaxies in our sample, we detect warm ionized gas in 28, which yields a global detection fraction of 38\pm6% down to a typical [O II] equivalent width limit of 2AA. MASSIVE fast rotators are more likely to have gas than MASSIVE slow rotators with detection fractions of 80\pm10% and 28\pm6%, respectively. The spatial extents span a wide range of radii (0.6 - 18.2 kpc; 0.1 - 4R_e), and the gas morphologies are diverse, with 17/28 = 61\pm9% being centrally concentrated, 8/28 = 29\pm9% exhibiting clear rotation out to several kpc, and 3/28 = 11\pm6% being extended but patchy. Three out of four fast rotators show kinematic alignment between the stars and gas, whereas the two slow rotators with robust kinematic measurements available exhibit kinematic misalignment. Our inferred warm ionized gas masses are roughly ~10^5M_sun. The emission line ratios and radial equivalent width profiles are generally consistent with excitation of the gas by the old underlying stellar population. We explore different gas origin scenarios for MASSIVE galaxies and find that a variety of physical processes are likely at play, including internal gas recycling, cooling out of the hot gaseous halo, and gas acquired via mergers.

We present spatially-resolved two-dimensional stellar kinematics for the 41 most massive early-type galaxies (MK <~ -25.7 mag, stellar mass M* >~ 10^11.8 Msun) of the volume-limited (D < 108 Mpc) MASSIVE survey. For each galaxy, we obtain high-quality spectra in the wavelength range of 3650 to 5850 A from the 246-fiber Mitchell integral-field spectrograph (IFS) at McDonald Observatory, covering a 107" x 107" field of view (often reaching 2 to 3 effective radii). We measure the 2-D spatial distribution of each galaxy's angular momentum (lambda and fast or slow rotator status), velocity dispersion (sigma), and higher-order non-Gaussian velocity features (Gauss-Hermite moments h3 to h6). Our sample contains a high fraction (~80% ) of slow and non-rotators with lambda <~ 0.2. When combined with the lower-mass ETGs in the ATLAS3D survey, we find the fraction of slow-rotators to increase dramatically with galaxy mass, reaching ~50% at MK ~ -25.5 mag and ~90% at MK <~ -26 mag. All of our fast rotators show a clear anti-correlation between h3 and V/sigma, and the slope of the anti-correlation is steeper in more round galaxies. The radial profiles of sigma show a clear luminosity and environmental dependence: the 12 most luminous galaxies in our sample (MK <~ -26 mag) are all brightest cluster/group galaxies (except NGC 4874) and all have rising or nearly flat sigma profiles, whereas five of the seven "isolated" galaxies are all fainter than MK = -25.8 mag and have falling sigma. All of our galaxies have positive average h4; the most luminous galaxies have average h4 ~ 0.05 while less luminous galaxies have a range of values between 0 and 0.05. Most of our galaxies show positive radial gradients in h4, and those galaxies also tend to have rising sigma profiles. We discuss the implications for the relationship among dynamical mass, sigma, h4, and velocity anisotropy for these massive galaxies.

Studies of the physical properties of local elliptical galaxies (e.g., gas temperatures, halo masses, stellar kinematics) are shedding new light on galaxy formation. Here we present the hot X-ray gas properties of 33 early-type systems within the MASSIVE galaxy survey sample that have archival Chandra X-ray observations. Through careful X-ray spectral modeling, we derive X-ray luminosities (L_X) and plasma temperatures (T_gas) for the diffuse gas components in these galaxies. We combine the MASSIVE sample with 41 galaxies from the ATLAS^3D survey to investigate the X-ray and optical properties of a statistically significant sample of nearby early-type galaxies across a wide-range of environments. We deduce that all early-type galaxies (independent of galaxy mass and rotational support) follow a universal scaling law such that L_X~T_gas^4.5. When X-ray measurements are performed consistently in apertures set by the galaxy stellar content, the wide-scale environment does not contribute to the intrinsic scatter (~0.5dex) within the scaling relation. We further demonstrate that the scatter in L_X around both K-band luminosity (L_K) and the galaxy stellar velocity dispersion is primarily driven by T_gas, with no clear trends with halo mass, radio power, or angular momentum of the stars. It is not trivial to tie the origin of the gas directly to either the stellar mass or the galaxy potential. Indeed, our data require a steeper relation between L_X, L_K, and sigma_e than predicted by standard mass-loss models. Finally, we find a statistically significant correlation between sigma_e and T_gas, suggesting that T_gas is set by the galaxy potential inside the optical effective radius. We conclude that within the inner-most 10-30kpc region, early-types maintain pressure-supported hot gas, with a minimum T_gas set by the virial temperature, but the majority show evidence for some additional heating.

Quasars are associated with and powered by the accretion of material onto massive black holes; the detection of highly luminous quasars with redshifts greater than z = 6 suggests that black holes of up to ten billion solar masses already existed 13 billion years ago. Two possible present-day dormant descendants of this population of active black holes have been found in the galaxies NGC 3842 and NGC 4889 at the centres of the Leo and Coma galaxy clusters, which together form the central region of the Great Wall - the largest local structure of galaxies. The most luminous quasars, however, are not confined to such high-density regions of the early Universe; yet dormant black holes of this high mass have not yet been found outside of modern-day rich clusters. Here we report observations of the stellar velocity distribution in the galaxy NGC 1600 - a relatively isolated elliptical galaxy near the centre of a galaxy group at a distance of 64 Mpc from Earth. We use orbit superposition models to determine that the black hole at the centre of NGC 1600 has a mass of 17 billion solar masses. The spatial distribution of stars near the centre of NGC 1600 is rather diffuse. We find that the region of depleted stellar density in the cores of massive elliptical galaxies extends over the same radius as the gravitational sphere of influence of the central black holes, and interpret this as the dynamical imprint of the black holes.

Samples of early-type galaxies show a correlation between stellar velocity dispersion and the stellar initial mass function (IMF) as inferred from gravity-sensitive absorption lines in the galaxies' central regions. To search for spatial variations in the IMF, we have observed two early-type galaxies with Keck/LRIS and measured radial gradients in the strengths of absorption features from 4000-5500 \AA $\,$ and 8000-10,000 \AA. We present spatially resolved measurements of the dwarf-sensitive spectral indices NaI (8190 \AA) and Wing-Ford FeH (9915 \AA), as well as indices for species of H, C$_2$, CN, Mg, Ca, TiO, and Fe. Our measurements show a metallicity gradient in both objects, and Mg/Fe consistent with a shallow gradient in \alpha-enhancement, matching widely observed trends for massive early-type galaxies. The NaI index and the CN$_1$ index at 4160 \AA $\,$ exhibit significantly steeper gradients, with a break at $r \sim 0.1 r_{\rm eff}$ ($r \sim 300$ pc). Inside this radius NaI strength increases sharply toward the galaxy center, consistent with a rapid central rise in [Na/Fe]. In contrast, the ratio of FeH to Fe index strength decreases toward the galaxy center. This behavior cannot be reproduced by a steepening IMF inside $0.1 r_{\rm eff}$ if the IMF is a single power law. While gradients in the mass function above $\sim 0.4 M_\odot$ may occur, exceptional care is required to disentangle these IMF variations from the extreme variations in individual element abundances near the galaxies' centers.

In this work we present CO(1-0) and CO(2-1) observations of a pilot sample of 15 early-type galaxies (ETGs) drawn from the MASSIVE galaxy survey, a volume-limited integral-field spectroscopic study of the most massive ETGs ($M_* >10^{11.5}M_\odot$) within 108 Mpc. These objects were selected because they showed signs of an interstellar medium and/or star formation. A large amount of gas ($>$2$\times$10$^8$ M$_{\odot}$) is present in 10 out of 15 objects, and these galaxies have gas fractions higher than expected based on extrapolation from lower mass samples. We tentatively interpret this as evidence that stellar mass loss and hot halo cooling may be starting to play a role in fuelling the most massive galaxies. These MASSIVE ETGs seem to have lower star-formation efficiencies (SFE=SFR/M$_{\rm H2}$) than spiral galaxies, but the SFEs derived are consistent with being drawn from the same distribution found in other lower mass ETG samples. This suggests that the SFE is not simply a function of stellar mass, but that local, internal processes are more important for regulating star formation. Finally we used the CO line profiles to investigate the high-mass end of the Tully-Fisher relation (TFR). We find that there is a break in the slope of the TFR for ETGs at high masses (consistent with previous studies). The strength of this break correlates with the stellar velocity dispersion of the host galaxies, suggesting it is caused by additional baryonic mass being present in the centre of massive ETGs. We speculate on the root cause of this change and its implications for galaxy formation theories.

We examine stellar population gradients in ~100 massive early type galaxies spanning 180 < sigma* < 370 km/s and M_K of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (Ma et al. 2014). Using integral-field spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population properties extending to beyond twice the effective radius (~2.5 R_e). Specifically, we examine trends in age, metallicity, and abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories and elemental yields. At a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high redshift) stellar age and [alpha/Fe] increase with increasing sigma* and depend only weakly on stellar mass, as we might expect if denser galaxies form their central cores earlier and faster. If we instead focus on 1-1.5 R_e, the trends in abundance and abundance ratio are washed out, as might be expected if the stars at large radius were accreted by smaller galaxies. Finally, we show that when controlling for \sigmastar, there are only very subtle differences in stellar population properties or gradients as a function of group richness; even at large radius internal properties matter more than environment in determining star formation history.

Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ~100 most massive early-type galaxies within a distance of 108 Mpc. This survey probes a stellar mass range M* > 10^{11.5} Msun and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results.

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).

New kinematic data and modeling efforts in the past few years have substantially expanded and revised dynamical measurements of black hole masses (Mbh) at the centers of nearby galaxies. Here we compile an updated sample of 72 black holes and their host galaxies, and present revised scaling relations between Mbh and stellar velocity dispersion (sigma), V-band luminosity (L), and bulge stellar mass (Mbulge), for different galaxy subsamples. Our best-fitting power law relations for the full galaxy sample are log(Mbh) = 8.32 + 5.64*log(sigma/200 kms), log(Mbh) = 9.23 + 1.11*log(L/10^{11} Lsun), and log(Mbh) = 8.46 + 1.05*log(Mbulge/10^{11} Msun). A log-quadratic fit to the Mbh-sigma relation with an additional term of beta_2*[log(sigma/200 kms)]^2 gives beta_2 = 1.68 +/- 1.82 and does not decrease the intrinsic scatter in Mbh. When the early- and late-type galaxies are fit separately, we obtain similar slopes of 5.20 and 5.06 for the Mbh-sigma relation but significantly different intercepts -- Mbh in early-type galaxies are about 2 times higher than in late types at a given sigma. Within early-type galaxies, our fits to Mbh(sigma) give Mbh that is about 2 times higher in galaxies with central core profiles than those with central power-law profiles. Our Mbh-L and Mbh-Mbulge relations for early-type galaxies are similar to those from earlier compilations. When the conventional quadrature method is used to determine the intrinsic scatter in Mbh, our dataset shows weak evidence for increased scatter at Mbulge < 10^{11} Msun or L_V < 10^{10.3} Lsun, while the scatter stays constant for 10^{11} < Mbulge < 10^{12.3} Msun and 10^{10.3} < L_V < 10^{11.5} Lsun. A Bayesian analysis indicates a larger sample of Mbh measurements would be needed to detect any statistically significant trend in the scatter with galaxy properties.

We present stellar kinematics and orbit superposition models for the central regions of four Brightest Cluster Galaxies (BCGs), based upon integral-field spectroscopy at Gemini, Keck, and McDonald Observatories. Our integral-field data span radii from < 100 pc to tens of kpc. We report black hole masses, M_BH, of 2.1 +/- 1.6 x 10^10 M_Sun for NGC 4889, 9.7 + 3.0 - 2.6 x 10^9 M_Sun for NGC 3842, and 1.3 + 0.5 - 0.4 x 10^9 M_Sun for NGC 7768. For NGC 2832 we report an upper limit of M_BH < 9 x 10^9 M_Sun. Stellar orbits near the center of each galaxy are tangentially biased, on comparable spatial scales to the galaxies' photometric cores. We find possible photometric and kinematic evidence for an eccentric torus of stars in NGC 4889, with a radius of nearly 1 kpc. We compare our measurements of M_BH to the predicted black hole masses from various fits to the relations between M_BH and stellar velocity dispersion, luminosity, or stellar mass. The black holes in NGC 4889 and NGC 3842 are significantly more massive than all dispersion-based predictions and most luminosity-based predictions. The black hole in NGC 7768 is consistent with a broader range of predictions.

We have designed an inquiry-based laboratory activity on transiting extrasolar planets for an introductory college-level astronomy class. The activity was designed with the intent of simultaneously teaching science process skills and factual content about transits and light curves. In the activity, groups of two to four students each formulate a specific science question and design and carry out an investigation using a table-top model of a star and orbiting planet. Each group then presents their findings to other students in their class. In a final presentation, the instructors integrate students' findings with a summary of how measured light curves indicate properties of planetary systems. The activity debuted at Hartnell College in November 2009 and has also been adapted for a lecture-based astronomy course at U.C. Santa Cruz. We present the results of student surveys before and after the astronomy course at Hartnell and discuss how well our activity promotes students' confidence and identity as scientists, relative to traditional lab activities.

As a part of our ongoing Volume-limited A-Star (VAST) adaptive optics survey, we have obtained observations of 26 binary systems with projected separations <100 AU, 13 of which have sufficient historical measurements to allow for refinement of their orbital elements. For each system with an estimated orbit, the dynamical system mass obtained was compared with the system mass estimated from mass-magnitude relations. Discrepancies between the dynamical and theoretical system mass can be explained by the presence of a previously unresolved spectroscopic component, or by a non-solar metallicity of the system. Using this approach to infer the presence of additional companions, a lower limit to the fraction of binaries, triples, and quadruples can be estimated as 39, 46, and 15 per cent, for systems with at least one companion within 100 AU. The fraction of multiple systems with three or more components shows a relative increase compared to the fraction for Solar-type primaries resolved in previous volume-limited surveys. The observations have also revealed a pair of potentially young ($<$100 Myr) M-dwarf companions, which would make an ideal benchmark for the theoretical models during the pre-Main Sequence contraction phase for M-dwarfs. In addition to those systems with orbit fits, we report 13 systems for which further orbital monitoring observations are required, 11 of which are newly resolved as a part of the VAST survey.

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.

We present the first direct measurement of the central black hole mass, M_BH, in NGC 6086, the Brightest Cluster Galaxy (BCG) in Abell 2162. Our investigation demonstrates for the first time that stellar dynamical measurements of M_BH in BCGs are possible beyond the nearest few galaxy clusters. We observed NGC 6086 with laser guide star adaptive optics and the integral-field spectrograph (IFS) OSIRIS at the W.M. Keck Observatory, and with the seeing-limited IFS GMOS-N at Gemini Observatory North. We combined the two IFS data sets with existing major-axis kinematics, and used axisymmetric stellar orbit models to determine M_BH and the R-band stellar mass-to-light ratio, M*/L_R. We find M_BH = 3.6(+1.7)(-1.1) x 10^9 M_Sun and M*/L_R = 4.6(+0.3)(-0.7) M_Sun/L_Sun (68% confidence), from models using the most massive dark matter halo allowed within the gravitational potential of the host cluster. Models fitting only IFS data confirm M_BH ~ 3 x 10^9 M_Sun and M*/L_R ~ 4 M_Sun/L_Sun, with weak dependence on the dark matter halo structure. When data out to 19 kpc are included, the unrealistic omission of dark matter causes the best-fit black hole mass to decrease dramatically, to 0.6 x 10^9 M_Sun, and the best-fit stellar mass-to-light ratio to increase to 6.7 M_Sun/L_Sun. The latter value is at further odds with stellar population studies favoring M*/L ~ 2 M_Sun/L_Sun,R. Biases from dark matter omission could extend to dynamical models of other galaxies with central stellar cores, and new measurements of M_BH from models with dark matter could steepen the empirical scaling relationships between black holes and their host galaxies.

While both society and astronomy have evolved greatly over the past fifty years, the academic institutions and incentives that shape our field have remained largely stagnant. As a result, the astronomical community is faced with several major challenges, including: (1) the training that we provide does not align with the skills that future astronomers will need, (2) the postdoctoral phase is becoming increasingly demanding and demoralizing, and (3) our jobs are increasingly unfriendly to families with children. Solving these problems will require conscious engineering of our profession. Fortunately, this Decadal Review offers the opportunity to revise outmoded practices to be more effective and equitable. The highest priority of the Subcommittee on the State of the Profession should be to recommend specific, funded activities that will ensure the field meets the challenges we describe.