Rebuilding the theory of isotope fractionation for evaporation of silicate melts under vacuum condition

Isotope effects are pivotal in understanding silicate melt evaporation and planetary accretion processes. Based on the Hertz–Knudsen equation, the current theory often fails to predict observed isotope fractionations of laboratory experiments due to its oversimplified assumptions. Here, we point out that the Hertz-Knudsen-equation-based theory is incomplete for silicate melt evaporation cases and can only be used for situations where the vaporized species is identical to the one in the melt. We propose a new model designed for silicate melt evaporation under vacuum. Our model considers multiple steps including mass transfer, chemical reaction, and nucleation. Our derivations reveal a kinetic isotopic fractionation factor（KIFF or α） αourmodel= [ m（ 1 species）/m（ 2 species）] 0.5, where m（species） is the mass of the reactant of reaction/nucleation-limiting step or species of diff usion-limiting step and superscript 1 and 2 represent light and heavy isotopes, respectively. This model can eff ectively reproduce most reported KIFFs of laboratory experiments for various elements, i.e., Mg, Si, K, Rb, Fe, Ca, and Ti. And, the KIFF-mixing model referring that an overall rate of evaporation can be determined by two steps jointly can account for the eff ects of low PH2 pressure, composition, and temperature. In addition, we find that chemical reactions, diff usion, and nucleation can control the overall rate of evaporation of silicate melts by using the fitting slope in ln（-ln f） versus ln（ t）. Notably, our model allows for the theoretical calculations of parameters like activation energy （Ea）, providing a novel approach to studying compositional and environmental eff ects on evaporation processes, and shedding light on the formation and evolution of the proto-solar and Earth-Moon systems.