Main-group metal cyclophane complexes with high coordination numbers

Density functional theory calculations using the PBE0-D3BJ hybrid functional have been employed to investigate the complexation of main-group metal-cations with [2.2.2]paracyclophane and deltaphane. Geometry optimization under symmetry constraints was performed to observe the mode of coordination that a metal-cation adopts when it resides inside the cyclophane cavity. Thermodynamic properties were investigated to note the trends of stability along a group of metals. To further investigate the bonding properties, Morokuma–Ziegler energy decomposition analysis, natural bond orbital analysis and Bader's analysis were employed. It was observed that most of the main-group metal complexes with cyclophanes prefer an η6η6η6 coordination mode where the metal-cation sits in the centre of the cyclophane cavity. There is an increased thermodynamic stability in [2.2.2]paracyclophane complexes compared to their deltaphane analogues while the reverse is true regarding the strength of coordination based on interaction energy.


Introduction
Cyclophanes consist of two or more aromatic rings connected through aliphatic bridges, to form a cyclic cavity. 1 The aromatic rings are characterized by a planar geometry but the strain imposed by the aliphatic components renders the geometry of the cyclophanes twisted. 2 There is a growing interest in the donor-acceptor complexes of cyclophanes with metals due to their various applications such as the development of ion-selective electrodes, 3 catalysis 4 and chelation. 5 Metal-chelating agents involving cyclophanes are potential candidates to be used in waste-water treatment since they are known for metal-ion scavenging. 3,6 Host molecules undergo conformational changes upon incorporation of the guest species in their macrocyclic cavity and these differences may be observed sometimes through uorescence. This allows recognition of targetted guest metal-cations. 7 Moreover, derivatives of metallacyclophane hosts have been reported as potential candidates for optical biomolecular recognition. 8 Luminescent metal-complexes of cyclophanes have been reported to have their potential applications in thin-lm nonoporous materials. 9 All these applications are due to the potential of the p-rich cyclophane cavity to host electrophilic or cationic guests by making inclusion complexes, despite the fact that some smaller cyclophanes are also known to form exclusion complexes. 10 The selectivity and sensitivity of cyclophanes can be improved as evident from the efforts in the past such as derivatization and functionalization using different coordination groups, 4,11,12 controlling the size of the cavity to capture guest species 13 and replacing the usual phenyl ring with heterocyclic aromatic rings such as imidazolium. 14 Owing to the variety of their potential applications, it is important to explore the bonding properties of different possible metallacyclophanes. In the current study, the two cyclophanes [2.2.2]paracyclophane (pCp) and deltaphane (Dp) given in Fig. 1 were selected for this purpose.
pCp was rst synthesized by Pierre and co-workers where they called it a p-prismand due to its p-rich prism-shaped cavity 15 and demonstrated its complexation with silver triate. Extending the concept and introducing an increased rigidity compared to that of pCp, Kang et al. reported the synthesis of Dp. 16 As the systematic name of the latter i.e. [2,2,2,2,2,2](1,2,4,5)cyclophane suggests, Dp has six ethano bridges compared to three in pCp. pCp and Dp consist of three phenyl rings each bridged through aliphatic chains. However, pCp is more exible as each aromatic ring is attached to two ethyl bridges unlike its attachment to four in case of Dp.
Due to the presence of phenyl rings, both these molecules contain p-rich cavities that can coordinate to the metal-cations to form stable complexes. It was shown earlier that the silver triate complex of pCp is much more stable compared to other such complexes with aromatic systems. 15 Similarly, the crystal structure of silver triate complex of Dp was also reported. 16 During the same period (mid to late 1980s), Schmidbaur and coworkers reported the groundbreaking complexes of pCp with some of the main-group metals. [17][18][19] Despite these contributions in synthesis of metallacyclophanes of pCp and Dp, a lot of effort is still required to explore the bonding properties of these complexes to make use of these promising p-donating ligands on industrial scale.
In this regard, we previously reported the bonding properties of coinage metal complexes of pCp and Dp. 20 Earlier, the groups of Frenking and Castro carried out a computational study on the coordination mode and bonding properties of inclusion complexes of Sn 2+ and Ag + with pCp. 21 They further extended the concept to computationally understand the role of formal charge of a cation in p-cation interactions by comparing the complexes of isoelectronic In + and Cd 2+ with pCp. 22 In further instances, Castro et al. investigated helicenes 23 and Dp 24 as potential p-donors to form various p-cation interactions through relativistic DFT approach. A review of the use of relativistic computational tools to study the structural and bonding properties of these p-cation interactions has been published recently 25 which shows that different avenues in this eld have been opened by the groups of Castro and Frenking over the last few years. The current work is aimed at investigating the coordination and bonding properties of the complexes of group 13 (Ga + , In + and Tl + ), 14 (Ge 2+ , Sn 2+ and Pb 2+ ) and 15 (As 3+ , Sb 3+ and Bi 3+ ) with the p-rich cavity of pCp and Dp using density functional theory (DFT).

Computational details
All calculations related to geometry optimization were performed with the Gaussian 09 suite of programs (Revision D.01). 26 The PBE0 hybrid functional 27 in conjunction with Grimme's empirical D3 correction with Becke-Johnston damping (D3BJ) 28 was used in combination with the def2-TZVP basis-set 29 of triple-z quality in all these calculations.
The optimized structures were further subject to the Morokuma-Ziegler Energy Decomposition Analysis (MZEDA) 30 that was carried out with the ADF2014 program. 31 The TZ2P (Slater Type Orbital) basis set 32 was employed along with the relativistic ZORA Hamiltonian. 33,34 MZEDA involves the decomposition of total energy (DE) of a molecule as: In eqn (1), DE int is the instantaneous interaction between the two molecular fragments and (DE 1 +DE 2 ) is the sum of their individual energies. The above equation implies that DE int is the difference between the total energy of a molecule and its fragments. DE int can be further subdivided as in eqn (2): Here, DE elstat is the energy due to electrostatic interactions (mostly attractive in nature) between the molecular fragments. DE Pauli is the repulsion term and it arises due to the electrons with same spin. DE orb indicates the interactions involving charge transfer polarization effects. Bader's analysis based on quantum theory of atoms in molecules (QTAIM) 35 given by Richard Bader was performed using Multiwfn soware. 36 NBO analysis was carried out with NBO 6.0 program 37 as interfaced with Gaussian09.

Structural features
Un-complexed cyclophanes. pCp possesses a rigid geometry and has an internal cavity with diameter 2.5Å. 39 The phenyl rings are bridged at the para positions through ethyl chains. Dp is even more rigid as the phenyl rings are doubly connected to each other at the ortho and meta positions through ethano bridges. The D 3 symmetric geometries were optimized at PBE0-def2TZVP level of DFT and have been characterized as minima on the potential energy surface. Additionally, the two have D 3h symmetric transition structures optimized previously at the same level. 20 Calculated structural data of both of these ligands is consistent with their crystal structural information 16,39 except the fact that gas-phase calculated bond lengths are sometimes longer than those in a crystal structure (solid-phase) due to crystal packing forces that are absent in a gas-phase calculation.
In pCp, the C-C bond lengths in bridging ethyl units and that for bridgeheads are 1.54Å and 1.51Å compared to the experimental mean bond lengths of 1.43Å and 1.52Å respectively. The calculated C-C bond lengths in the phenyl rings are equivalent at 1.39Å compared to the experimental bond lengths i.e. 1.37Å. Dp presents a different case with a small variation of bond lengths in phenyl rings. Here, the calculated C-C bond length of the two sides of each phenyl ring attached to ethyl bridges is 1.40Å compared to its experimental value of 1.390 (2) A while the other four bonds are calculated to be 1.39Å compared to the experimental 1.387(3)Å bond length.

pCp-M n+ and Dp-M n+ complexes
Inclusion complexes of pCp (pCp-M n+ ) and Dp (Dp-M n+ ) with nine metal-cations of interest were optimized in their D 3 and C 3 symmetry and conrmed as either minima or transition structures through vibrational analysis. In the former case, the group 13 complexes and pCp-As 3+ are C 3 symmetric minima while group 14 and the rest of group 15 pCp-M n+ complexes have D 3 symmetry in their ground state structures. In case of group 13 pCp-M + complexes and pCp-As 3+ , the D 3 symmetry was possible for the transition structures only. The D 3 symmetric pCp-M n+ complexes show h 6 h 6 h 6 metal-phenyl rings coordination. However, pCp-Ga + and pCp-As 3+ present a case with C 3 symmetry showing h 6 h 6 h 6 coordination while its counterparts i.e. pCp-In + and pCp-Tl + exhibit h 2 h 2 h 2 coordination. As for the Dp complexes, all the nine structures were optimized as minima Table 1 Key structural features of pCp-M n+ and Dp-M n+ . Exp. shows experimentally reported pCp-Ga + , pCp-In + , pCp-Ge 2+ , pCp-Sn 2+ and pCp-As 3+ . The computational results for pCp-In + and pCp-Sn 2+ reported earlier given in parentheses. Avg. C-M, Cent-M and C]C-M denote average metal-carbon distances, distance of metal cation from the center of the cavity, and average distance between C]C of phenyl rings (upper and lower in pCp-M n+ and sideways in Dp-M n+ ) from metal ion, respectively  in C 3 as well as D 3 symmetry. Dp-M n+ complexes with C 3 show h 1 h 1 h 1 coordination of a metal ion staying on the top of Dp cavity in each complex. D 3 symmetric complexes, on the other hand, exhibit h 6 h 6 h 6 coordination of the metal-cation present inside the Dp cavity. These results suggest that despite few exceptions, the pCp and Dp complexes of main-group metalcations under discussion differ from the transition metal complexes where the metal-cation tends to come out of the cavity 20,21 preferring peripheral coordination unlike current complexes where central (h 6 h 6 h 6 ) coordination is preferred.
Some key structural parameters of pCp-M n+ and Dp-M n+ of interest are given in Table 1. The experimental evidence is available only for pCp-Ga + , pCp-In + , pCp-Ge 2+ , pCp-Sn 2+ and pCp-As 3+ where it can be seen that computational results align well with the experimental ones. However, the distance of metal ion from the center of the pCp cavity in experimental and computational instances differs, which can be attributed to the proximity of corresponding counter anion(s) to the metal cation in experimentally reported complexes. This is in line with the computational ndings reported earlier in case of pCp-Sn 2+ . 18 The difference is highly pronounced in case of pCp-As 3+ where both the calculated conformations (C 3 and D 3 ) exhibit h 6 h 6 h 6 coordination mode, as mentioned earlier, with As 3+ in C 3 symmetry displaced by 0.266Å from the center of the host cavity. On the other hand, the experimentally reported pCp-As 3+ shows h 2 h 2 h 2 coordination where AsCl 3 coordinates from the top of the cavity. In addition to pCp-Sn 2+ , the computational ndings for pCp-In + have also been reported earlier 22 which, despite small numerical differences, show a great deal of similarity with the pCp-In + calculated in our case.
The structural parameters of C 3 and D 3 symmetric conformers of pCp-Ga + are identical except the displacement of Ga + from the center of the cavity by 0.466Å in C 3 conformer in contrast to its exactly central location in the D 3 counterpart. This is supported by identical DE int of both the conformers in Table 4. However, in case of pCp-In + and pCp-Tl + , the metal ion is located signicantly further from the center in C 3 symmetry while it is exactly in the center in D 3 analogue. In both of these complexes, cavity size increases upon inclusion of metal cation for h 6 h 6 h 6 coordination as evident from internal radius C]C-M (Table 1) while it shrinks when the metal cation approaches from the top of the pCp cavity for h 2 h 2 h 2 interaction. In line with these observations, C 3 conformers of pCp-In + (DE int ¼ À51.51 kcal mol À ) and pCp-Tl + (DE int ¼ À46.78 kcal mol À ) are energetically favourable by À4.73 kcal mol À and À5.2 kcal mol À than their D 3 (DE int ¼ À46.78 kcal mol À for pCp-In + and À41.37 kcal mol À for pCp-Tl + ) counterparts. Similarly, the D 3 symmetric pCp-As 3+ is less favourable by 8.47 kcal mol À than its C 3 analogue as evidenced by DE int (Table  4).
In case of Dp-M n+ , all the C 3 conformers exhibit h 1 h 1 h 1 coordination with the metal cation located on the top of the cavity where their D 3 analogues are perfectly h 6 h 6 h 6 coordinated. A comparison of DE int (Table 5) shows that in most of the cases, D 3 conformers of Dp-M n+ are energetically favourable compared to their C 3 counterparts except Dp-In + , Dp-Tl + and Dp-As 3+ where h 1 h 1 h 1 coordination is more favourable.

Thermodynamic parameters
The trends of thermodynamic stability were established by calculating the enthalpies of reaction for the complexes under discussion given by the reaction as given in eqn (3).
Here, Lg (ligand) denotes pCp or Dp as the case may be and M n+ represents the respective metal cation with n ¼ 1, 2, 3 for group 13, 14 and 15 respectively. LgM n+ shows the resulting cyclophane-metal complex. It can be seen in Table 2 that while moving from le to right in a period, thermodynamic feasibility of pCp or Dp complexes is enhanced as indicated by increasing exothermic enthalpy from group 13 through 14 to 15 in the same period. On the other hand, there is a trend of decreasing thermodynamic feasibility while moving from lighter to heavier elements in a group. These trends can be attributed to the size of ionic radii and formal charges on metal-cations. The smaller the ionic radius (as on the top of a group and the le side of a period), the greater the thermodynamic feasibility of the corresponding complex. Ionic radii of the metal cations under discussion are given in Table 3. Moreover, the role of formal charge has been found crucial in determining the strength of interaction and stability previously in such complexes. 22 Upon moving from group 13 to 15 in a period, the formal charge in our case increases from +1 to +3 and so does the exothermic enthalpy of reaction. The trends of thermodynamic feasibility can be correlated with the results obtained from EDA and are discussed in the next section.

Bonding properties
The nature of bonding in the inclusion complexes under discussion was carried out using MZEDA technique, Bader's and NBO analyses. The former decomposes total interaction energy into various energy terms and hence it explains the strength of different interactions. NBO and Bader's analyses predict the type of bonding.

Morokuma-Ziegler energy decomposition analysis
The results for MZEDA of pCp-M n+ complexes are given in Table  4. It is to be noted that DE orb indicates the strength of covalent interaction while DE ele shows the strength of electrostatic attraction. The ratio of DE orb to DE ele explains the relative importance of covalent and ionic interactions i.e. the greater the ratio, the higher the percentage of DE orb will be. The total steric repulsion present in a complex is depicted as DE Pauli . The sum of the above-mentioned three terms accounts for the instantaneous interaction energy DE int . However, it is advised in the literature to describe the overall interaction in terms of three separate quantities; DE Pauli , DE ele and DE orb . 43 Table 4 shows that in group 13 pCp-M + complexes, interaction energy is highest on the top of group and lowest in case of pCp-Tl + . This is in line with the trends of enthalpy of reaction ( Table 2) that depicts a decrease in thermodynamic feasibility down the group. It can also be argued that increasing cationic radius (Table 3) down the group makes the p-donation from cyclophane cavity to metal less convenient. It can be seen that the C 3 symmetric pCp-Ga + has an h 6 h 6 h 6 coordination with the three aromatic rings of pCp which facilitates an overall stronger interaction. The metal-cation in pCp-In + and pCp-Tl + is located further from the centre of the cavity building an h 2 h 2 h 2 coordination in each case, thus contributing to a comparatively lower interaction energy. A deeper insight into the EDA results (Table 4) shows that D 3 symmetric h 6 h 6 h 6 transition structures of group 13 experiences a greater repulsion than their C 3 symmetric complexes as indicated by DE Pauli . This is then compensated by a greater DE prep in case of the former compared to the latter. However, both the C 3 and D 3 conformers of pCp-Ga + are h 6 h 6 h 6 with Ga + locating a little further from the center in C 3 symmetric complex while DE int is identical for both. Moreover, coordination in pCp-Ga + is facilitated by a higher percentage of orbital interaction (DE orb ) than the electrostatic interaction (DE ele ). In case of pCp-In + and pCp-Tl + , the percentage of DE orb is less compared to that in pCp-Ga + and DE ele fraction substantially increases which ultimately accounts for a decrease in strength of coordination down the group as suggested by a decreasing DE int . The strength of both the attraction terms may be expressed as the ratio of DE orb to DE ele which is 1.2 for pCp-Ga + while it is 1.0 for each of pCp-In + and pCp-Tl + . The trends of EDA results for group 14 and 15 are identical to those of group 13.
It is important to note that Table 4 includes EDA results reported earlier for pCp-In +22 and pCp-Sn 2+ (ref. 21) calculated at BP86/TZ2P+ level. Although there is difference between the numerical values of those earlier and current studies which may be attributed to the different methodologies (BP86/TZ2P+ earlier versus PBE0/TZ2P current) used in both studies, the trends are identical. For instance, the h 2 h 2 h 2 -coordinated pCp-In + is favoured compared to its h 6 h 6 h 6 analogue based on their DE int in both the cases. Similarly, DE orb is a major attractive term in h 6 h 6 h 6 -coordinated pCp-Sn 2+ in earlier and current studies as evidenced by the percentages of their DE orb and DE ele ( Table 4).
The trends of EDA in the case of the Dp complexes under discussion (Table 5) are similar to that for pCp complexes. There occurs a decrease in overall interaction energy from top to bottom in a group. However, Dp-In + and Dp-Tl + exhibit relatively stronger electrostatic interactions compared to their pCp analogues. This is evident from the ratio of DE orb to DE ele that is 0.8 for both of these complexes while the same is 1.0 for both of their pCp counterparts. All Dp complexes demonstrate the presence of strong coordination as evident from the interaction energy. However, Table 5 shows that DE int of D 3 symmetric h 6 h 6 h 6 -coordinated Dp-M n+ complexes is comparatively higher than that of their C 3 symmetric h 1 h 1 h 1 -coordinated counterparts. This difference can be correlated to the greater contribution of covalent interactions (DE orb ) compared to that of Fig. 2 Molecular graphs of C 3 symmetric pCp-In + (representing all the complexes with C 3 symmetry) and D 3 symmetric pCp-Ge 2+ and Dp-As 3+ (representing all the D 3 symmetric complexes) calculated at PBE0-B3BJ/def2TZVP. BCPs are shown as orange-coloured dots.
This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 30796-30805 | 30801 electrostatic interactions (DE ele ) in overall DE int in case of D 3 complexes. On the other hand, C 3 symmetric complexes are characterized by higher contribution of DE ele than that of DE orb in the overall DE int which causes a comparatively weaker coordination in C 3 complexes than their D 3 analogues.
It can be deduced from Tables 4 and 5 that on moving from group 13 to 15 in a period, DE int increases to a great extent from one metal ion to the next (Table 4). This trend can be attributed to the formal charge of a cation as witnessed previously. 22 In our case, the formal charge is +1, +2 and +3 on group 13, 14, and 15 metal ions whereas DE int in a period increases in the order of group 13 < group 14 < group 15. Moreover, the percentage of DE orb successively increases and that of DE ele subsequently decreases from le to right in a period. This whole discussion suggests that an increase in formal charge of a metal ion strengthens its coordination with the cyclophane host characterized by an increasing DE orb which further accounts for an increase in overall DE int along a period.

QTAIM (Bader's) analysis
The molecular graphs were extracted from the Bader's analysis results that show the bond paths for all the electron pairs that would be expected for the host molecules i.e. CC and CH bonds. Additionally, there are bond critical paths (BCPs) connecting a metal-cation to the host (Fig. 2). For h 2 h 2 h 2 complexes i.e. pCp-In + and pCp-Tl + , there is single BCP between the metal and one carbon atom of each of the aromatic rings toward the surface of the cavity that coordinates with the metal-cation as in the case of pCp-In + in Fig. 2 where the coordination is on the top of the cavity. For all the pCp and Dp D 3 symmetric complexes there are two BCPs connecting the metal from the centre of the cyclophane cavity with two carbon atoms one each on top and bottom sides of the cavity as in the case of pCp-Ge 2+ and Dp-As 3+ (Fig. 2). The graphs for C 3 complexes are similar to each other and same is the case with D 3 complexes. This suggests a similar bonding in these complexes.
The results for Bader's analysis of all the complexes are given in Table 6 which may be used to classify the types of interaction in these complexes. For that, Popelier has recently devised a mechanism based on electron density (r) and its various other functions. 44 According to the set of these rules, a small r accompanied by a negative Laplacian L (or V 2 r > 0 since L ¼ ÀV 2 r) usually denotes a depletion of electron density along a BCP and is the characteristic of either closed-shell or donoracceptor interactions while a small r along with L approaching zero indicates a shared interaction. On the other hand, a large electron density and a positive L (or V 2 r < 0) shows that electron density is concentrated along a BCP and the interaction will usually be classied as covalent or intermediate. To further conrm the nature of interaction, some additional parameters such as local energy density and the ratio of kinetic energy (G) to r are also advised by Popelier. 44 A careful classication based on Table 6 suggests that the complexes under discussion possess shared interactions i.e. electron sharing is enabled between the cyclophane hosts and the cationic guests. These ndings are in line with the results of EDA where DE orb in most of the cases is a major contributor in total attraction energy compared to DE ele . All these complexes have a small r and its laplacian (L), a G/r less than 1 and a negative value of H except pCp-Tl + that has a positive H suggesting that it has interaction between a shared and a donoracceptor interaction. The EDA results of pCp-Tl + support this assumption since DE ele has a greater percentage than DE orb .

Molecular orbital analysis
The electronic structure of the pCp and Dp complexes under consideration was further investigated based on natural population analyses (Table 7). Quantitative molecular orbital (MO) diagrams of pCp-Ga + (C 3 h 6 h 6 h 6 ), pCp-In + (C 3 h 2 h 2 h 2 ), pCp-Sb 3+ (D 3 h 6 h 6 h 6 ) and Dp-Sn 2+ (D 3 h 6 h 6 h 6 ) are given in Fig. S2 in ESI † to examine bonding and antibonding interactions between the cation and ligand fragments of a complex. The p bonding and antibonding orbitals of a cyclophane ring correspond to those of an aromatic system and have been labelled as p 1 , p 2 and p 3 . In all the h 6 h 6 h 6 complexes where the metal cation is centrally located inside the ligand cavity, p 1 orbital of the Table 6 Electron density (r) and its Laplacian (L), ratio of kinetic energy to electron density (G/r) and local energy density (H) calculated through Bader's analysis of the complexes of interest. These QTAIM parameters may be used as criteria to characterize the type of interaction in the current metal-cyclophane complexes based on cyclophane ligand coordinates with ns orbital of metal cation (with n ¼ 4, 5, 6 for 4th, 5th and 6th row of elements respectively) giving rise to fully occupied bonding and antibonding p 1 ns interactions as in Fig. S2(a) † thus not taking part in overall bonding ('n' has the value 4, 5, 6 for 4th, 5th, and 6th-row elements respectively). However, np orbital of the metal cation can be bonded to p 2 and p 3 orbitals of the aromatic system such that p z interacts with p 2 based on orientation while p 3 has an equal chance of interaction with p x and p y orbitals of the metal cation. In pCp-In + and pCp-Tl + complexes, the bonding scheme is somewhat different as both have h 2 h 2 h 2 coordination mode. We propose in these two cases that only p 2 and p 3 of the cyclophane cavity coordinate with ns orbital and one of the np sub-orbitals, respectively leaving behind p 1 without any interaction as in Fig. S2 (b). † This is also supported by the comparison of the amount of ligand to metal charge transfer (LMCT) in these two complexes (Table 7) with h 6 h 6 h 6 -coordinated pCp-Ga + where LMCT is double the amount of that in its other two counterparts.
It can be seen in Fig. S2(a) † in MO diagram of pCp-Ga + that p 1 of pCp and 4s Ga + mix to give fully occupied bonding and antibonding MOs where bonding orbital has 15.87% contribution of Ga + while the remaining 84.13% come from pCp. The contribution of Ga + in antibonding MO, however, increases to 50.90%. Next, p 2 and 4p z mix with a contribution of 6.86% and 83.79% from Ga + in the resultant bonding and antibonding MOs, respectively. Moreover, p 3 mixes equally with 4p x and 4p y (as discussed above) to give bonding (1.93% from Ga + ) and antibonding (60.62% from Ga + ) interactions. In contrast, we can observe an overall lesser contribution of 5p orbital of In + in Fig. S2 (b) compared to 4p of Ga + in (a). † This is in with a higher LMCT in case of h 6 h 6 h 6 -coordinated pCp-Ga + compared to that in h 2 h 2 h 2 -coordinated pCp-In + (Table 7). A comparison of (a), (b), (c) and (d) in Fig. S2 † suggests that the greater the contribution of p orbital of corresponding metal cation, the larger the LMCT would be.
Comparison of pCp and Dp complexes pCp and Dp have p-rich cavities and can efficiently host the main-group metal-cations making inclusion complexes with them. However, both possess different structures that cause these p-prismands to behave somewhat differently from each other. Although pCp has a rigid geometry as stated earlier, it is still somewhat exible due to three ethano-bridges linking the three aromatic rings compared to Dp which has six aliphatic chains bridging the three phenyl rings. It readily undergoes conformational changes and tends to adjust its geometry accordingly to accommodate the metal-cations more conveniently. Hence, the pCp complexes are thermodynamically more feasible compared to the Dp complexes.
The comparison is not so simple when it comes to the strength of bonding interaction based on an overall interaction energy. Although, there is no denite trend followed while comparing pCp complexes with their Dp analogues, there are instances where Dp complexes have an increased coordination strength. For example, group 13 pCp-M + complexes (Table 4) have higher interaction energy than corresponding Dp-M + complexes (Table 5) which follows the same trend as in their thermodynamic feasibility (Table 2). However, the interaction energy in case of most of the group 14 and 15 pCp-M n+ is lower than their Dp analogues. This anomaly can be attributed to the collective effects of DE orb and DE prep . As a general trend, DE orb is greater while DE prep is smaller for Dp-M n+ complexes compared to pCp-M n+ complexes.

Conclusion
Quantum chemical calculations were employed to investigate the pCp and Dp complexes of main-group metals. Geometry optimization under symmetry constraints shows that an h 6 h 6 h 6 mode of coordination is preferred in most of the cases. These complexes present excellent examples of host-guest interactions unlike the previously reported transition metal complexes of cyclophanes where a metal-cation takes peripheral position on top of the cavity. Most of the complexes are minima in D 3 symmetry whereas group 13 pCp-M + and pCp-As 3+ have only C 3 symmetric minima. MZEDA shows based on an overall interaction energy that all the complexes of interest bear strong metal-cyclophane coordination. The thermodynamic stability of pCp complexes is higher than their Dp analogues which is as expected due the more exible and "adjustable" structure of pCp compared to Dp. However, the trends of coordination strength in both the cases are mixed. For example, group 13 pCp complexes have higher interaction energy than their Dp counterparts (a trend that is in line with their thermodynamic stability) while the majority of the rest of Dp complexes have an increased coordination strength compared to the corresponding pCp complexes (thus opposing the trend of thermodynamic stability). This can be correlated with the fact that generally DE prep is lower and DE orb is higher for Dp-M n+ complexes than corresponding pCp-M n+ except group 13 complexes which result in an overall higher DE int for Dp-M n+ than pCp-M n+ of group 14 and 15. NBO analysis provided the basis for explaining in detail the electronic structure of complexes. Different parameters of Bader's analysis suggest the shared nature of M-C interactions in all these complexes except pCp-Tl + that has a donor-acceptor type of interaction.

Conflicts of interest
The authors declare that they have no conict of interest.