Comparison and statistical evaluation of neutron activation methodologies for the determination of gold in copper concentrate

Gold was determined in copper concentrate by an instrumental neutron activation analysis technique using the relative method. The determination was repeated by radiochemical and chemical neutron activation analysis with the aid of an anion exchange strategy for the chemical separation of gold from the matrix (Fe, Cu, and Zn). Irradiation of copper concentrate samples was carried out in APSARA reactor (BARC, India) with a thermal neutron flux of about 10 cm 2 s . Analysis of variance established that the results obtained by instrumental, radiochemical and chemical methodologies were statistically indistinguishable. Among the three approaches, radiochemical activation showed superior detection limits, which was attributed to both the reduction in background and increase in sensitivity. The best measurement repeatability was observed in chemical activation methodology, compared to the other two. The rationale behind the improvement in measurement repeatability was the capability of chemical activation to process a larger mass of sample for each replicate, which resulted in improved counting statistics and reduction in sampling error with respect to gold. Combined uncertainty for all the three methodologies was evaluated through the bottom-up approach. Systematic evaluation of various uncertainty parameters showed that the major contributor to the combined standard uncertainty was the counting statistics during the instrumental approach and the chemical separation yield during radiochemical and chemical activation. The combined uncertainty for all the three approaches was less than the measurement repeatability indicating the inhomogeneous distribution of gold in copper concentrate. Calculations showed that the sampling constant for gold in copper concentrate was about 22 g.


Introduction
Copper occurs as oxides, carbonates or suldes in nature and sulde ores are the most common. 1 The copper ore is crushed and ground before it is concentrated to have a copper concentration between 20 and 40%. High grade copper sulde ores are concentrated using the froth otation process and the dried froth is known as copper concentrate. Copper concentrates are the raw materials for copper smelters and are used to produce high purity metals. They are composed of approximately equal parts of copper, iron and sulfur by mass. Besides these elements, they generally contain precious metals like gold (Au) in the range of 0.01 to 200 mg kg À1 . Determination of Au in copper concentrate is essential from the recovery point of view.
Methods for the determination of noble metals have been reviewed by Beamish and Van Loon. 2 Natural levels of noble metals are in the region of limits of detection for most of the current analytical techniques. 3 In order to achieve an interference-free determination at trace levels, it is necessary to separate and/or pre-concentrate noble metals. The most common method for determination of noble metals in geological samples is atomic absorption spectrometry with a re assay pre-concentration step. The re assay technique requires special facilities, equipment and a highly skilled analyst. Also, this technique is unsuitable for sulphide minerals. 4 Cyanidation is not a feasible separation route, since it is non-specic and most of the metal cyanides, including that of copper, are highly toxic. 5 Ion exchange methods have been extensively used for the separation of Au from different matrices. [6][7][8][9] Neutronactivation analysis is one of the most sensitive techniques for determination of noble elements. [10][11][12][13][14] Instrumental neutron activation analysis (INAA) has been successfully applied for the determination of noble elements in ores. However, it is necessary to separate noble metals from the matrix, for determination at a low level. Chemical separation of analytes from the matrix can be carried out either before or aer irradiation. Radiochemical neutron activation analysis (RNAA) involves chemical separation of the analyte (radionuclide) from the matrix, aer neutron irradiation. When pre-irradiation chemical separations are employed, the procedure is known as chemical neutron activation analysis (CNAA). 15 As the matrix is removed prior to irradiation, larger quantities of samples can be processed in CNAA, leading to preconcentration of the analyte. CNAA methodology has the additional incentive of minimal radiation exposure compared to RNAA, since the matrix separation is carried out prior to irradiation. The gamma ray spectra in RNAA and CNAA are simpler than those of INAA, owing to the removal of matrix activity. The reagent blank correction is not required in RNAA, as the separation is carried out aer irradiation.
The reliability of the results can be ensured by adopting quality control measures in each step of the quantication procedure. Any analytical result is complete only if it is accompanied by a statement of uncertainty in the measurement. 16 Measurement uncertainties can be estimated using statistical analysis (Type A evaluation) of a set of measurements or by using the existing/available information (Type B evaluation) about the measurement process. 17,18 For evaluating the uncertainty associated with the complete measurement process, either the bottom-up or the top-bottom approach can be adopted. 19 Uncertainty evaluation helps identifying the chief sources of uncertainty, which can be easily visualized from the bar diagram representation, showing the fractional contribution from each source. Once the major contributors are iden-tied, the analyst can pay attention to the key step and attempt to improve the quality of the results.
In the present work, Au was determined in copper concentrate by the INAA technique. The determination was repeated using RNAA and CNAA techniques and the results were compared. Detection limits for all the three methodologies were estimated and compared. The bottom-up approach of uncertainty evaluation was adopted for deriving the combined uncertainty of all the three methodologies. Various uncertainty sources of the measurement process were identied and estimated systematically. Measurement repeatability was compared with the combined uncertainty, in order to test for the heterogeneity of Au in copper concentrate. Statistical evaluation of the results is discussed.

Materials and methods
All reagents used were of analytical grade. Polyethylene sheets and glassware were cleaned with dilute nitric acid and deionized water. Analytical grade anion exchange resin Dowex 1X8 (50-100 mesh), Biorad Laboratories, USA was used for ionexchange separation.
The energy dispersive X-ray uorescence (EDXRF) spectrum of copper concentrate was acquired using a Jordan Valley EX3600M spectrometer, having a stabilized low power (50 watt) X-ray tube and a Rh-target for excitation.

Instrumental neutron activation analysis
Gold standard solution (100 mg L À1 ) was prepared by dissolving high purity Au wire in aqua-regia. Accurately weighed (150-200 mg) copper concentrate samples were heat-sealed in polyethylene. 50 mL of the Au standard solution was evaporated on lter paper placed on clean polyethylene, heat-sealed, and was used as the comparator. The dimensions of the sample and comparator were identical. The samples, comparator and blank (polyethylene) were triply-sealed in another polyethylene packet to make it waterproof. Irradiation was carried out for 7 hours at the F7 position (thermal neutron ux $10 12 cm À2 s À1 ) in APSARA, a swimming pool type nuclear reactor, BARC, Trombay, Mumbai, India.
Gamma ray spectrometric measurements were carried out using a High Purity Germanium (HPGe) detector, Eurisys Mesures, France (Relative efficiency: 45%, resolution: 1.95 keV at 1332 keV gamma ray of 60 Co), coupled to an 8k channel analyzer. Counting geometries for the blank, comparator and samples were identical during gamma ray measurements.

Radiochemical neutron activation analysis
Sample preparation and irradiation conditions were identical to those described in the preceding section (instrumental neutron activation analysis). Irradiated copper concentrate samples ($150-200 mg each) were transferred to polytetrauoroethylene (PTFE) beakers; small volumes of HF and HNO 3 were added and heated on a sand bath. During digestion of the samples, small volumes of HF and HNO 3 were added repeatedly for complete decomposition. Aer dissolution, solutions were evaporated to dryness; residues were dissolved in a minimum volume of concentrated HCl and evaporated to dryness on a water bath. Treatment with HCl was repeated. Finally the residues were added to 25 mL of 1 M HCl (clear green solution). Approximately 200 mg of Dowex 1X8 resin (pre-equilibrated with 1 M HCl and air-dried) was added to the sample solution and agitated on a mechanical shaker for 2 hours. The resin beads were ltered, washed with 0.1 M HCl, air-dried and heat sealed using polyethylene. Gamma ray spectrometric measurements of resin beads (loaded with gold) were carried out using a HPGe detector as described previously (instrumental neutron activation analysis).

Determination of chemical yield for the anion exchange separation of gold
Standard solution containing 25 mg of Au was evaporated on clean polyethylene, heat-sealed and was irradiated in the F7 position of an APSARA reactor for 7 hours. Aer irradiation, the polyethylene packet was opened, treated with concentrated HCl-HNO 3 mixture, made up to 10 mL in 1 M HCl and used as radiotracer solution.
About 500 mg of accurately weighed copper concentrate sample was taken in a PTFE beaker. The radiotracer solution of Au (250 mL) was added, the sample was decomposed with HF and HNO 3 and ion exchange separation was carried out as described previously (radiochemical neutron activation analysis). The resin beads were ltered and washed with 0.1 M HCl. Filtrate and washings (25 mL) were collected in a clean ask. Gold reference solution was prepared by adding 250 mL of Auradiotracer to 25 mL of deionized water. Gamma ray measurements of the ltrate and reference solution were carried out in identical geometry using a HPGe detector and count rates for 411.8 keV gamma rays of 198 Au were compared. The chemical yield for the anion exchange separation of Au was calculated using eqn (1).
where, cps ref and cps ltrate are count rates for the reference Au solution and ltrate, respectively.

Chemical neutron activation analysis
About 500 mg of accurately weighed copper concentrate sample was taken in a PTFE beaker. The sample was decomposed with a HF-HNO 3 mixture and ion exchange separation was carried out as described previously (radiochemical neutron activation analysis). Resin beads were ltered, washed with 0.1 M HCl, air dried and heat sealed using polyethylene. The gold comparator was prepared by drying the standard solution on lter paper having dimensions identical with polyethylene packets containing Au-loaded resin beads. The sample (Au-loaded resin) along with the comparator and blank (resin beads in polyethylene) was irradiated and radioactive assay was carried out as described previously (instrumental neutron activation analysis).

Calculation of concentration of Au in copper concentrate
The relative (comparator) method of NAA was used for calculating the concentration of Au. The mass of Au (m sam ) in the copper concentrate sample was obtained using eqn (2).
where, N is the net counts under the photopeak of 411.8 keV gamma rays of 198 Au corrected for decay time, m is the mass (mg), sam and std represent the sample and standard, respectively. The m sam is converted to the concentration (mg kg À1 ) by dividing with the sample mass (g). In the case of RNAA and CNAA, eqn (2) was modied, considering the chemical yield factor as shown in eqn (3).

Results and discussion
The EDXRF spectrum of the copper concentrate sample (Fig. 1) showed that Cu, Fe and S were the major constituents and Zn, Pb, Ti, Mo, Si, Sr and Ca were the minor constituents. The peak corresponding to the X-ray energy of Au (L III M V : 9.7 keV) was not visible in the EDXRF spectrum, due to the low concentration of Au as well as the absorption of characteristic X-rays of Au by the matrix (Cu K-edge ¼ 8.9 keV). 20 The nuclear reaction 197 Au(n,g) 198 Au was used for the determination of Au during NAA, and the relevant nuclear parameters are given in Table 1. 21 During INAA, since the activity due to the matrix ( 64 Cu) was high, both the sample and comparator were counted at an appropriate distance from the detector to keep the dead time below 5%. As a result, the observed count rate of 198 Au reduced signicantly. An alternative option was to cool/decay the samples for about 7 days, in order to bring down the activity of the matrix ( 64 Cu, t 1/2 ¼ 12.70 h) to a sufficiently low level. However, cooling would have resulted in the reduction of 198 Au activity (t 1/2 ¼ 64.67 h) also. These practical difficulties resulted in an inferior count rate ( 198 Au) and hence the poor detection limit.
NAA, in principle, is a non-destructive analysis technique. However, when elements having a high neutron absorption cross-section are the major/minor constituents, whose activation products have appreciable half lives, gamma rays from the matrix complicate the measurements by contributing to the background and dead time of the counting system. These may lead to inferior detection limits. In such cases, it is essential to remove the matrix by resorting to chemical methods (RNAA and CNAA).
The most important factor to be considered during RNAA and CNAA is the yield of the separation process, as no process has 100% efficiency. The chemical yield of a separation process is determined in order to correct for the possible losses of the analyte. 22 The radiotracer method is an efficient and convenient tool for the determination of the chemical separation yield for various metal ions. 23 Ion exchange procedures have been reported for the separation of noble elements prior to NAA. 3,9 Noble elements form anionic complexes with moderate molarities of HCl, which are readily retained on an anion exchanger. The distribution coef-cient (k d ) for Au, Cu and Fe for the anion exchanger from HCl solution has been reported by Kraus and Nelson. 24 The distribution coefficient for Au(III) is of the order of 10 4 at lower acidities (<3 M HCl) and decreases with increasing acid concentration. Cu(II) is not sorbed below 2 M HCl whereas the distribution coefficient for Fe(III) is negligible below 1 M HCl. Therefore anion exchange separation of Au from Cu and Fe was carried out using 1 M HCl during RNAA and CNAA. The ion exchange separation yield for Au from copper concentrate was found to be 98.6 AE 1.1% (n ¼ 5) using the radiotracer 198 Au.
The cardinal effects of radiochemical separations are improved counting statistics and reduced interference. 22 Greenberg, 25 Kucera and Zeisler 26 have reported the importance of radiochemical separations in improving the accuracy, uncertainty and detection limits of NAA. The detection limit depends on the irradiation, decay and counting conditions. It also depends on the interference situation including the ambient background, the Compton continuum from high energy gamma rays, as well as any gamma ray interference from such factors as the blank from pre-irradiation treatment and packing materials. 23 In RNAA and CNAA, due to the reduction in count rates (matrix activity), samples could be counted closer to the detector. In addition, the post-irradiation cooling step was not required in RNAA and CNAA which was necessary during INAA, for reducing the dead time of the counting system. The higher count rates for the analyte obtained in RNAA and CNAA were due to improved counting geometries as well as the obviation of the decay time.
The results obtained by all the three NAA methodologies are in good agreement ( Table 2). Values obtained for Au during RNAA and CNAA were corrected for the chemical yield.

Analysis of variance of the results
Analysis of variance (ANOVA) was carried out to test the statistical parity of the Au values obtained by INAA, RNAA and CNAA methodologies. Results of ANOVA are tabulated in Table 3. The experimentally derived F-ratio (F (2,12),experimental ¼ 0.94) was less than the critical F-ratio (F (2,12),critical ¼ 3.89), 27 which established the statistical equality of the results obtained by all the three methodologies. Table 2 that the measurement repeatability for INAA and RNAA is inferior compared to that of CNAA. Generally, the measurement repeatability in activation analysis is governed by both the counting statistics and sample homogeneity characteristics. Since the count rates ( 198 Au) were sufficiently high for samples and standards, during all the three methodologies, the inhomogeneous distribution of Au in the copper concentrate matrix could be the key rationale for the inferior repeatability in INAA and RNAA. The better measurement repeatability in the case of CNAA could be due to the higher sample mass.

Limit of detection
The limit of detection (LOD) was calculated as recommended by IUPAC and ACS and approximating normal distribution for counting statistics. 28,29 Table 2 shows the calculated detection limits for INAA, RNAA and CNAA. Among all the three methodologies, RNAA showed a superior detection limit, which was the combined outcome of the reduction in the background (matrix activity) and improved sensitivity (higher net counts due to better counting geometry). The detection limit of CNAA was slightly inferior to that of RNAA because of higher background (due to less cooling time compared to RNAA and 82 Br from the resin matrix). The detection limit was poor in the case of INAA because of high background (matrix activity) and low sensitivity (source to detector distance was large to keep the dead time of the counting system low).
For qualitative comparison, gamma ray spectra recorded during INAA, RNAA and CNAA methodologies are shown in Fig. 2. The counting geometry and counting time in all the three cases were the same whereas the cooling time was different. It is vivid from Fig. 2 that the background is minimal in the region of interest (around 411.8 keV) in RNAA, compared to the other two.

Evaluation of the combined and expanded uncertainties
The most important characteristic of a result obtained from any measurement is its uncertainty. However, most of the analytical  results are accompanied by an indication of repeatability alone. ISO/IEC recommends that the laboratories should express measurement uncertainty. 16 The rst step in the uncertainty evaluation of an analytical system 30 is the identication of all probable sources. Further, one should quantify the contribution from each source to the combined uncertainty, followed by the proper combination of standard uncertainties. INAA is a primary method with well known sources of uncertainties. 31,32 Although there are large number of parameters involved in activation analysis, the use of the relative method results in cancellation of nearly all parameters, except the spectral peak area and mass of the standard (eqn (2)). The major source of uncertainty in all the methodologies of NAA is the counting statistics along with the chemical yield determination in the cases of RNAA and CNAA. 33 The relevant sources of uncertainty in the present NAA determination of Au could be identied as (i) the determination of mass of the sample and standard (ii) gamma ray spectrometric measurements, and are shown in the cause and effect diagram (Fig. 3). The uncertainty associated with weighing and volumetric operations was evaluated as described in the Eurachem/CITAC Guide. 34 The contributors to the combined uncertainty during gamma spectrometry could be identied as the counting statistics, counting geometry and dead time. Effects of counting geometry on the measurement uncertainty were brought down to negligible levels, by matching the dimensions of the sample and comparator and maintaining an identical distance from the detector in each set of measurement. The dead time was minimized (<5%) by optimizing the counting geometries in each analysis.
Combined standard uncertainty for RNAA and CNAA was derived by propagating the standard uncertainty from each source, using eqn (4), based on the standard method of error propagation. 34 where, u represents the standard uncertainty associated with the measurement of the respective component, C is the concentration, m sam is the mass of the sample, m Au is the mass of the Au comparator, N is the net area under the photopeak of 411.8 keV gamma ray of 198 Au corrected for decay time and % yield is the ion exchange separation yield. In INAA, the last term in eqn (3) was not included in the combined uncertainty calculation, as the chemical separation step was not involved in this approach. The second term under the square root in eqn (3) is expanded as eqn (5).   View Article Online u m Au m Au where, u represents the standard uncertainty associated with the measurement of the respective component, m Au is the mass of the Au comparator, w Au represents the mass of Au used for preparing the stock solution of the comparator, V 1 is the volume of the stock solution of Au, P 1 is the volume of stock solution pipetted for preparing the working standard solution of Au, V 2 is the total volume of the working standard solution of Au and P 2 is the volume of the working standard of Au, pipetted for irradiation. The contributions from individual components towards the combined uncertainty of the complete measurement process were evaluated via the bottom-up approach and the results are depicted in Fig. 4. It is vivid from Fig. 4 that the major share in combined uncertainty in both RNAA and CNAA determinations arises from the chemical yield determination of the ion exchange separation process of Au from the matrix, whereas in INAA it comes up from the counting statistics. The expanded uncertainties were calculated at 95% condence level.
In order to have a quick selection among the three methodologies for determination of Au in the copper concentrate sample, the sample mass has been plotted against the Au concentration as shown in Fig. 5. The plot has been prepared considering (i) uncertainty in the measurement; (ii) dead time of the HPGe detector (#10%); (iii) maximum expected concentration of Au in copper concentrate (6%); (iv) maximum sample   mass for irradiation; (v) resin size (0.2 g), limit of quantitation and sample mass. It is important to note that this diagram is valid for given experimental conditions which include neutron ux in the irradiation position, efficiency of the HPGe detector and homogeneity of the copper concentrate sample. The measurement repeatability (precision) and the combined uncertainties of INAA, RNAA and CNAA methodologies are given in Table 2. It is instructive to compare measurement repeatability with uncertainty obtained from Type A evaluation. 34 If the measurement repeatability is higher than the uncertainty from Type A evaluation, it could indicate that either the measurement process is not fully understood or there is a signicant error in sub-sampling (indicative of the degree of heterogeneity). The total variance, S 2 as determined from a set of replicate samples of approximately equal mass, is composed of the variance of the analytical system, S a 2 and the sampling variance from the heterogeneity, S s 2 , as given by eqn (6).
Taking measurement repeatability as S 2 and combined uncertainty as S a 2 , the sampling variance can be calculated using eqn (7). Ingamells 35 introduced a term "sampling constant", K s (in unit of mass) where, w is the sample mass and S s 2 is the sampling variance. K s is numerically equal to the subsample mass necessary to ensure a relative sampling error of 1% at 68% condence level. Quantitative estimation of the sampling error is carried out by repetitive determination of analytes in the solid (sample) using an analytical technique having high precision. 36 The analytical uncertainty should be well understood and sufficiently small. The NAA method is precise, accurate and meets the above mentioned criteria. 31 An accurate measurement in activation analysis is the one whose result is indistinguishable from the true value within the interval set by counting statistics. Gold was determined in the Au-lm standard by INAA using a conventional counting system and standard procedure to obtain accuracy comparable to counting statistics i.e., well below 1%. 37 The INAA technique was used to assess whether the observed variability was due to the analytical system alone, or sampling error also. 38 Calculations showed that the sampling constant, K s , was $22 g for Au in copper concentrate in the present NAA determinations. It is interesting to note that the sampling constant was the highest for Au among all the elements in IAEA-396. 36 Kontas also has reported a sample mass of 30 g as the representative analytical portion for determination of Au in ore reference materials. 39 The representativeness of the analytical portion (laboratory subsample) is debatable, especially in the case of rocks and minerals 36 and more so in the case of gold ores or gold bearing concentrate. Typically, gold occurs as native grains in these materials and it is very difficult to homogenize such samples. Improvement in the precision of Au values on increasing the sample mass has been reported in the literature. 39 The combined uncertainties for INAA, RNAA and CNAA methodologies were practically identical, as shown in Table 2. The smaller combined uncertainties, compared to measurement repeatability, were indicative of the in-homogeneous distribution of Au in the copper concentrate sample.

Conclusion
Gold was determined in copper concentrate by INAA, RNAA and CNAA methodologies using the relative method. The results obtained by all the three methodologies were in very good agreement. The statistical indistinguishability of the results was evidenced by ANOVA. INAA methodology can be used for determination of gold in copper concentrate as it is simple to perform, even though the turn-around time is high. The detection limits are signicantly lower in the case of RNAA and CNAA compared to INAA and are more suitable for determination of gold at a very low level provided that the uncertainty permits. The observed better measurement repeatability of CNAA was due to the large sample mass and improved counting geometry (higher count rate). RNAA stands rst in terms of the detection limit, which was the combined outcome of the reduction in the background, increase in the sensitivity and the inherent freedom from the reagent blank.
Step-wise, systematic uncertainty budgeting could throw light upon the chief component towards the combined uncertainty. It was the counting statistics in the case of the instrumental approach and the chemical separation yield for the radiochemical and chemical activation analysis methodologies. NAA could be used for establishing the heterogeneity of Au in copper concentrate samples, as the combined uncertainty was smaller than the measurement repeatability.