We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at. Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG = 54 ± 6 M e . These results establish the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index −1.95 ± 0.51 (−1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics. Unified Astronomy Thesaurus concepts: Astrophysical black holes (98); Black holes (162); Gravitational waves (678); Particle astrophysics (96); Gravitational wave astronomy (675); Nuclear astrophysics (1129); Stellar evolution (1599); Stellar populations (1622); Population III stars (1285); Helium burning (716)