The Role of Cold Gas in Low-Level Supermassive Black Hole Activity

Similar documents
The Correlation Between Supermassive Black Hole Mass and the Structure of Ellipticals and Bulges

A Supermassive Black Hole in the Dwarf Starburst Galaxy Henize Amy Reines Einstein Fellow National Radio Astronomy Observatory

Supermassive Black Hole Ac3vity in Field and Cluster Early Type Galaxies

Exploiting HI surveys with stacking. The HI content of massive galaxies from ALFALFA and GASS

Active Galaxies & Quasars

Part 2. Hot gas halos and SMBHs in optically faint ellipticals. Part 3. After Chandra?

Radio emission from galaxies in the Bootes Voids

Demographics of radio galaxies nearby and at z~0.55. Are radio galaxies signposts to black-hole mergers?

Lecture Outlines. Chapter 24. Astronomy Today 8th Edition Chaisson/McMillan Pearson Education, Inc.

The AGN / host galaxy connection in nearby galaxies.

There are three main ways to derive q 0 :

A zoo of transient sources. (c)2017 van Putten 1

Astro2010 Science White Paper: Tracing the Mass Buildup of Supermassive Black Holes and their Host Galaxies

Observing the Formation of Dense Stellar Nuclei at Low and High Redshift (?) Roderik Overzier Max-Planck-Institute for Astrophysics

Chapter 15 2/19/2014. Lecture Outline Hubble s Galaxy Classification. Normal and Active Galaxies Hubble s Galaxy Classification

Luminous radio-loud AGN: triggering and (positive?) feedback

Thus Far. Intro / Some Definitions Hubble Classification Components of Galaxies. Specific Galaxy Types Star Formation Clusters of Galaxies

Lecture Two: Galaxy Morphology:

GALAXIES. Edmund Hodges-Kluck Andrew Ptak

Black Holes in Hibernation

Quasars and Active Galactic Nuclei (AGN)

Probing the Origin of Supermassive Black Hole Seeds with Nearby Dwarf Galaxies. Amy Reines Einstein Fellow NRAO Charlottesville

Surface Photometry Quantitative description of galaxy morphology. Hubble Sequence Qualitative description of galaxy morphology

The Milky Way Galaxy. Some thoughts. How big is it? What does it look like? How did it end up this way? What is it made up of?

Journal Club Presentation on The BIMA Survey of Nearby Galaxies. I. The Radial Distribution of CO Emission in Spiral Galaxies by Regan et al.

THE GAS MASS AND STAR FORMATION RATE

Revealing new optically-emitting extragalactic Supernova Remnants

GALAXIES. I. Morphologies and classification 2. Successes of Hubble scheme 3. Problems with Hubble scheme 4. Galaxies in other wavelengths

Coevolution (Or Not) of Supermassive Black Holes and Galaxies

Debate on the toroidal structures around hidden- vs non hidden-blr of AGNs

What HI tells us about Galaxies and Cosmology

Feeding the Beast. Chris Impey (University of Arizona)

Black Holes and Active Galactic Nuclei

Quasars ASTR 2120 Sarazin. Quintuple Gravitational Lens Quasar

Part two of a year-long introduction to astrophysics:

Ultraluminous X-ray Sources in Nearby Galaxies

The Classification of Galaxies

Multi-wavelength Surveys for AGN & AGN Variability. Vicki Sarajedini University of Florida

Hubble Space Telescope ultraviolet spectroscopy of blazars: emission lines properties and black hole masses. E. Pian, R. Falomo, A.

MASSIVE BLACK HOLES AMY REINES IN NEARBY DWARF GALAXIES HUBBLE FELLOW NATIONAL OPTICAL ASTRONOMY OBSERVATROY

What HI tells us about Galaxies and Cosmology

Galaxies. Need a (physically) meaningful way of describing the relevant properties of a galaxy.

Quasars: Back to the Infant Universe

Galaxies. Galaxy Diversity. Galaxies, AGN and Quasars. Physics 113 Goderya

LeMMINGs the emerlin radio legacy survey of nearby galaxies Ranieri D. Baldi

New Results on the AGN Content of Galaxy Clusters

FORMATION OF SUPERMASSIVE BLACK HOLES Nestor M. Lasso Cabrera

Active Galactic Nuclei

ROSAT Roentgen Satellite. Chandra X-ray Observatory

Quasars and AGN. What are quasars and how do they differ from galaxies? What powers AGN s. Jets and outflows from QSOs and AGNs

Chapter 19 Galaxies. Hubble Ultra Deep Field: Each dot is a galaxy of stars. More distant, further into the past. halo

A galaxy is a self-gravitating system composed of an interstellar medium, stars, and dark matter.

Lecture 19: Galaxies. Astronomy 111

XRBs in Star-forming Galaxies X-ray Luminosity Functions Ultra-Luminous X-ray Sources

X-ray emission from star-forming galaxies

An analogy. "Galaxies" can be compared to "cities" What would you like to know about cities? What would you need to be able to answer these questions?

AGN Feedback In an Isolated Elliptical Galaxy

The cosmic distance scale

Citation for published version (APA): Fathi, K. (2004). Dynamics and morphology in the inner regions of spiral galaxies Groningen: s.n.

Galaxies with Active Nuclei. Active Galactic Nuclei Seyfert Galaxies Radio Galaxies Quasars Supermassive Black Holes

ASTRO504 Extragalactic Astronomy. 2. Classification

The parsec scale of. ac-ve galac-c nuclei. Mar Mezcua. International Max Planck Research School for Astronomy and Astrophysics

A100H Exploring the Universe: Quasars, Dark Matter, Dark Energy. Martin D. Weinberg UMass Astronomy

ACTIVE GALACTIC NUCLEI: optical spectroscopy. From AGN classification to Black Hole mass estimation

X-raying galactic feedback in nearby disk galaxies. Q. Daniel Wang University of Massachusetts

Fundamental Planes and Galaxy Formation

AGN Selec)on Techniques. Kris)n Kulas Astro 278 Winter 2012

Physics of Galaxies 2016 Exercises with solutions batch I

TEMA 3. Host Galaxies & Environment

Connection between AGN and Star Formation activities, or not?

Galaxy Morphology. - a description of the structure of galaxies

The Bright Side of the X-ray Sky The XMM-Newton Bright Survey. R. Della Ceca. INAF Osservatorio Astronomico di Brera,Milan

Scaling Correlations Among Central Massive Objects and Their Host Galaxies

Quasars in the SDSS. Rich Kron NGC June 2006 START CI-Team: Variable Quasars Research Workshop Yerkes Observatory

Our Galaxy. We are located in the disk of our galaxy and this is why the disk appears as a band of stars across the sky.

Active Galactic Nuclei (AGNs): A type of AGNs: Quasars. Whatever is powering these QSO s must be very small!!

The Rela(on Between Bars and AGN locally and at high redshi8

AST Cosmology and extragalactic astronomy. Lecture 20. Black Holes Part II

Halo Gas Velocities Using Multi-slit Spectroscopy

Major Review: A very dense article" Dawes Review 4: Spiral Structures in Disc Galaxies; C. Dobbs and J Baba arxiv "

Dual and Binary MBHs and AGN: Connecting Dynamics and Accretion

Frequency of Seyfert Type Transitions in a Sample of 102 Local Active Galactic Nuclei

Gas and Stellar Dynamical Black Hole Mass Measurements:

Radio Quiet AGN: Black Hole & Host Galaxy Proper;es

Radio infrared correlation for galaxies: from today's instruments to SKA

The Milky Way. Mass of the Galaxy, Part 2. Mass of the Galaxy, Part 1. Phys1403 Stars and Galaxies Instructor: Dr. Goderya

Detecting the cold neutral gas in young radio galaxies. James Allison Triggering Mechanisms for AGN

Introduction to the Sloan Survey

Astronomy 1 Fall 2016

PERSPECTIVES of HIGH ENERGY NEUTRINO ASTRONOMY. Paolo Lipari Vulcano 27 may 2006

Luminosity Functions of Planetary Nebulae & Globular Clusters. By Azmain Nisak ASTR 8400

Hubble s Law and the Cosmic Distance Scale

In a dense region all roads lead to a black Hole (Rees 1984 ARAA) Deriving the Mass of SuperMassive Black Holes

Active Galactic Nuclei - Zoology

ALFALFA status Jan Martha Haynes Cornell University UAT13 Jan 2013

A new catalogue of ultraluminous X-ray sources (and more!)

Survey of dusty AGNs based on the mid-infrared all-sky survey catalog. Shinki Oyabu (Nagoya University) & MSAGN team

ASTRON 449: Stellar (Galactic) Dynamics. Fall 2014

Results from the Chandra Deep Field North

A PRELIMINARY CLASSIFICATION SCHEME FOR THE CENTRAL REGIONS OF LATE-TYPE GALAXIES

Transcription:

Macalester Journal of Physics and Astronomy Volume 3 Issue 1 Spring 2015 Article 1 May 2015 The Role of Cold Gas in Low-Level Supermassive Black Hole Activity Erik D. Alfvin Macalester College, ealfvin@macalester.edu Abstract The nature of the relationship between low-level supermassive black hole activity and galactic cold gas, if any, is currently unclear. It has been hypothesized that feedback may heat or expel gas and quench star formation; alternatively, central black holes may feed at higher rates (either directly or as a secondary effect from stellar winds) in gas-rich galaxies. We use a combination of radio data from the on-going ALFALFA survey and from the literature, along with archival X-ray flux measurements from the Chandra X-ray observatory, to investigate this potential relationship. We construct a sample of 136 late-type galaxies, with MB < 18 out to 50 Mpc, that have both HI masses and sensitive X-ray coverage. Of these, 76 host a nuclear X-ray source, a 56% detection fraction. There is a highly significant correlation between LX and Mstar with a slope of 1.5±0.2, and a tentative correlation (significant at the 2.5σ level) between LX and MHI. However, a joint fit to LX as a function of both Mstar and MHI finds no significant dependence on MHI, and similarly the residuals of LX LX(Mstar) show no trend with MHI. We conclude that the galaxy-wide cold gas content in these spirals does not strongly influence their low-level supermassive black hole activity. Follow this and additional works at: http://digitalcommons.macalester.edu/mjpa Part of the External Galaxies Commons, Physical Processes Commons, and the Physics Commons Recommended Citation Alfvin, Erik D. (2015) "The Role of Cold Gas in Low-Level Supermassive Black Hole Activity," Macalester Journal of Physics and Astronomy: Vol. 3 : Iss. 1, Article 1. Available at: http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 This Capstone is brought to you for free and open access by the Physics and Astronomy Department at DigitalCommons@Macalester College. It has been accepted for inclusion in Macalester Journal of Physics and Astronomy by an authorized editor of DigitalCommons@Macalester College. For more information, please contact scholarpub@macalester.edu.

The Role of Cold Gas in Low-Level Supermassive Black Hole Activity Cover Page Footnote We gratefully acknowledge the assistance of M. Haynes, R. Giovanelli, and R. Koopmann with the preliminary α70 ALFALFA catalog, and we thank the UAT for support of this work. This capstone is available in Macalester Journal of Physics and Astronomy: http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1

Alfvin: Cold Gas and Black Hole Activity Draft version April 16, 2015 Preprint typeset using L A TEX style emulateapj v. 11/10/09 THE ROLE OF COLD GAS IN LOW-LEVEL SUPERMASSIVE BLACK HOLE ACTIVITY Erik Alfvin, Brendan P. Miller, Martha Haynes, Elena Gallo, Riccardo Giovanelli, Rebecca A. Koopmann Draft version April 16, 2015 Abstract The nature of the relationship between low-level supermassive black hole activity and galactic cold gas, if any, is currently unclear. It has been hypothesized that feedback may heat or expel gas (and consequently quench star formation); alternatively, central black holes may feed at higher rates (either directly or as a secondary effect from stellar winds) in gas-rich galaxies. We use a combination of radio data from the on-going ALFALFA survey and from the literature, along with archival X-ray flux measurements from the Chandra X-ray observatory, to investigate this potential relationship. We construct a sample of 136 late-type galaxies, with M B < 18 out to 50 Mpc, that have both HI masses and sensitive X-ray coverage. Of these, 76 host a nuclear X-ray source, a 56% detection fraction. There is a highly significant correlation between L X and M star with a slope of 1.5±0.2, and a tentative correlation (significant at the 2.5σ level) between L X and M HI. However, a joint fit to L X as a function of both M star and M HI finds no significant dependence on M HI, and similarly the residuals of L X L X (M star ) show no trend with M HI. We conclude that the galaxy-wide cold gas content in these spirals does not strongly correlate with their low-level supermassive black hole activity. Subject headings: galaxies black holes; accretion; cold gas 1. INTRODUCTION Nearly every large galaxy (i.e., M star > 10 10 M ) appears to contain a supermassive black hole (SMBH) at its center, with M BH correlating with the bulge stellar velocity dispersion σ, optical luminosity, or stellar mass (e.g., Gültekin et al. 2009; McConnell & Ma 2013, and references therein). Some accrete dust and gas while others do not, or accrete at different rates. The galactic processes that drive accretion are complex and gas movement within disks is not yet well understood. However, black hole mass is well studied and there are useful trends for this measurement in different types of galaxies. Compared to bulge-dominated early-type galaxies, where the bulge accounts for a significant amount of the galaxies luminosity and stellar mass, late-type galaxies follow a M BH σ relation with a similar slope (black hole mass positively correlating with greater velocity dispersion) but a lower intercept. The SMBHs hosted in spirals are therefore on average less massive by a factor of 2 than in ellipticals at a given σ (McConnell & Ma 2013). Late-type galaxies also offer disk tracers of SMBH mass, such as spiral arm pitch angle (e.g., Davis et al. 2014, and references therein), and it is now clear that there are bulgeless galaxies with SMBHs (Reines et al. 2011; Kormendy & Ho 2013; Bizzocchi et al. 2014; Satyapal et al. 2014). This suggests that black hole growth may be tied to galactic gas, and not just to bulge properties. The relationship between SMBH growth and activity and host galaxy evolution is complex, even in late type galaxies for which major mergers are not a triggering event. SMBH feedback may heat or expel gas (shutting down star formation), or cold gas reservoirs could mutually fuel star formation and SMBH accretion (the latter perhaps indirectly, from stellar winds). Or, particularly in smaller galaxies with smaller SMBHs, timescale mismatches and a limited energy budget may causally decouple these processes (e.g., Kormendy & Ho 2013, and references therein). Here, we explore the black hole accretion - cold gas mass relationship through investigating how nuclear X-ray luminosities are linked to HI masses in nearby spirals. Our HI masses are primarily drawn from the ALFALA survey. ALFALFA is a blind HI survey using the Arecibo L-band Feed Array (ALFA) at the Arecibo telescope to scan a portion of the 0 < δ < 36 sky at frequencies surrounding the HI 1420 MHz line (1335 1435 MHz). The ALFALFA survey has a limiting two-pass sensitivity of 1.8 mjy per beam and a half power beam width of 3.3 x3.8 for each of seven feeds (Giovanelli et al. 2005). The α40 preliminary catalog (40% coverage complete) is presented in Haynes et al. (2011). Here, we use the new α70 catalog (Haynes et al., in preparation; 70% coverage complete) to obtain HI masses for a sample of local spiral galaxies. We supplement the ALFALFA data with literature HI masses taken from the HyperLeda database 1 (Makarov et al. 2014) For X-ray luminosities we use archival Chandra data. High angular resolution Chandra X-ray observations provide a proven method of identifying low-level supermassive black hole activity in early-type galaxies (e.g., Gallo et al. 2010; Pellegrini 2010; Zhang et al.; Miller et al. 2012). In contrast to ellipticals, which do not typically contain large amounts of cold gas (particularly in clusters; e.g., Grossi et al. 2009 and references therein), spirals are often gas-rich. A gas-rich environment provides a complicating contamination possibility from highmass X-ray binaries (which can be present near sites of recent star formation), but here too Chandra generally has sufficient resolution to identify SMBH activity (e.g., O Sullivan et al. 2014; Tzanavaris et al. 2014; Mathur et al. 2010). We develop a statistical technique to quantify potential X-ray binary contamination and account for it in our linear regression modeling. 1 http://leda.univ-lyon1.fr/ Published by DigitalCommons@Macalester College, 2015 1

Macalester Journal of Physics and Astronomy, Vol. 3, Iss. 1 [2015], Art. 1 2 Figure 1. Left: Full band X-ray Chandra image of NGC 4501, showing nuclear X-ray source as well as a few XRBs. Right: DSS optical image at matched angular scale. Bottom: ALFALFA HI spectrum, showing double-peaked emission at the redshift of NGC 4501. 2. SAMPLE SELECTION AND CHARACTERISTICS Here we describe the optical selection, measurement of HI masses, and determination of X-ray coverage for our sample. X-ray, optical, and ALFALFA HI data are shown for reference for NGC 4501, the most luminous spiral in our sample, in Figure 1. 2.1. Optical selection of parent sample We selected a volume-limited parent sample of latetype galaxies from the HyperLeda database. We required each source to be a spiral galaxy with morphological type between t = 1 (Sa) and t = 7 (Sd), and imposed a luminosity limit of M B < 18. We additionally restrict consideration to spirals with inclination between the polar axis and the line-of-sight < 70 (i.e., not edge-on) to reduce the complicating effects of internal extinction. We initially select all late-type galaxies with a distance modulus less than 34.8 (i.e., within 90 Mpc). The HyperLeda SQL query that we used to construct our parent sample is: select objname, al2000, de2000, t, modz, mabs, bvtc, m21c where modz <= 34.8 and mabs <= -18 and t >= 1 and t <= 7 and incl < 70 order by mabs There are 6491 galaxies that meet these criteria; the full table is available as a text file. The distance in Mpc is calculated from the redshift distance modulus. The absolute magnitude is given under M B. The B V color is calculated from the difference between the B and V optical magnitudes corrected for galactic extinction. The stellar mass in solar units is calculated using the relations given in Bell et al. (2003) as: M star = 1.737(B V ) 0.942 + 0.4(B B) (1) Here the absolute magnitude for the sun is taken to be B = 5.515. 2 (2) 2.2. HI masses Nearly all ALFALFA HI detected galaxies have an optical counterpart (Haynes et al. 2011). The ALFALFA catalog uses source codes to classify objects. Code 1 objects have a high signal-to-noise ratio (S/N > 6.5) and are high confidence detections. Objects with a source code of 2 have a moderate signal-to-noise ratio (4 < S/N < 6.5) with an optical counterpart at the same redshift, while entries with a source code of 3 also have a moderate signal-to-noise ratio but without an optical counterpart; deeper follow-up radio observations indicate that code 2 objects are generally reliable while code 3 events are generally spurious. Code 4 and 5 objects have lower signal to noise ratios or RFI contamination and are not reliable sources. Finally, code 9 objects are likely high velocity clouds within the Milky Way. We used ALFALFA sources with code 1 and 2 only. To get HI data, we matched the resulting HyperLeda source positions to those in α70 with a search radius of 1 arc-min. We retrieved HI fluxes from both α70 and HyperLeda literature when possible. The HI mass was calculated from the relationship M HI = 2.36 10 5 d 2 MpcS HI (3) and is expressed as a logarithm. Our HI sample contains 3966 spirals that have either ALFALFA or HyperLeda literature HI measurements. In Figure 3 (left panel), we plot HI masses derived from each catalog for 780 sources with data in both. For all sources we used redshift distances from HyperLeda. The agreement between catalogs is excellent; the difference has median, mean, and standard deviation of 0.003, 0.008, and 0.158 dex. A text file is available with the optical properties and HI fluxes calculated from either literature Hyper- Leda values or listed in the α70 catalog for these 3966 spirals. 2.3. X-ray luminosities and limits Chandra is required to determine nuclear X-ray luminosities for our sample of local late-type galaxies due to its high angular resolution, which helps avoid contamination from X-ray binaries that may be present throughout the galaxy including near the (projected) center. Archival observations from ROSAT or XMM- Newton do not provide sufficient angular resolution and consequently are not considered. The Chandra Source Catalog (CSC; Evans et al. 2010) is a list of X-ray detected point sources from archived Chandra observations. The CSC contains many parameters including X-ray fluxes for each source. Where an X-ray source is not present at the specified coordinates, the catalog provides upper limit sensitivities, which we include in our analysis. The CSC includes imaging spectroscopy (ACIS-I or ACIS-S without gratings) data that is public as of 2011. The HyperLeda source positions were matched to CSC X-ray sources with a search radius of 2 arc-seconds to 2 http://sites.google.com/site/mamajeksstarnotes/ basic-astronomical-data-for-the-sun http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 2

Alfvin: Cold Gas and Black Hole Activity 3 optical images in Figure 4. Figure 2. Histogram of Chandra X-ray coverage (open) and nuclear X-ray source detection (filled) as a function of distance. We impose a cut of d < 50 Mpc to maintain a high detection fraction. obtain only nuclear matches. The X-ray luminosities were calculated from CSC full-band source fluxes using HyperLeda distance moduli for each source as simply L X = 4πd 2 f X (i.e., no k-correction is required for these low-redshift objects) and are also expressed as logarithms. For sources that had Chandra coverage but were not detected in the CSC, we retrieved the sensitivities at those positions. The sensitivities were then used as upper limits in the fitting. We also obtained nuclear X-ray luminosities from the catalog compiled by Liu (2011). The agreement between the CSC and Liu luminosities is also generally good (Figure 3, right), with offset median/mean/standard deviation for 35 objects of 0.11/ 0.10/0.52. To verify the accuracy of the X-ray flux measurements, we reduced 15 sources by hand using the Chandra Interactive Analysis of Observations (CIAO) software. After the observations were obtained from archival Chandra data, we identified the central X-ray source and extracted the flux from a circular aperture, then subtracted the background. Our hand-measured fluxes closely matched those given in the catalog search. Our initial optical selection went out to a distance of 90 Mpc. We examined the detection fraction of X-ray nuclear sources as a function of distance (Figure 2) to see how deep we could reasonably make our sample. A detection refers to a spiral in this Hyperleda sample that is also classified as a Chandra source. We noticed that there is a significant decline in the detection fraction beyond 50 Mpc and therefore imposed d < 50 Mpc for further analysis of X-ray properties. The resulting optical sample contains spiral galaxies with B-band optical magnitudes < 18 within 50 Mpc. Our sample of spirals with Chandra coverage and HI data consists of 136 galaxies, of which 76 (56%) are detected in X-rays. Table 1 lists the optical properties, HI mass, and nuclear X-ray luminosity or limit for each galaxy. A subset of our galaxies are shown with SDSS 3. CORRELATION ANALYSIS We illustrate the relationship between stellar mass and cold gas fraction (of the baryonic mass), here calculated as M HI /(M HI + M star ), in Figure 5. The best-fit linear regression slope, calculated robustly against outliers, is 0.5 ± 0.02. Our data therefore show the expected trend that optically brighter and more massive spirals possess more HI gas in absolute terms (because the slope is > 1) but have lower gas fractions (because the slope is < 0). The correlation is extremely significant with p < 0.001. This trend has been quantified using α40 ALFALFA data by Huang et al. (2012, 2014; see also Maddox et al. 2015) and we here illustrate it with the new α70 catalog for nearby spirals. We also investigated the relationship between M HI and M B. First we fit M HI versus M B, then we fit M B versus M HI. Here we calculated the best-fit trend using the bisector slope since it is not clear which direction the causal relationship, if any, should go (Isobe et al. 1990). We find a strong correlation, but this is not surprising since M B is used to calculate M star, so this is another perspective on the gas fraction trend noted above. 3.1. X-ray binary contamination Our main caveat is that it is difficult to reduce the chance of high-mass X-ray binary contamination in latetype galaxies, even with the outstanding angular resolution of Chandra. The X-ray surface brightness profiles, and the independent classification of Liu (2011), suggest that at least our most luminous sources are indeed dominated by emission produced by the SMBH. First, we examined X-ray surface brightness profiles from several deep Chandra observations of relatively nearby galaxies (Figure 6) to map X-ray contamination as a function of radial distance. These profiles suggested that we are able to identify nuclear sources above the radial X-ray binary contamination out to about 50 Mpc. This does not guarantee that the central source is associated with the SMBH, but the likelihood increases with L X > 40, which many of our galaxies surpass. For each galaxy, we estimate potential X-ray binary contamination within the Chandra aperture. Low-mass X-ray binaries (LMXBs) have a luminosity function that scales with stellar mass (Gilfanov 2004), whereas highmass X-ray binaries (HMXBs) have a luminosity function that scales with star formation rate (Grimm et al. 2003). This is because LMXBs are fed by Roche-lobe overflow from their low-mass long-lived companion star whereas HMXBs also accrete from the winds of their high-mass short-lived companion star. We calculate the stellar mass within the 2 aperture from which the X-ray luminosity was extracted by treating the light distribution within the central region of each galaxy as following a simple Sersic profile, using the formalism given by Ciotti & Bertin (1999). The effective radius for this inner Sersic component is fixed to a typical value and the index is taken as n = 3.0. Then the overall X-ray luminosity from LMXBs within the X-ray aperture is estimated by integrating the luminosity function given in Gilfanov (2004, their Table 3, using the parameters for all types which are similar to the values for late types) Published by DigitalCommons@Macalester College, 2015 3

Macalester Journal of Physics and Astronomy, Vol. 3, Iss. 1 [2015], Art. 1 4 Figure 3. Left: HI masses from Hyperleda against new ALFALFA HI measurements. Right: X-ray luminosity from Liu (2011) against values from the Chandra Source Catalog. The red line indicates equal measurement values. L X = 9.05 10 28 M star + 1.62 10 39 SF R, (4) with M star and SF R in units of M and M /yr, respectively. We check that this empirical prescription provides similar results to integrating the LMXB and HMXB luminosity functions separately. We are in the process of calculating the contamination probabilities from LMXBs or HMXBs for all sources with a nuclear X-ray source detected by Chandra. Because most of the nuclear X-ray luminosities are greater than those typically reached by LMXBs or even HMXBs, and because the Chandra PSF out to our adopted distance limit only encloses the very central region of each galaxy, these contamination probabilities are expected to be generally small. However, they will be fully taken into account in our correlation analysis, as described below. Figure 4. SDSS optical cutouts of 20 representative galaxies in our sample with ALFALFA HI measurements. All show clear spiral morphology. normalized to the stellar mass. Finally the LMXB X-ray luminosity is compared to the measured L X to determine the likelihood of LMXB contamination. The overall X-ray luminosity from HMXBs within the X-ray aperture is estimated by integrating the luminosity function given in Grimm et al. (2003, their equation 6) normalized to the star formation rate. Finally the HMXB X-ray luminosity is compared to the measured L X to determine the likelihood of HXMB contamination. A simplified empirical recipe for these steps is given by Lehmer et al. (2008, 2010), who parameterize the total X-ray luminosity in binaries as a function of both M star and SFR. They find 3.2. X-ray fitting In Figure 7, we plot nuclear X-ray luminosity versus steller mass and HI mass. Chandra detections and upper limits are indicated with plus symbols and arrows, respectively. We use the Bayesian linear regression code of Kelly (2007; implemented in IDL as linmix err) to determine the correlation slope and intercept, taking both measurement uncertainties and the X-ray upper limits into account. For completeness, we also show the regression results for the galaxies with nuclear X-ray detections only. We report the preferred model parameters as the median of 10000 draws from the posterior distribution, with uncertainties corresponding to 1σ from the 16th and 84th percentiles. For the stellar mass versus nuclear X-ray luminosity, fit as log L X 39 = α + β (log M star 10.5) (5) the best-fit intercept, slope, and intrinsic scatter are α = 0.06 ± 0.12, β = 1.47 ± 0.23, and σ 0 = 1.11 ± 0.11. While the intrinsic scatter is large, there is a highly http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 4

Alfvin: Cold Gas and Black Hole Activity 5 Figure 5. Comparison of HI to stellar mass, with fractional M HI/(M HI + M star as a function of M star. There is a decline in the gas fraction of late-type galaxies with increasing stellar mass (e.g., Huang et al. 2012; Maddox et al. 2015). The solid red line shows the robust best-fit linear regression model. Figure 6. X-ray surface brightness profiles for three example galaxies. NGC 1316 (Fornax A) is a complex quasi-lenticular with radio jets and a weak nuclear X-ray source, excluded from our sample. NGC 4579 is an Sb and NGC 4594 is an Sa spiral; both are also Seyfert 2 galaxies and contain a nuclear X-ray source powered by low-level supermassive black hole activity. The vertical line shows 1 at 50 Mpc. Published by DigitalCommons@Macalester College, 2015 5

Macalester Journal of Physics and Astronomy, Vol. 3, Iss. 1 [2015], Art. 1 6 significant correlation between the variables, in the sense that optically-brighter and more massive spiral galaxies also have greater L X values. For the HI mass versus nuclear X-ray luminosity, fit as log L X 39 = α + β (log M HI 9.5) (6) the best-fit intercept, slope, and intrinsic scatter are α = 0.15 ± 0.15, β = 0.82 ± 0.30, and σ 0 = 1.35 ± 0.14. The intrinsic scatter is quite large, but there is tentative evidence for a positive correlation between HI mass and nuclear X-ray luminosity; the slope is greater than zero at the 2.5σ level. However, we then fit L X as a function of both M star and M HI and found that the dependence upon M HI in this joint fit has a slope of 0.01 ± 0.28, consistent with zero. Further, the residuals of L X L X (M star ) show no trend with M HI. The tentative relationship between L X and M HI appears to be entirely driven by M star. 4. DISCUSSION AND CONCLUSIONS 4.1. Multiwavelength activity indicators We were able to use indicators of a number of wavelengths to distinguish black hole over non-black hole activity. X-ray, optical and radio emissions can be of use to probe this. First, several galaxies are known Seyferts or have variable X-ray emission, which again improves the chances that we are probing SMBH activity rather than high-mass X-ray binary contamination. Seyferts have strong optical features of black hole activity that include strong emission lines as well as absorption lines from the host galaxy. Seyferts 1 and 2 are thought to be the same objects seen from different vantage points. Seyfert 1s have broadened emission lines while Seyfert 2s do not because of our viewing angle. Table 1 shows optical indicators of black hole activity for selected sources under notes. 1.4 GHz measurements would indicate that there is significant radio emission from the source. This likely means that the source is a black hole and not a high-mass X-ray binary since black holes produce strong continuum radio emission. This data will be collected in ongoing work. The average X-ray luminosities of the galaxies identified as Seyferts are larger than those for which no optical classification is given in NED (taking X-ray limits into account, using ASURV; Isobe et al. 1990). This may suggest that non-seyferts accrete at a slower rate then Seyferts. On the other hand, the gas fractions of the Seyfert galaxies are similar to those of the non-seyferts at matched M star values. This also suggests that accretion rate is not dependent on gas fraction alone. from that in a matched comparison sample of galaxies, or depend on the O III/M BH accretion indicator (Fabello et al. 2011). Our results are consistent with this finding, to the extent that X-ray luminosities provide a complementary (if not entirely overlapping) view of nuclear activity. The lack of a strong dependence of L X upon M HI in our analysis suggests that the galaxy-wide cold gas does not directly influence low-level supermassive black hole activity in these spirals. We also investigated the relationship between the presence of bar features (barred spirals) and the HI and X- ray luminosities. A bar is believed to funnel gas down to the central regions of the galaxy by decreasing its angular momentum, so we might expect barred galaxies to have greater nuclear L X values at a given HI mass. However, we did not find any significant difference in the average nuclear X-ray luminosities of late-type galaxies with versus without bars. This suggests that bars may not channel a significant amount of gas directly into the black hole or that there is some kind of balance between feeding and accretion. 4.3. Future work Future work could use optical activity, X-ray variability and hardness ratios to increase the confidence of black hole identification. Formalization of 1.4 GHz continuum radio measurements could also be used to increase the confidence of black hole presence. A larger sample of spirals spanning a greater dynamic range in HI mass distribution could test for more subtle dependencies of SMBH activity on cold gas. 4.2. SFR, cold gas, and SMBH activity Our fits test for a relationship between black hole accretion rate and HI content. Since cold gas is available for star formation, we might expect a positive trend based on the correlation between star formation rate and active galactic nuclei fraction found by Rafferty et al. (2011). On the other hand, Fabello et al. (2011a, 2012) find that bulge properties and SMBH accretion rate are not linked to HI gas in blue galaxies (although they do appear to correlate in red galaxies). Specifically, the gas fraction in galaxies that host active galactic nuclei does not differ http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 6

Alfvin: Cold Gas and Black Hole Activity 7 Figure 7. Top: X-ray luminosity versus stellar mass (left) and versus HI mass (right). The dashed lines are fits to detections only (plus symbols), while the solid lines are fits including the X-ray upper limits (down arrows). Measurement uncertainties of 0.2 dex and 0.1 dex were assumed for L X and both M star and M HI, respectively. Linear regression was performed using the Bayesian methodology of Kelly (2007). Bottom: Residuals after subtracting the best-fit relations from the top panels. For reference, zero is marked with a solid horizontal red line. Published by DigitalCommons@Macalester College, 2015 7

Macalester Journal of Physics and Astronomy, Vol. 3, Iss. 1 [2015], Art. 1 8 REFERENCES Bizzocchi, L., Filho, M. E., Leonardo, E., et al. 2014, ApJ, 782, 22 Davis, B. L., Berrier, J. C., Johns, L., et al. 2014, ApJ, 789, 124 Evans, I. N., Primini, F. A., Glotfelty, K. J., et al. 2010, ApJS, 189, 37 Isobe, T., Feigelson, E. D., Akritas, M. G., et al. 1990, ApJ, 364, 104 Fabello, S., Kauffmann, G., Catinella, B., et al. 2011, MNRAS, 416, 1739 Fabello, S., Kauffmann, G., Catinella, B., et al. 2012, MNRAS, 427, 2841 Gallo, E., Treu, T., Marshall, P. J., et al. 2010, ApJ, 714, 25 Giovanelli, R., Haynes, M. P., Kent, B. R., et al. 2005, AJ, 130, 2598 Grossi, M., di Serego Alighieri, S., Giovanardi, C., et al. 2009, A&A, 498, 407 Gültekin, K., Richstone, D. O., Gebhardt, K., et al. 2009, ApJ, 698, 198 Huang, S., Haynes, M. P., Giovanelli, R., & Brinchmann, J. 2012, ApJ, 756, 113 Huang, S., Haynes, M. P., Giovanelli, R., et al. 2014, ApJ, 793, 40 Haynes, M. P., Giovanelli, R., Martin, A. M., et al. 2011, AJ, 142, 170 Jenkins, L. P., Brandt, W. N., Colbert, E. J. M., et al. 2011, ApJ, 734, 33 Table 1 Sample Properties Kelly, B. C. 2007, ApJ, 665, 1489 Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511 Maddox, N., Hess, K. M., Obreschkow, D., Jarvis, M. J., & Blyth, S.-L. 2015, MNRAS, 447, 1610 Makarov, D., Prugniel, P., Terekhova, N., Courtois, H., & Vauglin, I. 2014, A&A, 570, AA13 Mathur, S., Ghosh, H., Fiore, F., & Ferrarese, L. 2010, X-ray Astronomy 2009; Present Status, Multi-Wavelength Approach and Future Perspectives, 1248, 241 McConnell, N. J., & Ma, C.-P. 2013, ApJ, 764, 184 Miller, B., Gallo, E., Treu, T., & Woo, J.-H. 2012, ApJ, 747, 57 Oosterloo, T., Morganti, R., Crocker, A., et al. 2010, MNRAS, 409, 500 O Sullivan, E., Zezas, A., Vrtilek, J. M., et al. 2014, arxiv:1407.7543 Rafferty, D. A., Brandt, W. N., Alexander, D. M., et al. 2011, ApJ, 742, 3 Reines, A. E., Sivakoff, G. R., Johnson, K. E., & Brogan, C. L. 2011, Nature, 470, 66 Satyapal, S., Secrest, N. J., McAlpine, W., et al. 2014, ApJ, 784, 113 Shankar, F. 2009, NAR, 53, 57 Simmons, B. D., Lintott, C., Schawinski, K., et al. 2013, MNRAS, 429, 2199 Tzanavaris, P., Gallagher, S. C., Hornschemeier, A. E., et al. 2014, ApJS, 212, 9 Name RA Dec Dist M B M HI M star L X Notes (deg) (deg) (Mpc) (M ) (M ) (erg/s) NGC4501 187.996950 14.420140 34.7 23.1 9.9 11.6 40.1 Sy NGC4254 184.706550 14.416410 36.3 22.6 10.4 11.2 <39.3 NGC4594 189.997650 11.623010 16.1 22.3 8.8 11.7 40.6 Sy NGC4321 185.728200 15.821890 24.8 22.1 9.8 11.3 39.4 B NGC4450 187.123350 17.084990 29.9 21.9 9.0 11.3 40.6 NGC4303 185.478900 4.474180 23.6 21.9 10.0 10.9 39.4 SyB NGC5194 202.469550 47.195150 10.9 21.9 9.6 11.0 39.5 Sy NGC4725 192.610950 25.500910 20.1 21.8 10.1 11.2 39.5 SyB NGC1365 53.401650 36.140510 20.7 21.8 10.0 11.0 40.6 SyB NGC4579 189.431400 11.818000 23.6 21.7 9.1 11.3 41.4 SyB UGC01378 29.080050 73.282780 46.3 21.7 10.2 10.9 <39.7 B NGC7331 339.267300 34.415620 14.9 21.7 9.9 11.2 39.1 NGC3256 156.963750 43.903760 38.2 21.6 9.8 10.7 40.8 NGC2276 111.804450 85.751390 39.3 21.6 9.8 10.6 39.0 B NGC4258 184.740000 47.303880 10.3 21.6 9.9 10.9 40.0 SyB NGC5728 220.599600 17.253010 40.5 21.5 9.5 11.2 41.1 SyB NGC1068 40.669650 0.013240 16.1 21.5 9.2 11.1 41.6 Sy NGC0988 38.866050 9.356420 20.7 21.4 9.4 10.4 <38.8 B NGC1097 41.579550 30.274970 16.1 21.3 9.8 11.0 40.8 B IC2560 154.077750 33.563870 40.0 21.3 10.2 10.8 40.8 SyB NGC3627 170.062500 12.990990 12.0 21.3 9.1 10.9 38.7 B NGC5055 198.955500 42.029220 10.9 21.3 9.8 10.9 38.9 NGC5457 210.802500 54.349060 8.0 21.2 10.0 10.5 38.3 B NGC5427 210.858750 6.030590 39.3 21.2 10.2 10.7 39.7 Sy NGC3521 166.452450 0.035780 12.0 21.2 9.9 11.0 38.8 B NGC4559 188.990400 27.959680 14.3 21.2 10.1 10.4 39.1 B NGC4666 191.286000 0.461860 23.0 21.1 9.7 10.8 39.4 B NGC6643 274.943100 74.568390 25.9 21.1 9.6 10.6 <38.7 NGC3631 170.262000 53.169940 20.1 21.0 9.5 10.6 39.4 NGC5643 218.169600 44.174480 15.5 21.0 9.3 10.6 40.1 SyB NGC4736 192.720750 41.119990 8.0 21.0 8.8 10.9 39.2 Sy NGC6764 287.068200 50.933190 39.3 21.0 9.6 10.6 40.1 B NGC2841 140.510250 50.976830 12.6 21.0 9.7 11.0 39.1 NGC3368 161.690400 11.819810 13.7 21.0 9.5 11.0 39.1 B NGC5033 198.364650 36.593710 16.1 20.9 9.9 10.5 41.1 Sy NGC5350 208.340100 40.363940 37.0 20.9 9.8 10.9 39.9 B NGC2782 138.521100 40.113770 40.5 20.9 9.6 10.7 40.6 B NGC4826 194.182500 21.681950 8.6 20.9 8.8 10.9 <38.5 Sy NGC5970 234.625050 12.186110 30.6 20.9 9.7 10.7 <39.1 B NGC6500 268.999050 18.338290 46.3 20.9 9.8 10.6 40.6 NGC7727 354.974400 12.292840 25.9 20.8 8.4 11.0 38.6 B NGC4689 191.939850 13.762720 25.4 20.8 9.1 10.6 <39.0 NGC3568 167.702400 37.447790 33.0 20.8 9.9 10.5 39.5 B NGC3556 167.878950 55.674150 13.7 20.7 9.7 10.6 <37.4 B NGC0628 24.174000 15.783320 10.3 20.7 10.0 10.4 38.2 NGC7582 349.596750 42.369990 20.7 20.7 9.5 10.7 40.7 SyB NGC5678 218.023350 57.921420 31.8 20.7 9.6 10.5 39.3 B NGC4654 190.985850 13.126570 16.7 20.7 9.5 10.4 39.0 B http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 8

Alfvin: Cold Gas and Black Hole Activity 9 Table 1 Sample Properties NGC3898 177.313650 56.084030 20.7 20.7 9.6 11.0 39.4 NGC0470 19.936800 3.409870 34.0 20.6 9.6 10.7 <41.0 NGC7314 338.941650 26.050500 19.5 20.6 9.4 10.4 <38.4 SyB NGC5656 217.606200 35.320910 48.8 20.5 9.5 10.4 <38.0 NGC7552 349.044900 42.584280 20.7 20.5 9.5 10.6 39.5 B NGC7714 354.058800 2.155020 40.5 20.5 9.9 10.2 41.0 B NGC4750 192.530100 72.874470 27.7 20.5 9.2 10.4 40.4 NGC3169 153.562050 3.466360 18.4 20.4 9.8 10.8 40.5 Sy ESO386 039 224.105550 37.600920 39.3 20.4 9.5 10.4 <39.8 B NGC5954 233.646000 15.199390 29.9 20.4 9.3 10.4 39.0 B NGC2139 90.283350 23.672410 24.2 20.3 9.7 9.9 <38.1 B NGC7716 354.131100 0.297290 37.0 20.3 9.7 10.5 <39.7 B NGC3395 162.458850 32.982830 25.9 20.2 9.5 9.8 <38.6 B NGC4414 186.612600 31.223390 13.7 20.2 9.3 10.7 39.2 NGC4274 184.960650 29.614290 16.1 20.2 8.6 10.8 <39.0 B NGC0925 36.820350 33.579000 9.7 20.2 9.6 10.1 38.5 B NGC4151 182.635950 39.405790 17.8 20.2 9.4 10.5 41.5 SyB NGC3718 173.145150 53.067870 18.4 20.1 9.9 10.6 40.8 B NGC3310 159.689850 53.501740 17.8 20.1 9.5 9.9 40.1 B NGC1367 53.755350 24.933230 19.0 20.1 9.6 10.7 40.3 B NGC3938 178.205700 44.120800 14.9 20.1 9.5 10.1 <38.5 NGC2906 143.026050 8.441590 33.4 20.0 9.1 10.2 <39.2 NGC1073 40.918800 1.375700 17.2 20.0 9.7 10.0 39.1 B NGC4470 187.407450 7.823900 35.2 20.0 9.3 10.2 <39.6 NGC3887 176.769150 16.854710 16.7 20.0 9.3 10.2 <38.6 B NGC3184 154.570800 41.424360 11.5 20.0 9.3 10.2 38.2 B NGC3683 171.882600 56.876990 28.2 20.0 9.6 10.2 39.6 B NGC3982 179.117250 55.124930 20.1 19.9 9.2 10.2 39.3 SyB NGC4647 190.885500 11.582210 21.9 19.9 8.9 10.2 <39.1 B NGC1055 40.437300 0.443180 13.7 19.8 9.7 10.4 <38.4 B NGC2993 146.451300 14.368380 32.4 19.8 9.5 9.8 40.2 NGC7424 344.326500 41.070640 11.5 19.8 9.7 9.9 <37.6 B NGC3344 160.879650 24.922150 10.3 19.8 9.6 10.1 <38.7 B IC5325 352.180950 41.333390 19.5 19.8 8.9 10.1 <38.7 B NGC3389 162.116400 12.533130 20.1 19.7 9.4 9.8 39.1 NGC4651 190.927650 16.393390 13.7 19.7 9.4 10.1 <39.2 NGC5347 208.324200 33.490800 37.5 19.7 9.6 10.4 40.5 SyB NGC5068 199.728900 21.039050 9.2 19.7 9.3 10.1 <38.0 B NGC4498 187.914900 16.852790 23.6 19.7 9.2 10.1 <39.3 B NGC0255 11.947050 11.468730 22.4 19.6 9.6 9.8 <38.9 B NGC3314A 159.303450 27.683760 40.0 19.5 9.3 10.1 <39.5 NGC1703 73.217250 59.742190 18.4 19.5 9.2 9.9 <38.8 B NGC4772 193.371600 2.168330 16.1 19.4 8.9 10.5 39.9 NGC0685 26.928300 52.761790 16.1 19.4 9.4 9.8 <38.6 B NGC7741 355.976850 26.075470 13.2 19.4 9.2 9.8 <38.6 B NGC3066 150.545400 72.125250 34.7 19.4 9.2 10.0 39.5 B IC0239 39.115950 38.969030 15.5 19.4 9.6 10.1 <38.6 B NGC4102 181.597800 52.710970 16.1 19.4 8.8 10.0 40.3 B NGC5774 223.427100 3.582530 24.2 19.4 9.8 9.9 <38.0 B NGC1640 70.560600 20.434740 20.7 19.3 8.9 10.2 39.8 B NGC5240 203.979900 35.588250 35.8 19.3 9.0 10.0 <38.6 B NGC4639 190.718250 13.257040 16.1 19.2 9.1 10.1 41.0 SyB IC0396 74.495700 68.323520 16.7 19.2 8.3 9.9 38.9 NGC4548 188.860200 14.496010 9.2 19.2 8.3 10.2 39.3 B NGC0672 26.975100 27.432200 8.0 19.2 9.3 9.7 <38.3 B NGC1341 51.993300 37.150060 23.6 19.1 8.6 9.7 <38.9 B NGC4904 195.244350 0.027520 18.4 19.1 9.1 9.9 <39.0 B NGC4900 195.163350 2.501010 14.9 19.1 9.0 9.8 38.5 B NGC3507 165.855750 18.135990 15.5 19.1 9.1 9.9 39.1 B NGC4492 187.748850 8.077690 26.6 19.1 8.0 9.9 <39.0 PGC052935 222.387300 10.164170 27.0 19.0 8.5 9.8 <38.0 NGC4355 186.727650 0.877560 31.8 19.0 8.8 9.8 38.9 SyB NGC3185 154.410750 21.688320 19.5 18.9 8.5 10.1 39.3 SyB ESO184 060 290.672250 54.585270 41.1 18.8 9.2 9.8 <38.9 B NGC0991 38.885850 7.155020 21.3 18.8 9.2 9.8 <38.8 B NGC4448 187.064400 28.620310 12.6 18.8 7.9 10.3 <38.8 B NGC5585 214.950750 56.729440 8.6 18.7 9.2 9.4 <38.0 B UGC09235 216.175650 35.267150 47.0 18.7 8.8 9.7 <39.8 B NGC1058 40.875450 37.341100 9.7 18.7 9.1 9.6 38.2 Sy NGC1493 59.364450 46.210900 12.0 18.7 9.0 9.6 38.8 B NGC2541 123.667200 49.061580 10.9 18.6 9.4 9.3 <38.6 NGC4411B 186.696750 8.884500 19.5 18.6 9.3 9.7 <38.8 B UGC11466 295.744800 45.298110 16.1 18.6 9.8 9.7 <38.9 NGC1637 70.367400 2.858090 9.2 18.5 9.0 9.7 38.6 B NGC0949 37.702950 37.136590 10.9 18.5 8.6 9.5 <38.5 NGC4411A 186.625050 8.871670 20.1 18.5 9.2 9.8 <38.8 B UGC08041 193.802700 0.116680 20.1 18.5 9.3 9.6 <38.8 B Published by DigitalCommons@Macalester College, 2015 9

Macalester Journal of Physics and Astronomy, Vol. 3, Iss. 1 [2015], Art. 1 10 Table 1 Sample Properties UGC05707 157.809300 43.137410 43.5 18.4 9.6 9.6 <39.4 B NGC5474 211.256400 53.662240 8.0 18.4 9.1 9.4 <37.2 NGC7320 339.014250 33.948160 14.3 18.4 8.5 9.4 39.8 NGC4136 182.323800 29.927670 11.5 18.4 9.1 9.6 38.0 B UGC06930 179.323050 49.283720 14.9 18.4 9.2 9.6 <38.5 B IC5332 353.614500 36.101450 8.0 18.3 9.2 9.6 <37.0 NGC2082 85.462950 64.301080 13.7 18.3 8.6 9.4 <38.4 B NGC4713 192.491100 5.311290 10.9 18.2 9.1 9.3 38.2 B NGC3913 177.662250 55.353900 17.8 18.1 8.9 9.4 <38.7 NGC4020 179.736150 30.412470 13.7 18.1 8.7 9.5 <39.0 B NGC2500 120.471600 50.737230 10.3 18.1 8.8 9.4 <38.5 B Note. HI mass, stellar mass, and X-ray luminosity are expressed as logarithms. Notes are Sy for Seyfert and B for barred, from NED. http://digitalcommons.macalester.edu/mjpa/vol3/iss1/1 10