Cation recognition by benzene sandwich compounds - a DFT perspective

Cation-π interactions between alkali, alkaline earth and ammonium cations and sandwich compounds of benzene and the cyclopentadienyl (Cp) anion were studied using quantum chemical CCSD(T)/CBS and DFT (B3LYP/def2-TZVP) calculations. The results show significantly stronger interactions of sandwich compounds with respect to (uncoordinated) benzene. Moreover, very strong cation-π interactions of Cp sandwich compounds are furthermore surpassed by cation-π interactions of benzene sandwich compounds, which are capable of reaching a remarkable interaction energy value of -196.8 kcal mol-1 (Mg2+/W(benzene)2). While there are only small variations of interaction energy values for sandwich compounds of different transition metals (3d metals < 4d < 5d), cation-π interactions progressively become stronger in the following order: (uncoordinated) benzene < Cp sandwich < benzene sandwich. Aside from interaction energies, the cation-π interactions can be assessed by means of their influence on the geometries of sandwich compounds, which are found to strongly correlate with the strength of cation-π interactions. These results emphasize sandwich compounds, particularly those containing C6 aromatic rings, as promising candidates for new receptors for common metal cations.

A practical post-Hartree-Fock approach describing open-shell metal cluster-support interactions. Application to Cu3 adsorption on benzene/coronene

Current advances in synthesizing and characterizing atomically precise monodisperse metal clusters (AMCs) at the subnanometer scale have opened up fascinating possibilities in designing new heterogeneous (photo)catalysts as well as functional interfaces between AMCs and biologically relevant molecules. Understanding the nature of AMC-support interactions at molecular-level is essential for optimizing (photo)catalysts performance and designing novel ones with improved properties. Møller-Plesset second-order perturbation theory (MP2) is one of the most cost-efficient single-reference post-Hartree-Fock wave-function-based theories that can be applied to AMC-support interactions considering adequate molecular models of the support, and thus complementing state-of-the-art dispersion-corrected density functional theory. However, the resulting AMC-support interaction is typically overestimated with the MP2 method and must be corrected. The coupled MP2 (MP2C) scheme replacing the uncoupled Hartree-Fock dispersion energy by a coupled dispersion contribution, has been proven to describe accurately van-der-Waals (vdW)-dominated interactions between closed-shell AMCs and carbon-based supports. In this work, the accuracy of a MP2C-based scheme is evaluated in modelling open-shell AMC-cluster interactions that imply charge transfer or other strong attractive energy contributions beyond vdW forces. For this purpose, we consider the interaction of Cu3 with molecular models of graphene of increasing size (benzene and coronene). In this way, it is shown that subchemical precision (within 0.1 kcal mol-1) is achieved with the modified MP2C scheme, using the explicitly correlated coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)-F12] as a benchmark method. It is also revealed that the energy difference between uncoupled and coupled dispersion terms closely follows benchmark values of the repulsive intramonomer correlation contribution. The proposed open-shell MP2C-based approach is expected to be of general applicability to open-shell atomic or molecular species interacting with coronene for regions of the potential landscape where single-reference electronic structure descriptions suffice.

Ion Coupling, Bonding, and Transfer in Narrow Carbon Nanotubes

Narrow carbon nanotubes (nCNT) are unique mimics of biological channels with water-ion selectivity attractive for applications such as water purification and osmotic energy harvesting, yet their understanding is still incomplete. Here, an ab initio computation is employed to develop the full picture of ion transfer in nCNT including specificity and coupling between ions. The thermodynamic costs of ion transfer are computed for single ions and ion pairs and used to evaluate different local coupling scenarios including strong (pairing) and weak (free-ion) coupling as well as "electroneutrality breakdown" (EB), possible for cations only due to their chemisorption-like interaction with nCNT. The results also indicate that nCNT behaves as a highly polarizable metal-like shell, which eliminates the dielectric energy when CNT accommodates coupled cation and anion. This allows facile computation and comparison of the full transfer costs, including translation entropy, for different ions in different coupling modes to identify the dominant regime. EB transfer appears most favorable for K+, while anions strongly favor transfer as pairs, except for chloride which favors weak coupling and, at neutral pH, transfers as a trace ion coupled to both cation and OH-. The results demonstrate that, in general, observed ion permeation and conduction in nCNT, especially for anions, reflect a complex ion-specific and composition-dependent interplay between different ions.

Single-Molecule Force Spectroscopy Reveals Cation-π Interactions in Aqueous Media Are Highly Affected by Cation Dehydration

Cation-π interactions underlie many important processes in biology and materials science. However, experimental investigations of cation-π interactions in aqueous media remain challenging. Here, we studied the cation-π binding strength and mechanism by pulling two hydrophobic polymers with distinct cation binding properties, i.e., poly-pentafluorostyrene and polystyrene, in aqueous media using single-molecule force spectroscopy and nuclear magnetic resonance measurement. We found that the interaction strengths linearly depend on the cation concentrations, following the order of Li^{+}<NH_{4}^{+}<Na^{+}<K^{+}. The binding energies are 0.03-0.23  kJ mol^{-1} M^{-1}. This order is distinct from the strength of cation-π interactions in gas phase and may be caused by the different dehydration ability of the cations. Taken together, our method provides a unique perspective to investigate cation-π interactions under physiologically relevant conditions.

Contrasting the Mechanism of H2 Activation by Monomeric and Potassium-Stabilized Dimeric AlI Complexes: Do Potassium Atoms Exert any Cooperative Effect?

Aluminyl anions are low-valent, anionic, and carbenoid aluminum species commonly found stabilized with potassium cations from the reaction of Al-halogen precursors and alkali compounds. These systems are very reactive toward the activation of σ-bonds and in reactions with electrophiles. Various research groups have detected that the potassium atoms play a stabilization role via electrostatic and cation ⋯ π interactions with nearby (aromatic)-carbocyclic rings from both the ligand and from the reaction with unsaturated substrates. Since stabilizing K⋯H bonds are witnessed in the activation of this class of molecules, we aim to unveil the role of these metals in the activation of the smaller and less polarizable H₂ molecule, together with a comprehensive characterization of the reaction mechanism. In this work, the activation of H₂ utilizing a NON-xanthene-Al dimer, [K{Al(NON)}]₂ (D) and monomeric, [Al(NON)]^- (M) complexes are studied using density functional theory and high-level coupled-cluster theory to reveal the potential role of K⁺ atoms during the activation of this gas. Furthermore, we aim to reveal whether D is more reactive than M (or vice versa), or if complicity between the two monomer units exits within the D complex toward the activation of H₂ . The results suggest that activation energies using the dimeric and monomeric complexes were found to be very close (around 33 kcal mol^-1 ). However, a partition of activation energies unveiled that the nature of the energy barriers for the monomeric and dimeric complexes are inherently different. The former is dominated by a more substantial distortion of the reactants (and increased interaction energies between them). Interestingly, during the oxidative addition, the distortion of the Al complex is minimal, while H₂ distorts the most, usually over 0.77 Δ E d i s t ≠ . Overall, it is found here that electrostatic and induction energies between the complexes and H₂ are the main stabilizing components up to the respective transition states. The results suggest that the K⁺ atoms act as stabilizers of the dimeric structure, and their cooperative role on the reaction mechanism may be negligible, acting as mere spectators in the activation of H₂ . Cooperation between the two monomers in D is lacking, and therefore the subsequent activation of H₂ is wholly disengaged.

The effect of cation-π interactions on the stability and electronic properties of anticancer drug Altretamine: a theoretical study

The present work utilizes density functional theory (DFT) calculations to study the influence of cation-π interactions on the electronic properties of the complexes formed by Altretamine [2,4,6-tris(dimethylamino)-1,3,5-triazine], an anticancer drug, with mono- and divalent (Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺ and Ca²⁺) metal cations. The structures were optimized with the M06-2X method and the 6-311++G(d,p) basis set in the gas phase and in solution. The theory of `Atoms in Molecules' (AIM) was applied to study the nature of the interactions by calculating the electron density ρ(r) and its Laplacian at the bond critical points. The charge-transfer process during complexation was evaluated using natural bond orbital (NBO) analysis. The results of DFT calculations demonstrate that the strongest/weakest interactions belong to Be²⁺/K⁺ complexes. There are good correlations between the achieved densities and the amounts of charge transfer with the interaction energies. Finally, the stability and reactivity of the cation-π interactions can be determined by quantum chemical computation based on the molecular orbital (MO) theory.

Curvature and size effects hinder halogen bonds with extended π systems

The curvature of aromatic systems strengthens the interaction by the concave face while it weakens by the convex one. Parallel structures are favoured over halogen bonded ones.

Substituent effects in the so-called cationπ interaction of benzene and its boron-nitrogen doped analogues: overlooked role of σ-skeleton

Despite massive efforts to pinpoint the substituent effects in the so-called cationπ systems, no consensus has been yet reached on how substituents exercise their effects in the interaction of the aromatic molecule with the metal ion. The π-polarization (the Hunter model) and the direct local effect (the Wheeler-Houk model) are two lines of thought applied to this problem, but the justification of both approaches is based on insufficiently proven assumptions and approximations. In order to shed more light on this issue we propose a new approach which enables us to gauge directly the energetic trends resulting from the interaction of the ring with the cation. In our method we add one more partitioning level to the interacting quantum atoms (IQA) scheme and decompose the IQA interaction energies into contributions resulting from σ and π electron densities of the aromatic ring. The new approach, which is named partitioned-IQA, abbreviated as p-IQA, has been applied to complexes of derivatives of benzene or azaborine interacting with a sodium cation. The p-IQA approach reveals that in these systems both σ and π electronic moieties are polarized. Interestingly, for the majority of cases the σ-polarization outweighs the π one, contrary to the Hunter model. However, the Wheeler-Houk model is not precise, either, since the σ-polarization shows some degree of non-locality. In addition, the substituents are found to have a negligible influence on the ring orbital-overlapping capability, i.e. the covalency. Therefore, the substituent effect in the cationπ interaction is a nonlocal classical effect, indicating that neither Hunter model nor Wheeler-Houk model is able to fully describe all the aspects of the substituent effects. The p-IQA conclusions for the considered systems have been compared with the results from the functional-group SAPT (F-SAPT) method. We believe that the presented partitioning in the IQA framework will provide a deeper insight into the substituent effects in the cationπ interactions, which is beyond the σ-π atomic charge population separation.

On the Nature of σ-σ, σ-π, and π-π Stacking in Extended Systems

Stacking interactions have been evaluated, employing computational methods, in dimers formed by analogous aliphatic and aromatic species of increasing size. Changes in stability as the systems become larger are mostly controlled by the balance of increasing repulsion and dispersion contributions, while electrostatics plays a secondary but relevant role. The interaction energy increases as the size of the system grows, but it does much faster in π-π dimers than in σ-π complexes and more remarkably than in σ-σ dimers. The main factor behind the larger stability of aromatic dimers compared to complexes containing aliphatic molecules is related to changes in the properties of the aromatic systems due to electron delocalization leading to larger dispersion contributions. Besides, an extra stabilization in π-π complexes is due to the softening of the repulsive wall in aromatic species that allows the molecules to come closer.

Ion-Migration Inhibition by the Cation-π Interaction in Perovskite Materials for Efficient and Stable Perovskite Solar Cells

Migration of ions can lead to photoinduced phase separation, degradation, and current-voltage hysteresis in perovskite solar cells (PSCs), and has become a serious drawback for the organic-inorganic hybrid perovskite materials (OIPs). Here, the inhibition of ion migration is realized by the supramolecular cation-π interaction between aromatic rubrene and organic cations in OIPs. The energy of the cation-π interaction between rubrene and perovskite is found to be as strong as 1.5 eV, which is enough to immobilize the organic cations in OIPs; this will thus will lead to the obvious reduction of defects in perovskite films and outstanding stability in devices. By employing the cation-immobilized OIPs to fabricate perovskite solar cells (PSCs), a champion efficiency of 20.86% and certified efficiency of 20.80% with negligible hysteresis are acquired. In addition, the long-term stability of cation-immobilized PSCs is improved definitely (98% of the initial efficiency after 720 h operation), which is assigned to the inhibition of ionic diffusions in cation-immobilized OIPs. This cation-π interaction between cations and the supramolecular π system enhances the stability and the performance of PSCs efficiently and would be a potential universal approach to get the more stable perovskite devices.

Cation-π interactions in protein-ligand binding: theory and data-mining reveal different roles for lysine and arginine

The interactions of neutral aromatic ligands with cationic arginine, histidine and lysine amino acid residues have been studied with ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of X-ray structures.

We have studied the cation–π interactions of neutral aromatic ligands with the cationic amino acid residues arginine, histidine and lysine using ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of all available Protein Data Bank (PDB) X-ray structures. Quantum chemical potential energy surfaces (PES) for these interactions were obtained at the DLPNO-CCSD(T) level of theory and compared against the empirical distribution of 2012 unique protein–ligand cation–π interactions found in X-ray crystal structures. We created a workflow to extract these structures from the PDB, filtering by interaction type and residue pK_a. The gas phase cation–π interaction of lysine is the strongest by more than 10 kcal mol^–1, but the empirical distribution of 582 X-ray structures lies away from the minimum on the interaction PES. In contrast, 1381 structures involving arginine match the underlying calculated PES with good agreement. SAPT analysis revealed that underlying differences in the balance of electrostatic and dispersion contributions are responsible for this behavior in the context of the protein environment. The lysine–arene interaction, dominated by electrostatics, is greatly weakened by a surrounding dielectric medium and causes it to become essentially negligible in strength and without a well-defined equilibrium separation. The arginine–arene interaction involves a near equal mix of dispersion and electrostatic attraction, which is weakened to a much smaller degree by the surrounding medium. Our results account for the paucity of cation–π interactions involving lysine, even though this is a more common residue than arginine. Aromatic ligands are most likely to interact with cationic arginine residues as this interaction is stronger than for lysine in higher polarity surroundings.

Structural and energetic study of cation–π–cation interactions in proteins

Statistical and energetic analysis of cation–π–cation motifs in protein structures suggests a potential stabilizing role in the protein fold.

A theoretical study of complexes formed between cations and curved aromatic systems: electrostatics does not always control cation–π interaction

Cation–π interactions in curved aromatic systems are not controlled by electrostatics; induction and dispersion dominate in most cases studied.

Effect of the cation on the stability of cation–glyme complexes and their interactions with the [TFSA]− anion

The interactions of glymes with alkali or alkaline earth metal cations depend strongly on the metal cations.

Improving the Force Field Description of Tyrosine-Choline Cation-π Interactions: QM Investigation of Phenol-N(Me)4+ Interactions

Cation-π interactions between tyrosine amino acids and compounds containing N,N,N-trimethylethanolammonium (N(CH3)3) are involved in the recognition of histone tails by chromodomains and in the recognition of phosphatidylcholine (PC) phospholipids by membrane-binding proteins. Yet, the lack of explicit polarization or charge transfer effects in molecular mechanics force fields raises questions about the reliability of the representation of these interactions in biomolecular simulations. Here, we investigate the nature of phenol-tetramethylammonium (TMA) interactions using quantum mechanical (QM) calculations, which we also use to evaluate the accuracy of the additive CHARMM36 and Drude polarizable force fields in modeling tyrosine-choline interactions. We show that the potential energy surface (PES) obtained using SAPT2+/aug-cc-pVDZ compares well with the large basis-set CCSD(T) PES when TMA approaches the phenol ring perpendicularly. Furthermore, the SAPT energy decomposition reveals comparable contributions from electrostatics and dispersion in phenol-TMA interactions. We then compared the SAPT2+/aug-cc-pVDZ PES obtained along various approach directions to the corresponding PES obtained with CHARMM, and we show that the force field accurately reproduces the minimum distances while the interaction energies are underestimated. The use of the Drude polarizable force field significantly improves the interaction energies but decreases the agreement on distances at energy minima. The best agreement between force field and QM PES is obtained by modifying the Lennard-Jones terms for atom pairs involved in the phenol-TMA cation-π interactions. This is further shown to improve the correlation between the occupancy of tyrosine-choline cation-π interactions obtained from molecular dynamics simulations of a bilayer-bound bacterial phospholipase and experimental affinity data of the wild-type protein and selected mutants.

Theoretical insight into the interaction between SnX2 (X = H, F, Cl, Br, I) and benzene

For a series of five model complexes composed of a singlet SnX₂ molecule (X = H, F, Cl, Br, I) and a benzene molecule, the first-principles calculations of their energetics and the analysis of their electron density topology have been performed. The CCSD(T)/CBS interaction energy between SnX₂ and C₆H₆ fall into the range between −10.0 and −11.2 kcal/mol, which indicates that the complexes are rather weakly bound. The relevant role of electrostatic and dispersion contributions to the interaction energy between SnX₂ and C₆H₆ is highlighted in the results obtained from the symmetry-adapted perturbation theory (SAPT). The electron density topological analysis has been carried out using the quantum theory of atoms in molecules (QTAIM) and the noncovalent interactions (NCI) visualization index. Both QTAIM and NCI prove the closed-shell, noncovalent and attractive character of the interaction. A very small charge transfer from C₆H₆ to SnX₂ has been detected. The formation of the five complexes is accompanied by the electron density deformations that are spatially restricted mostly to the region around the Sn atom and its adjacent C atom. The results presented in this work shed some light on the nature of the interactions associated with crystalline structural motifs involving low-valent tin complexed with neutral aryl rings.

The n → π* interaction: a rapidly emerging non-covalent interaction

This perspective describes the current status of a recently discovered non-covalent interaction named as the n → π* interaction, which is very weak and counterintuitive in nature. In this review, we have provided a brief overview of the widespread presence of this interaction in biomacromolecules, small biomolecules and materials, as well as the physical nature of this interaction explored using various experimental and theoretical techniques. It has been found that this interaction is equally important to other non-covalent interactions for the stability and specific structures of biomolecules and materials. An in-depth understanding of this interaction can help in designing more efficient functional materials as well as drugs. The review also provides a future outlook in terms of exploring the detailed functional role of this interaction in biological processes and its direct spectroscopic evidence, which other commonly known non-covalent interactions (conventional hydrogen bonding, π-hydrogen bonding, π-stacking, etc.) have.

Is Aerogen-π Interaction Capable of Initiating the Noncovalent Chemistry of Group 18?

The interactions between atoms of noble gases and π systems are generally considered as van der Waals interaction, which have not attracted attention yet. Herein, we present high-level ab initio calculations to show the unexpected noncovalent interaction between a covalently bonded noble gas atom and a delocalized aromatic π electron using XeO3⋅benzene as the prototype. The CCSD(T)/CBS reference data show its strength amounting to -10.2 kcal mol(-1), comparable to a typical H-bond or an anion-π interaction. The energy decomposition analysis reveals that the aerogen-π interaction is favored by the electrostatic interaction (27.7%), the induction (13.4%), and the dispersion (21.6%). This interaction may prompt us to consider the noncovalent chemistry of aerogen derivatives in the near future.

On the interaction between the imidazolium cation and aromatic amino acids. A computational study

Phe, Tyr and Trp form parallel complexes with cation⋯π interactions. His complexes are the strongest, but without making contact with the aromatic cloud.

Effect of stepwise microhydration on the guanidinium···π interaction

The characteristics of the interaction of microhydrated guanidinium cation with the aromatic moieties present in the aromatic amino acids side chains have been studied by means of computational methods. The most stable minima found for non-hydrated complexes correspond in all cases to structures with guanidinium oriented toward the ring and interacting by means of N-H···π hydrogen bonds. The interaction becomes stronger when going from benzene (-14 kcal mol⁻¹) to phenol (-17 kcal mol⁻¹) to indole (-21 kcal mol⁻¹). These complexes are held together mainly by electrostatics, but with important contributions from induction and dispersion. The presence of a small number of water molecules significantly affects the characteristics of the complexes. Hydrogen bonds formed by water with the cation, another water molecule, or the aromatic units become more and more similar in intensity as water molecules are included in the complex, leading to a great variety of minima with similar stability but showing very different structural patterns. The behavior is similar with the three aromatic units, the differences in stability mainly being a consequence of the different strength of the cation···π contact.

Predicting the cation–π binding of substituted benzenes: energy decomposition calculations and the development of a cation–π substituent constant

This work proposes a new substituent constant,Π⁺, to describe cation–π binding using computational methods at the MP2(full)/6-311++G** level of theory with Symmetry Adapted Perturbation Theory (SAPT) calculations on selected cation–π complexes.

Microsolvation of the potassium ion with aromatic rings: comparison between hexafluorobenzene and benzene

We employ a recently developed methodology to study structural and energetic properties of the first solvation shells of the potassium ion in nonpolar environments due to aromatic rings, which is important to understand the selectivity of several biochemical phenomena. Our evolutionary algorithm is used in the global optimization study of clusters formed of K(+) solvated with hexafluorobenzene (HFBz) molecules. The global intermolecular interaction for these clusters has been decomposed in HFBz-HFBz and in K(+)-HFBz contributions, using a potential model based on different decompositions of the molecular polarizability of hexafluorobenzene. Putative global minimum structures of microsolvation clusters up to 21 hexafluorobenzene molecules were obtained and compared with the analogous K(+)-benzene clusters reported in our previous work (J. Phys. Chem. A 2012, 116, 4947-4956). We have found that both K(+)-(Bz)n and K(+)-(HFBz)n clusters show a strong magic number around the closure of the first solvation shell. Nonetheless, all K(+)-benzene clusters have essentially the same first solvation shell geometry with four solvent molecules around the ion, whereas the corresponding one for K(+)-(HFBz)n is completed with nine HFBz species, and its structural motif varies as n increases. This is attributed to the ion-solvent interaction that has a larger magnitude for K(+)-Bz than in the case of K(+)-HFBz. In addition, the ability of having more HFBz than Bz molecules around K(+) in the first solvation shell is intimately related to the inversion in the sign of the quadrupole moment of the two solvent species, which leads to a distinct ion-solvent geometry of approach.

Accurate prediction of cation-π interaction energy using substituent effects

Substituent effects on cation-π interactions have been quantified using a variety of Φ-X···M(+) complexes where Φ, X, and M(+) are the π-system, substituent, and cation, respectively. The cation-π interaction energy, E(M(+)), showed a strong linear correlation with the molecular electrostatic potential (MESP) based measure of the substituent effect, ΔV(min) (the difference between the MESP minimum (V(min)) on the π-region of a substituted system and the corresponding unsubstituted system). This linear relationship is E(M(+)) = C(M(+))(ΔV(min)) + E(M(+))' where C(M(+)) is the reaction constant and E(M(+))' is the cation-π interaction energy of the unsubstituted complex. This relationship is similar to the Hammett equation and its first term yields the substituent contribution of the cation-π interaction energy. Further, a linear correlation between C(M(+))() and E(M(+))()' has been established, which facilitates the prediction of C(M(+)) for unknown cations. Thus, a prediction of E(M(+)) for any Φ-X···M(+) complex is achieved by knowing the values of E(M(+))' and ΔV(min). The generality of the equation is tested for a variety of cations (Li(+), Na(+), K(+), Mg(+), BeCl(+), MgCl(+), CaCl(+), TiCl(3)(+), CrCl(2)(+), NiCl(+), Cu(+), ZnCl(+), NH(4)(+), CH(3)NH(3)(+), N(CH(3))(4)(+), C(NH(2))(3)(+)), substituents (N(CH(3))(2), NH(2), OCH(3), CH(3), OH, H, SCH(3), SH, CCH, F, Cl, COOH, CHO, CF(3), CN, NO(2)), and a large number of π-systems. The tested systems also include multiple substituted π-systems, viz. ethylene, acetylene, hexa-1,3,5-triene, benzene, naphthalene, indole, pyrrole, phenylalanine, tryptophan, tyrosine, azulene, pyrene, [6]-cyclacene, and corannulene and found that E(M)(+) follows the additivity of substituent effects. Further, the substituent effects on cationic sandwich complexes of the type C(6)H(6)···M(+)···C(6)H(5)X have been assessed and found that E(M(+)) can be predicted with 97.7% accuracy using the values of E(M(+))' and ΔV(min). All the Φ-X···M(+) systems showed good agreement between the calculated and predicted E(M(+))() values, suggesting that the ΔV(min) approach to substituent effect is accurate and useful for predicting the interactive behavior of substituted π-systems with cations.

Alkali-ion microsolvation with benzene molecules

The target of this investigation is to characterize by a recently developed methodology, the main features of the first solvation shells of alkaline ions in nonpolar environments due to aromatic rings, which is of crucial relevance to understand the selectivity of several biochemical phenomena. We employ an evolutionary algorithm to obtain putative global minima of clusters formed with alkali-ions (M(+)) solvated with n benzene (Bz) molecules, i.e., M(+)-(Bz)(n). The global intermolecular interaction has been decomposed in Bz-Bz and in M(+)-Bz contributions, using a potential model based on different decompositions of the molecular polarizability of benzene. Specifically, we have studied the microsolvation of Na(+), K(+), and Cs(+) with benzene molecules. Microsolvation clusters up to n = 21 benzene molecules are involved in this work and the achieved global minimum structures are reported and discussed in detail. We observe that the number of benzene molecules allocated in the first solvation shell increases with the size of the cation, showing three molecules for Na(+) and four for both K(+) and Cs(+). The structure of this solvation shell keeps approximately unchanged as more benzene molecules are added to the cluster, which is independent of the ion. Particularly stable structures, so-called "magic numbers", arise for various nuclearities of the three alkali-ions. Strong "magic numbers" appear at n = 2, 3, and 4 for Na(+), K(+), and Cs(+), respectively. In addition, another set of weaker "magic numbers" (three per alkali-ion) are reported for larger nuclearities.