There are three areas in ICP emission spectrometry for which it would be beneficial to provide accurate simulation. The first is the generation of the physical environment which establishes the bulk properties of the plasma. The second is a group of processes which result in the transformation of analyte aerosol into atomic vapor; i.e. desolvation, vaporization and atomization. The third is generation of observed analyte emission signals, i.e. excitation and ionization of analyte atoms.
There are two steps involved in the computer simulation. The first is development of a suitable set of mathematical equations for description of all relevant processes. The second is the development of an algorithm for solution of the equations. The simulation of each of the three areas mentioned above is discussed in further detail below.
A recent comparison of results for a dry argon ICP from this model with experimental results obtained using Thomson and Rayleigh scattering methods shows agreement to within 20% for most areas of the plasma. In particular, the electron temperature and electron number density show very good agreement (10%) with experiment. The gas temperature shows fairly good agreement but has an apparent systematic error of approximately +1000K. Modifications to the model and examination of data collection and processing methods are underway to determine the cause of the discrepancy.
The droplet/particle model currently describes the loss of solvent (desolvation) and the subsequent vaporization of a solute particle, atomization of the solute vapor and partial ionization of the resulting analyte atoms. Single aerosol-droplet trajectories as well as multi-aerosol droplet profiles can be simulated for a variety of ICP operating conditions. Separate calculations have been carried out for pure solvent droplets, for pure solute particles, and for aerosol droplets containing one or a mixture of two solutes. The results of these studies are outlined below.
The effect of droplet size, ICP central-channel flow rate and applied rf power on single- droplet trajectories has been determined. The time (tcd) or height (hcd) of complete desolvation of a droplet is quadratically related to its diameter. This is because the desolvation process is controlled by heat transfer. In addition, the tcd is related roughly linearly to the central channel flow rate, and inversely to the applied rf power. These dependences reflect the changes in gas temperature that occur with changing flow rate and applied power. Multi-droplet profiles have been calculated for droplet distributions from several commercial nebulizer/spray-chamber combinations. These results indicate that, in the absence of any effects of the droplets on plasma parameters, the aerosol from a typical sample introduction system is almost completely desolvated at 30 mm above the load coil (ALC). Droplet histories obtained using this droplet model have recently been compared with those obtained by experiment. Experiments show that 30 Ám diameter droplets are completely desolvated at a height of 18 mm ALC, whereas the theory indicates that a 20 Ám droplet is desolvated at 18 mm ALC. Work is underway to determine the origin of the difference in results.
The effect of particle size, ICP central-channel flow rate and applied rf power on single- particle trajectories has also been determined for several solutes. The time (tcv) or height (hcv) of complete vaporization is linearly related to the particle diameter. The hcv is not related in any obvious way to the thermal properties of the solute, however, it is related to the central channel flow rate and applied power in the same fashion as was the hcd for droplet desolvation.
The changes in particle vaporization that occur on addition of a second solute to the sample have been examined using the same desolvation and vaporization model. Properties of the mixed solute are calculated from the individual solute properties and from the mole fractions of each solute. In order to compare directly with experimental results, radially resolved profiles have been calculated by including diffusion and flow of the analyte away from the central channel. With these inclusions, analyte number densities as a function of height and radial distance in the plasma have been simulated and compared to those obtained by experiment. Absolute values of number density obtained from the simulation are similar to those obtained by experiment; however simulated ion to atom number density ratios are in general lower than experimental values. Good correlation between simulated and experimentally observed changes in analyte densities as a result of the added element was obtained for Ca in the presence of Zn, Li and Ba. For some of the other mixed solute pairs (Ca with Na, Ag and K), however, there are notable differences between simulation and experiment. There are several possible reasons for the differences. First, the simulation package currently assumes that no changes in plasma properties take place on addition of an aerosol; experiments have been carried out which are in disagreement with this assumption. Second, the properties of the mixed solute may not be well represented by a linear combination of the individual properties weighted by their mole fractions. And, third, atomization and ionization are assumed to be the only chemical changes taking place for the analyte solute. Further, these processes are assumed to be at equilibrium and to be instantaneous, whereas they are likely under kinetic control. Any or all of these assumptions may be incorrect and could result in differences between simulated and experimental analyte number density profiles.
The ultimate goal of the ICP simulation work in this group is to develop a single model which accurately describes all of the three areas of ICP emission for any desired set of operating conditions. In order to achieve this goal, several changes must be made to the existing model. First, the aerosol-diffusion model must take into account radial movement of the analyte due to convection, i.e. forced movement. Second, the analyte excitation/ionization model must be expanded to include processes involving the added element. The complete aerosol model should then be integrated with the fundamental model so that reciprocal interactions could be examined. Reciprocal interactions refer to the effect of droplet desolvation, particle vaporization, atomization, ionization and excitation on the plasma. This software could be used to simulate, for example, spatial analyte emission profiles under a variety of operating conditions and in the presence of a variety of concomitant species. The simulated profiles could then be compared with experimental profiles to gain an improved fundamental understanding of the ICP.
 S. V. Patankar, "Numerical Heat Transfer and Fluid Flow",Hemisphere Publishing Corporation, New York, NY.
 P. Yang, R. M. Barnes, J. Mostaghimi and M. I. Boulos, Spectrochim. Acta Part B, 44B, 657-666 (1989).