Exploration of quantitative kinetic models for the evaporation of silicate melts in vacuum and in hydrogen

被引:35
作者
Alexander, CMO [1 ]
机构
[1] Carnegie Inst Sci, Dept Terr Magnetism, Washington, DC 20015 USA
关键词
D O I
10.1111/j.1945-5100.2001.tb01870.x
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Two basic approaches (pure component reference (PCR) and equilibrium reference (EQR)) to modeling silicate melt evaporation are explored. The PCR model calculates the maximum possible evaporation rates of the pure oxides from their equilibrium vapor pressures and rescales these rates according to the activities of the oxides in the silicate melts and the melt densities. The EQR model calculates the maximum possible evaporation rates based on the equilibrium vapor pressures of the melts. Differences between the calculated and experimentally determined evaporation rates are accounted for with evaporation (alpha (evap)) coefficients that are only dependent on temperature. Two versions of the PCR model, Cases 1 and 2, are explored to try to resolve apparently contradictory conclusions about the composition of the evaporating species based on Mg and Si isotope fractionation during evaporation (species are not in thermodynamic equilibrium proportions) and direct measurements of gas species in Langmuir experiments (species are in roughly equilibrium proportions). The Case 2 and EQR models cannot explain the observed isotope fractionations unless evaporation occurred under non-Rayleigh conditions, either because there was significant recondensation during the experiments or because diffusion was playing a limiting role. Whether or not the role of diffusion is included, the PCR and EQR models are able to reproduce the elemental results of evaporation experiments of "chondritic" melts from temperatures of 1700 to 2000 degreesC, and up to mass losses of about 95%. However, the models underestimate absolute evaporation rates in very Ca- and Al-rich melts. This may reflect errors in the model used to estimate oxide activities. The EQR model can only reproduce the observed evaporation behavior of Na if, unlike the other oxides, its alpha (evap) coefficient is close to unity. Based on available diffusion data, diffusion is not slow enough in "chondritic" or forsteritic melts to explain the isotopic fractionations of Mg and O in the evaporation experiments, but it may play a role in limiting Si isotope fractionation. Provided recondensation was not a significant factor in the experiments, at present PCR Case 1 appears to be the best model if both the Langmuir and the isotopic fractionation experiments are to be explained.
引用
收藏
页码:255 / 283
页数:29
相关论文
共 76 条
[1]  
Alexander C. M. O'D., 1996, CHONDRULES PROTOPLAN, P233
[2]   TRACE-ELEMENT DISTRIBUTIONS WITHIN ORDINARY CHONDRITE CHONDRULES - IMPLICATIONS FOR CHONDRULE FORMATION CONDITIONS AND PRECURSORS [J].
ALEXANDER, CMO .
GEOCHIMICA ET COSMOCHIMICA ACTA, 1994, 58 (16) :3451-3467
[3]  
Alexander CMO, 1998, METEORIT PLANET SCI, V33, pA9
[4]   The lack of potassium-isotopic fractionation in Bishunpur chondrules [J].
Alexander, CMO ;
Grossman, JN ;
Wang, J ;
Zanda, B ;
Bourot-Denise, M ;
Hewins, RH .
METEORITICS & PLANETARY SCIENCE, 2000, 35 (04) :859-868
[5]  
[Anonymous], THERMODYNAMIC PROPER
[6]  
ASH RD, 1998, LUNAR PLANET SCI, V29
[7]   MASS DEPENDENCE OF THE SELF-DIFFUSION COEFFICIENTS IN 2 EQUIMOLAR BINARY-LIQUID LENNARD-JONES SYSTEMS DETERMINED THROUGH MOLECULAR-DYNAMICS SIMULATION [J].
BEARMAN, RJ ;
JOLLY, DL .
MOLECULAR PHYSICS, 1981, 44 (03) :665-675
[8]  
Berman R.G., 1983, THESIS U BRIT COLUMB
[9]   OXYGEN DIFFUSION IN AN FE-RICH BASALT MELT [J].
CANIL, D ;
MUEHLENBACHS, K .
GEOCHIMICA ET COSMOCHIMICA ACTA, 1990, 54 (11) :2947-2951
[10]   HIGH-TEMPERATURE PROPERTIES OF SILICATE LIQUIDS - APPLICATIONS TO EQUILIBRATION AND ASCENT OF BASIC MAGMA [J].
CARMICHAEL, ISE ;
NICHOLLS, J ;
SPERA, FJ ;
WOOD, BJ ;
NELSON, SA .
PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1977, 286 (1336) :373-431