logo Norris Scientific Norris Scientific » Expet Tools » MetalAct
 

HomeAbout this Metal Activity Calculator

Purpose

This website calculates the activities of all elements (for which we have defined interaction parameters) in a liquid solution of molten iron at a given temperature that contains some mix of other elements. The temperature and composition of interest are defined in a data file and uploaded to this server for computation.

Attribution

This website is based on the previous work[1] of Jon Wade and Bernard Wood with implementation in PHP by Ashley Norris.

In 2016 the site was moved to Ashley's personal website norris.org.au as a result of changes made to the Oxford Earth Sciences web server that left it incompatible with MetalAct. In Jan 2024 the Australian domain registrar revoked the norris.org.au domain name resulting in the site being moved to norrisau.org. To ensure you can always access the site, please bookmark the launch page here:

http://www.earth.ox.ac.uk/~expet/metalact/

Results from this website can be cited in the following way:

MetalAct Website, http://www.earth.ox.ac.uk/~expet/metalact/, Wood, BJ;Wade, J;Norris, A, [Date of Access].

We also ask that you cite the original paper that this website is based upon:

Wade, J & Wood, BJ, "Core formation and the oxidation state of the Earth", EARTH AND PLANETARY SCIENCE LETTERS, DOI:10.1016/j.epsl.2005.05.017, 2005 [DOI Lookup]

Data Sources

Note (8th June 2012): We have noticed a correction to be made to Table 4 in Wade, Wood and Tuff (2012), "Metal/silicate partitioning of Mo and W at high pressures and temperatures: Evidence for late accretion of sulphur to the Earth". (γW0) was misprinted as "1" and should read "3", which is the value used on this website since it was created.

We apologise for any inconvenience this may have caused. If you believe that this may affect your research please do not hesitate to contact us to discuss this matter further.

Regards - Oxford Experimental Petrology Group.

The interaction parameters in liquid iron alloys are sourced from Table 2.1, p280-293 in [2] and have been converted to the εij form using Equation (20), p276, in the same text. Any data derived at a temperature other than the reference temperature (1873K) is adjusted to the reference temperature using Equation (17) from [2].

There have been a number of revisions to these values, as listed here:

  • εSiSi was revised using data from [4].

  • εSiNi, εSiCo, εSiV, εSiCr, εSiNb, εSiMo, εSiW were revised using data from [7].

The gamma zero (γi0) values for each element in liquid iron are taken from Table 1, p278 in [2] and for convenience the natural logarithm is computed. These values are then adjusted to the reference temperature (1873K) using Equation (16) from [2].

There have been a number of revisions, and additions, to these values, as listed here:

  • γH0, γN0, γO0, γP0, γS0, γGa0, γGe0, γHf0, γRe0 were estimated by Wade and Wood.

  • γSi0, γNi0, γCo0, γW0, γMo0, γCr0, γV0, γNb0 were revised using p479 from [7].

Method

The algorithm is implemented in PHP and computed in double precision by the web server at the time that the data file is uploaded. For convenience results are temporarily retained by the server, but it is encouraged that users should download the results as soon as possible.

The algorithm follows these steps:

  1. The entire data file is parsed for correctness. If errors are detected then these are reported to the user and the algorithm aborted.
  2. From the samples listed in the data file a list of unique temperatures is generated. For each unique temperature a table of interaction parameters (εij) and gamma zeros (γi0) are computed using Equations (17) and (16), respectively, from [2].
  3. Then, for each sample composition:

    1. The composition is converted to mole fraction (if required).
    2. The composition is normalised.
    3. The lnγ1 term is computed using Equation [23] from [3].
    4. For each element (i) the lnγi term is computed using Equation [24] from [3].
    5. If desired, the theoretical saturation of C is computed using the method described on p604 of [5] and the computation is then repeated from (ii.) using this adjusted composition. Convergence is typically achieved in 4 to 5 iterations.
    6. The activity coefficients (γi) are reported to the user and can be downloaded as a CSV file.

References

[1] Wade, J & Wood, BJ, "Core formation and the oxidation state of the Earth", EARTH AND PLANETARY SCIENCE LETTERS, DOI:10.1016/j.epsl.2005.05.017, 2005 [DOI Lookup]

[2] Gordon and Breach, "Steelmaking Data Sourcebook", Gordon and Breach Science Publishers, ISBN:2-88124-153-0, 1988 [ISBN Lookup]

[3] MA, Zhongting, "Thermodynamic description for concentrated metallic solutions using interaction parameters", Metallurgical and Materials Transactions B, DOI:10.1007/s11663-001-0011-0, 2001 [DOI Lookup]

[4] Bouchard and Bale, "Simultaneous optimization of thermochemical data for liquid iron alloys containing C, N, Ti, Si, Mn, S, and P", Metallurgical and Materials Transactions B, DOI:10.1007/BF02653863, 1995 [DOI Lookup]

[5] Wood, BJ, "Carbon in the core", Earth and Planetary Science Letters, DOI:10.1016/0012-821X(93)90105-I, 1993 [DOI Lookup]

[6] Malhotra, Chang, Schlesinger, "Liquid solution thermodynamics in the alumina-saturated Fe-C-Pb system", ISIJ INTERNATIONAL, DOI:10.2355/isijinternational.38.1, 1998 [DOI Lookup]

[7] Tuff, J, Wood, BJ, Wade, J, "The effect of Si on metal-silicate partitioning of siderophile elements and implications for the conditions of core formation", GEOCHIMICA ET COSMOCHIMICA, DOI:10.1016/j.gca.2010.10.027, 2011 [DOI Lookup]