Sarah J.
Pike
a,
Ennio
Lavagnini
a,
Lisa M.
Varley
b,
Joanne L.
Cook
c and
Christopher A.
Hunter
*a
aDepartment of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK. E-mail: herchelsmith.orgchem@ch.cam.ac.uk
bDepartment of Chemistry, University of Sheffield, Sheffield, S3 7HF, UK
cUnilever R&D Port Sunlight, Quarry Road East, Bebington, Wirral CH63 3JW, UK
First published on 16th May 2019
UV/Vis absorption and NMR spectroscopy titrations have been used to investigate the formation of complexes between cations and neutral H-bond acceptors in organic solvents. Complexes formed by two different H-bond acceptors with fifteen different cations were studied in acetone and in acetonitrile. The effects of water and ion pairing with the counter-anion were found to be negligible in the two polar solvents employed for this study. The data were used to determine self-consistent H-bond donor parameters (α) for a series of organic and inorganic cations; guanidinium, primary, tertiary and quaternary ammonium, imidazolium, methylpyridinium, lithium, sodium, potassium, rubidium and caesium. The results demonstrate the transferability of α parameters for cations between different solvents and different H-bond acceptor partners, allowing reliable prediction of cation recognition properties in different environments. Lithium and protonated nitrogen cations form the most stable complexes, but the α parameter is only 5.0, which is similar to the neutral H-bond donor 3-trifluoromethyl,4-nitrophenol (α = 5.1). Quaternary ammonium is the weakest H-bond donor investigated with an α value of 2.7, which is comparable to an alcohol. The α parameters for alkali metal cations decrease down the group from 5.0 (Li+) to 3.5 (Cs+).
Quantitative scales that describe the H-bond acceptor (HBA) and H-bond donor (HBD) properties of a wide variety of neutral organic functional groups have been developed by Abraham.11 Experimentally determined association constants (K) for simple complexes that form a single H-bond in non-polar solvents are the basis upon which the scales were established.11–14 To develop a universal H-bonding scale, Hunter extended this method to address the influence of solvent on solution phase equilibria between H-bonded solutes. The solvent competition model illustrated in Fig. 1 treats solution phase H-bonding interactions as an equilibrium between pairwise contacts between solvent and solute.15
The H-bond donor and acceptor parameters for the solute (α and β) and the solvent (αs and βs) can be used in eqn (1) to predict the Gibbs free energy change (ΔG°) for formation of a H-bonded complex in any solvent.16
ΔG°/kJ mol−1 = −(α − αs)(β − βs) + 6 | (1) |
Experimentally measured association constants for H-bond formation (K) can be used in conjunction with eqn (1) to determine H-bond parameters for solutes or solvents.15–20 For example, the α value of a HBD can be determined through the rearrangement of eqn (1) to give eqn (2), when the β, αs and βs parameters are known.
α = αs + (RTlnK + 6)/(β − βs) | (2) |
This approach has been employed to quantify the H-bond properties of neutral organic functional groups,21–23 and we have recently demonstrated that charged species can be placed on the same H-bonding scale through the determination of β values for a range of anions.19 Here, the same methodology is used to place cations on the H-bond scale through the experimental determination of α values.
To date, there have been few systematic studies of non-covalent interactions with cations. Marcus investigated the solvation properties of cations using linear solvation energy relationships.24 The thermodynamics of phase transfer of ions from water to organic solvents was used to determine empirical parameters to describe the solvation properties of ions.25 Gilkerson used conductance experiments and UV/Vis absorption titrations to study interactions between cations and HBAs in organic solvents.26–30 The relative HBD strength of the cations was found to decrease in the order protonated amine > alkali metal cation > quaternary ammonium.28a Conductance experiments were also used to measure association constants (K) for 1:1 complexes formed between triethylammonium picrate and a series of pyridines in nitrobenzene. A linear correlation was found between the value of logK and the corresponding Hammett parameter for the substituent on the H-bond acceptor.31,32
In this paper, we report experiments that establish α H-bond donor parameters for a number of inorganic and organic monovalent cations by using the results of titration experiments conducted in two different solvents with two different HBAs in conjunction with eqn (2). The metal ions used in these experiments are not H-bond donors, but as we will show, it is possible to provide a quantitative description of the stabilities of the complexes formed with H-bond acceptors using the same α parameter scale that is used to describe the non-covalent interaction properties of H-bond donors.
Scheme 1 (a) H-Bond acceptors and (b) cations (ArF = 3,5-bis(trifluoromethyl)phenyl, R = 2-ethylhexyl, 12 is the tris(1,2-dimethoxyethane) adduct). |
Twenty different salts were used in this study providing variation in both the counter-anion and the cation (Scheme 1b). Weakly coordinating anions, tetraphenylborate, hexafluorophosphate and tetrafluoroborate, were used as counterions to minimize competing interactions with the formation of H-bonded complexes. The influence of the counterion on the stabilities of the complexes was investigated by studying Na+ salts with six different anions. The influence of the water content of the solvent on the association constant was also studied for a subset of the cation complexes. Salts 8–11 were prepared from the corresponding hydrochloride salt using an anion exchange reaction with NaBPh4 in water.26,35,36 The other salts (12–26) are commercially available.
A number of additional titration experiments were carried out to extend the range of cations, but the results proved unsuitable for accurate determination of α parameters. For example, association constants were too high to be measured for the iodide salts of Ca2+ and Mg2+ and the tetrafluoroborate salts of Pd2+, Fe2+ and Cu+. The tetrafluoroborate salt of methyl(diphenyl)sulfonium caused methylation of 1 and did not bind strongly enough to 2 to measure association constants. The hexafluorophosphate salt of the trityl cation caused disappearance of the UV/Vis absorption of 1 and gave an anomalously large change in the 31P chemical shift of 2, suggesting that the interactions are covalent in nature.
Two different solvents were used for the titration experiments: acetonitrile and acetone. Acetonitrile is a stronger H-bond donor but significantly weaker H-bond acceptor than acetone,15 and the two solvents have different dielectric constants (37.5 for acetonitrile, and 20.7 for acetone). The original neutral H-bond scales were established based on association constants measured in carbon tetrachloride and 1,1,1-trichloroethane,11–13 so in order to establish transferability between solvents the complexes formed by 1 and 2 with a series of neutral HBDs (3–7 in Scheme 1b) were studied in four solvents, carbon tetrachloride, chloroform, acetonitrile and acetone.21,22 The neutral HBDs (3–7) all have a UV/Vis absorption band that is sensitive to H-bond formation, so for the titrations with 2, these compounds were used as hosts in UV/Vis absorption titration experiments.
HBD | HBA | MeCN | Acetone | CHCl3 | CCl4 |
---|---|---|---|---|---|
a Average of at least two titrations. Errors are quoted at the 95% confidence limit. b Association constant too high to be measured using UV/Vis spectroscopy. c Poor solubility in carbon tetrachloride prevented the acquisition of data. d Average of K values obtained from both 31P NMR and UV/Vis spectroscopy data spectroscopy titrations. e | |||||
3 | 1 | 260 ± 63 | 340 ± 6 | 260 ± 55 | —c |
4 | 1 | 200 ± 50 | 860 ± 350 | 1200 ± 500 | —c |
5 | 1 | 3100 ± 650 | 2000 ± 500 | 4000 ± 700 | —c |
3 | 2 | 20 ± 1 | 26 ± 1 | 47 ± 24d | 1860 ± 226d |
4 | 2 | 6 ± 1 | 11 ± 1 | 18 ± 6e | 1280 ± 226d |
5 | 2 | 63 ± 18d | 64 ± 45d | 217 ± 66d | 11775 ± 7280d |
6 | 2 | 106 ± 11d | 134 ± 18 | 1410 ± 23d | 120300 ± 61376d |
7 | 2 | 283 ± 48d | 259 ± 80 | 5400 ± 42 | —b |
Some of the relevant H-bond parameters are available in the literature,15–21 and the other H-bond parameters were optimized to obtain the best fit of the experimental association constants to eqn (1). The resulting H-bond parameters are shown in Table 2. Fig. 2 shows that eqn (1) provides an excellent description of the association constants shown in Table 1, if the H-bond parameters in Table 2 are used.
Fig. 2 Comparison of experimental free energies of complexation with values calculated using eqn (1) for complexes formed by 1 and 2 with neutral HBDs in carbon tetrachloride, chloroform, acetonitrile and acetone. The line represents . |
Anion | Cation | HBD | HBA/solvent | |||
---|---|---|---|---|---|---|
1 | 2 | |||||
MeCN | Acetone | MeCN | Acetone | |||
a Average of at least two titrations. Errors are quoted at the 95% confidence limit. Unless otherwise stated greater than 50% saturation of the binding isotherm was achieved. b Association constant was too high to be measured using UV/Vis spectroscopy. c Protonation of the H-bond acceptor was observed upon addition of guest. d Li salt is LiBPh4 tris(1,2-dimethyoxyethane). e Association constant was too low to be measured, because the salt was not sufficiently soluble to obtain 50% saturation of the binding isotherm. f Poor fit to a 1:1 binding isotherm. | ||||||
BPh4− | Guanidinium | 8 | 33800 ± 2000 | —b | 198 ± 61 | 225 ± 40 |
BPh4− | 2-Ethylhexyl ammonium | 9 | —c | —c | 147 ± 25 | 123 ± 34 |
BPh4− | Triethyl ammonium | 10 | —c | —c | 85 ± 6 | 203 ± 50 |
BPh4− | N-Methyl imidazolium | 11 | —c | —c | 102 ± 12 | 156 ± 18 |
BPh4− | Li+ | 12 | 27400 ± 5900 | —b | —b | — |
BPh4− | Na+ | 13 | 390 ± 51 | 1700 ± 160 | —e | 73 ± 5 |
BPh4− | K+ | 14 | —e | 270 ± 75 | —e | 13 ± 3 |
BPh4− | Rb+ | 15 | —e | 220 ± 50 | —e | —e |
BPh4− | Cs+ | 16 | —e | 210 ± 35 | —e | —e |
BF4− | Ag+ | 17 | —f | —b | 8 ± 1 | 966 ± 40 |
PF6− | +I(4-tBuPh)2 | 18 | 38300 ± 1100 | 48000 ± 3000 | 62 ± 3 | 83 ± 12 |
PF6− | MOIM+ | 19 | 23 ± 9 | 200 ± 60 | —e | —e |
PF6− | N-Butyl-4-methyl pyridinium | 20 | 12 ± 2 | 180 ± 40 | —e | —e |
PF6− | Tetra(n-butyl) ammonium | 21 | —e | 10 ± 2 | —e | —e |
PF6− | Na+ | 22 | 391 ± 35 | 1700 ± 140 | —e | 55 ± 13 |
BF4− | Na+ | 23 | 300 ± 90 | 1900 ± 100 | 57 ± 6 | 57 ± 2 |
BArF4− | Na+ | 24 | 400 ± 160 | 1400 ± 400 | 53 ± 3 | 72 ± 2 |
N(SO2CF3)2− | Na+ | 25 | 410 ± 30 | 1500 ± 200 | 56 ± 5 | 78 ± 12 |
I− | Na+ | 26 | 320 ± 120 | 1000 ± 300 | 45 ± 1 | 59 ± 6 |
The association constants measured for the complexes span four orders of magnitude (Table 3). The relative polarities of the solvents and solutes determine the stabilities (Table 2). The complexes formed with 1 are generally more stable than the corresponding complexes formed with 2, which is in agreement with the H-bond parameters of the two acceptors (β = 14.0 and 10.7 respectively). The differences in association constant are most significant in acetone, where the complexes formed with 1 are up to three orders of magnitude more stable than corresponding complexes formed with 2. For the weaker HBA 2, the association constants for the complexes formed with 14–16, 19 and 21–23 were too low to be measured in acetonitrile. For stronger HBA 1, the association constants for the complexes formed with 6–8, 12 and 17 were too high to be measured in acetone.
The presence of traces of water might perturb the measured values of the association constants reported in Table 3. The water content of acetonitrile and acetone used in these experiments was determined to be 0.02% in both cases.38 To establish whether these quantities of water can influence the stabilities of the complexes and hence perturb the measured association constants, water was deliberately added to the stock solutions of both solvents, and the titration experiments were repeated for a subset of the complexes. The influence of added water on the association constants measured in acetonitrile and in acetone is illustrated in Fig. 4. Large quantities of water are required to have a significant effect in either solvent, so we conclude that the traces of water present in the titration experiments do not affect the values of the association constants in Table 3.
Anion | Cation | HBD | HBD/solvent | |||
---|---|---|---|---|---|---|
1 | 2 | |||||
MeCN | Acetone | MeCN | Acetone | |||
a Errors at the 95% confidence limit. b No experimental data available. | ||||||
BPh4− | Guanidinium+ | 8 | 5.0 ± 0.1 | —b | 4.9 ± 0.1 | 5.1 ± 0.1 |
BPh4− | 2-Ethylhexyl ammonium | 9 | —b | —b | 4.8 ± 0.1 | 4.8 ± 0.2 |
BPh4− | Triethyl ammonium | 10 | —b | —b | 4.5 ± 0.1 | 5.0 ± 0.1 |
BPh4− | N-Methyl imidazolium+ | 11 | —b | —b | 4.6 ± 0.1 | 4.9 ± 0.1 |
BPh4− | Li+ | 12 | 5.0 ± 0.2 | —b | —b | —b |
BPh4− | Na+ | 13 | 3.8 ± 0.1 | 4.1 ± 0.1 | —b | 4.5 ± 0.1 |
BPh4− | K+ | 14 | —b | 3.6 ± 0.1 | —b | 3.7 ± 0.1 |
BPh4− | Rb+ | 15 | —b | 3.5 ± 0.1 | —b | —b |
BPh4− | Cs+ | 16 | —b | 3.5 ± 0.1 | —b | —b |
BF4− | Ag+ | 17 | —b | —b | 3.5 ± 0.1 | 5.8 ± 0.1 |
PF6− | +I(4-tBuPh)2 | 18 | 5.1 ± 0.1 | 5.1 ± 0.1 | 4.4 ± 0.1 | 4.6 ± 0.1 |
PF6− | MOIM+ | 19 | 3.0 ± 0.4 | 3.5 ± 0.1 | —b | —b |
PF6− | N-Butyl-4-methyl pyridinium | 20 | 2.9 ± 0.1 | 3.5 ± 0.1 | —b | —b |
PF6− | Tetra(n-butyl) ammonium | 21 | —b | 2.6 ± 0.1 | —b | —b |
PF6− | Na+ | 22 | 3.8 ± 0.1 | 4.1 ± 0.1 | —b | 4.4 ± 0.2 |
BF4− | Na+ | 23 | 3.8 ± 0.1 | 4.2 ± 0.1 | 4.4 ± 0.1 | 4.4 ± 0.1 |
BArF− | Na+ | 24 | 3.8 ± 0.2 | 4.1 ± 0.3 | 4.3 ± 0.1 | 4.5 ± 0.1 |
N(SO2CF3)2− | Na+ | 25 | 3.8 ± 0.1 | 4.1 ± 0.2 | 4.4 ± 0.1 | 4.6 ± 0.2 |
I− | Na+ | 26 | 3.8 ± 0.2 | 4.0 ± 0.2 | 4.3 ± 0.1 | 4.4 ± 0.1 |
The type of ion pair and/or aggregate formed by salts in solution is highly dependent on both the nature of the counterion and the polarity of the solvent.39–41 The dielectric constant of a solvent is known to be inversely proportional to the association constant for the ion pairing of the salt42 and in solvents with high dielectric constants the presence of loose ion pairs (solvent-shared ion pairs wherein oppositely charged species are separated by one layer of solvent molecules and solvent-separated ion pairs which have more than one layer of solvent separating the oppositely charged species) are known to exist.40 Contact ion pairs which involve direct interaction between oppositely charged species and higher aggregates of ions are generally found in solvents of low dielectric constant and at high salt concentrations.41 Acetonitrile and acetone have high dielectric constants of 37.5 and 20.7 respectively.42 In both solvents loose ion pairs are mostly likely to dominate, either as solvent-separated or solvent-shared, whilst the presence of contact ion-pairs (and aggregates) would not be expected to be significant.40
To establish whether ion pairing (or aggregation) of the salt has a significant effect on the H-bond parameters reported in Table 4, six different counter-anions were used for the sodium cation; tetraphenylborate (13), hexafluorophosphate (23), tetrafluoroborate (24), tetrakis(3,5-bis(trifluoromethyl)phenyl)borate (25), trifluoromethanesulfonimide (26) and iodide (27) (Scheme 1). While most of the anions are weakly coordinating, iodide is a relatively strongly coordinating counterion. The HBA parameters (β) of three of the counter-anions have previously been determined:19a hexafluorophosphate (7.0) trifluoromethanesulfonimide (7.3) and iodide (8.9).19a Accordingly, the influence of the strength of the H-bond acceptor parameter of the counterion on the interaction of Na+ with HBAs 1 and 2 was investigated. Although there are some variations in the value of α measured for Na+ with different HBAs and solvents, the values for different counterions are practically identical in all cases. This result suggests that the values reported in Table 4 are not significantly perturbed by interactions with the counter-anion.
Table 4 shows that the analysis described here results in a range of different α parameters for the same cation. For example, the 22 independent measurements of the value α for Na+ fall between 3.8 and 4.6. This variation between individual measurements of α for the same cation reflects the magnitude of the errors associated with the approach. However, it is possible to describe the all of the experimental data by using a single representative value of α for each cation. Fig. 5 shows the result of optimising a single α parameter for each cation to fit the experimental association constants in Table 3. The calculated free energies of complexation in Fig. 5 agree well with the experimental data. The rmsd between the experimental and calculated values is 1.5 kJ mol−1, which provides an indicator of the error associated with the spread of values of the individual measurements of α for each cation in Table 4. The optimised cation α values are reported in Table 5. These α parameters can be used to describe the interactions of cations with different HBAs in different solvent environments and with different counterions.
Fig. 5 Comparison of experimental free energies of complexation with values calculated using eqn (1) for 1:1 complexes formed with cations. The line represents . |
Cation | α |
---|---|
Guanidinium | 5.0 |
2-Ethylhexyl ammonium | 4.8 |
Triethyl ammonium | 4.8 |
N-Methyl imidazolium | 4.7 |
Li+ | 5.0 |
Na+ | 4.1 |
K+ | 3.6 |
Rb+ | 3.5 |
Cs+ | 3.5 |
+I(4-tBuPh)2 | 4.9 |
MOIM+ | 3.3 |
N-Butyl-4-methyl pyridinium | 3.5 |
Tetra(n-butyl)ammonium | 2.7 |
The α parameters measured for the cations are illustrated graphically in Fig. 6. Lithium and guanidinium form the most stable complexes, but the value of α is only 5.0, which is similar to the most polar neutral HBDs, e.g. 4-nitro,3-trifluoromethyl-phenol (α = 5.1). The HBD properties of protonated nitrogen donors are all very similar and are comparable to hexafluoropropan-2-ol (α = 4.7–4.8).15 The α parameters for group 1 metal cations decrease in strength upon progression down the group. Sodium has properties similar to phosphoric acid (α = 4.1), and potassium, rubidium and caesium are similar to a carboxylic acid (α = 3.6).15 For the organic cations, changing the nature of donor from the NH in protonated amines to the CH in quaternary ammonium cations reduces α from 4.7–5.0 to 2.7–3.5. Tetrabutyl ammonium is the weakest HBD (α = 2.7), but the HBD properties are still significant, comparable to those of an alcohol (α = 2.7).15
Measurement of the α parameters for cations does not extend the α scale beyond the values measured for neutral functional groups. The highest α values are found for lithium and protonated nitrogens, with values similar to the strongest neutral H-bond donors (α ≈ 5).30 The variation in α values measured for alkali metal cations shows that upon progression down the group, the interaction properties of the cation are reduced. In organic cations, the nature of the donor heavy atom is important, so protonated amines have significantly higher HBD parameters than the CH donor in quaternary ammonium cations. These new H-donor parameters for cations will be useful in the development of our understanding and prediction of the behaviour of charged species in organic solution with applications in the design of new supramolecular systems.
BArF | Tetrakis(3,5-bis(trifluoromethyl)phenyl)borate |
br | Broad |
Bu | Butyl |
Bz | Benzoyl |
d | Doublet |
DMF | Dimethylformamide |
DMSO | Dimethylsulfoxide |
Et | Ethyl |
HBA | H-bond acceptor |
HBD | H-bond donor |
HMPA | Hexamethylphosphoramide |
LFER | Linear free energy relationship |
m | Multiplet |
Me | Methyl |
m.p. | Melting point |
MOIM | 1-Methyl,-3-octyl,imidazolium |
Ph | Phenyl |
Pip | Piperidine |
q | Quartet |
R | 2-Ethylhexyl |
s | Singlet |
t | Triplet |
TBA | Tetrabutylammonium |
THF | Tetrahydrofuran |
Footnote |
† Electronic supplementary information (ESI) available: Titration data and details of fitting the binding isotherm. See DOI: 10.1039/c9sc00721k |
This journal is © The Royal Society of Chemistry 2019 |