Abstract

The study of the evolution of the conduction band in dense gases and supercritical fluids near the critical point has been complicated by a lack of precise experimental measurements. Both photoemission from an electrode immersed in the fluid and field ionization of a molecule doped into the fluid have been used to probe solvent density effects on the energy of an excess electron as a function of fluid number density and temperature. In this perspective, we present recent experimental results that show a strong critical point effect on the minimum conduction band energy near the critical density and temperature of a fluid. We also discuss the recent development of a new theoretical model that advances our understanding of the density and temperature dependence of the conduction band minimum in near critical point fluids.

1. Introduction

The critical point of a dense gas represents a state of matter under extreme conditions. (In argon, e.g., the critical point occurs at °C and  bar.) In studies of the evolution of the conduction band in insulators—from the low density gas to the liquid to the solid—the unique properties of the critical point have only recently been investigated, as will be described below. Our fundamental understanding of the behavior of electrons in an insulator follows from experimental and theoretical studies on conduction in solid state systems, where translational symmetry allows one to predict accurately the shape and electronic characteristics of both the valence and conduction bands. The energy of an electron at the bottom of the lowest conduction band is the minimum energy of the conduction band relative to the vacuum level. In comparison to this picture, a free electron injected into a perturber gas has an energy that can be predicted based on single scattering interactions between the electron and individual perturbers. However, as the density of a perturber gas increases towards fluid densities, the electron/perturber interaction cannot be treated as a single scattering problem. Unlike in a solid, however, there is no long range translational symmetry. Since the density of the fluid can be adjusted dramatically using both pressure and temperature, understanding the minimum energy of the conduction band of an electron in a fluid of density is a complex problem.

2. Minimum Conduction Band Energy in Dense Gases

Two general methods [16] are available to measure the minimum energy of the conduction band in a fluid, namely, photoemission from an electrode placed in the fluid and photoionization of a dopant perturbed by the fluid. Photoemission allows for the direct measurement of if the photoemission threshold of the metal in a vacuum and in the fluid can be accurately determined. However, contamination of the metal surface and coating of the electrode surface by the fluid make accurate correction of photoemission data problematic (cf. Figure 1(a)). The second method involves the measurement of the photoionization spectrum of a dopant perturbed by fluid. The exact ionization energy of a dopant in a perturbing fluid can be difficult to determine, because the perturbing fluid broadens and shifts the ionization threshold as the density increases. This problem can be overcome by measuring the photoionization spectrum of the dopant at two different electric field strengths, since an electric field ionizes high- Rydberg states converging to the ionization threshold. Thus, taking the difference between two photoionization spectra measured with different applied electric fields (after intensity normalizing the spectra to remove the effects of secondary ionization) yields a peak that represents the Rydberg states ionized by the high electric field but not by the low field . The energy shift of the field ionization peak is given by [1, 2, 4, 713] where is a dopant-dependent field ionization constant. The zero-field dopant ionization energy is, therefore, [1, 2, 4, 713] where is the energy of the maximum of the dopant field ionization peak at the perturber density . Although dopant field ionization allows one to determine the dopant ionization threshold accurately, this ionization energy is not the minimum energy of the conduction band of the perturbing fluid. The relationship between the perturber-shifted dopant ionization energy and is [1, 2, 4, 713] where is the ionization energy of the unperturbed dopant, and is the ensemble-averaged dopant core/perturber polarization energy, which can be calculated using standard statistical mechanical techniques [1, 2, 713]. Figure 1(b) shows an example of obtained from (3). Clearly, the results from (3) are consonant with those obtained from photoemission measurements. In fact, determined from the experimental has less scatter, thus making field ionization of high- dopant Rydberg states an ideal technique for the study of critical point effects (cf. Figures 1(b) and 2).

3. Near Critical Point Fluids

is that temperature above which distinct gas and liquid phases do not exist. Above , the single phase system is generally referred to as a supercritical fluid (SCF) if the pressure is above the critical pressure. In an SCF, the fluid density, viscosity, solubility, and diffusivity can be tailored by adjusting the pressure and temperature. Thus, SCFs constitute controllable, tunable solvents that are currently being employed in the synthesis of various nanoparticles [1416] having small size distributions with shorter preparation and purification times. The bulk controllability of SCFs has also been used in the pharmaceutical industry [17, 18] to create high structural purity organic crystals. The lack of a detailed understanding of the intermolecular interactions in SCFs, however, hampers modeling the more complex reaction dynamics involved in the synthesis of nanomaterials and pharmaceutical compounds in an SCF solvent.

The myriad of factors that affect chemical reactions in near critical point fluids, where the correlation length is essentially unbounded, make the modeling and understanding of chemical processes in these fluids troublesome [1921]. The difficulty of isolating process control factors in near critical point fluids dictates the study of simple model systems, and therefore we have focused on the energy of the quasifree electron in near critical point fluids. The minimum energy of the conduction band, obtained from (3), is shown as a function of reduced number density (with being the critical number density) for Ar in Figure 1(b) and for Kr, Xe, CH4, and C2H6 in Figure 2. Clearly, there is little or no temperature dependence in along noncritical isotherms (solid markers). However, does vary considerably near the critical temperature and density of the perturber fluid. The variation of in near critical point SCFs may well have a bearing on the changes in product distribution and reactivity [22] observed in chemical reactions carried out in such systems. Understanding why changes in a near critical point SCF requires a new model for conduction in a fluid, however.

4. Local Wigner-Seitz Model [1, 2, 713]

In the local Wigner-Seitz model, which is a recent modification of the original Springett, Jortner, and Cohen model [23], is given by where is the ensemble-averaged electron/perturber polarization energy, is the zero-point kinetic energy of the quasifree electron, and ( Boltzmann's constant) is the thermal energy of the quasifree electron. The average electron/perturber polarization energy is determined using standard statistical mechanical techniques [2, 7]: where is the perturber/perturber radial distribution function and is the electron/perturber interaction potential. Any critical point behavior in arises from variations in near the critical density and temperature.

The zero-point kinetic energy of the quasifree electron is obtained from solving the Schrödinger equation for the quasifree electron in a dense perturber. In (6), is a short-ranged potential that accounts for local dynamic polarization of a perturber by the optical electron, is the mass of the electron, and is the reduced Planck constant. The local potential satisfies an average translational symmetry where is the interaction range. At any density, the minimum distance between a quasifree electron with low kinetic energy and a single perturber is given by the absolute value of the scattering length . At high perturber number densities, under the assumption that interactions in the first solvent shell dominate the dynamics of the problem, the maximum distance for a short-ranged interaction is one-half of the spacing between two perturbers in the first solvent shell. This maximum spacing, otherwise known as the local Wigner-Seitz radius, is where is the maximum of . Thus, the interaction range for the local short-ranged potential is . By applying these conditions to the asymptotic solutions of (6), the zero-point kinetic energy of the quasifree electron is In (9), is the phase shift induced by the short-ranged potential, and the density dependence arises from . Thus, any critical point behavior in will result from a change in near the critical density along the critical isotherm due to density fluctuations, as illustrated in Figure 3(a) for Ar. Similarly, Figure 3(b) presents the maximum local density in the first solvent shell (i.e., ) of methane plotted as a function of the reduced perturber number density . Clearly, the enhancement of the local density at the critical temperature extends to number densities around the critical density, and not just to the critical density itself. (In fact, the critical isotherm turning points that bound the saddle point in the perturber phase diagram coincide with the number densities that delimit the deviations of from the noncritical isotherm values.)

Figure 4 shows the average electron/perturber polarization energy and the zero-point kinetic energy for methane at noncritical temperatures and on an isotherm near the critical isotherm of methane. Although shows a small critical point effect (Figure 4(a)), the critical point variation in is dominated by the effects on the zero-point kinetic energy of the quasifree electron (Figure 4(b)), which is dependent on the local density within the first solvent shell through the local Wigner-Seitz radius. Thus, as the local density in the first solvent shell increases, the interaction range for the electron and a single perturber decreases. This decrease in interaction box size increases near the critical point of the fluid, thus leading to the observed critical point effect.

The lines in Figures 1(b) and 2, which are calculated from the local Wigner-Seitz model described above, accurately fit the experimental data across the entire density range both for noncritical temperatures and on an isotherm near the critical isotherm. In conclusion, the development of dopant field ionization as a technique to extract the minimum energy of the conduction band, , with minimal scatter led to the discovery of a strong critical point effect. This discovery, in turn, has resulted in the development of a new model for electron conduction in dense fluids that represents an order of magnitude improvement over previous theoretical treatments [1, 2, 713].

Acknowledgments

The University of Wisconsin Synchrotron Radiation Center is funded by NSF DMR-0537588. This work was supported by grants from the Petroleum Research Fund (45728-B6), from the Professional Staff Congress-City University of New York (60074-34 35), and from the Louisiana Board of Regents Support Fund (LEQSF(2006-09)-RD-A33).