The Propagation of Uncertainties in Stellar Population Synthesis Modeling. III. Model Calibration, Comparison, and Evaluation

Stellar population synthesis (SPS) provides the link between the stellar and dust content of galaxies and their observed spectral energy distributions. In the present work, we perform a comprehensive calibration of our own flexible SPS (FSPS) model against a suite of data. These data include ultraviolet, optical, and near-IR photometry, surface brightness fluctuations, and integrated spectra of star clusters in the Magellanic Clouds (MCs), M87, M31, and the Milky Way (MW), and photometry and spectral indices of both quiescent and post-starburst galaxies at z {approx} 0. Several public SPS models are intercompared, including the models of Bruzual and Charlot (BC03), Maraston (M05), and FSPS. The relative strengths and weaknesses of these models are evaluated, with the following conclusions: (1) the FSPS and BC03 models compare favorably with MC data at all ages, whereas M05 colors are too red and the age dependence is incorrect; (2) all models yield similar optical and near-IR colors for old metal-poor systems, and yet they all provide poor fits to the integrated J - K and V - K colors of both MW and M31 star clusters; (3) FSPS is able to fit all of the ultraviolet data because both the post-asymptotic giant branch (post-AGB) and horizontal branch evolutionary phases are handled flexibly, while the BC03 and M05 models fail in the far-UV, and both far- and near-UV, respectively; (4) all models predict ugr colors too red, D{sub n}4000 strengths too strong, and Hdelta{sub A} strengths too weak compared to massive red sequence galaxies, under the assumption that such galaxies are composed solely of old metal-rich stars; and (5) FSPS and, to a lesser extent, BC03 can reproduce the optical and near-IR colors of post-starburst galaxies, while M05 cannot. Reasons for these discrepancies are explored. The failure at predicting the ugr colors, D{sub n}4000, and Hdelta{sub A} strengths can be explained by some combination of a minority population of metal-poor stars, young stars, blue straggler and/or blue horizontal branch (HB) stars, but not by appealing to inadequacies in either theoretical stellar atmospheres or canonical evolutionary phases (e.g., the main-sequence turnoff). The different model predictions in the near-IR for intermediate age systems are due to different treatments of the thermally pulsating asymptotic giant branch stellar evolutionary phase. We emphasize that due to a lack of calibrating star cluster data in regions of the metallicity-age plane r