2.2       Ion sources

 

As noted above, the traditional start of mass spectrometric analysis is to heat the sample under vacuum, leading to thermal ionisation. With the exception of the rare gases, thermal ionisation is normally achieved by loading a solid deposit of the sample onto a metal filament, which can then be subjected to resistive heating. However, the manner by which the sample is deposited and then heated has a major effect on the efficiency of the analysis. Therefore these procedures will be discussed in some detail.

 

 

2.2.1    Thermal ionisation

 

For some elements such as Sr, stable emission of metal ions is achieved from a salt deposited directly onto a single metal filament (Fig 2.5a), usually tantalum (Ta). The loading procedure involves evaporating the salt solution onto the filament before insertion into the vacuum system. The sample is often loaded in phosphoric acid, which seems to a) displace all other anion species to yield a uniform salt composition, b) destroy organic residues (such as ion exchange resin) mixed with the sample, and c) glue the sample to the filament. During mass spectrometric analysis, the filament current is raised by means of a stabilised power supply to yield a temperature where simultaneous volatilisation and ionisation of the sample occurs.

Fig. 2.5. The arrangement of filament ribbons on commonly used single- and triple-filament bead assemblies. Note that only one side filament is shown attached to the ‘triple’ bead.

 

            However, for many elements, stable volatilisation and ionisation of metal species does not occur at the same temperature. This problem was noted by Ingram and Chupka (1953), who first proposed the use of multiple source filaments (Fig. 2.5b). In this configuration, one or more filaments bearing the sample load can be heated to the optimum temperature for stable volatilisation, while another hotter filament can be used to ionise the atomic cloud by bombarding it with electrons.

 

            This method is particularly effective for REE analysis, where the sample is usually loaded onto one or both of the Ta side filaments of a triple-filament bead. These are held at a moderate temperature (ca. 1400 oC) where REE volatilisation is most stable. The centre filament (usually Re) is held at a much higher temperature (ca. 2000 oC), which promotes ionisation of the metal vapour. To some extent the ratio of metal to oxide species can be controlled by the centre filament temperature, which may help to suppress REE isobaric interferences. The properties of the REE under such conditions vary from light to heavy rare earths. La and Ce tend to form oxides unless extremely high centre filament temperatures are used, while heavier REE tend to form the metal species (Hooker et al., 1975; Thirlwall, 1982).

 

            Uranium and thorium may also be analysed by the triple-filament technique. Again, the temperature of the centre filament controls the metal/oxide ratio of the emitted ions (Li et al., 1989). The triple-filament method was also used in the first successful analysis of Hf (Patchett and Tatsumoto, 1980). However, Hf analysis is now exclusively performed by ICP-MS (section 2.5.2).

 

            An alternative to TIMS analysis with multiple filaments is to use special conditions to control the evaporation)ionisation behaviour from a single filament. For example, in the case of Pb, the sample is usually loaded on a rhenium filament in a suspension of silica gel (Cameron et al., 1969). This is thought to form a blanket over the sample which effectively retards Pb volatilisation so that the filament can be raised to a higher temperature (where Pb fractionation is more reproducible) without burning off the sample uncontrollably.

 

            Early Nd isotope determinations (Lugmair et al., 1975; DePaolo and Wasserburg, 1976) were made using NdO+ ions, whose emission from a single filament source was promoted by relatively high oxygen pressures in the mass spectrometer source housing. Some workers have continued to use this method, rather than the more popular multiple-filament method, since it can yield higher efficiency. Oxygen may be bled into the source in minute amounts to increase oxide emission. Alternatively, loading with silica gel may achieve the same objective without degrading source vacuum (Thirlwall, 1991a). Uranium may also be analysed as the oxide by mixing with TaO2 powder on a tungsten filament.

 

            A different approach to single filament U and Th analysis is to use graphite to promote the formation of metal ions (Edwards et al., 1987). In this procedure it is critical to prevent oxidation during sample loading by avoiding oxidising acids and by maintaining temperatures below the point of visible glowing. Noble et al. (1989) applied a similar technique to Nd isotope analysis. The use of platinised graphite was argued to give greater thermal stability to the reducing agent.

 

2.2.2    Plasma source mass spectrometry

 

Technically, plasma source mass spectrometry is a kind of thermal ionisation mass spectrometry, because it relies on heating in a gas plasma to achieve ionisation of the sample. However, for practical purposes, the term ‘TIMS’ is usually restricted to the case where thermal ionisation achieved by a heated filament under vacuum. Therefore, plasma source mass spectrometry is regarded as a distinct technique.

 

            Plasma source mass spectrometry was invented in the late 1970s when the inductively coupled plasma (ICP) was first attached to a mass spectrometer (MS) to produce the ICP-MS (Houk et al., 1980). The ICP source consists of a plasma torch made by using a radio frequency (RF) generator to induce intense eddy currents in a stream of ionised argon gas. The RF generator transmits about a kilo-watt of power into the plasma, raising its temperature to about 5000 oC and causing very efficient ionisation of most elements (Houk, 1986). Furthermore, the extreme temperature of the plasma ensures that nearly all of these ions are monatomic ions, ideal for mass spectrometric analysis.

 

            When the ICP-MS was first conceived, it was found most convenient to use a quadrupole mass spectrometer (Dawson, 1976) as the analyser (Fig. 2.6). This analyser has four rods which are arranged parallel to the path of the ion beam. When alternating voltages with variable frequencies are applied to these rods, the analyser can be tuned to allow a certain mass of ions to pass through to the collector, while all other masses are de-focussed. By rapidly changing the electrical frequencies on the rods, the quadrupole analyser can be made to rapidly ‘scan’ the mass spectrum from mass zero to 250, allowing the detector (usually an ion multiplier) to measure the relative abundance of every mass in the ‘spectrum’. A detailed review is given by Potts (1987).

Fig. 2.6. Schematic illustration of an ICP-MS instrument with quadrupole analyser.

 

            The main technical breakthrough in the development of the ICP-MS was the physical feat of actually feeding a plasma at 5000 oC and atmospheric pressure into a mass spectrometer whose analyser pressure is ca. 108 times lower (10-5 mbar). This was achieved by firing the plasma at a two-stage water-cooled orifice, with continuous pumping of the intermediate space by a mechanical pump (Fig. 2.6). Subsequent technical developments in ICP-MS over the next 20 years have mainly involved greatly increased efficiency in the sampling of the plasma by the mass spectrometer. This has allowed ICP-MS to reach remarkable sensitivity, with detection limits as low as parts per trillion (pico grams per gram). However, the precision of isotope ratio measurements by conventional ICP-MS is limited by instabilities in the plasma. By rapidly scanning the ‘mass spectrum’, the quadrupole analyser can achieve precision of around 1%, but this is a practical limit for single collector analysis. This level of precision was found to be useful in early development of the Re-Os method (section 8.1), but is not useful for most radiogenic isotopic systems. Hence, to apply ICP-MS to these systems it was necessary to introduce multiple collection techniques to cancel out the instability of the source.

 

            To perform multiple collection mass spectrometry (MC-MS) with a plasma ion source it was necessary to link the ICP source with a magnetic sector mass spectrometer of the type normally used in TIMS instruments (section 2.3). This involved two main technical challenges. The first problem is that magnetic sector mass spectrometry requires a large accelerating voltage to raise ions to the high velocities where magnetic separation is efficient. However, in order to keep the whole analyser assembly at electrical ground, the ion source must be at up to 8000 v positive. For MC-ICP-MS, this means isolating the plasma at up to 8000 v positive, with all of the attendant engineering problems. The second problem is that magnetic sector mass spectrometry requires an ion beam with a very small range of ion energy, whereas the ICP source generates ions with a relatively large energy range. This was overcome in the first MC-ICP-MS instruments by using both quadrupole and electrostatic pre-filters to smooth the energy distribution in the ion beam before it entered the magnetic sector (Fig. 2.7). Subsequent instruments have used a variety of other filtering devices to achieve similar objectives.

 

Fig. 2.7. Simplified plan view of the VG Elemental Plasma 54 instrument. After Halliday et al. (1998).

 

 

2.2.3    Mass fractionation

 

The process of volatilisation and ionisation during mass spectrometry requires the breaking of chemical bonds, but the strength of these bonds is mass dependent. Therefore, excitation of the sample leads to mass-dependent fractionation, which can be understood by approximating the chemical bond between two atoms as a harmonic oscillator.

 

            The energy of a molecule (or part of an ionic lattice) decreases with decreasing temperature, but at absolute zero it has a certain finite value called the zero point energy, equal to 0.5 h< (where h is Plank’s constant and < is the vibrational frequency). A bond involving the light isotope of an element has a higher vibration frequency and hence a higher zero point energy than one involving a heavier isotope, as illustrated in Fig. 2.8. The difference in bond energies diminishes as temperature rises, but still persists. Because the potential energy well of the bond involving the lighter isotope is always shallower than for the heavier, the bond with the lighter isotope is more readily broken. Hence it is preferentially released from the hot filament, causing isotopic fractionation.

Fig. 2.8. Schematic diagram of potential energy against bond length for a hypothetical molecule made of two isotopes, based on the ‘harmonic oscillator’ model.

 

            In a plasma source mass spectrometer (section 2.2.2), fresh sample is continually fed into the plasma torch. Hence, mass fractionation produces a fairly constant (but large) discrepancy between the isotopic composition of the solid sample and the ion cloud. In contrast, solid source TIMS analysis produces smaller fractionation effects, but the continual process of fractionation starts to ‘use up’ the lighter isotope on the filament so that the isotopic composition of the sample gets progressively heavier (the ‘reservoir effect’). Eberhardt et al. (1964) showed that this process follows a Rayleigh fractionation law (Fig. 2.9). The magnitude of this effect could yield totally unacceptable errors of up to 1% in measured isotope ratio. However, for elements with two or more non-radiogenic isotopes, an internal normalisation for such mass-dependent fractionation can be performed.

 

Fig. 2.9. Effect of within-run fractionation, over time, on a sample of natural rubidium undergoing isotopic analysis. Points are observed ratios; dashed line schematically indicates actual composition of Rb on the filament. Data from Eberhardt et al. (1964).

 

            In the case of strontium, the fractionation of 87Sr/86Sr can be monitored using the 88Sr/86Sr ratio, since 88Sr and 86Sr are both non-radiogenic (i.e. produced only by nucleosynthetic processes in stars). The ratio 86Sr/88Sr is constant throughout the Earth and is taken to be 0.1194 by international convention. This value cannot be measured absolutely, but was originally estimated from the average beam composition half-way through very many TIMS runs. The deviation of observed 86Sr/88Sr from 0.1194 at each point through the run is divided by the difference between the two masses ()mass = 2.003) in order to calculate a fractionation factor per mass unit:

 

                        (86Sr/88Sr)obs

                        )))))))))     !  1

                          0.1194

            F  =   ))))))))))))))))                                                          [2.1]

                                    )mass

 

This fractionation factor can then be used to correct the observed (raw) 87Sr/86Sr ratio, for which )mass = 1.003:

 

            (87Sr)               (87Sr)

            ())))  =          ()))) @ ( 1  +  F )mass)                                    [2.2]

            (86Sr)true            (86Sr)obs

 

This has the effect of improving the within-run precision of the 87Sr/86Sr ratio from ca. 1% to better than 0.01%. Neodymium metal analyses are similarly normalised for fractionation (Fig. 2.10) using an internationally agreed value of 146Nd/144Nd = 0.7219 (O’Nions et al., 1979). However, Nd oxide analyses are normalised to different values (Wasserburg et al., 1981) which are incompatible with the Nd metal normalising value.

Fig. 2.10. Plot of raw 146Nd/144Nd ratios and fractionation-corrected 145Nd/144Nd ratios (normalised to 146Nd/144Nd = 0.7219) for a single mass spectrometer run. Each point is a mean of 10 scans of the mass spectrum, while horizontal lines are grand means. After Noble (pers. comm.).

 

            The fractionation correction described above is usually called the linear law, but the power law (Wasserburg et al., 1981; Thirlwall, 1991b) is effectively identical. Both of these laws assume that fractionation is proportional to mass difference only, and is independent of the absolute masses of the fractionating species. In other words, fractionation per mass unit is constant. However, this is an approximation to the real evaporation process, where fractionation per mass unit must vary inversely with the absolute masses of the evaporating species. Russell et al. (1978) first observed a break)down of the linear law in isotopic analysis of the ‘light’ element Ca. To remedy this, they introduced an ‘exponential’ law, where the fractionation factor depends also on the mass of the evaporating species. This gave a better fit to Ca isotope data than the linear law (Fig. 2.11).

Fig. 2.11. Plot of measured/true 44Ca/48Ca versus 40Ca/44Ca ratios showing fit of linear and exponential fractionation laws to typical data from two runs. After Russell et al. (1978).

 

            These problems are much less severe for Sr and Nd isotope analysis because of their heavier masses. However, Thirlwall (1991b) found small deviations from linear law behaviour in a large data set of Sr standard analyses. This is revealed by a correlation between normalised 87Sr/86Sr and average observed 86Sr/88Sr ratios for complete runs (Fig. 2.12). Thirlwall found that he could eliminate the correlation by retrospectively applying an exponential law correction to the data. This is described as follows:

 

  | (87Sr/86Sr)norm            |ln(86/88)              | (86Sr/88Sr)obs               |ln(87/86)

  | )))))))))          |           =          | )))))))))            |                       [2.3]

  | (87Sr/86Sr)corr             |                       |    .1194                      |

 

where ‘norm’ and ‘corr’ refer to products of linear normalisation and exponential correction, and where 86, 87 and 88 are the actual masses of evaporating ions.

 

Fig. 2.12. Plot of fractionation-corrected 87Sr/86Sr ratios against mean un-corrected 86Sr/88Sr ratios for measurements of the SRM987 standard over a period of several months. Error bar indicates average within-run precision. After Thirlwall (1991b).

 

            This model assumes that strontium evaporates from the filament as the species ‘Sr’. On the other hand, Habfast (1983) suggested that a strontium sample on a rhenium filament might evaporate as a species such as SrReO4 rather than atomic Sr. Hence the ‘apparent mass’ of 88Sr which should be used in the exponential correction would be ca. 330 rather than 88. Under such conditions the exponential law might actually produce a worse fit to the true fractionation behaviour of the sample than the linear law. However, Thirlwall’s data suggest that evaporation from a tantalum filament does occur as the metal species, so the exponential model is an improvement over the linear model.

 

            At the limits of analytical precision, fractionation errors have also been observed with the exponential law. However, Thirlwall (1991b) argued that between-bead precision of Nd isotope analyses could be improved to a level comparable with within-run precision by using a secondary normalisation against 142Nd/144Nd. This is based on an observed correlation between fractionation-corrected 143Nd/144Nd and 142Nd/144Nd ratios from standard runs over a period of months (Fig. 2.13). Subsequent work on Nd fractionation during MC (multi-collector) - ICP-MS analysis of Nd showed similar behaviour (Vance and Thirlwall, 2002). However, because mass fractionation in the plasma source is an order of magnitude larger, the effect is much bigger (Fig. 2.13).

Fig. 2.13. Plot of (exponential law) fractionation-normalised 143Nd/144Nd against 142Nd/144Nd for analyses of a laboratory standard over a period of several months using normalisation to 146Nd/144Nd = 0.7219. ( " ) = TIMS data; ( ! ) = MC-ICP-MS data. After Vance and Thirlwall (2002).

 

            Vance and Thirlwall showed that this problem could be solved using two alternative approaches. One is to apply a post-analysis correction to 143Nd/144Nd ratios using the information from the 142Nd/144Nd ratio, as in Fig. 2.13. An alternative approach is to choose masses for the fractionation monitor that bracket the masses to be corrected. For example, the mean mass of 143Nd and 144Nd is 143.5, which is the same as the mean mass of 142Nd and 145Nd. Therefore, 142Nd/145Nd represents the ideal monitor ratio for fractionation correction of 143Nd/144Nd ratios. This approach is only possible if Ce interferences on 142Nd are low, and is therefore not feasible for TIMS analysis using the popular ‘reverse phase’ Nd separation (section 2.1.2). However, for most TIMS analysis the fractionation problem can be avoided by collecting data near a 146Nd/144Nd ratio of 0.7219. For Nd analysis by MC-ICP-MS it is necessary to take account of this problem in order to achieve reproducibilities comparable with  TIMS. However, Ce isobaric interferences can be accurately corrected in MC-ICP-MS analysis, so 142Nd can be used as one of the fractionation monitors.

 

            Internal fractionation correction is only possible where there are two or more isotopes present in a constant ratio. This is not the case for Pb isotope analysis by TIMS, or in the isotope dilution analysis of Rb, so an external correction must be used. This depends on achieving uniform fractionation behaviour between standards and different samples, so that an across-the-board correction can be made to all runs. In the case of Pb, the use of a silica gel blanket on the filament achieves this objective by reducing the magnitude of fractionation processes drastically, from a previous between-lab. variation of ca. 3% to a present variation of ca. 0.3% . This improvement is partly due to the higher filament temperatures possible using silica gel. Because bond energy levels become closer with increasing temperature, the magnitude of isotope fractionation falls with increasing temperature. However, for this technique to work well, the Pb sample has to be well purified during the chemical separation (section 2.1.4). In the analysis of uranium, fractionation effects may similarly be reduced by running at high temperature as the oxide. Alternatively, analysis as the metal ion produces larger but relatively consistent degrees of fractionation, which can be corrected by comparison with standard runs.

 

            Where the relevant nuclides exist, double spiking with two artificial isotopes may be used to apply an internal fractionation correction to elements such as Pb with only one natural non-radiogenic isotope (section 2.4.2). In MC-ICP-MS analysis the two isotopes that form this double spike do not even have to be the same element as the isotope ratios being corrected. For example, the two isotopes of thallium can be used to correct Pb isotope data (section 2.5.4).

 

 

References

 

Arden, J. W. and Gale N. H. (1974). New electrochemical technique for the separation of lead at trace levels from natural silicates. Anal. Chem. 46, 2)9.

 

Aston, F. W. (1919). A positive ray spectrograph. Phil. Mag. (Series 6), 38, 707––14.

Aston, F. W. (1927). The constitution of ordinary lead. Nature 120, 224.

 

Barovich, K. M., Beard, B. L., Cappel, J. B., Johnson, C. M., Kyser, T. K. and Morgan, B. E. (1995). A chemical method for hafnium separation from high-Ti whole-rock and zircon samples. Chem. Geol. (Isot. Geosci. Sect.) 121, 303–8.

 

Blichert-Toft, J., Chauvel, C. and Albarede, F. (1997). Separation of Hf and Lu for high-precision isotope analysis by magnetic sector-multiple collector ICP-MS. Contrib. Mineral. Petrol. 127, 248–60.

 

Brooks, C., Wendt, I. and Harre, W. (1968). A two-error regression treatment and its application to Rb)Sr and initial Sr87/Sr86 ratios of younger Variscan granitic rocks from the Schwarzwald massif, Southwest Germany. J. Geophys. Res. 73, 6071)84.

 

Brooks, C., Hart, S. R. and Wendt, I. (1972). Realistic use of two-error regression treatments as applied to rubidium-strontium data. Rev. Geophys. Space Phys. 10, 551)77.

 

Cameron, A. E., Smith, D. H. and Walker, R. L. (1969). Mass spectrometry of nanogram-size samples of lead. Anal. Chem. 41, 525)6.

 

Cassidy, R. M. and Chauvel, C. (1989). Modern liquid chromatographic techniques for the separation of Nd and Sr for isotopic analyses. Chem. Geol. 74, 189)200.

 

Catanzaro, E. J. and Kulp, J. L. (1964). Discordant zircons from the Little Butte (Montana), Beartooth (Montana) and Santa Catalina (Arizona) Mountains. Geochim. Cosmochim. Acta 28, 87)124.

 

Chen, J. H. and Wasserburg, G. J. (1981). Isotopic determination of uranium in picomole and sub-picomole quantities. Anal. Chem. 53, 2060)7.

 

Christensen, J. N., Halliday, A. N., Godfrey, L. V., Hein, J. R. and Rea, D. K. (1997). Climate and ocean dynamics and the lead isotopic records in Pacific ferromanganese crusts. Science 277, 913–8.

 

Clayton, D. D. (1978). On strontium isotopic anomalies and odd-A p-process abundances. Astrophys. J. 224, L93)5.

 

Collerson, K. D., Kamber, B. S. and Schoenberg, R. (2002). Applications of accurate, high-precision Pb isotope ratio measurement by multi-collector ICP-MS. Chem. Geol. 188, 65–83.

 

Compston, W. and Oversby, V. M. (1969). Lead isotopic analysis using a double spike. J. Geophys. Res. 74, 4338)48.

 

Cotte, M. (1938). Recherches sur l’optique electronique. Ann. Physique 10, 333)405.

 

Crock, J. G., Lichte, F. E. and Wildeman, T. R. (1984). The group separation of the rare-earth elements and yttrium from geologic materials by cation-exchange chromatography. Chem. Geol. 45, 149)63.

 

Croudace, I. W. (1980). A possible error source in silicate wet-chemistry caused by insoluble fluorides. Chem. Geol. 31, 153)5.

 

Crumpler, T. B. and Yoe, J. H. (1940). Chemical Computations and Errors, Wiley, pp. 189)90.

 

Daly, N. R. (1960). Scintillation type mass spectrometer ion detector. Rev. Sci. Instrum. 31, 264)7.

 

David, K., Birch, J. L., Telouk, P. and Allegre, C. J. (1999). Application of isotope dilution for precise measurement of Zr/Hf and 176Hf/177Hf ratios by mass spectrometry (ID–TIMS/ID–MC–ICP–MS). Chem. Geol. 157, 1–12.

 

Davis, D. W. (1982). Optimum linear regression and error estimation applied to U)Pb data. Can. J. Earth Sci. 19, 2141)9.

 

Dawson, P. H. (1976). (Ed.), Quadrupole Mass Spectrometry and its Applications. Elsevier, 349 p.

 

DeBievre, P. J. and Debus, G. H. (1965). Precision mass spectrometric isotope dilution analysis. Nucl. Instrum. Meth. 32, 224)8.

 

Dempster, A. J. (1918). A new method of positive ray analysis. Phys. Rev. 11, 316-24.

 

DePaolo, D. J. and Wasserburg, G. J. (1976). Nd isotopic variations and petrogenetic models. Geophys. Res. Lett. 3, 249)52.

 

Dickin, A. P., Jones, N. W., Thirlwall, M. and Thompson, R. N. (1987). A Ce/Nd isotope study of crustal contamination processes affecting Palaeocene magmas in Skye, NW Scotland. Contrib. Mineral. Petrol. 96, 455)64.

 

Dodson, M. H. (1978). A linear method for second-degree interpolation in cyclical data collection. J. Phys. E (Sci. Instrum.) 11, 296.

 

Dodson, M. H. (1963). A theoretical study of the use of internal standards for precise isotopic analysis by the surface ionisation technique: Part I - General first-order algebraic solutions. J. Sci. Instrum. 40, 289)95.

 

Dosso, L. and Murthy, V. R. (1980). A Nd isotopic study of the Kerguelen islands: inferences on enriched oceanic mantle sources. Earth Planet. Sci. Lett. 48, 268)76.

 

Eberhardt, A., Delwiche, R. and Geiss, J. (1964). Isotopic effects in single filament thermal ion sources. Z. Natur. 19a, 736)40.

 

Edwards, R. L., Chen, J. H. and Wasserburg, G. J. (1987). 238U)234U)230Th)232Th systematics and the precise measurement of time over the past 500,000 years. Earth Planet. Sci. Lett. 81, 175)92.

 

Eugster, O., Tera, F., Burnett, D. S. and Wasserburg, G. J. (1970). The isotopic composition of gadolinium and neutron capture effects in some meteorites. J. Geophys. Res. 75, 2753)68.

 

Faul, H. (1966). Ages of Rocks, Planets, and Stars. McGraw-Hill, 109 p.

 

Gale, N. H. (1970). A solution in closed form for lead isotopic analysis using a double spike. Chem. Geol. 6, 305)10.

 

Galer, S. J. G. (1999). Optimal double and triple spiking for high precision lead isotopic measurement. Chem. Geol. 157, 255)74.

 

Habfast, K. (1983). Fractionation in the thermal ionization source. Int. J. Mass Spectrom. Ion Phys. 51, 165)89.

 

Halliday, A. N., Lee, D.-C., Christensen, J. N., Rehkamper, M., Yi, W., Luo, X., Hall, C. M., Ballentine, C. J., Pettke, T. and Stirling, C. (1998). Applications of multiple collector- ICPMS to cosmochemistry, geochemistry, and paleoclimatology. Geochim. Cosmochim. Acta 62, 919)40.

 

Hamelin, B., Manhes, G., Albarede, F. and Allegre, C. J. (1985). Precise lead isotope measurements by the double spike technique: a reconsideration. Geochim. Cosmochim. Acta 49, 173)82.

 

Hooker, P., O’Nions, R. K. and Pankhurst, R. J. (1975). Determination of rare-earth elements in U.S.G.S. standard rocks by mixed-solvent ion exchange and mass spectrometric isotope dilution. Chem. Geol. 16, 189)96.

Horn, I., Rudnick, R. L. and McDonough, W. F. (2000). Precise elemental and isotope ratio determination by simultaneous solution nebulization and laser ablation-ICP-MS: application to U–Pb geochronology. Chem. Geol. 164, 281–301. (Erratum = vol. 167, 405–25.)

 

Houk, R. S. (1986). Mass spectrometry of inductively coupled plasmas. Anal. Chem. 58, 97A–105A.

 

Houk, R. S., Fassel, V. A., Flesch, G. D., Svec, H. J., Gray, A. L. and Taylor, C. E. (1980). Inductively coupled argon plasma for mass spectrometric determination of trace elements. Anal. Chem. 52, 2283–9.

 

Ingram, M. G. and Chupka, W. A. (1953). Surface ionisation source using multiple filaments. Rev. Sci. Instrum. 24, 518)20.

 

Kalsbeek, F. and Hansen, M. (1989). Statistical analysis of Rb)Sr isotope data by the ‘bootstrap’ method. Chem. Geol. (Isot. Geosci. Section) 73, 289)97.

 

Krogh, T. E. (1973). A low contamination method for hydrothermal decomposition of zircon and extraction of U and Pb for isotopic age determination. Geochim. Cosmochim. Acta 37, 485)94.

 

Kuritani, T. and Nakamura, E. (2002). Precise isotope analysis of nanogram-level Pb for natural rock samples without use of double spikes. Chem. Geol. 186, 31–43.

 

Kurz, E. A. (1979). Channel electron multipliers. Amer. Lab. 11 (3), 67)74.

 

Lee D.-C. and Halliday, A. N. (1995). Precise determinations of the isotopic compositions and atomic weights of molybdenum, tellurium, tin and tungsten using ICP source magnetic sector multiple collector mass spectrometry. Int. J. Mass Spec. Ion Process. 146/147, 35–46.

 

Li, W. X., Lundberg, J., Dickin, A. P., Ford, D. C., Schwarcz, H. P., McNutt, R. H. and Williams, D. (1989). High-precision mass spectrometric uranium-series dating of cave deposits and implications for paleoclimate studies. Nature 339, 534)6.

 

Luais, B., Telouk, P. and Albarede, F. (1997). Precise and accurate neodymium isotopic measurements by plasma-source mass spectrometry. Geochim. Cosmochim. Acta 61, 4847)54.

 

Ludwig, K. R. (1980). Calculation of uncertainties of U)Pb isotope data. Earth Planet. Sci. Lett. 46, 212)20.

 

Ludwig, K. R. (1997). Isoplot. Program and documentation, version 2.95. Revised edition of U.S. Geol. Surv. Open-File report 91-445 (kludwig@bgc.org).

 

Lugmair, G. W., Scheinin, N. B. and Marti, K. (1975). Search for extinct 146Sm, 1. The isotopic abundance of 142Nd in the Juvinas meteorite. Earth Planet. Sci. Lett. 27, 79)84.

 

Luo, X., Rehkamper, M., Lee, D.-C. And Halliday, A. N. (1997). High precision 230Th/232Th and 234U/238U measurements using energy-filtered ICP magnetic sector multiple collector mass spectrometry. Int. J. Mass Spec. Ion Process. 171, 105–17.

 

McIntyre, G. A., Brooks, A. C. Compston, W and Turek, A. (1966). The statistical assessment of Rb)Sr isochrons. J. Geophys. Res. 71, 5459)68.

 

Manton, W. I. (1988). Separation of Pb from young zircons by single-bead ion exchange. Chem. Geol. (Isot. Geosci. Section) 73, 147)52.

 

Noble, S. R., Lightfoot, P. C. and Scharer, U. (1989). A new method for single-filament isotopic analysis of Nd using in situ reduction. Chem. Geol. (Isot. Geosci. Section) 79, 15)19.

 

Nier, A. O. (1940). A mass spectrometer for routine isotope abundance measurements. Rev. Sci. Instrum. 11, 212)16.

 

O’Nions, R. K., Hamilton, P. J. and Evensen, N. M. (1977). Variations in 143Nd/144Nd and 87Sr/86Sr ratios in oceanic basalts. Earth Planet. Sci. Lett. 34, 13)22.

 

O’Nions, R. K., Carter, S. R., Evensen, N. M. and Hamilton P. J. (1979). Geochemical and cosmochemical applications of Nd isotope analysis. Ann. Rev. Earth Planet. Sci. 7, 11)38.

 

Papanastassiou, D. A., Huneke, J. C., Esat, T. M. and Wasserburg, G. J. (1978). Pandora’s box of the nuclides. In: Lunar Planet. Sci. IX, Lunar Planet. Sci. Inst., Houston, Texas, pp. 859)61 (abstract).

 

Parrish, R. R. (1987). An improved micro-capsule for zircon dissolution in U)Pb geochronology. Chem. Geol. (Isot. Geosci. Section) 66, 99)102.

 

Parrish, R. R. and Krogh, T. E. (1987). Synthesis and purification of 205Pb for U)Pb geochronology. Chem. Geol. (Isot. Geosci. Section) 66, 103)10.

 

Patchett, P. J. (1980). Sr isotopic fractionation in Ca)Al inclusions from the Allende meteorite. Nature 283, 438)41.

 

Patchett, P. J. and Tatsumoto, M. (1980). A routine high-precision method for Lu)Hf isotope geochemistry and chronology. Contrib. Mineral. Petrol. 75, 263)7.

 

Pietruszka, A. J., Carlson, R. W. and Hauri, E. H. (2002). Precise and accurate measurement of 226Ra–230Th–238U disequilibria in volcanic rocks using plasma ionization multicollector mass spectrometry. Chem. Geol. 188, 171–91.

 

Potts, P. J. (1987). Handbook of Silicate Rock Analysis. Blackie, 622 p.

 

Powell, R., Hergt, J. and Woodhead, J. (2002). Improving isochron calculations with robust statistics and the bootstrap. Chem. Geol. 185, 191–204.

 

Powell, R., Woodhead, J. and Hergt, J. (1998). Uncertainties on lead isotope analyses: deconvolution in the double-spike method. Chem. Geol. 148, 95–104.

 

Richard, P., Shimizu, N. and Allegre, C. J. (1976). 143Nd/146Nd, a natural tracer: an application to oceanic basalts. Earth Planet. Sci. Lett. 31, 269)78.

 

Roddick, J. C., Loveridge, W. D. and Parrish, R. R. (1987). Precise U/Pb dating of zircon at the sub-nanogram Pb level. Chem. Geol. (Isot. Geosci. Section) 66, 111)21.

 

Russell, W. A., Papanastassiou, D. A. and Tombrello, T. A. (1978). Ca isotope fractionation on the Earth and other solar system materials. Geochim. Cosmochim. Acta 42, 1075)90.

 

Shen, C.-C., Edwards, R. L., Cheng, H., Dorale, J. A., Thomas, R. B., Moran, S. B., Weinstein, S. E. and Edmonds, H. N. (2002). Uranium and thorium isotopic and concentration measurements by magnetic sector inductively coupled plasma mass spectrometry. Chem. Geol. 185, 165–78.

 

Smythe, W. R. and Mattauch, J. (1932). A new mass spectrometer. Phys. Rev. 40, 429-

 

Tanaka, T. and Masuda, A. (1982). The La)Ce geochronometer: a new dating method, Nature 300, 515)18.

 

Thirlwall, M. F. (1982). A triple-filament method for rapid and precise analysis of rare-earth elements by isotope dilution. Chem. Geol. 35, 155)66.

 

Thirlwall, M. F. (1991a). High-precision multicollector isotopic analysis of low levels of Nd as oxide. Chem. Geol. (Isot. Geosci. Section) 94, 13)22.

 

Thirlwall, M. F. (1991b). Long-term reproducibility of multicollector Sr and Nd isotope ratio analysis. Chem. Geol. (Isot. Geosci. Section) 94, 85)104.

 

Thirlwall, M. F. (2000). Inter-laboratory and other errors in Pb isotope analyses investigated using a 207Pb204Pb double spike. Chem. Geol. 163, 299)322.

 

Thirlwall, M. F. (2002). Multicollector ICP-MS analysis of Pb isotopes using a 207Pb–204Pb double spike demonstrates up to 400 ppm/amu systematic errors in Tl-normalization. Chem. Geol. 184, 255–79.

 

Thirlwall, M. F. and Walder, A. J. (1995). In situ hafnium isotope ratio analysis of zircon by inductively coupled plasma multiple collector mass spectrometry. Chem. Geol. 122, 241–7.

 

Thompson J.J. (1913). Rays of Positive Electricity and their Application to Chemical Analysis, Longmans, Green and Co. Ltd.

 

Titterington, D. M. and Halliday, A. N. (1979). On the fitting of parallel isochrons and the method of maximum likelihood. Chem. Geol. 26, 183)95.

 

Todt, W., Cliff, R. A., Hanser, A. and Hofmann, A. W. (1996). Evaluation of a 202Pb –  205Pb double spike for high-precision lead isotopic analysis. In: Basu, A. and Hart, S. R. (Eds.) Earth Processes: Reading the Isotopic Code. Geophys. Monograph 95, American Geophysical Union, pp. 429–37.

 

Vance, D. and Thirlwall, M. (2002). An assessment of mass discrimination in MC-ICPMS using Nd isotopes. Chem. Geol. 185, 227–40.

 

Walder, A. J., Abell, I. D., Freedman, P. A. and Platzner, I. (1993a). Lead isotopic ratio measurement of NIST 610 glass by laser ablation-inductively coupled plasma-mass spectrometry. Spectrochim. Acta 48B, 397–402.

 

Walder, A. J. and Furuta, N. (1993). High precision lead isotope ratio measurement by inductively coupled plasma multiple collector mass spectrometry. Anal. Sci. 9, 675-80.

 

Walder, A. J., Platzner, I. And Freedman, P. A. (1993b). Isotope ratio measurement of lead, neodymium and neodymium–samarium mixtures, hafnium and hafnium–lutetium mixtures with a double focussing multiple collector inductively coupled plasma mass spectrometer. J. Anal. Atomic. Spectrom. 8, 19-23.

 

Wasserburg, G. J., Jacobsen, S. B., DePaolo, D. J. McCulloch, M. T. and Wen, T. (1981). Precise determination of Sm/Nd ratios, Sm and Nd isotopic abundances in standard solutions. Geochim. Cosmochim. Acta 45, 2311)23.

 

Wendt, I. and Carl, C. (1991). The statistical distribution of the mean squared weighted deviation. Chem. Geol. (Isot. Geosci. Section) 86, 275)85.

 

White, W. M., Albarede, F. and Telouk, P. (2000). High-precision analysis of Pb isotope ratios by multi-collector ICP-MS. Chem. Geol. 167, 257)270.

 

York, D. (1966). Least-squares fitting of a straight line. Can. J. Phys. 44, 1079)86.

 

York, D. (1967). The best isochron. Earth Planet. Sci. Lett. 2, 479)82.

 

York, D. (1969). Least-squares fitting of a straight line with correlated errors. Earth Planet. Sci. Lett. 5, 320)4.