Journal of Chemistry

Volume 2013 (2013), Article ID 318015, 9 pages

http://dx.doi.org/10.1155/2013/318015

Research Article

Conductometric Studies of Thermodynamics of 1,10-Didecyl-1,10-diaza-18-crown-6 Complexes with , , , , , , and Ions in Acetonitrile, Methanol, and Ethanol Solutions

Mehdi Taghdiri, Mahmood Payehghadr, Reza Behjatmanesh-Ardakani, Sirous Movahhedi, and Farzaneh Nourifard

Department of Chemistry, Payame Noor University, Tehran 19395-3697, Iran

Received 28 June 2012; Accepted 1 October 2012

Academic Editor: Pedro Avila-Pérez

Copyright © 2013 Mehdi Taghdiri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The complexation reactions between 1,10-didecyl-1,10-diaza-18-crown-6 (DD18C6) and , , , , , and ions were studied conductometrically in acetonitrile, methanol, and ethanol at various temperatures. The formation constants of the resulting 1 : 1 complexes were calculated from the computer fitting of the molar conductance-mole ratio data at different temperatures. In the solvents studied, the stability of the resulting complexes varied in the order acetonitrile ≫ methanol > ethanol. The enthalpy and entropy changes of the complexation reactions were evaluated from the temperature dependence of formation constants. The versus plot of thermodynamic data obtained shows a linear correlation indicating the existence of enthalpy-entropy compensation in the complexation reactions. The ab initio studies calculated at B3LYP/6-31G level of theory indicate that the binding energy of complexes decreases with increasing cation size in the gas phase. The experimental data obtained in the solvents show different trends so that the stability constant of Li⁺ and Mg²⁺ complexes is lower than other alkali and alkaline earth ions, respectively due to the solvation of these small ions in solution.

1. Introduction

Since the first discovery of crown ethers by Pedersen [1], the studies of these ligands and their metal ion complexes have become a very popular field of research [2–4], mainly due to their similarities in many respects to the naturally occurring ionophores [5]. Meanwhile, these ligands have demonstrated a high potential for use in many chemical and industrial processes, where cation selectivity and/or solubility in nonpolar solvents are required. Thus, due to the importance of selectivity and stability of the crown ether complexes with different cations in a wide variety of practical applications [6, 7], extensive amount of work in the thermodynamic aspects of the corresponding complexation reactions has also been reported in the literature [8–10].

The stability constants for the 1 : 1 complexes of the alkali and alkaline earth metal ions with 1,10-diaza-18-crown-6 (D18C6) have extensively been determined in various pure solvents [8, 10]. It was of interest to us to study the influence of the cation size (and nature), solvent properties, and bulky substituent in the ligand on the interaction of metal ions with crown ethers.

In this work, we report a conductometric study of the stoichiometry, stability, and thermodynamic parameters of the 1,10-didecyl-1,10-diaza-18-crown-6 (DD18C6) complexes with Li⁺, Na⁺, K⁺, Cs⁺, , Mg²⁺, and Ca²⁺ ions in acetonitrile (AN), methanol (MeOH), and ethanol (EtOH) solvents. It should be noted that, while AN, MeOH, and EtOH have about the close dielectric constants (35.94, 32.66, and 24.55, resp.), they possess quite different donor numbers (DN; i.e., DN = 14.1 for AN, DN = 19.0 for MeOH, and DN = 31.5 for EtOH) [11].

2. Experimental

2.1. Reagents

Reagent-grade perchlorate salts of lithium, sodium, potassium, ammonium, and magnesium, nitrate salt of cesium, dried AN (H₂O < 0.005%), MeOH, EtOH, and DD18C6 (all from Merck) were of the highest purity available and used as received. Ca(ClO₄)₂ was prepared by treating CaCO₃ with small excess of 3 M perchloric acid, followed by evaporation to dryness, recrystallization three times from deionized water, and drying at 120°C. Although we did not encounter any problem, it should be noted that dry perchlorate salts and perchlorate salts in organic solvents are potentially explosive and should be handled only in small quantities with appropriate precautions.

2.2. Conductometric Titration

Conductance measurements were carried out with a Metrohm 712 conductivity meter. A dip-type conductivity cell made of platinum black was used. The cell constant at the different temperatures used was determined by conductivity measurements of a 0.010 M solution of analytical-grade KCl (Merck) in triply distilled deionized water. The specific conductance of this solution at various temperatures has been reported in the literature [12]. In all measurements, the cell was thermostated at the desired temperature ±0.1°C using a Haake D1 thermostat-circulator water bath.

In a typical experiment, 10 mL of the desired metal ion (5.0 × 10−5 M) was placed in the titration cell, thermostated to the desired temperature, and the conductance of solution was measured. Then, a known amount of a concentrated crown ether solution was added in a stepwise manner using a calibrated micropipette. The conductance of the solution was measured after each addition. The ligand solution was continually added until the desired ligand to cation mole ratio was achieved.

2.3. Calculation of Formation Constant

The formation constants, , and the limiting molar conductances, , of the resulting 1 : 1 complexes between DD18C6 and the univalent and divalent cations used, in different solvents and at various temperatures, were calculated by fitting the observed molar conductance, , at varying [DD18C6]/[] mole ratios to a previously derived equation [13–15] which express the as a function of the free and complexed metal ions. A nonlinear least squares curve fitting using Microsoft Excel Solver (version 11.0) was applied for the evaluation of formation constant and limiting molar conductance of the resulting 1 : 1 complexes.

Chi-square statistic was used to evaluate the fitness of equation to the experimental data. The Chi-square test statistic is the sum of the squares of the differences between the experimental data and data obtained by calculating from (1). The equivalent mathematical statement is:

where and are the molar conductance obtained by experimental data and calculating from (2), respectively

where is the stability constant of complex, and are the analytical concentrations of DD18C6 and M⁺ ion, respectively, and and are the respective molar conductance of the free and complexed M⁺ ion. If calculated data are similar to the experimental data, will be a small number ; if they are different, will be a large number .

2.4. Computational Method

GAMESS R3 2010 Linux code [16] was used for all Density Functional Theory (DFT) calculations. All calculated results belong to the B3LYP/6-31G level of theory. The compilation of the code adopted with Intel’s ifort compiler of fortran in conjunction with MKL mathematical library (all are freely accessible from Intel site with non-commercial license). CentOS linux distribution was used as operating system on Intel SMP computer with 16 processors and 32 GB of RAM. All input z-matrixes were prepared with HyperChem program [17], and then Avogadro [18] was used to convert the “pdb” format to “inp” file. The output “log” files were also converted to the “pdb” format to able to visualize the optimized structures.

3. Results and Discussion

In order to evaluate the influence of adding DD18C6 on the molar conductance of Li⁺, Na⁺, K⁺, Cs⁺, , Mg²⁺, and Ca²⁺ ions in AN, MeOH, and EtOH, the molar conductance at a constant salt concentration (5.0 × 10−5 M) was monitored while increasing the crown ether concentration at various temperatures. Some of the resulting molar conductances versus DD18C6/cation mole ratio plots are shown in Figures 1 and 2.

(a)

(b)

(c)

Figure 1: Molar conductance (S cm² mol⁻¹) versus [DD18C6]/[] for various DD18C6- systems at 25°C in AN (a), MeOH (b), and EtOH (c). Solid lines represent the fitted curves.

(a)

(b)

(c)

(d)

Figure 2: Molar conductance (S cm² mol⁻¹) versus [DD18C6]/[] plots in AN at various temperatures. The cations are (a) Na⁺, (b) K⁺, (c) Cs⁺, and (d) Mg²⁺. Solid lines represent the fitted curves.

As it is seen, in AN, there is a gradual decrease in the molar conductance with an increase in the crown ether concentration, which begins to level off at a mole ratio greater than one (except Cs⁺). Such a conductance behavior is indicative of the formation of fairly stable 1 : 1 complexes in AN.

The formation constants of all DD18C6- complexes in different solvents at various temperatures, obtained by computer fitting of the molar conductance-mole ratio data and Chi square numbers (), are listed in Table 1. The computer fits of the molar ratio data are shown as solid lines in Figures 1 and 2. The small numbers of Chi square indicate that our assumption of 1 : 1 stoichiometry seems reasonable in the light of the fair agreement between the observed and calculated molar conductances. It should be noted that, in the procedure of calculation of formation constants, the association between cation and counter ions was considered negligible, under the highly dilute experimental conditions used (5.0 × 10−5 M). Furthermore, since the concentration of crown ether was kept below 1.0 × 10−3 M during the conductometric titrations, the corrections for viscosity changes were also neglected.

Table 1: Formation constants for different M-DD18C6 complexes in acetonitrile, methanol, and ethanol solvents.

Cationa
Solvent

15°C

25°C

35°C

45°C

AN
6.70
0.277
6.49
0.383
6.18
0.117
6.01
0.004

Li⁺ (0.74)
MeOH
3.04
0.006
3.69
0.009
3.18
0.009
2.83
0.167

EtOH
1.86
0.141
1.67
0.003
1.54
0.126
1.27
0.027

AN
6.87
0.153
6.50
0.104
6.23
0.072
6.00
0.075

Na⁺ (1.02)
MeOH
3.85
0.006
3.75
0.005
3.58
0.018
3.54
0.021

EtOH
2.22
0.049
2.10
0.025
1.71
0.009
1.45
0.003

AN
7.00
0.012
6.50
0.014
6.00
0.011
5.62
0.012

K⁺ (1.38)
MeOH
3.63
0.011
3.07
0.007
2.40
0.010
2.28
0.010

EtOH
2.26
0.012
2.09
0.006
1.72
0.003
1.60
0.004

AN
6.48
0.016
6.20
0.135
6.07
0.154
6.00
0.004

(1.48)
MeOH
3.43
0.005
2.91
0.054
2.63
0.068
1.83
0.203

EtOH
2.17
0.080
2.01
0.062
1.65
0.061
1.32
0.093

AN
2.95
0.008
2.61
0.033
2.31
0.030
1.62
0.007

Cs⁺ (1.70)
MeOH
3.73
0.031
3.34
0.056
3.15
0.024
2.34
0.016

EtOH
2.32
0.023
2.13
0.010
1.85
0.010
1.71
0.006

AN
5.97
0.040
6.04
0.031
6.04
0.031
5.70
0.019

Mg²⁺ (0.72)
MeOH
4.02
0.021
4.05
0.047
4.09
0.046
4.17
0.026

EtOH
1.98
0.017
1.87
0.010
1.77
0.003
1.48
0.013

AN
5.85
0.125
5.80
0.097
5.75
0.093
5.69
0.089

Ca²⁺ (1.00)
MeOH
4.27
0.059
—b
1.02
—b
0.86
—b
0.64

EtOH
3.09
0.125
2.65
0.045
2.16
0.105
1.83
0.220

aThe values in parenthesis are the ionic sizes in [15]

bInsufficient fit.

In order to have a better understanding of the thermodynamics of complexation reactions of Li⁺, Na⁺, K⁺, Cs⁺, , Mg²⁺, and Ca²⁺ ions with DD18C6, it is useful to investigate the enthalpic and entropic contributions to these reactions. The and of the complexation reactions in different solvents were evaluated from the temperature dependence of the formation constants by applying a linear least-squares analysis according to the Van’t Hoff equation. The Van’t Hoff plots of log  versus are shown in Figure 3. The enthalpies and entropies of complexation were determined in the usual manner from the slopes and intercepts of the plots, respectively, and the results are listed in Table 2.

Table 2: Thermodynamic parameters for different -DD18C6 complexes in acetonitrile, methanol, and ethanol solvents.

Cation
Solvent
° (kJ mol⁻¹)
° (J mol⁻¹ K⁻¹)
° (kJ mol⁻¹)

AN

Li⁺
MeOHa

EtOH

AN

Na⁺
MeOHb

EtOH

AN

K⁺
MeOHc

EtOH

AN

MeOHb

EtOH

ANc

Cs+
MeOHb

EtOH

ANd
—d
—d
—d

Mg²⁺
MeOH

EtOHc

Ca²⁺
AN

EtOH

aWithout considering at 15°C

bWithout considering at 35°C

cWithout considering at 45°C

dHigh uncertainty.

(a)

(b)

(c)

Figure 3: Van’t Hoff plots for the 1 : 1 complexation of various DD18C6- systems in AN (a), MeOH (b), and EtOH (c).

It is interesting to note that the formation constants obtained in this study at 25°C in AN and MeOH are in satisfactory arrangement with the previously reported values for D18C6 (Table 3).

Table 3: Comparison of stability constants for complexes of 1,10-diaza-18-crown-6 (D18C6) and DD18C6 with Li⁺, Na⁺, K⁺, Cs⁺, Mg²⁺, and Ca²⁺ ions in AN and MeOH solvents at 25°C.

Ligand
Cation
Solvent
Method

Reference

D18C6
Li⁺
AN
NMR
4.39
[19]

DD18C6

Conductometry
6.49
This work

D18C6
Na⁺
AN
NMR
4.49
[20]

Conductometry
4.30
[21]

DD18C6

Conductometry
6.50
This work

D18C6

MeOH
Potentiometry
1.50
[22]

DD18C6

Conductometry
3.75
This work

D18C6
K⁺
AN
NMR
4.35
[20]

Conductometry
4.32
[21]

DD18C6

Conductometry
6.50
This work

D18C6

MeOH
Potentiometry
1.80
[22]

DD18C6

Conductometry
3.07
This work

D18C6
Cs⁺
AN
NMR
2.25
[20]

Conductometry
2.48
[21]

DD18C6

Conductometry
2.61
This work

D18C6
Mg²⁺
MeOH
Potentiometry
3.81
[10]

Spectrophotometry
3.40
[23]

DD18C6

Conductometry
4.05
This work

D18C6
Ca²⁺
MeOH
Calorimetry
3.87
[24]

DD18C6

Conductometry
4.27a
This work

aAt 15°C.

In the case of complexation of macrocyclic ligands, there are at least three factors which can make significant contributions to the stability of their metal ion complexes: (i) the cation size, (ii) the ionic solvation of the charged species involved, and (iii) conformations of the free and complexed crown ethers. The data given in Table 1 clearly illustrate the fundamental role of the solvent properties in the -DD18C6 complexation reactions studied. In the case of all cations except Cs⁺, the stability of the resulting complexes with DD18C6 varied in the order AN ≫ MeOH > EtOH. It is well known that the solvating ability of the solvent, as expressed by the Gutmann donor number [11], plays an important role in different complexation reactions [9, 10]. There is actually an inverse relationship between the stabilities of the complexes and the solvating abilities of the solvents. AN has a lower donicity () than MeOH () and EtOH (), therefore, shows less competition with the crown ether for mentioned ions; thus, it is not unexpected to observe the highest stability of the complexes in AN and the lowest in EtOH. In the case of DD18C6-Cs⁺, the interaction of the large size central metal ion with solvent is stronger for MeOH than for AN due to higher donicity; therefore, the DD18C6-Cs⁺ complex will be stabilized by MeOH solvation and formation constant will be larger.

Table 2 shows that, as expected, for -DD18C6 systems studied, the thermodynamic data vary significantly with the solvent properties. However, the observed increase (or decrease, depending on the nature of the metal ion) in value will be compensated by an increase (or decrease) in the corresponding value. The existence of such a compensating effect (Figure 4) between and values has been frequently reported for a variety of metal-ligand systems [25, 26]. This effect causes the overall change in the value of the complex to be smaller than might be expected from the change in either or independently.

Figure 4: Plot of (kJ mol⁻¹) versus (kJ mol⁻¹) for 1 : 1 complexation of Li⁺, Na⁺, K⁺, Cs⁺, Mg²⁺, Ca²⁺, and ions in AN, MeOH, and EtOH solvents ().

The formation constants for -DD18C6 and for -D18C6 have been compared in Table 3. The increased formation constant of -DD18C6 as compared to that of -D18C6 must be ascribed to replacement of the solvation layer(s) around cation by the bulky DD18C6 ligand and also to electron donating of the decyl groups, which increase the electron-donor ability of the nitrogen atoms of the ring.

To obtain more information about the structures of complexes, and to investigate that whether ions form complexes in the gas phase or not, the quantum chemical calculations were conducted. The structures were full optimized with no symmetry or any constraints about bond lengths, bond angles, or bond torsions. Five conformations were tested for DD18C6. The conformation with the lowest gas phase energy was chosen for further studies. Figure 5 shows the optimized geometries of both free and complexes of DD18C6 with Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ ions. The figure shows that for K⁺-DD18C6 complex, all donor atoms are in one direction, and the ion is placed out of the plane of the cavity, while for Li⁺- and Na⁺-DD18C6, the cavity is twisted, and the ions are placed inside the cavity. The structure of K⁺-DD18C6 is similar to a boat structure, while for Li+ and Na⁺ the structures are similar to a chair. K⁺ interact with 5 donor atoms in its boat structure, while Li+ and Na⁺ interact with four oxygen atoms.

Figure 5: The optimized geometries of free DD18C6 (a) and its complexes with Li⁺ (b), Na⁺ (c), K⁺ (d), Mg²⁺ (e), and Ca²⁺ (f). Hydrogen atoms have not been shown for clarity.

Table 4 shows the electronic and binding energies of the species. The data show that there is a strong interaction between each of the studied ions with DD18C6 in the gas phase. From the table, it can be concluded that in the gas phase the binding energy of K⁺ with DD18C6 is less than the binding energies of Li+ and Na⁺ with DD18C6. It seems that this less stability (even with 5 interactions) could be related to its boat structure compared to the two other chair structures. Moreover, the surface charge density (charge number/surface area) is smaller for K⁺ than Li+ and Na⁺ and hence the interaction of K⁺ with crown is weaker. Based on the ionic radius of Li⁺, Na⁺, and K⁺ (Table 1) and the cavity of 18C6 with a radius of 1.3–1.6 Å [15], it is clear that the potassium ion is placed out of the cavity plane. In comparison with the above alkaline ions, Mg²⁺ and Ca2+ ions have higher charge density; therefore, it can be concluded that their binding energy in the gas phase should be larger (see the data in Table 4). With enlarging the ionic radius from Mg²⁺ to Ca²⁺, the crown cavity in the complex form tend to be more planar. Mg²⁺-DD18C6 has more binding energy than Ca²⁺-DD18C6.

Table 4: Calculated electronic energies of different species and binding energies of complexes formed by interaction between DD18C6 and Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺. All calculations have been done with B3LYP/6-31G.

Compounds
Electronic energy (a.u.)
Binding energy (a.u.)
Binding energy

(kcal/mol)a

DD18C6
−1668.315091
—
—

Li⁺
−7.277405391
—
—

Na⁺
−162.0485932
—
—

K⁺
−599.6912799
—
—

Mg²⁺
−199.2015568
—
—

Ca²⁺
−676.8353899
—
—

DD18C6-Li⁺
−1675.786365
−0.193868538
−121.6543491

DD18C6-Na⁺
−1830.512579
−0.148894724
−93.43285381

DD18C6-K⁺
−2268.111697
−0.105325935
−66.093025

DD18C6-Mg²⁺
−1868.047261
−0.530613766
−332.965179

DD18C6-Ca²⁺
−2345.535887
−0.385406425
−241.8461934

aBinding energy () is defined as: = (complex) (ligand) (ion). 1 a.u. is equal to 627.5095 kcal/mol.

As it is seen in Table 4, the binding energy decreases monotonically with increasing the size of metal cation. The experimental data in solution (Table 1) shows different trend (Li⁺ has the weaker binding energy than Na⁺ and K⁺ whereas Mg²⁺ has the weaker binding energy than Ca²⁺). This behavior has also been reported by other researchers [27–30]. The most probable reason for this inconsistency is due to the solvation process of ions. Small ions such as Li⁺ and Mg²⁺ have high solvation free energy that causes them to be more solvalized in the solution, and have a less interaction with free ligand.

4. Conclusion

From the conductomeric results obtained on the thermodynamics of complexation of DD18C6 with some alkali and alkaline earth metal ions in acetonitrile, methanol, and ethanol solvents and quantum chemical calculations, it can be concluded the following.

(1)The solvents illustrate the fundamental role in the -DD18C6 complexation reactions. The stability of the resulting complexes (except Cs⁺) varied in the order acetonitrile ≫ methanol > ethanol. (2)Although the enthalpy and entropy changes are strongly solvent dependent, the observed increase (or decrease) in value will be compensated by an increase (or decrease) in the corresponding value.(3)The quantum chemical calculations in the gas phase indicate that the binding energy of complexes decreases with increasing cation size.

References

C. J. Pedersen, “Cyclic polyethers and their complexes with metal salts,” Journal of the American Chemical Society, vol. 89, no. 26, pp. 7017–7036, 1967.

S. R. Cooper, Compounds: Toward Future Applications, Wiley-VCH, 1992.

L. F. Lindoy, “New super-and supramolecular receptor systems-cages, chains, squares, and dendrimers incorporating macrocycles as structural elements,” Journal of the Iranian Chemical Society, vol. 1, pp. 1–9, 2004.

K. Gloe, Macrocyclic Chemistry: Current Trends and Future Perspectives, Kluwer Academic, 2005.

K. A. Rubinson, “Chemistry and nerve conduction,” Journal of Chemical Education, vol. 54, no. 6, pp. 345–347, 1977.

P. Bühlmann, E. Pretsch, and E. Bakker, “Carrier-based ion-selective electrodes and bulk optodes. 2. Ionophores for potentiometric and optical sensors,” Chemical Reviews, vol. 98, no. 4, pp. 1593–1687, 1998.

R. M. Izatt, G. A. Clark, J. S. Bradshaw, J. D. Lamb, and J. J. Christensen, “Macrocycle-facilitated transport of ions in liquid membrane systems,” Separation and Purification Methods, vol. 15, no. 1, pp. 21–72, 1986.

R. M. Izatt, J. S. Bradshaw, S. A. Nielsen, J. D. Lamb, J. J. Christensen, and D. Sen, “Thermodynamic and kinetic data for cation-macrocycle interaction,” Chemical Reviews, vol. 85, no. 4, pp. 271–339, 1985.

R. M. Izatt, K. Pawlak, J. S. Bradshaw, and R. L. Bruening, “Thermodynamic and kinetic data for macrocycle interaction with cations, anions, and neutral molecules,” Chemical Reviews, vol. 95, no. 7, pp. 2529–2586, 1995.

R. M. Izatt, K. Pawlak, J. S. Bradshaw, and R. L. Bruening, “Thermodynamic and kinetic data for macrocycle interaction with cations and anions,” Chemical Reviews, vol. 91, no. 8, pp. 1721–2085, 1991.

V. Gutmann, The Donor-Acceptor Approach to Molecular Interactions, Plenum Press, New York, NY, USA, 1978.

Y. C. Wu and W. F. Koch, “Absolute determination of electrolytic conductivity for primary standard KCl solutions from 0 to ,” Journal of Solution Chemistry, vol. 20, no. 4, pp. 391–401, 1991.

D. P. Zollinger, E. Bulten, A. Christenhusz, M. Bos, and W. E. van der Linden, “Computerized conductometric determination of stability constants of complexes of crown ethers with alkali metal salts and with neutral molecules in polar solvents,” Analytica Chimica Acta, vol. 198, pp. 207–222, 1987.

S. Katsuta, Y. Ito, and Y. Takeda, “Stabilities in nitromethane of alkali metal ion complexes with dibenzo-18-crown-6 and dibenzo-24-crown-8 and their transfer from nitromethane to other polar solvents,” Inorganica Chimica Acta, vol. 357, no. 2, pp. 541–547, 2004.

F. A. Christy and P. S. Shrivastav, “Conductometric studies on cation-crown ether complexes: a review,” Critical Reviews in Analytical Chemistry, vol. 41, no. 3, pp. 236–269, 2011.

M. Schmidt, K. Baldridge, J. Boatz et al., “General atomic and molecular electronic structure system,” Journal of Computational Chemistry, vol. 14, pp. 1347–1363, 1993.

HyperChem, Professional 7.5, Hypercube, Gainesville, Fla, USA.

Avogadro: an open-source molecular builder and visualization tool. version 1.03., 2011

M. Shamsipur and A. I. Popov, “Multinuclear Nmr study of some alkali and thallium complexes with 1,10-diaza-18-crown-6 in nonaqueous solutions,” Inorganica Chimica Acta, vol. 43, pp. 243–247, 1980.

R. D. Boss and A. I. Popov, “Competitive NMR study of the complexation of several cations with 18-crown-6, 1,10-diaza-18-crown-6, and cryptand-2,2,2 in nonaqueous solutions,” Inorganic Chemistry, vol. 25, no. 11, pp. 1747–1750, 1986.

S. Kulstad and L. Å. Malmsten, “Diaza-crown ethers-II. Stability constants of a cyclic and a bicyclic diaza-polyether with some alkali cations in acetonitrile,” Journal of Inorganic and Nuclear Chemistry, vol. 42, no. 4, pp. 573–578, 1980.

V. J. Gatto, K. A. Arnold, A. M. Viscariello, S. R. Miller, and G. W. Gokel, “Novel synthetic access to 15- and 18-membered ring diaza-bibracchial lariat ethers (BiBLEs) and a study of sidearm-macroring cooperativity in cation binding,” Tetrahedron Letters, vol. 27, no. 3, pp. 327–330, 1986.

S. Kashanian and M. Shamsipur, “Complex formation of alkaline earth cations with benzo-15-crown-5 and some 18-crowns in methanol, dimethylformamide and dimethyl sulfoxide solutions,” Inorganica Chimica Acta, vol. 155, no. 2, pp. 203–206, 1989.

H. J. Buschmann, “Complex formation of alkaline-earth cations with crown ethers and cryptands in methanol solutions,” Journal of Solution Chemistry, vol. 15, no. 6, pp. 453–461, 1986.

Y. Inoue, F. Amano, N. Okada et al., “Thermodynamics of solvent extraction of metal picrates with crown ethers: enthalpy-entropy compensation. Part 1. Stoicheiometric 1:1 complexation,” Journal of the Chemical Society, Perkin Transactions 2, no. 7, pp. 1239–1246, 1990.

E. Grunwald and C. Steel, “Solvent reorganization and thermodynamic enthalpy-entropy compensation,” Journal of the American Chemical Society, vol. 117, no. 21, pp. 5687–5692, 1995.

E. D. Glendening and D. Feller, “An ab initio investigation of the structure and alkaline earth divalent cation selectivity of 18-crown-6,” Journal of the American Chemical Society, vol. 118, no. 25, pp. 6052–6059, 1996.

M. B. More, D. Ray, and P. B. Armentrout, “Intrinsic affinities of alkali cations for 15-crown-5 and 18-crown-6: bond dissociation energies of gas-phase M⁺-crown ether complexes,” Journal of the American Chemical Society, vol. 121, no. 2, pp. 417–423, 1999.

D. Ray, D. Feller, M. B. More, E. D. Glendening, and P. B. Armentrout, “Cation-ether complexes in the gas phase: bond dissociation energies and equilibrium structures of Li⁺(1,2-dimethoxyethane)x, and 2, and Li⁺ (12-crown-4),” Journal of Physical Chemistry, vol. 100, no. 40, pp. 16116–16125, 1996.

A. L. Sargent, B. J. Mosley, and J. W. Sibert, “A theoretical investigation on the Wurster's crown analogue of 18-crown-6,” Journal of Physical Chemistry A, vol. 110, no. 10, pp. 3826–3837, 2006.