In "Interfacial structures and acidity constants (pKa) {{editor, shouldn' the a be a subscript?}} of goethite from first principles molecular dynamics siluations," authors Y. Zang, X. Lui, J. Cheng, and X. Lu (Zhang et al. 2021) apply First Principles molecular dynamics (FPMD, also called Density Functional Theory MD, DFT/MD, or ab initio MD, AIMD) to evaluate the complete set of acidy constants (pKa) of the hydroxyl groups on the most prominent goethite crystal facets. The pKa of these OH and OH+ (Sprik 2000) groups are compared with available data from the multisite complexation (MUSIC) model traditionally used to estimate pKa on mineral surfaces. The authors have presented eloquent rational for the importance and implications of understanding goethite acidity constants in room temperature geochemistry settings. Here I focus on the computational aspects, the strengthes of FPMD, and its possibilities. pK calculations represent a significant triumph of the FPMD free energy formulation. Many FPMD simulations in the literature are conducted in a non-equilibrium mode and are "numerical This is the peer-reviewed, final accepted version for American Mineralogist, published by the Mineralogical Society of America. The published version is subject to change. Cite as Authors (Year) Title. American Mineralogist, in press.