Active galactic nuclei with massive black holes have closer galactic neighbors

Context. The large-scale environments of active galactic nuclei (AGNs) reveal important information on the growth and evolution of supermassive black holes (SMBHs). Previous AGN clustering measurements using two-point correlation functions have hinted that AGNs with massive black holes preferentially reside in denser cosmic regions than AGNs with less massive SMBHs. At the same time, little to no dependence on the accretion rate has been found; however, the significance of such trends has been limited. Aims. Here, we present kth-nearest-neighbor (kNN) statistics of 2MASS galaxies around AGNs from the Swift/BAT AGN Spectroscopic survey. These statistics have been shown to contribute additional higher order clustering information on the cosmic density field. Methods. By calculating the distances to the nearest seven galaxy neighbors in angular separation to each AGN within two redshift ranges (0.01 < z < 0.03 and 0.03 < z < 0.06), we compared their cumulative distribution functions to that of a randomly distributed sample to demonstrate the sensitivity of this method to the clustering of AGNs. We also split the AGNs into bins of bolometric luminosity, black hole mass, and Eddington ratio (while controlling for redshift) to search for trends between kNN statistics and fundamental AGN properties. Results. We find that AGNs with massive SMBHs have significantly closer neighbors than AGNs with less massive SMBHs (at the 99.98% confidence level), especially in our lower redshift range. We find less significant trends with luminosity or Eddington ratio. By comparing our results to empirical SMBH-galaxy-halo models implemented in N-body simulations, we show that small-scale kNN trends with black hole mass may go beyond stellar mass dependencies. Conclusions. This suggests that massive SMBHs in the local universe reside in more massive dark matter halos and denser regions of the cosmic web, which may indicate that environment is important for the growth of SMBHs, bolstering prior conclusions based on correlation functions.