Published in J. Environ. Qual. 33:276-284 (2004).
© ASA, CSSA, SSSA
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Organic Compounds in the Environment
Characterization of Cation
Interactions in Aqueous Solution Using Deuterium Nuclear Magnetic Resonance Spectroscopy
Dongqiang Zhu*,a,e,
Bruce E. Herbertb,
Mark A. Schlautmanc,d and
Elizabeth R. Carrawayd
a Department of Civil Engineering, Texas A&M University, College Station, TX 77840
b Department of Geology and Geophysics, Texas A&M University, College Station, TX 77843
c Department of Agricultural and Biological Engineering, Clemson University, Clemson, SC 29634-0357
d Department of Environmental Toxicology and the Clemson Institute of Environmental Toxicology, Clemson University, Pendleton, SC 29670
e Connecticut Agricultural Experiment Station, 123 Huntington Street, New Haven, CT 06504
* Corresponding author (Don.Zhu{at}po.state.ct.us).
Received for publication April 23, 2003.
 |
ABSTRACT
|
|---|
Chemical interactions of aromatic organic contaminants control their fate, transport, and toxicity in the environment. Recent molecular modeling studies have suggested that strong interactions can occur between the
electrons of aromatic molecules and metal cations in aqueous solutions and/or on mineral surfaces, and that such interactions may be important in some environmental systems. However, spectroscopic evidence for these so-called cation
interactions has been extremely limited to date. In this paper, cation
interactions in aqueous salt solutions were characterized via 2H nuclear magnetic resonance (NMR) spinlattice relaxation times (T1) and calculations of molecular correlation times (
c) for a series of perdeuterated (d6benzene) benzenecation complexes. The T1 values for d6benzene decreased with increasing concentrations of LiCl, NaCl, KCl, RbCl, CsCl, and AgNO3, with the largest effects observed in the AgNO3 and CsCl solutions. Upon normalizing
c values by solution viscosity effects, an overall affinity trend of Ag+ >> Cs+ > K+ > Rb+ > Na+ > Li+ was derived for the d6benzenecation complexes. The ability of Ag+ to complex d6benzene was significantly reduced upon addition of NH3, which strongly coordinates Ag+ at high pH. Results with d6benzene, d8naphthalene, d2dichloromethane, and d12cyclohexane in 0.1 M methanolic salt solutions confirmed that spinlattice relaxation rates are characterizing cation
interactions. The relatively strong cation
bonding observed between Ag+ and aromatic hydrocarbons probably results from covalent interactions between the aromatic
electrons and the d orbitals of Ag+, in addition to the normal electrostatic interaction.
Abbreviations: NMR, nuclear magnetic resonance PAH, polycyclic aromatic hydrocarbon T1, spinlattice relaxation time
c, molecular correlation time
 |
INTRODUCTION
|
|---|
POLYCYCLIC AROMATIC HYDROCARBONS (PAHs) are a group of ubiquitous hydrophobic organic contaminants (HOCs) that continue to receive considerable attention because of their toxicity, persistence, and extensive distribution in the environment (e.g., Means et al., 1980; Kubicki et al., 1999; Schwarzenbach et al., 2003). In general, the chemical state in which PAHs are present in the environment determines their fate, transport, and overall environmental effects (e.g., Schnoor, 1996; Kubicki et al., 1999; Schwarzenbach et al., 2003). At many contaminated sites, a complex suite of organic chemicals and inorganic elements are commonly present together (United States Department of Energy, 1992; National Research Council, 1994; Spiro and Stigliani, 1996). For example, polluted biogeochemical systems often contain PAHs, other HOCs (e.g., chlorinated hydrocarbons, pesticides), and heavy metals in addition to all naturally occurring chemical species (e.g., alkali and alkaline earth metals, trace metals, anions, organic matter, mineral surfaces). Characterization of the complex interactions between PAHs and cations in aqueous solution will therefore aid our understanding of the environmental fate and risk posed by these ubiquitous contaminants.
Relatively strong noncovalent bonds probably form between benzene and base cations in aqueous solution due to electrostatic attractions between the permanent quadrupole of benzene and the positively charged cations (Kumpf and Dougherty, 1993). The importance of these and other so-called cation
interactions has been recognized for some time in biological processes, including bimolecular recognition, proteinligand binding, and the selectivity of K+ within K+ channels in cell membranes (Miller, 1991; Dougherty, 1996; Ma and Dougherty, 1997; Gokel et al., 2001). The importance of cation
interactions from the perspective of environmental science and engineering, however, has gone largely unexplored. For example, it is reasonable to expect that cation
interactions may affect the distribution of PAHs between water and mineral surfaces in low organic carbon environments (Kubicki et al., 1999; Zhu et al., 2004). Cation
interactions may also influence the relative toxicities of PAHs and heavy metals that co-exist in aqueous solution, because their association with one another would probably alter their transport across cell membranes. However, few toxicity studies of co-contaminant mixtures of heavy metals and PAHs have been performed to date (e.g., van den Hurk et al., 1998a, 1998b, 2000; Babu et al., 2001), and none of these has tried to elucidate the potential importance of cation
effects on toxicity.
Spectroscopic evidence for cation
interactions in various organic solvents and in the gas or solid phase has been widely reported in the literature over the past several decades (e.g., Baddiel et al., 1965; Sunner et al., 1981; Munakata et al., 2000; Gokel et al., 2001 and references therein). However, spectroscopic evidence for cation
interactions in aqueous solution has been extremely limited. Instead, much of the supporting evidence for cation
interactions in aqueous solution has come from theoretical studies. For example, formation of aqueous-phase cation
complexes between benzene and base cations was first evaluated by Kumpf and Dougherty (1993) using an electrostatic model. Based on their calculations, the ordering K+ > Rb+ >> Na+, Li+ was predicted for the overall binding energies of cationbenzene complexes when the cation desolvation energies were considered. More recently, a molecular modeling study of benzene and selected PAHs in the aqueous and adsorbed states identified the cation
interaction as a potentially important environmental process (Kubicki et al., 1999). Microscopic characterization of weak solute interactions in aqueous solution has been technically challenging. Routine spectroscopic techniques such as 1H-NMR and Fourier-transform infrared (FTIR) are not able to provide useful information on cation
interactions in aqueous solution because of the relative weakness of the interactions and the insensitivity of these two spectroscopies. Therefore, no significant NMR chemical shifts or appearance of new absorbance bands are observed for aqueous-phase cation
complexes using these two conventional spectroscopic techniques. The strongest spectroscopic evidence to date for cation
interactions in aqueous solution comes from circular dichroism studies of synthetic cyclophane hosts, which contained clustered benzene rings able to compete with water molecules for cations via formation of cation
bonds (Shepodd et al., 1986; Forman et al., 1995). However, complete interpretation of circular dichroism spectra is indirect and complicated when detecting chirooptical properties of such molecules. Furthermore, the method is not applicable to studies of compounds such as benzene, which have rigid molecular structures and no chirooptical properties.
Spinlattice relaxation times (T1) can be used to obtain detailed information on molecular interactions (Smith, 1983). Deuterium NMR has the advantage of being sensitive to relatively weak solute interactions, and is well suited for NMR relaxation studies because deuterium relaxation is dominated by the quadrupole relaxation mechanism (Smith, 1983). This allows for a relatively simple correlation between the relaxation rate and molecular correlation time (
c). For example, complexation of CD3CN by various cations was elucidated based on 14N- and 2H-NMR T1 values (Bopp, 1967). Noncovalent interactions between perdeuterated monoaromatic compounds and geosorbents have been characterized using measured T1 and T2 (spinspin relaxation times) values (Nanny and Maza, 2001; Zhu et al., 2003, 2004). Because the cation
interaction affects the
c of benzene and other aromatic compounds, the interaction can be characterized by quantifying the T1 of perdeuterated compounds, which is inversely related to
c.
It is widely assumed that entropic effects control the distribution of aromatic molecules in biogeochemical systems (Schwarzenbach et al., 2003). We hypothesize, however, that specific covalent and noncovalent interactions such as cation
interactions may be important in influencing aromatic molecule distributions and reactivities in the environment. The specific objectives of this investigation were to (i) test the utility of 2H-NMR relaxation measurements for probing the interactions between aromatic molecules and cations in solution and (ii) evaluate the influence of different cations on the formation of cation
interactions. We used 2H-NMR to measure T1 values of d6benzene and other perdeuterated probes in a variety of aqueous and methanolic electrolyte solutions. We then calculated viscosity-normalized
c values from the T1 values to separate the cation
interactions from normal bulk solution viscosity effects.
 |
MATERIALS AND METHODS
|
|---|
Materials
The I A group monovalent cations (Li+, Na+, K+, Rb+, Cs+) varying in softness (i.e., positively correlated to cation radius) and a representative soft transition metal (Ag+) were selected for studying their cation
interactions with benzene. Perdeuterated compounds including d6benzene (Cambridge Isotope Laboratories, Andover, MA), d8naphthalene (Aldrich Chemical, St. Louis, MO), d2dichloromethane (Aldrich Chemical), d12cyclohexane (Aldrich Chemical), d4methanol (Aldrich Chemical), and d2water (Aldrich Chemical) were used as received, as were LiCl (Fisher Scientific, Hampton, NH), NaCl (EM Science, Gibbstown, NJ), KCl (EM Science), RbCl (Aldrich Chemical), CsCl (Fisher Scientific), NaNO3 (EM Science), AgNO3 (Fisher Scientific), and 1.0 M NH4OH (Fisher Scientific). All solutions were prepared using high performance liquid chromatography (HPLC)-grade methanol (Fisher Scientific) and/or double-distilled water.
Sample Preparation
Aqueous solutions of LiCl, NaCl, KCl, RbCl, CsCl, and AgNO3 ranging from 0 to 3 M were prepared. Diammonium-silver solutions ranging from 0.1 to 0.5 M were prepared by mixing and diluting AgNO3 and 1.0 M NH4OH at a mole ratio of 1:2 (AgNO3 and NH4OH). Methanolic solutions of 0.1 M NaCl, CsCl, NaNO3, and AgNO3 were prepared by shaking the solutions vigorously. The pH of all aqueous solutions except those containing Ag+ was adjusted to 7 using 0.5 M HCl and/or 0.5 M NaOH unless otherwise indicated. After transferring 5 mL of the solution of interest to a 10-mm broadband, 17.8-cm (7-in) NMR tube (Wilmad, Buena, NJ), 5 µL (0.1% v/v, 0.0110.025 M for different compounds) perdeuterated solute or 0.0136 g (0.020 M) d8naphthalene was added and then the tube was immediately sealed with silicone stoppers (Fisher Scientific). These samples were kept in the dark at ambient temperature and shaken overnight (812 h) to allow full equilibration before NMR analysis.
Deuterium Nuclear Magnetic Resonance Experiments
Most NMR spectra were recorded at ambient temperature (19.0 ± 0.5°C) using a Varian (Palo Alto, CA) XL-200 FT-NMR spectrometer operating at 30.7 MHz with a 10-mm broadband probe. A Varian 400 FT-NMR spectrometer was used to run RbCl samples at 61.3 MHz with a 5-mm switchable probe. After locking and shimming on a reference sample containing pure d2water, the spectrometer was run unlocked following a standard inversion-recovery pulse sequence after carefully shimming for the NMR signal of interest. A recycle delay of at least 5T1 was applied, and at least 10 transients were used to give a good signal-to-noise ratio. Spectra were processed with an exponential multiplication corresponding to line broadening of 2 Hz. Chemical shifts were referenced to the natural abundance of deuterium in water or to d4methanol spiked into the methanolic salt samples. The value for T1 was calculated by an exponential regression with an associated error reported by the instrument.
 |
RESULTS AND DISCUSSION
|
|---|
Experimental Nuclear Magnetic Resonance Results
No detectable chemical shift (i.e., >0.05 ppm) was observed for any deuterated probe in the different aqueous or methanolic solutions. For mobile isotropic fluids in the extreme narrowing case, the spinlattice relaxation rate (1/T1) of deuterium can be expressed as (Abragam, 1961; Mantsch et al., 1977):
 | [1] |
where
c is the molecular rotational correlation time,
is the asymmetry parameter of the electric field gradient and is very close to zero for a C2H fragment (Mantsch et al., 1977), and e2Qq/ℏ is the quadrupolar coupling constant in rad s1, equal to 196.5 ± 1.3 kHz for benzene (Diehl and Khetrapal, 1969), 160 ± 3 kHz for dichloromethane (Glasel, 1969), and 165 ± 2 kHz (assumed uncertainty) for methanol (Glasel, 1969). As shown in Eq. [1], an inverse relationship exists between
c and the spinlattice relaxation time T1. Figure 1a
shows T1 values of d6benzene for LiCl, NaCl, KCl, RbCl, CsCl, and AgNO3 aqueous solutions. Viscosity can also affect the solute correlation time, which in turn affects the measured T1 of that solute. Assuming the simplest case of spherical molecules in a viscous medium, the relationship between the dynamic viscosity (
) of a solution and
c is given by the StokesEinstein hydrodynamic formula (Glasel, 1969)
 | [2] |
where f is a microviscosity factor that accounts for influence of a noncontinuous medium on the diffusion of molecules, k is the Boltzmann constant, and T is the absolute temperature. In general, the value of f is set equal to 1/6 or 1/12 (Mantsch et al., 1977). The molecular volume is calculated from pure liquid properties using 0.74 (M/
N), assuming hexagonal close-packed spherical molecules where N is Avogadro's number, M is the molecular weight, and
is the density. By comparing Eq. [1] and [2], it can be seen that T1 and
are inversely related. This inverse relationship is particularly evident when the nuclear spin relaxation is dominated by molecular rotation.
Dynamic viscosities of various aqueous electrolyte solutions are given in Table 1. A decrease in T1, which indicates an increase in
c, is observed with increasing salt concentration for all salts, with the two largest variations occurring in AgNO3 and CsCl solutions. The relationship between d6benzene
c values and salt concentration cannot be fully explained by solution viscosity effects. For example, at a salt concentration of 3.0 M, the LiCl solution has the highest viscosity while the CsCl solution has the lowest viscosity. Therefore, a higher T1 value of d6benzene in the CsCl versus the LiCl solution would be expected if
c were only affected by viscosity. This can be further demonstrated by comparing T1 values of d6benzene in the NaCl and AgNO3 solutions. Although the viscosity of the NaCl solutions is always higher than that for AgNO3 solutions of equal salt concentration, T1 values in AgNO3 solutions are substantially lower than in NaCl solutions, a result that would not be observed if solution viscosity alone dominated the T1 of d6benzene. These results suggest that interactions between cations and benzene, characterized as cation
interactions, are responsible for the observed T1 difference between different electrolyte solutions.
To isolate the nonviscosity effects on
c values, a normalized molecular rotational correlation time can be defined by dividing the
c value calculated from Eq. [1] by the
c value calculated from Eq. [2] (Fig. 1b). In using Eq. [2], we set the value of f to be 1/12. The strength of cation
interactions for d6benzene is positively correlated to the normalized
c, which has the sequence Ag+ >> Cs+ > K+ > Rb+ > Na+ > Li+. This result is consistent with calculations from the electrostatic model in literature (Kumpf and Dougherty, 1993), where a binding energy sequence of K+ > Rb+ >> Na+, Li+ was found between benzene and cations in aqueous solution (see below). It is interesting to note that the cation
interaction for Cs+ is stronger than that for K+ and Rb+, indicating that non-electrostatic forces play an important role in the cation
interaction between benzene and Cs+. This most likely can be attributed to the relatively high polarizability of Cs+ (Israelachvili, 1991). For this particular case, an electrostatic model treating the cation as a point charge cannot correctly describe the cation
interaction. Figure 1b also shows that normalized
c values for benzene decreased slightly with increasing LiCl and NaCl concentrations, whereas increasing trends were observed for all other electrolyte solutions. The decreasing trend for LiCl and NaCl is unexplained, because we expected a constant normalized
c value if a cation exhibits no interaction with benzene. However, our results for LiCl and NaCl most likely can be explained by shortcomings of the simple StokesEinstein model (Eq. [2]) used to correct for viscosity effects. For example, we selected a value of 1/12 for f in Eq. [2]; as noted above, researchers generally select either a value of 1/6 or 1/12 when using this equation but there is no reason to believe that f should be constant with increasing salt concentration. In addition, the assumption that the planar benzene molecule can be adequately represented by a simple spherical geometry in Eq. [2] is questionable. Finally, one could argue that the appropriate volume (or diameter) to use in Eq. [2] should be one that considers both complexed and uncomplexed benzene molecules. From a practical standpoint, the decreasing trends observed in Fig. 1b for LiCl and NaCl suggest that the other cations are similarly affected by these shortcomings in our viscosity correction, and that their interactions with benzene are, in fact, actually stronger than shown.
Compared with the other cations tested, Ag+ exhibited a much stronger capability to decrease T1 (Fig. 1a), and increase the normalized
c (Fig. 1b) of d6benzene in aqueous solution, indicating a relatively strong cation
interaction between Ag+ and benzene. A strong affinity between Ag+ and aromatic hydrocarbons was reported for many years (Andrews and Keefer, 1949) in aqueous solubility studies. The investigators observed that benzene and phenanthrene solubilities increased by about 4 and 10 times, respectively, in 1 M AgNO3 versus KNO3 solutions. Andrews and Keefer postulated the presence of two water-soluble complexes, AgAr+ and Ag2Ar+, to explain their solubility data but offered no other independent evidence to corroborate this hypothesis. Recently, studies have shown that Ag+ has a very high affinity for some spherical hydrocarbons containing several benzene rings that form a cavity to host Ag+ (Gross et al., 1995; Munakata et al., 2000). However, no spectroscopic evidence of the strong interaction between Ag+ and benzene in aqueous solution has been reported previously.
Figure 1c shows that the normalized
c of d2dichloromethane increases with increasing concentrations of CsCl and AgNO3, but not NaCl. The binding energy sequence Cs+ > Ag+ > Na+ can be derived for dichloromethanecation interactions by comparing the normalized
c values. Considering the relatively small polarity of dichloromethane compared with water, dichloromethanecation interactions should have much smaller interaction energies than the corresponding cation hydration energies. Thus, water can easily out-compete dichloromethane for cations in aqueous solution, and the strongest electrostatic interaction between cations and dichloromethane would be expected for the cation with the lowest hydration energy (i.e., 302 kJ mol1 for Cs+, versus 496 and 422 kJ mol1 for Ag+ and Na+, respectively; Burgess, 1988). Because both Cs+ and dichloromethane have relatively high polarizabilities (Israelachvili, 1991), non-electrostatic forces between Cs+ and dichloromethane due to mutual polarization can also be important. Thus, the dichloromethaneCs+ complex has the highest binding energy. Similarly, Ag+ interacts with dichloromethane more strongly than Na+ because the former has a larger radius and is more polarizable, even though it has an overall higher energy of hydration than does Na+.
Figure 2
shows normalized
c values of perdeuterated benzene, dichloromethane, and methanol (CD3 group) in CsCl aqueous solutions. It can be seen that benzene interacts more strongly with Cs+ than does dichloromethane; a relationship is also observed with the other cations (Fig. 1). This is because the intensity of benzenecation interactions is higher than that for dichloromethanecation interactions. Figure 2 also shows that the interaction between methanol and Cs+ is very weak. One reason for this is that the "desolvation penalty" for breaking hydrogen bonds between methanol and water is very high. In addition, because methanol is much less polarizable than benzene and dichloromethane, the contribution from mutual polarization to methanolCs+ interactions is negligible. Therefore, little interaction between methanol and cations in aqueous solution would be expected.

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 2. Normalized molecular correlation time ( c) values for d6benzene, d2dichloromethane, and d4methanol (CD3 group) in aqueous solution as a function of CsCl concentration.
|
|
To gain further spectroscopic evidence for the cation
interaction between Ag+ and benzene, T1 values of d6benzene in 1:2 (mol mol1) AgNO3NH4OH solutions were measured (Fig. 3a)
. Normalized
c values were calculated for AgNO3NH4OH solutions assuming that their viscosities were equivalent to the corresponding AgNO3 solutions. This assumption was based on independent measurements of d2dichloromethane T1 values in AgNO3 and AgNO3NH4OH solutions, which were indistinguishable at any given Ag(I) concentration (data not shown). As can be seen in Fig. 3b, normalized
c values decreased upon addition of NH4OH. Thus, the cation
interaction between Ag+ and benzene is impaired due to coordination of the Ag+ by NH3; note that for the conditions used here, Ag+ is strongly coordinated by two NH3 molecules [i.e., pH above 9.9 results in over 96% Ag
+2 chemical speciation; Smith and Martell, 1981; Allison et al., 1991]. The fact that different pH conditions (i.e., pH of approximately 5 to 6 and approximately 10 for AgNO3 and AgNO3NH4OH solutions, respectively) were used for the two studies shown in Fig. 3 is not the reason for their different T1 values, as verified by the negligible effect of pH on T1 of d6benzene (Fig. 4)
. Recent studies have suggested that cation
interactions may exist between benzene and transition-metal-complexed cations such as Co
3+6 (Zaric, 1999). Therefore, although the cation
interaction may still exist here between Ag
+2 and benzene, the interaction must be much weaker than that between benzene and Ag+; however, it is likely the reason why the normalized
c does decrease slightly as the concentration of Ag
+2 increases (Fig. 3b). A comparison between Fig. 1b and 3b shows that the normalized
c values for d6benzene and Ag
+2 fall between those for K+ and Rb+, indicating that the intensity of the cation
interaction between benzene and Ag
+2 is decreased to the level of Rb+ and K+benzene complexes when Ag+ is strongly coordinated by NH3.

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 3. Relaxation of d6benzene in aqueous solution as a function of Ag(I) concentration. (a) Spinlattice relaxation times (T1). (b) Normalized molecular correlation times ( c).
|
|

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 4. Spinlattice relaxation times (T1) of d6benzene in 0.1 M NaCl aqueous solution as a function of pH.
|
|
In a final test of cation
interactions for additional probes, T1 values of d6benzene, d8naphthalene, d2dichloromethane, and d12cyclohexane were measured in 0.1 M methanolic salt solutions (Fig. 5)
. The purpose of using methanol as the solvent was to circumvent the relatively low sensitivity of NMR spectroscopy by enhancing the solubility of these probes. As expected, T1 values of d6benzene and d8naphthalene were lower in the presence of Ag+ than for the other cations. However, for d2dichloromethane and d12cyclohexane, which are incapable of forming cation
bonds, no significant difference in T1 was observed among the different cations. The result for cyclohexane is consistent with the computational result that cyclohexane is not a better cation binder than benzene even though it is more polarizable than benzene (Dougherty, 1996). It also provides further experimental evidence that the electrostatic component is very important in cation
interactions. However, although the electrostatic contribution dominates the cation
interaction for base cations with relatively simple electron configurations such as K+ (Kumpf and Dougherty, 1993; Dougherty, 1996), the interaction between the Ag+ d orbitals and benzene apparently plays an important role in the cation
interaction, based on the experimental results obtained here. It can be further concluded that for transition metals with d electron configurations, the cation
interaction cannot be completely accounted for if only electrostatic effects are considered. In this case, it is not appropriate to view the cation
interaction as only an electrostatic bond. Examining the ionic and covalent indices of the cations tested in this study (Table 2), it is clear that cations favoring cation
interactions are soft cations with a relatively high covalent index and a low ionic index (e.g., Ag+, Cs+). Viewed from the Lewis acidbase concept, benzene, a very soft base, can have a relatively strong interaction with soft cations (i.e., soft acids) in aqueous solution. Based on a more theoretical model, Caldwell and Kollman (1995) pointed out that the cation
interaction can be more exactly quantified when electrostatic effects and polarization and dispersion contributions are all considered. For soft cations with relatively large ionic radii such as Ag+ and Cs+, the contribution of polarization and dispersion of the cation may be very important to the overall cation
interaction.

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 5. Spinlattice relaxation times (T1) for perdeuterated solutes in 0.1 M methanolic salt solutions: (a) d6benzene, (b) d8naphthalene (white bars for peaks and black bars for ß peaks), (c) d12cyclohexane, and (d) d2dichloromethane.
|
|
Quantitative Comparison with Literature Data
To confirm that our NMR relaxation measurements were quantitatively, as well as qualitatively, characterizing cation
interactions, we compared them with results available from the literature. A previous study calculated the fraction (A) of humic-associated phenol from the 13C T1 relaxation times using the formula (Bortiatynski et al., 1997):
 | [3] |
where T1F and T1A are the relaxation times for the free and associated species, respectively, and T1 is the observed value averaged on these two. Similarly, the fraction of benzene associated with cations can be calculated by assuming 1:1 M+benzene complexation (where M is a cation). To remove viscosity effects from our calculation, normalized
c values should be used instead of relaxation times:
 | [4] |
where
cF and
cA are the normalized
c values for the free and associated benzene, respectively. The term
cF can be calculated from the observed T1 value for d6benzene in pure water, and
cA can be obtained from model fitting of the binding constant at different salt concentrations. Using this approach, the binding constant for a 1:1 Ag+benzene complex was determined to be 2.26 ± 0.36 M1, which is very close to the value of 2.41 M1 determined by Andrews and Keefer in their solubility enhancement study (Andrews and Keefer, 1949). A 2:1 Ag+benzene complex with a binding constant of 0.212 M2 was also reported by those investigators, but was neglected here because it is relatively less important.
Using Monte Carlo simulations based on statistical perturbation theory, Kumpf and Dougherty (1993) calculated the free energy (
G) needed to exchange K+ with Li+, Na+, and Rb+ in M+benzene complexes in aqueous solution. As shown in Fig. 6
, a very good correlation is observed between the difference in normalized
c values (
c) for the maximum salt concentrations minus no salt and Kumpf and Dougherty's 
G values for the different cations. This result indicates that the change in normalized
c values observed with increasing salt concentration is a direct measure of the cation
binding energy. It should be noted that higher binding constants would be expected for PAHs compared with benzene because their higher
electron densities favor cation
interactions.
Environmental Implications
2H NMR relaxation appears to be a useful technique for characterizing cation
interactions in aqueous solution, and we have recently extended its application to include mineral sorbent systems (Zhu et al., 2004). One shortcoming of the technique is that high aqueous salt concentrations had to be used here to observe the cation
interactions in solution because of their relative weakness and the experimental limitations of NMR (i.e., high detection limit restricts probing solutes [e.g., PAHs] having higher
-electron densities but low aqueous solubilities than benzene; relatively low sensitivity compared with other spectroscopic techniques such as fluorescence quenching [J.-H. Lee, personal communication, 2003]). However, despite these limitations, the results obtained in the present study are important for environmental systems. For example, the present aqueous solution results provide a baseline against which other, more environmentally relevant systems can be compared (Zhu et al., 2004). In addition, equilibrium constants that are determined experimentally in these high salt solutions can be extrapolated down to more realistic environmental solution conditions, or even infinite solution conditions, merely by using an appropriate activity correction model. It is also important to note that some environmentally relevant aqueous systems do contain relatively high concentrations of salts (e.g., brines, brackish underground water, estuaries, oceans). It is also well known that minerals can serve to concentrate cations at their solidwater interfaces because of their negative surface charges. For the latter item, this can happen by actual cation adsorption on mineral surfaces and/or by electrostatic attraction of cations into the electric double layer region. For example, Stumm and Morgan (1981) demonstrated with a simple calculation based on GuoyChapman theory that a low surface charge clay such as kaolinite will increase the concentration of a monovalent cation in the interfacial region by two orders of magnitude over its bulk concentration. The concentrating effects of cations in the electric double layer region by negatively charged surfaces is even higher for bivalent cations and/or for higher-charged mineral surfaces (e.g., vermiculite). The concentrating effect of cations by adsorption on mineral surfaces, especially those having high cation exchange capacities (CEC), can be very significant. For example, the CEC can be as high as 1200 to 1500 mmol kg1 for vermiculite (Bohn et al., 1985) and 400 to 500 mmol kg1 for zeolite (Dixon and Weed, 1989). Other negatively charged materials that can collect and concentrate cations include ionized natural organic matter (NOM) and anionic micellessurfactants that have been proposed for use in soilaquifer remediation schemes. Because these interfacial regions have much smaller volume compared with the bulk aqueous phase, the cations collected within become highly concentrated.
Development of "smart" materials that can efficiently and selectively bind the chemical species of interest is the current trend in modern system design. Cation
interactions have been used as important processes in some engineered treatment systems. For example, it has been reported that calix arenes can selectively remove certain cations such as Ag+ and Cs+ from aqueous solution through the formation of cation
complexes (Ikeda and Shinkais, 1994, 1997; Bocchi et al., 1995). The relatively strong cation
bonding capability of Ag+ and Cu+ has also been used to prepare cation-charged minerals as adsorbents for effective olefin separation (Takahashi et al., 2001; Padin and Yang, 2002). More recently, photooxidation of PAHs via the formation of cationPAH complexes with Nb, V, and Ti in the gas phase using laser vaporization has been reported (Duncan et al., 2001; Foster et al., 2001).
Cation
interactions are also likely to be important in some environmental toxicity and bioavailability studies. For example, intracellular and extracellular cation concentrations are high enough (e.g., 0.14 M intracellular K+; Alberts et al., 1994) to allow for relatively intensive cation
interactions with PAH molecules, which may affect their transport across cell membranes and thereafter metabolism reactions. Despite the small thermodynamic gradient driving force, long-term evolutionary changes have enabled biological systems to selectively and efficiently extract cations from the aqueous phase (i.e., 1000:1 K+ over Na+ selectivity by K+ channels; Kumpf and Dougherty, 1993; Ma and Dougherty, 1997) using the underlying mechanism of cation
interactions. Because transition metals are expected to have even higher cation
interaction intensities compared with base cations (e.g., K+), they can be more efficiently extracted and enriched in biological systems by ion channels. However, further studies will be necessary to better understand the potential influence of cation
interactions on issues related to environmental toxicology and bioavailability.
 |
CONCLUSIONS
|
|---|
Cation
interactions between benzene and cations in aqueous solution were characterized via 2H NMR relaxation times and calculations of molecular correlation times. Soft base cations and soft transition metals (e.g., Cs+, Ag+) have stronger cation
interactions compared with hard base cations (e.g., Li+, Na+). Results with d6benzene, d8naphthalene, d2dichloromethane, and d12cyclohexane in 0.1 M methanolic salt solutions confirmed that spinlattice relaxation rates are characterizing cation
interactions. Overall findings from this study provide a methodology for investigating and quantifying cation
interactions and present preliminary experimental evidence suggesting that these interactions are likely to be important in some environmental systems.
 |
ACKNOWLEDGMENTS
|
|---|
We gratefully acknowledge the constructive comments of Dr. Linda Lee and two other anonymous reviewers and their suggestions for improving this paper. The authors also gratefully acknowledge financial support from the National Science Foundation (CTS-0096053) and the gift of perdeuterated aromatic hydrocarbons from Cambridge Isotope Laboratories.
 |
REFERENCES
|
|---|
- Abragam, A. 1961. The principles of nuclear magnetism. Oxford Univ. Press, London.
- Alberts, B., D. Bray, J. Lewis, M. Raff, K. Roberts, and J.D. Watson. 1994. Molecular biology of the cell. 3rd ed. Garland Publ., New York.
- Allison, J.D., D.S. Brown, and K.J. Novo-Gradac. 1991. MINTEQA2/PRODEFA2, A geochemical assessment model for environmental systems. Version 3.0 user's manual. EPA/600/3-91/021. USEPA, Office of Res. and Development, Washington, DC.
- Andrews, L.J., and R.M. Keefer. 1949. Cation complexes of compounds containing carbon-carbon double bonds. 4. The argentation of aromatic hydrocarbons. J. Am. Chem. Soc. 71:36443647.
- Anonymous. 1983. Handbook of chemistry and physics. CRS Press, Boca Raton, FL.
- Babu, T.S., J.B. Marder, S. Tripuranthakam, D.G. Dixon, and B.M. Greenberg. 2001. Synergistic effects of a photooxidized polycyclic aromatic hydrocarbon and copper on photosynthesis and plant growth: Evidence that in vivo formation of reactive oxygen species is a mechanism of copper toxicity. Environ. Toxicol. Chem. 20:13511358.[Web of Science][Medline]
- Baddiel, C.B., M.J. Tait, and G.J. Janz. 1965. Nonaqueous silver nitrate solutions. Spectral studies in acetonitrile. J. Phys. Chem. 69:36343638.
- Bocchi, C., M. Careri, A. Casnati, and G. Mori. 1995. Electivity of calix[4]arene crown-6 for cesium ion in ISEEffect of the conformation. Anal. Chem. 67:42344238.
- Bohn, H.L., B.L. McNeal, and G.A. O'Connor. 1985. Soil chemistry. 2nd ed. John Wiley & Sons, New York.
- Bopp, T. 1967. Magnetic resonance studies of anisotropic molecular rotation in liquid acetonitrile-d3. J. Chem. Phys. 47:36213626.
- Bortiatynski, J.M., P.G. Hatcher, and R.D. Minard. 1997. The development of 13C labeling and 13C NMR spectroscopy techniques to study the interaction of pollutants with humic substances. p. 2650. In M.A. Nanny, R.A. Minear, and J.A. Leenheer (ed.) Nuclear magnetic resonance spectroscopy in environmental chemistry. Oxford Univ. Press, New York.
- Burgess, J. 1988. Ions in solution: Basic principles of chemical interactions. Halsted Press, New York.
- Caldwell, J.W., and P.A. Kollman. 1995. Cation-
interactionsNonadditive effects are critical in their accurate representation. J. Am. Chem. Soc. 117:41774178.
- Diehl, P., and C.L. Khetrapal. 1969. Study of isotopic effects and determination of quadrupole coupling constant from nuclear magnetic resonance spectra of monodeuterobenzene in a nematic phase. Can. J. Chem. 47:14111416.
- Dixon, J.B., and S.B. Weed (ed.) 1989. Minerals in soil environments. SSSA Book Ser. 1. 2nd ed. SSSA, Madison, WI.
- Dougherty, D.A. 1996. Cation-
interactions in chemistry and biology: A new view of benzene, Phe, Tyr, and Trp. Science (Washington, DC) 271:163168.[Abstract]
- Duncan, M.A., A.M. Knight, Y. Negishi, S. Nagao, K. Judai, A. Nakajima, and K. Kaya. 2001. Photoelectron spectroscopy of V-x(coronene)(y) and Ti-x(coronene)(y) anions. J. Phys. Chem. A 105:1009310097.
- Forman, J.E., R.E. Barrans, and D.A. Dougherty. 1995. Circular-dichroism studies of molecular recognition with cyclophane hosts in aqueous-media. J. Am. Chem. Soc. 117:92139228.
- Foster, N.R., J.W. Buchanan, N.D. Flynn, and M.A. Duncan. 2001. Ring destruction and carbide formation in niobium-PAH complexes. Chem. Phys. Lett. 341:476482.
- Glasel, J.A. 1969. Deuteron magnetic relaxation times and molecular and intermolecular motion in some organic liquids. J. Am. Chem. Soc. 91:45694571.
- Gokel, G.W., L.J. Barbour, S.L. De Wall, and E.S. Meadows. 2001. Macrocyclic polyethers as probes to assess and understand alkali metal cation-
interactions. Coord. Chem. Rev. 222:127154.
- Gross, J., G. Harder, F. Vögtle, H. Stephan, and K. Gloe. 1995. C60H60 and C54H48: Silver ion extraction with new concave hydrocarbons. Angew. Chem. Int. Ed. Engl. 34:481484.
- Ikeda, A., and S. Shinkais. 1994. On the origin of high ionophoricity of 1,3-altenate calix[4]arenes
-donor participation in complexation of cations and evidence for metal-tunneling through the calix[4]arene cavity. J. Am. Chem. Soc. 116:31023110.
- Ikeda, A., and S. Shinkais. 1997. Novel cavity design using calix[n]arene skeletons: Toward molecular recognition and metal binding. Chem. Rev. 97:17131734.[Web of Science][Medline]
- Israelachvili, J.N. 1991. Intermolecular and surface forces. 2nd ed. Academic Press, San Diego, CA.
- Kubicki, J.D., C.A. Blake, and S.E. Apitz. 1999. Molecular models of benzene and selected polycyclic aromatic hydrocarbons in the aqueous and adsorbed states. Environ. Toxicol. Chem. 18:16561662.
- Kumpf, R.A., and D.A. Dougherty. 1993. A mechanism for ion selectivity in potassium channels: Computational studies of cation-
interactions. Science (Washington, DC) 261:17081710.[Abstract/Free Full Text]
- Lobo, V.M.M., and L.J. Quaresma. 1989. Handbook of electrolyte solutions. Elsevier Sci. Publ., New York.
- Ma, J.C., and D.A. Dougherty. 1997. The cation-
interaction. Chem. Rev. 97:13031324.[Web of Science][Medline]
- Mantsch, H.H., H. Saitô, and I.C.P. Smith. 1977. Deuterium magnetic resonance, applications in chemistry, physics and biology. Prog. Nucl. Magn. Reson. Spectrosc. 11:211271.
- Means, J.C., S.G. Wood, J.J. Hassett, and W.L. Banwart. 1980. Sorption of polynuclear aromatic hydrocarbons by sediments and soils. Environ. Sci. Technol. 14:15241528.
- Miller, C. 1991. 1990: Annus mirabilis of potassium channels. Science (Washington, DC) 252:10921096.[Free Full Text]
- Munakata, M., L.P. Wu, and G.L. Ning. 2000. A new type of multilayer system-silver(I) complexes of polycyclic aromatic compounds. Coord. Chem. Rev. 198:171203.
- Nanny, M.A., and J.P. Maza. 2001. Noncovalent interactions between monoaromatic compounds and dissolved humic acids: A deuterium NMR T1 relaxation study. Environ. Sci. Technol. 35:379384.[Medline]
- National Research Council. 1994. Alternatives for ground water cleanup. Natl. Academy Press, Washington, DC.
- Nieboer, E., and D.H.S. Richardson. 1980. The replacement of the non-descript term heavy-metals by a biologically and chemically significant classification of metal-ions. Environ. Pollut. Ser. B 1:326.
- Padin, J., and R.T. Yang. 2002. New sorbents for olefin/paraffin separations by adsorption via-complexation: Synthesis and effects of substrates. Chem. Eng. Sci. 55:26072616.
- Schnoor, J.L. 1996. Environmental modeling: Fate and transport of pollutants in water, air, and soil. John Wiley & Sons, New York.
- Schwarzenbach, R.P., P.M. Gschwend, and D.M. Imboden. 2003. Environmental organic chemistry. 2nd ed. Wiley-Interscience, New York.
- Shepodd, T.J., M.A. Petti, and D.A. Dougherty. 1986. Tight, oriented binding of an aliphatic guest by a new class of water-soluble molecules with hydrophobic binding-sites. J. Am. Chem. Soc. 108:60856087.
- Smith, I.P. 1983. NMR of newly accessible nuclei. Chemically and biochemically important elements. Academic Press, New York.
- Smith, M., and A.E. Martell. 1981. Critical stability constants. Vol. 4. Inorganic complexes. Plenum Press, New York.
- Spiro, T.G., and W.M. Stigliani. 1996. Chemistry of the environment. Prentice Hall, Upper Saddle River, NJ.
- Stumm, W., and J.J. Morgan. 1981. Aquatic chemistry. 2nd ed. Wiley Interscience, New York.
- Sunner, J., K. Nishizawa, and P. Kebarle. 1981. Ion-solvent molecule interactions in the gas phase. The potassium ion and benzene. J. Phys. Chem. 85:18141820.
- Takahashi, A., R.T. Yang, C.L. Munson, and D. Chinn. 2001. Cu(I)-Y-zeolite as a superior adsorbent for diene/olefin separation. Langmuir 17:84058413.
- United States Department of Energy. 1992. Chemical contaminants on DOE lands and selection of contaminant mixtures for subsurface science research. Dep. of Energy Rep. DOE/ER-05471. USDOE, Washington, DC.
- Van den Hurk, P., M. Faisal, and M.H. Roberts. 1998a. Interaction of cadmium and benzo[a]pyrene in mummichog (Fundulus heteroclitus): Effects on acute mortality. Mar. Environ. Res. 46:525528.
- Van den Hurk, P., M. Faisal, and M.H. Roberts. 2000. Interactive effects of cadmium and benzo[a]pyrene on metallothionein induction in mummichog (Fundulus heteroclitus). Mar. Environ. Res. 50:8387.[Web of Science][Medline]
- Van den Hurk, P., M.H. Roberts, and M. Faisal. 1998b. Interaction of cadmium and benzo[a]pyrene in mummichog (Fundulus heteroclitus): Biotransformation in isolated hepatocytes. Mar. Environ. Res. 46:529532.
- Zaric, S.D. 1999. Cation-
interaction with transition-metal complex as cation. Chem. Phys. Lett. 311:7780.
- Zhu, D., B.E. Herbert, and M.A. Schlautman. 2003. Molecular-level investigation of monoaromatic compound sorption to suspended soil particles by deuterium NMR. J. Environ. Qual. 32:232239.[Abstract/Free Full Text]
- Zhu, D., B.E. Herbert, M.A. Schlautman, and E.R. Carraway. 2004. Cation-
bonding: A new perspective on the sorption of polycyclic aromatic hydrocarbons to hydrated mineral surfaces. J. Environ. Qual. (in review).
This article has been cited by other articles:

|
 |

|
 |
 
D. Zhu, B. E. Herbert, M. A. Schlautman, E. R. Carraway, and J. Hur
Cation-{pi} Bonding: A New Perspective on the Sorption of Polycyclic Aromatic Hydrocarbons to Mineral Surfaces
J. Environ. Qual.,
July 1, 2004;
33(4):
1322 - 1330.
[Abstract]
[Full Text]
[PDF]
|
 |
|