Byline: Mirjam Kan Parker (1), Axel Liebscher (2,3), Dirk Frei (3,4), Jelle Sijl (1), Wim Westrenen (1), Jon Blundy (5), Gerhard Franz (2) Keywords: Orthopyroxene; Partition coefficient; Experimental determination; Computer simulations; Substitution mechanisms Abstract: Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420degC in the system [CaO--MgO--Al.sub.2]O.sub.3--SiO.sub.2 and subsystem CaO--MgO--SiO.sub.2. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu high field strength: Zr, Nb, Hf, Ta, Ti transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from D.sub.La.sup.Opx - melt 0.0005 to D.sub.Lu.sup.Opx - melt 0.109 , D values for highly charged elements vary from D.sub.Th.sup.Opx - melt 0.0026 through D.sub.Nb.sup.Opx - melt 0.0033 and D.sub.U.sup.Opx - melt 0.0066 to D.sub.Ti.sup.Opx - melt 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite ([Mg.sub.2][Si.sub.2]O.sub.6) and ferrosilite ([Fe.sub.2][Si.sub.2]O.sub.6). Calculated solution energies show that R.sup.2+ cations are more soluble in Opx than R.sup.3+ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R.sup.3+ cations by coupled substitution with Li.sup.+ on the M1 site that is energetically favoured over coupled substitution involving Al--Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r .sub.0, maximum partition coefficients D .sub.0, and apparent Young's moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r .sub.0 values for R.sup.3+ substitutions are 0.66--0.67 A for M1 and 0.82--0.87 A for M2. Simulations for enstatite result in r .sub.0 = 0.71--0.73 A for M1 and ~0.79--0.87 A for M2. Ferrosilite r .sub.0 values are systematically larger by ~0.05 A for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in r.sub.0.sup.M2 in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models. Author Affiliation: (1) Faculty of Earth and Life Science, VU University Amsterdam, De Boelelaan 1085, 1081 HV, Amsterdam, The Netherlands (2) Institut fur Angewandte Geowissenschaften, Technische Universitat Berlin, Ernst-Reuter-Platz 1, 10587, Berlin, Germany (3) GFZ German Research Centre for Geosciences, Telegrafenberg, 14473, Potsdam, Germany (4) Geological Survey of Denmark and Greenland (GEUS), Aster Voldgade 10, 1350, Copenhagen K, Denmark (5) CETSEI, Department of Earth Sciences, University of Bristol, Wills Memorial Building, Queens Road, Bristol, BS8 1RJ, United Kingdom Article History: Registration Date: 04/08/2009 Received Date: 20/04/2009 Accepted Date: 04/08/2009 Online Date: 22/08/2009 Article note: Communicated by J. Hoefs.