Short communication Extraction and DFT Study on the Complexation of K+ with a Hexaarylbenzene - Based Polyaromatic Receptor Emanuel Makrllk,1'* Petr Toman,2 Petr Vanura3 and Rajendra Rathore4 1 Faculty of Applied Sciences, University of West Bohemia, Husova 11, 306 14 Pilsen, Czech Republic 2 Institute of Macromolecular Chemistry, Academy of Sciences of the Czech Republic, Heyrovskeho sq. 2, 162 06 Prague 6, Czech Republic 3 Department of Analytical Chemistry, Institute of Chemical Technology, Technicka 5, 166 28 Prague 6,Czech Republic 4 Department of Chemistry, Marquette University, P.O.Box 1881, Milwaukee, Wisconsin 53201 - 1881, U.S.A. * Corresponding author: E-mail: makrlik@centrum.cz Received: 29-04-2010 Abstract From extraction experiments and /-activity measurements, the exchange extraction constant corresponding to the equilibrium K+ (aq) + 1'Cs+(nb)sipky ^ 1-K+ (nb) + Cs+(aq) taking part in the two-phase water-nitrobenzene system (1 = hexaarylbenzene - based polyaromatic receptor; aq = aqueous phase, nb = nitrobenzene phase) was evaluated as log Kex (K+, 1-Cs+) = -1.0 ± 0.1. Further, the stability constant of the hexaarylbenzene - based polyaromatic receptor • K+ complex (abbrev. 1^K+) in nitrobenzene saturated with water was calculated for a temperature of 25 °C: log Pnb (1^K+) = 5.1 ± 0.2. By using quantum mechanical DFT calculations, the most probable structure of the 1^K+ complex species was solved. In this complex having C3 symmetry, the cation K+ synergistically interacts with the polar ethereal oxygen fence and with the central hydrophobic benzene bottom via cation - n interaction. Finally, the calculated binding energy of the resulting complex 1K is -282.8 kJ mol-1, confirming the relatively high stability of the considered cationic complex species. Keywords: Hexaarylbenzene-based polyaromatic receptor, potassium cation, complexation, extraction and stability constants, water-nitrobenzene system, DFT, complex structure 1. Introduction Recently, hexaarylbenzene (HAB) derivatives attract a great attention because of their unique propeller-shaped structure and potential application in molecular electronics and nanotechnology. It has been previously described by employing NMR spectroscopy and X-ray crystallography that HAB - based receptor (abbrev. 1; see Scheme 1) binds a single potassium cation because it interacts both with the polar ethereal fence and with the central benzene ring via cation - n interaction.1 Cation - n interaction is a well-established phenomenon in gas phase, as well as in solid state,2 and is known to play an important role in the stabilization of tertiary structures of various proteins.3 Scheme 1. Structural formula of a hexaarylbenzene (HAB)-based polyaromatic receptor (abbrev. 1). Dicarbollylcobaltate anion4 and some of its halogen derivatives are very useful reagents for the extraction of various metal cations (especially Cs+, Sr2+, Ba2+, Eu3+ and Am3+) from aqueous solutions into a polar organic phase, both under laboratory conditions for purely theoretical or analytical purposes,5-35 and on the technological scale for the separation of some high-activity isotopes in the reprocessing of spent nuclear fuel and acidic radioactive waste.36-38 In the current work, the stability constant of the HAB - based receptor K+ complex species (1-K+) in nitrobenzene saturated with water was determined. Moreover, applying quantum mechanical DFT calculations, the most probable structure of the mentioned cationic complex was derived. 2. Experimental Preparation of the electroneutral HAB - based receptor 1 (see Scheme 1) is described elsewhere.1 Cesium 3,3'-commo-bis(undecahydro-1,2-dicarba-3-cobalta-clo-so-dodecabor)ate, Cs+{[n-(3)-B9C5H11]2Co(III)}- (called also cesium dicarbollylcobaltate, abbrev. CsDCC), was synthesized by means of the method published by Hawthorne et al.39 The other chemicals used (Lachema, Brno, Czech Republic) were of reagent grade purity. The radionuclide 137Cs+ (Techsnaveksport, Russia) was of standard radiochemical purity. The extraction experiments were carried out in 10 cm3 glass test-tubes with polyethylene stoppers: 2 cm3 of an aqueous solution of KCl of the concentration in the range from 1 x 10-3 to 3 x 10-3 mol dm-3 and microa-mounts of 137Cs+ were added to 2 cm3 of a nitrobenzene solution of 1 and CsDCC, the concentrations of which varied also from 1 x 10-3 to 3 x 10-3 mol dm-3 (in all experiments, the initial concentration of 1 in nitrobenzene, Ci1n,nb, was always equal to the initial concentration of CsDCC in this medium, CiCns,nDbCC). The test-tubes filled with the solutions were shaken for 2 hours at 25 ± 1 °C, using a laboratory shaker. Then the phases were separated by centrifu-gation. Afterwards, 1 cm3 samples were taken from each phase and their /-activities were measured using a welltype NaI(T1) scintillation detector connected to a /-analy-zer NK 350 (Gamma, Budapest, Hungary). The equilibrium distribution ratios of cesium, DCs, were determined as the ratios of the measured radioactivities of 137Cs+ in the nitrobenzene and aqueous samples. 3. Results and Discussion Regarding the results of previous papers,4,40 the two-phase water-KCl-nitrobenzene-CsDCC extraction system can be described by the following equilibrium with the corresponding exchange extraction constant Kex (K+, Cs+); aq and nb denote the presence of the species in the aqueous and nitrobenzene phases, respectively. For the constant K (K+, Cs+) one can write40 logKex(K+,Cs+) = logK' -logK' (2) where KK+ and KCs+ are the individual extraction constants for K+ and Cs+, respectively, in the water-nitrobenzene system.40 Knowing the values log KK+= -4.140 and log KCs+ = -2.7,40 the exchange extraction constant Kex (K+, Cs+) was simply calculated from Eq. (2) as log Kex (K+, Cs+) = -1.4. Previous results41-44 indicated that the two-phase water-KCl-nitrobenzene-1 (1 = HAB - based receptor)-CsDCC extraction system (see Experimental), chosen for determination of the stability constant of the complex 1*K+ in water-saturated nitrobenzene, can be characterized by the main chemical equilibrium K+(aq) + 1-Cs+(nb) 1-K+(nb) + Cs+(aq) (3) with the respective equilibrium extraction constant Ke (K+, 1-Cs+): Kex(K+,l-Cs+) = [l-K+]nb[Cs+]a [K+l[l.Cs+]n (4) It is necessary to emphasize that 1 is a considerably hydrophobic receptor, practically present in the nitrobenzene phase only, where it forms - with K+ and Cs+ - the relatively stable complexes 1-K+ and 1*Cs+ in the nitrobenzene phase. Taking into account the conditions of elec-troneutrality in the organic and aqueous phases of the system under study, the mass balances of the univalent cations studied at equal volumes of the nitrobenzene and aqueous phases, as well as the measured distribution ratio of cesium, DCs = [1*Cs+]nb / [Cs+]aq, combined with Eq. (4), we obtain the final expression for Kex (K+, 1*Cs+) in the form Kex(K+, 1 • Cs+) = 1 C" D n + D Cs V Cs * KCl CsDCC (5) K+(aq) + Cs+(nb) K+(nb) + Cs+(aq) (1) where CiKnC,nlb is the initial concentration of KCl in aqueous phase of the system under consideration. In this study, from the extraction experiments and /-activity measurements (see Experimental) by means of Eq. (5), the following value of the constant Kex (K+, 1*Cs+) was determined: log Kex (K+, 1-Cs+) = -1.0 ±0.1. Furthermore, with respect to previous results,41-44 for the extraction constants Kex (K+, Cs+) and Kex (K+, 1*Cs+) defined above, as well as for the stability constants of the complexes 1-K+ and 1*Cs+ in nitrobenzene saturated with water, denoted by Pnb (1*K+) and Pnb (1*Cs+), respectively, one gets log pnb (1-K+) = log pnb (1-Cs+) + log K (K+, 1-Cs+) - log Kex (K+, Cs+) (6) Using the constants log Kex (K+, Cs+) and log Kex (K+, 1-Cs+) given above, the value log Pnb (1-Cs+) = 4.7 ± 0.145 and applying Eq. (6), we obtain the stability constant of the 1-K+ complex in water-saturated nitrobenzene at 25 °C as log Pnb (1-K+) = 5.1 ± 0.2. In this context it should be noted that the stability constant of the complex species KL+, where L is 18-crown-6, in nitrobenzene saturated with water is log Pnb (KL+) = 9.5.46 This means that in the mentioned nitrobenzene medium, the stability of the 1-K+ complex under study is substantially lower than that of the cationic complex species KL+ (L = 18-crown-6). The quantum mechanical calculations were carried out at the density functional level of theory (DFT, B3LYP functional) using the Gaussian 03 suite of programs.47 The 6 - 31G(d) basis set was used and the optimizations were unconstrained. Although a possible influence of a polar solvent on the detailed structures of 1 and the 1-K+ complex species could be imagined, our quantum mechanical calculations in similar cases, performed in an analogous way, showed very good agreement of experiment with theory.48-53 In the model calculations, we optimized the molecular geometries of the parent HAB - based receptor 1 and its complex with K+. The optimized structure of a free receptor 1 having C3 symmetry, involving a bowl - shaped cavity, which is comprised of an aromatic bottom (i. e. central benzene ring) and an ethereal fence formed by six oxygen atoms from the peripheral aryl groups, is illustrated in Fig. 1. The depth of the cavity, i. e. the distance between the mean plane of the aromatic bottom and that of the ethereal fence, is 2.15 A (1 A = 0.1 nm); the diameter of this cavity in 1 is 6.19 A (see Fig. 1). The structure obtained by the full DFT optimization Figure 1: Two projections of the DFT optimized structure of a free HAB - based receptor 1 [B3LYP/6-31G(d)] (hydrogen atoms omitted for clarity). Figure 2: Two projections of the DFT optimized structure of the 1-K+ complex [B3LYP/6-31G(d)] (hydrogen atoms omitted for clarity). of the cationic complex species 1-K+ is depicted in Fig. 2. In this complex having also C3 symmetry, the cation K+ synergistically interacts with the hydrophilic polar ethereal oxygen fence (the corresponding K+ ••• O distances are gradually 3.01, 3.10, 3.01, 3.10, 3.01 and 3.10 A) and with the central hydrophobic benzene bottom via cation -n interaction (the distance between the mean plane of the bottom benzene ring and the K+ cation is 2.73 A, as pictured in Fig. 2). At this point it should be noted that the depth of the cavity in the 1-K+ complex is about 2.26 A, while its diameter is roughly 6.04 A. This means that the formation of the complex species 1-K+ results in the small tapering of the mentioned cavity and at the same time, in its getting longer. Finally, the interaction energy, E(int), of the 1-K+ complex [calculated as the difference between electronic energies of 1-K+ and isolated K+ and 1 species, E(int) = E(1-K+) - E(K+) - E(1)] was found to be -282.8 kJ mol-1,which confirms quite unambiguously the formation of the considered cationic complex 1-K+. 4. Acknowledgements This work was supported by the Czech Ministry of Education, Youth and Sports, Projects MSM 4977751303 and MSM 6046137307. The computer time at the MetaCen-trum (Project MSM 6383917201), as well as at the Institute of Physics (computer Luna/Apollo), Academy of Sciences of the Czech Republic, is gratefully acknowledged. 5. References 1. R. Shukla, S. V. Lindeman, R. Rathore, J. Am. Chem. Soc. 2006, 128, 5328-5329. 2. J. C. Ma, D. A. Dougherty, Chem. Rev. 1997, 97, 1303-1324. 3. K. Sakurai, T. Mizuno, H. Hiroaki, K. Gohda, J. Oku, T. Ta-naka, Angew. Chem., Int. Ed. Engl. 2005, 44, 6180-6183. 4. E. Makrlik, P. Vanura, Talanta 1985, 32, 423-429. 5. E. Makrlik, P. Vanura, P. Selucky, J. Solution Chem. 2009, 38, 1129-1138. 6. Z. Valentova, E. Makrlik, Acta Chim. Slov. 2007, 54, 175-178. 7. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2008, 55, 430-433. 8. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2008, 55, 223-227. 9. E. Makrlik, J. Budka, P. Vanura, Acta Chim. Slov. 2009, 56, 278-281. 10. E. Makrlik, P. Vanura, P. Selucky, V. A. Babain, I. V. Smir-nov, Acta Chim. Slov. 2009, 56, 718-722. 11. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2009, 56, 475-479. 12. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2009, 56, 973-976. 13. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2010, 57, 470-474. 14. E. Makrlik, P. Vanura, P. Selucky, Acta Chim. Slov. 2010, 57, 485-490. 15. E. Makrlik, P. Vanura, Z Phys. Chem. 2007, 221, 881-886. 16. E. Makrlik, P. Vanura, Z Phys. Chem. 2009, 223, 247-252. 17. E. Makrlik, P. Vanura, P. Selucky, Z Phys. Chem. 2009, 223, 253-261. 18. E. Makrlik, J. Dybal, P. Vanura, Z Phys. Chem. 2009, 223, 713-718. 19. E. Makrlik, J. Budka, P. Vanura, J. Dybal, Monatsh. Chem. 2008, 139, 1349-1351. 20. E. Makrlik, P. Vanura, P. Selucky, Monatsh. Chem. 2008, 139, 597-600. 21. E. Makrlik, J. Budka, P. Vanura, P. Selucky, Monatsh. Chem. 2009, 140, 157-160. 22. E. Makrlik, J. Dybal, P. Vanura, Monatsh Chem. 2009, 140, 251-254. 23. E. Makrlik, P. Vanura, P. Selucky, J. Hálová, J. Radioanal. Nucl. Chem. 2007, 274, 625-629. 24. E. Makrlik, P. Vanura, J. Radioanal. Nucl. Chem. 2008, 275, 217-220. 25. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2008, 275, 229-232. 26. E. Makrlik, P. Vanura, J. Budka, J. Radioanal. Nucl. Chem. 2008, 275, 463-466. 27. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2008, 275, 309-312. 28. E. Makrlik, P. Vanura, J. Radioanal. Nucl. Chem. 2008, 275, 673-675. 29. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2009, 279, 137-142. 30. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2009, 279, 287-291. 31. J. Dybal, E. Makrlik, P. Vanura, J. Radioanal. Nucl. Chem. 2009, 279, 553-559. 32. E. Makrlik, P. Vanura, P. Selucky, V. A. Babain, I. V. Smir-nov, J. Radioanal. Nucl. Chem. 2009, 279, 743-747. 33. E. Makrlik, P. Vanura, Z. Sedláková, J. Radioanal. Nucl. Chem. 2009, 280, 607-611. 34. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2009, 281, 547-551. 35. E. Makrlik, P. Vanura, P. Selucky, J. Radioanal. Nucl. Chem. 2009, 281, 633-638. 36. J. D. Law, K. N. Brewer, R. S. Herbst, T. A. Todd, D. J. Wood, Waste Management 1999, 19, 27-37. 37. V. N. Romanovskiy, I.V. Smirnov, V. A. Babain, T. A. Todd, R. S. Herbst, J. D. Law, K. N. Brewer, Solvent Extr. Ion Exch. 2001, 19, 1-21. 38. J. D. Law, R. S. Herbst, T. A. Todd, V. N. Romanovskiy, V. A. Babain, V. M. Esimantovskiy, I. V. Smirnov, B. N. Zait-sev, Solvent Extr. Ion Exch. 2001, 19, 23-36. 39. M. F. Hawthorne, D. C. Young, T. D. Andrews, D. V. Howe, R. L. Pilling, A. D. Pitts, M. Reintjes, L. F. Warren, P. A. Wegner, J. Am. Chem. Soc. 1968, 90, 879-896. 40. J. Rais, Collect. Czech. Chem. Commun. 1971, 36, 3253-3262. 41. M. Dankova, E. Makrlik, P. Vanura, J. Radioanal. Nucl. Chem. 1997, 227, 251-253. 42. E. Makrlik, J. Halova, M. Kyrs, Collect. Czech. Chem. Commun. 1984, 49, 39-44. 43. E. Makrlik, P. Vanura, J. Radioanal. Nucl. Chem. 1996, 214, 339-346. 44. E. Makrlik, P. Vanura, Monatsh. Chem. 2006, 137, 157-161. 45. E. Makrlik, P. Vanura, P. Selucky, unpublished results. 46. E. Makrlik, P. Vanura, M. Dankova, J. Radioanal. Nucl. Chem. 1999, 240, 579-583. 47. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vre-ven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. Moroku- ma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Zakr-zewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Fore-sman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cio-slowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-La-ham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, M. W. Wong, C. Gonzalez, J. A. Pople, Gaussian 03, Revision C. 02, Gaussian, Inc., Wallingford, CT, 2004. 48. J. Kriz, J. Dybal, E. Makrlik, Biopolymers 2006, 82, 536548. 49. J. Kriz, J. Dybal, E. Makrlik, P. Vanura, J. Lang, Supramol. Chem. 2007, 19, 419-424. 50. J. Kriz, J. Dybal, E. Makrlik, P. Vanura, Supramol. Chem. 2008, 20, 387-395. 51. J. Krriz, J. Dybal, E. Makrlik, J. Budka, P. Varnura, Supramol. Chem. 2008, 20, 487-494. 52. J. Kriz, J. Dybal, E. Makrlik, J. Budka, J. Phys. Chem. A 2008, 112, 10236-10243. 53. J. Kriz, J. Dybal, E. Makrlik, J. Budka, P. Vanura, J. Phys. Chem. A 2009, 113, 5896-5905. Povzetek S pomočjo ekstrakcije in meritev /-aktivnosti smo določili konstanto ekstrakcije log Kex (K+, 1.Cs+) = -1.0 ± 0.1, za ravnotežje K+ (aq) + 1.Cs+(nb) 1.K+ (nb) + Cs+(aq), ki ga dosežemo v dvofaznem sistemu voda-nitrobenzen (1 = heksaaril benzen-poliaromatski receptor; aq = vodna faza, nb = faza nitrobenzena). Pri 25 °C smo določili tudi konstanto stabilnosti K+ kompleksa (1.K+) v nitrobenzenu nasičenem z vodo, log Pnb (1.K+) = 5.1 ± 0.2. Z uporabo DFT kvantomehan-skega računa smo določili najbolj verjetno strukturo 1.K+ kompleksa, ki predvideva C3 simetrijo, K+ kation pa siner-gistično interagira s polarno etersko skupino in centralnim hidrofobnim benzenovim obročem preko kation - n interakcije. Izračunali smo tudi vezno energijo nastalega kompleksa, 1.K+, ki znaša -282.8 kJ mol-1 in potrjuje njegovo relativno visoko stabilnost.