“…To do this, we applied Perple_X thermodynamic modeling and free energy minimization software (Connolly, 2009(Connolly, , 2005 to model major-element chemistry and major-phase equilibria of a cooling and differentiating magma using the internally consistent thermodynamic database of Holland and Powell (2011), solid solution models based on the thermodynamic data of Tajčmanová et al (2009), Green et al (2016), and Jennings and Holland (2015), and the melt model of Holland et al (2018). After running major-element calculations with Perple_X, we assessed accessory phase stability using saturation models for sphene (Ayers et al, 2018), zircon (Boehnke et al, 2013), apatite (Harrison and Watson, 1984;Bea et al, 1992), monazite (Montel, 1993), and xenotime (Rusiecka and Baker, 2019), though the latter did not saturate. Finally, we partitioned trace elements among all major and accessory phases using partition coefficients compiled and averaged from the literature (including via the GERM database, https:// earthref.org /GERM) in the registered Julia package StatGeochem.jl (doi.org/10.17605/osf.io/tjhmw), which was also used to automate our Perple_X calculations (code available in File 4, see footnote 1).…”