The Fes Content of Sphalerite Along The Chalcopyrite-: Pyrite-Bornite Sulfur Fugacity Buffer
The Fes Content of Sphalerite Along The Chalcopyrite-: Pyrite-Bornite Sulfur Fugacity Buffer
Abstract
Concentrationsof 0.19 ñ 0.02 and 0.55 ñ 0.08 mole percent FeS have been measured
in sphalerite crystallized at 1 kbar in equilibrium with the chalcopyrite•pyrite-bornite
sulfur fugacity buffer at temperaturesof 271ø and 395øC, respectively. These data yield
a relationshipof log C----1353/T-I-1.77, where C is measuredin mole percent FeS,
T is in øK, and the pressureis 1 kb. The effect of pressureon C is evaluated. Hence-
forth, indications of equilibrium or disequilibrium crystallization and estimates of tem-
perature of formation may be made for the assemblagesphalerite-chalcopyrite-pyrite-
bornite at temperaturesin the range 150ø to 450øC.
The compositionof bornire in equilibrium with chalcopyriteand pyrite is found to be
Cus.24Fe0.s0S4.00at 271øC and Cus.02Fe0.0aS4.00 at 395øC. At thesetemperatureschalco-
pyrite and pyrite in the buffer assemblagedo not deviate measurablyfrom stoichiom-
etry. Bornite and chalcopyriteare found to incorporate0.03 to 0.18 and 0.00 to 0.43
atomicpercentZn in the range 197ø to 503øC.
Sat. indicatessolutionsaturatedat 25øC. Concentrationof NHd in runs 265 and 266 was 7.7 gm/10 mi.
Runs with consecutive numbers and identical run duration are from same reaction vessel.
Run carried out in cold-seal vessel and contained PbS.
this way were: a seriesof (Zn,Fe)S solidsolutions; tion texturesmay haveresultedfrom the relatively
FeS2; and an intimate, fine-grainedbuffer mixture slow quench. No reactionbetweenthe sulfidesand
containing roughly equal portions by weight of the platinumtubewasnoted.
CuFeS2, FeSa, and CusFeS4. Consumptionof the Crystallizationproductswere analyzed with an
FeSa componentof the buffer mixture was common ARL EMX-SM electronmicroprobeusing an ac-
in early experiments,so experimentsreported in celeratingvoltageof 15 KV and a specimen current
Table 1 typicallyused300 mg of the cpy-py-bnmix- of 2 x 104 ampereson brass. Lines analyzedand
ture, 75 mg of additionalFeSa, and 50 to 75 mg of crystals used for each element are as follows: Cu,
(Zn,Fe)S. Runs were startedwith compositions of Fe, and Zn, Ka, LiF; and S, Ka, ADP. Standards
(Zn,Fe)S that bracketed anticipated equilibrium used were: chemically analyzed, stoichiometric
compositions, so that equilibriumwas approached bornitefromButte,Montana,andchalcopyrite from
from both directions. O'okiep,Namaqualand, SouthAfrica; syntheticZnS
Reactantsand solutions (Table 1) were loaded containing0.00, 1.00, 2.00, and 5.00 mole percent
into 8 mm O.D. platinumtubes,which were flushed FeS; syntheticCuFeS•containing 0.50 weightper-
with hydrogen and flame-weldedwhile partly im- centZn; synethetic FeS containing2.00 weightper-
mersed in ice water. Groups of six sampletubes cent Cu; and a syntheticbornite of composition
were placedinto medium-sized(22.9 cm x 6.4 cm) Cus.44Fe0.82S4.oo. Reportedanalysesof sphalerites
pressurevesselswith modifiedBridgman-typeseals. represent 10-second countsat 30 to 60 discretepoints
Vesselswere loweredinto vertical,gradientfurnaces, on 5 to 10 homogeneous grains;for (Cu,Fe)-sulfide
such that nutrient and buffer phasessettledto the phases8 to 20 pointswere analyzed. All analyses
lower (hotter) end of the platinumtube and crystal- reportedare basedon directcomparison of counting
lization productsformedat the upper (cooler) end rates for unknowns and standards of similar com-
of the tube. All experimentswere conductedat 1 position, with appropriate correctionsfor background
kb pressureto promotesolutiontransportvia a rela- and instrument drift calculated
by hand. For anal-
tively dense medium. Temperature uncertainties ysesof sulfides,it is very importantto use either a
reportedin Table 1 includecalibrationof the thermo- rigorous, theoreticalcorrectionscheme,or a low
coupleand recorder,short- and long-term furnace acceleration voltage(Desborough et al., 1971), or
fluctuations,and an estimatedcorrectionfor separa- close standards. For example,at 15 KV using
tion of thermocoupleand sample (see Czamanske CusFeS4 as a standard,38.2 ratherthan34.9 weight
and Rye, 1974). Runs were quenchedby plunging percentS was obtainedfor CuFeS• by direct com-
the vesselsat pressureinto a bucketfilled with ice parison of count rates.
and water. Quenchesto room temperaturetook Many of the (Cu,Fe)-sulfide productsdisplayed
severalminutes,and some (Cu,Fe)-sulfide exsolu- exsolutiontextures,apparentlyrelatedto quenching
.1330 GERALD K. CZAMANSKE
503øC (Table 3). Considerableadditional study TAmE 3. Zinc Contents (in atomic percent)of Hydrothermal
Bornite and ChalcopyriteCrystallizedin Association
will be required to clarify relationshipsin the Cu- with Pyrite and Sphalerite
Fe-Zn-S systemin the temperaturerange 450ø to
535øC. T, øC Bornitc Chalcopyrite
0.44 +
-
0.15+.02
-Io
-12
-15
MONOCLINIC
PYRRHOTITE
/,,,/
/ / <q"ø•x
-
2O0 250 275 300 400 500 600 700 800
TEMPERATURE, øC
Fro. 1. Relationshipsbetweenthe compositionof sphalerite in mole percent FeS (light, solid lines) and pertinent uni-
variant equilibria, as a function of fugacity of sulfur and temperature at 1 bar. Several isoplethsfor avesin the area of
interest are shown as dotted lines. Measured values for mole percent FeS in sphalerite at 271ø and 395øC have been cor-
rected for the pressurechangefrom 1 kbar to 1 bar and are indicated. Methods of compilationdiscussedin text.
in the pyrrhotitefield, basedon their equation(2). (3) Isoplethsfor aveswithin the pyrite field have
The Scott-Barnesunivariant curve for pyrite • been locatedas follows: a) at 475øC, the approxi-
pyrrhotitehas beenchosenover that of Toulmin and mate crossoverpoint of the Toulmin-Barton and
Barton (1964) for threereasons:a) with the study Scott-Barnespyrite-pyrrhotite univarlant curves,
of Toulmin and Barton in hand, Scott and Barnes aveshas been calculated from the relation 2FeS -I- S•.
locateda univariant curve which is consistent,within •-2FeS•., using the expression•xG*=-70,740-I-
statederror limits, with both setsof data; b) inde- 69.04T (Barton and Skinner, 1967, table 7.2).
pendentstudiesby Giletti et al. (1968) and Schnee- Based on the calculatedspacings,isoplethswere
berg (1973) have eachsupportedthe slightlyflatter constructedparallel to the Scott-Barnesunivariant
curve (in log fs2-1/T space) of Scott and Barnes; curve. This procedureis not entirely arbitrary be-
and c) data from the presentstudyrelate better to causethe constancyof the iron content of sphal-
the Scott-Barnesunivariant curve. (2) The 'uni- erite in equilibriumwith pyrite and pyrrhotite,from
variant curve for the reactionchalcopyrite + S•.• 250ø to 550øC,stronglysuggests that the isopleths
bornite+ pyrite has beenlocatedon the basisof the of FeS activityin the pyrite field mustparallelthe
electrochemicalmeasurements of Schneeberg(1973), pyrite-pyrrhotite'univariant curve. (Minor inac-
which fit the equationlog fs•- 12.560- 1.1067x curaciesthat may exist-becausea free energy ex-
10b(i/T). Schneeberg's univariantcurve was se- pressionhas not been derived for the Scott-Barnes
lected because it is based on more extensive and univariant curve are considered to be well within
inherently more precise data than the univariant the overall uncertaintiesof this treatment.) (4) A
curve of Barton and Toulmin (1964); the two uni- Henry's law constantof 2.5 for sphaleriteof low
variant curvesare parallel and that of Barton and FeS content is obtained from Barton and Toulmin
Toulmin lies at valuesof rs2only 0.2 log unit lower. (1966) by drawinga straightline throughthe origin
THE FeS CONTENT OF SPHALERITE 1333
of their figure 14 tangent to their array of data. rite-pyrrhotitesolvusdid not agreewith measure-
Basedon this value, sphaleritecompositional points mentunlessvolumechangesin the coexistingphases
were located at 475 øC from the calculated ares values. as a functionof temperatureand pressurewere con-
Compositional isoplethswere constructed parallelto sidered. Sphaleritecompositions measured in this
the ares isopleths(and to the pyrite-pyrrhotiteuni~ studyagreereasonably with existingdata whenun-
variant curve) as requiredfor a Henry's law con- correctedfor pressure,and well when pressurecor-
stant insensitiveto compositionand temperature. rected;however,it seemsadvisableto apply correc-
(5) A field for monoclinicpyrrhotite has been tionsto high pressurewith caution.
shown schematically below its apparent upper It is appropriateto discussthe implications
of the
stabilitylimit at 248øC (Kissin and Scott, 1972). presentstudy for the sphaleritefrom the Iimori
(6) The two data points representingthe results mine,Japan,studied by Yui andCzamanske (1971),
of this study (see Table 1) have beenpressurecor- as this is at presentthe only knownapp!ication of
rected to 1 bar (see below) and have been placed data from the Cu-Fe-Zn-S systemto geothermome-
on the chalcopyrite•-bornite + pyrite univariant try. Reanalysis of the Iimori sphalerite
duringthe
curveof Schneeberg
(1973) at 271ø and 395øC. courseof the presentstudysuggests
that.0.053-
The effectof pressureuponsphaleritecomposition 0.009 molepercentFeS is probablya bettermeasure
has received attention from Barton and Toulmin of its FeS content than the value 0.044- 0.009
(1964), ScottandBarnes(1971), andScott(1973). given by Yui and Czamanske(1971). According
Along the pyrite-pyrrhotiteunivariant curve, in- to the presentstudy,0.0.53mole percentFeS in
creasedpressuremarkedlydecreases N•-•} at con- ZnS shouldcQrrespond to a formationtemperature
stant temperature. In contrast,along the chalco- of ,•170øC, if the sphaleriteis in equilibriumwith
pyrite• bornite+ pyrite univariantcurve,increased cpy-py-bnat 1 kb. Although this temperatureap-
pressureincreasesN•-V&slightly at constanttem- pearslow in comparison with publishedestimatesof
perature. Using the molar volume data of Robie conditionsduring ore formation(see Yui and Cza-
et al. (1966) for bornite and pyrite, Hall and manske, 1971), some studentsof Sanbagawameta-
Stewart (1973) for chalcopyrite,and Barton and morphismwill find it quite realistic(Yui, oral com-
Toulmin(1966) for •rvessph, AV• is calculated
to mun., 1973). The possibleeffectbf pressureon the
be----•-0.5 calfbar for the reaction5CuFeS2•- calculatedtemperatureis marked for the Iimori
CusFeS4+ 2FeS2+ [FeS](•n ,pha•erite). It can be sphaleritebecausepressuresof 5-7 kb have been
shown that postulatedduringformationof the Iimori ore. This
pressurecorrectionwould call for even lower tem-
-- AVrdP -- 2 RT d In peraturesof equilibration.
Assumingthat AVr is constant,partial integration
Acknowledgments
givesthe desiredequation,
Critical reviews by G. M. Anderson,P. B. Bar-
-/xVr (P2- P1) = 2.303 RT [log(aees)r.o ton, Jr., L. J. Cabri, J. T. Nash, and S. D. Scott
- log(a•,es)r, ]
have materiallyimprovedthe clarity and contentof
At 1 bar (the pressurefor whichthe frameworkof this manuscript,and I haveprofitedfrom numerous
Figure 1 is constructed),ares will be •80 percent stimulatingdiscussions with C. L. Christ.
of its valueat 1 kb, the pressureof experimentation. [J. S. GEOI•OGICAL SURVEY
Values for the mole fraction of FeS in sphaleriteat 345 M•DDLEFIEL•) ROAD
1 bar are.thereforeestimatedat 0.15 -- 0.02 (271øC) M•/•o PAR•C,C^t.tFOR•tA94025
and 0.44-0.07 (395øC); these values fit quite February 28, May 7, 1974
well into the framework of existing data (Fig. 1).
REFERENCES
Becauseares= (y•-•}) (N•"•), where y•.o• rep-
resentsthe Henry'slaw constantand 'NTsph
the mole
ß x FeS Barton, P. B., Jr., 1970, Sulfide petrology: Mineralog. Soc.
America Spec. Paper 3, 187-198.
fraction of FeS in ZnS, another way of considering 1973, Solid solutionsin the system Cu-Fe-S: Part 1:
the new data is to derive the Henry's law constant The Cu-S and CuFe-S joins: Eco•v. Gzor..,v. 68, p. 455-
which relates the measured ß, res values to the con- 465.
and Skinner, B. J., 1967, Sulfide mineral stabilities,
toured a•os values. The calculatedvalues of 2.5 ----- in Barnes, H. L., ed., Geochemistry and hydrothermal
0.3 and 2.6-+ 0.4 for 271ø and 395øC, respectively, or• deposits: New York, Holt, Rinehart and Winston, p.
236-333.
are in goodagreementwith the data of Barton and Barton, P. B., Jr., and Toulmin, P., III, 1964, Experimental
Toulmin(1966, Fig. 14) for ],res at low FeS con- determination of the reaction chalcopyrite + sulfur =
tents. pyrite + bornite from 350ø to 500øC: Ecoa. GzoI.., v.
59, p. 747-752.
Scott (1973) found that calculationof the pres- 1966, Phase relations involving sphalerite in the Fe-
sure effecton sphaleritecompositionalong the py- Zn-S system: Ecoa. GzoI.., v. 61, p. 815-849.
1334 GERALD K. CZztMANSKE
Boorman, R. S., 1967, Sobsolidusstudiesin the ZnS-FeS- data to problemsof ore d•positsII. Mineral assemblages
FeSu system: Ecoa. GEoI.., v. 62, p. 614-631. and the compositionof ore-forming fluids: Eco•q. GF.oI..,
. Sutherland, J. K., and Chernyshev, L. V., 1971, New v. 60, p. 1101-1166.
data on the sphalerite-pyrrhotite-pyrite solvus: Eco•q. Kissin, S. A., and Scott, S. D., 1972, Phase relations of
GEo•..,v. 66, p. 670-673. intermediate pyrrhotite [abs.]: Geol. Soc. America Ab-
Brett, R., 1964,Experimentaldata from the systemCu-Fe-S stracts with Programs for 1972, p. 562, also in Eco•q.
and their bearing on exsolution textures in ores: Ecoa. GEoI..,v. 67, p. 1007.
Gr.oI.., v. 59, p. 1241-1269. Kullerud, G., 1953, The FeS-ZnS system, a geological
Cabri, L. J., 1973, New data on phaserelations in the Cu- thermometer: Norsk. geol. tidsskr., v. 32, p. 61-147.
Fe-S system: Ecoa. G•oI.., v. 68, p. 443-454. Robie, R. A., Bethke, P.M., Toulmin, M. S., and Edwards,
Cameron, E. N., and Van Rensburg,W. C. J., 1965, Chemi- J. L., 1966, X-ray crystallographic data, densities, and
cal-mechanical polishing of ores: Eco•q. G,v. oI.., v. 60, p. molar volumes of minerals, in Clark, S. P., Jr., ed., Hand-
630-632. book of physical constants,rev. ed.: Geol. Soc. America
Chernyshev, L. V., and Anfilogov, V. N., 1967, Subsolidus Mem. 97, p. 27-73.
phase relations in the ZnS-FeS-FeSu system: Ecoa. G•oI.., Schneeberg,E. P., 1973, Sulfur fugacity measurementswith
v. 63, p. 841-844. the electrochemical cell AgIAglIAg•+xS, fs•: Ecoa. G•oI..,
Czamanske, G. K., and Rye, R. O., 1974, Experimentally v. 68, p. 507-517.
determined sulfur isotope fractionations between sphal- Scott, S. D., 1973, Experimental calibration of the sphal-
erite and galena in the temperature range 600ø to 275øC: erite geobarometer: Ecoa. Gini.., v. 68, p. 466-474.
Ecoa. G•oI.., v. 69, p. 17-25. and Barnes, H. L., 1971, Sphalerite geothermometry
Desborough,G. A., Heidel, R. H., and Czamanske,G. K., and geobarometry: Ecoa. G•oI.., v. 66, p. 653-669.
1971, Improved electron microprobe analysis at low Toulmin, P., III and Barton, P. B., Jr., 1964, A thermo-
operating voltage: II. Sulfur: Am. Mineralogist, v. 56, dynamic study of pyrite and pyrrhotite: Geochim. et
p. 2136-2141. Cosmochim.Acta., v. 28, p. 641-671.
Giletti, B. J., Yund, R. A., and Lin, T. J., 1968, Sulfur Yui, S., and Czamanske, G. K., 1971, Iron content of
vapor pressure of pyrite-pyrrhotite [abs.]: Ecoa. G•o•-., v. sphalerite from the Iimori mine, Japan: Soc. Mining
63, p. 702. GeologistsJapan, Spec. Issue 3, p. 277-279.
Hall, S. R., and Stewart, J. M., 1973, The crystal struc- Yund, R. A., and Kullerud, G., 1960, The Cu-Fe-S system:
ture refinement of chalcopyrite, CuFeSu: Acta Crystallo- Carnegie Inst. Washington Yearbook 59, p. 111-114.
graphica,B29, p. 579-585. 1966, Thermal stability of assemblagesin the Cu-Fe-S
Holland, H. D., 1965, Some applicationsof thermochemical system: Jour. Petrology, v. 7, p. 454-488.