“…While the electron density, n e , is expected to be much lower than for hydrogen, n H , it has larger collisional de-excitation rates. Electron collision rates are available for S I from Tayal (2004) and Fe II from Bautista et al (2015, and Priv Comm), S+H collision rates can be estimated using the O+H rates of Lique et al (2018), and Fe+H rates are given by Hollenbach & McKee (1989). Using the thermodynamic models of Rodgers & Glassgold (1991), Harper et al (2001), or Harper et al (2017b) (where the gas temperatures are reduced close to the star) both lines are formed predominantly above the critical densities (i.e., thermalized) when T gas > 500 K. If Fe II and S I and are dominant ionization states, then the ratio of [Fe II] 25.99 µm and [S I] 25.25 µm fluxes, would be ≤ 4.5 in the optically thin limit, and lower, ∼ 1.5, when allowing for optical depth effects.…”