C70 Fullerene Cage as a Novel Catalyst for Efficient Proton Transfer Reactions between Small Molecules: A Theoretical study

When acids are supplied with an excess electron (or placed in an Ar or the more polarizable N2 matrix) in the presence of species such as NH3, the formation of ion-pairs is a likely outcome. Using density functional theory and first-principles calculations, however, we show that, without supplying an external electron or an electric field, or introducing photo-excitation and -ionization, a single molecule of HCl or HBr in the presence of a single molecule of water inside a C70 fullerene cage is susceptible to cleavage of the σ-bond of the Brønsted-Lowry acid into X− and H+ ions, with concomitant transfer of the proton along the reaction coordinate. This leads to the formation of an X−···+HOH2 (X = Cl, Br) conjugate acid-base ion-pair, similar to the structure in water of a Zundel ion. This process is unlikely to occur in other fullerene derivatives in the presence of H2O without significantly affecting the geometry of the carbon cage, suggesting that the interior of C70 is an ideal catalytic platform for proton transfer reactions and the design of related novel materials. By contrast, when a single molecule of HF is reacted with a single molecule of H2O inside the C70 cage, partial proton transfers from HF to H2O is an immediate consequence, as recently observed experimentally. The geometrical, energetic, electron density, orbital, optoelectronic and vibrational characteristics supporting these observations are presented. In contrast with the views that have been advanced in several recent studies, we show that the encaged species experiences significant non-covalent interaction with the interior of the cage. We also show that the inability of current experiments to detect many infrared active vibrational bands of the endo species in these systems is likely to be a consequence of the substantial electrostatic screening effect of the cage.


Results and Discussion
Geometries and nature of potential energy surface. Figure 1a and d show the energy-minimized geometries of isolated H 2 O···HF and H 2 O···HF@C 70 , respectively, obtained with Gaussian 09 26 ; details of the intermolecular bond distances and angles are listed in Table S1.
There are significant differences between the intermolecular geometries of these two systems. On passing from isolated H 2 O···HF to H 2 O···HF@C 70 , the O···H hydrogen bond distance is significantly contracted by 0.200 Å, the H-F covalent bond distance is elongated by 0.023 Å and the O-H bond in H 2 O is elongated by 0.005 Å. The contraction of the O···H intermolecular hydrogen bond suggests that it is strengthened, whereas the H-F and O-H bonds are weakened inside C 70 . The latter observations concur with the IUPAC recommendations for identifying hydrogen bonding 27 . The contraction of the O···H hydrogen bond leads to the immediate conclusion that the bond gains stability upon encapsulation, possibly indicative of a gain in partial covalent character. This is in agreement with experiment ( Fig. S1) and with MP2 calculations, in which, the latter shows that the O···F (2.496 Å) distance decreases by 5.5%, whereas the H-F bond distance increases by 2.1% 16 . Other studies reported similar geometrical changes for the H 2 O···H 2 O and HF···HF dimers inside C 70 28 .
As found for isolated H 2 O···HF, the isolated H 2 O···HCl and H 2 O···HBr dimers are also hydrogen bonded (Fig. 1e,f). Their encapsulation inside the C 70 cage, however, caused the bonds in HCl and HBr to stretch sufficiently such that the resulting geometries (Fig. 1b,c) are distinctly different from that of the H 2 O···HF dimer in C 70 . This has required a reorganization energy of −6.76 and −9.34 kcal mol −1 for these dimers, respectively, that is significantly larger than the energy of −0.89 kcal mol −1 for the H 2 O···HF dimer in C 70 . The relative change in the free energies associated with the geometric reorganization of these three dimers is 2.16, 6.92 and 9.73 kcal mol −1 , respectively.
The stretching of HCl and HBr bonds in C 70 is presumably driven by dissociative electron density transfer from H to the Cl and Br atoms, respectively. This has caused cleavage of the bonds in HCl and HBr, and the formation of H + and X − ion pairs in the preliminary step. Although the critical external electric field required for the dissociation of the σ-bonds in isolated HCl and HBr molecules is 510 and 462 MV cm −1 , respectively, solvation of the proton and the halide anion by water molecules is known to substantially lower the critical electric field by about 300 MV cm −1 29 . This is the likely scenario inside the cage as the HX bond breaking in the presence of the H 2 O molecule is due to the encapsulation efficiency of the C 70 cage.
The dissociation of the acids (HCl and HBr) did not enable the H + ion to migrate towards the interior carbon surface of the cage. Rather, the proton moved apart from X − along the reaction (hydrogen-bond) coordinate, a process that involves a simple proton motion between the X − and H 2 O. The bridging proton resides at an intermediate position between X − and O at 1.405 and 1.563 Å from the Cl − and Br − ions, respectively, and was found at 1.283 and 1.221 Å from O(H 2 ); the latter distances suggest the formation of the hydronium ion (H 3 O + ). Clearly, the structure of the bridging proton for these contact ion-pairs, H 2 O-H + ···Cl − and H 2 OH + ···Br − , is similar to the structure of the water Zundel ion, (H 2 )O-H-O(H 2 ) + 30 . This is facilitated by the near linearity of the double charge-assisted hydrogen bonds (∠O-H + ···X − = 172.9° for H 2 O-H + ···Cl − and 164.8° for H 2 O-H + ···Br − in C 70 ) and the short distance between the H + and X − ion pairs (r(H + ···X − ) 1.405 Å for H 2 O-H + ···Cl − and 1.563 Å for H 2 O-H + ···Br − ). For comparison, a detail of the nature of the geometry of isolated H 2 O···HX (X = F, Cl, Br) dimers in their first excited and in their anionic ground states is given in TEXT S1 of the Supplementary Information.
The proton transfer from the acids HX was made possible because the cage substantially minimized the potential barrier between HX and H 2 O (Fig. 2). The proton transfer is essentially a one-dimensional process in these systems. The chemical process involved for H 2 O-H + ···X − ion-pair formation is analogous to that previously described by Gutberlet et al. 1 . Specifically, it was shown in that study that the successive aggregation of  Table S1 for details). The smaller red spheres at the center of five-membered rings of C 70 represent the centroids.   31 .
The formation of the hydronium anion inside the C 70 cage can be understood from the difference between the proton coordinate distances, Δ. This parameter is positive when the H atom is bound to the halogen atom, and negative when the H atom is detached and moves away 32 . For instance, the HCl and HBr bond distances in the monomers are 1.294 and 1.434 Å, respectively (Fig. 1i,j), whereas those in the isolated dimers H 2 O···HCl and H 2 O···HBr are 1.337 and 1.479 Å, respectively (Fig. 1e,f). These show that isolated H 2 O···HF, H 2 O···HCl and H 2 O···HBr dimers are associated with Δ > 0 (Tables S1 and S2), meaning the halides are covalently bound to the H atoms in these dimers. The same is also true for the H 2 O···HF dimer in H 2 O···HF@C 70 (Table S1).
The H 3 O + ···Cl − and H 3 O + ···Br − ion-pairs, on the other hand, are characterized by negative Δ values (Table S2). This signifies that the halide anions are indeed separated from the H + and the latter is bound covalently to H 2 O molecule in these dimers inside C 70 .
Stabilization energies. We estimated the binding energies of the H 2 O···HF dimer in the absence and in the presence of the C 70 cage at various levels of theory, including [CCSD(T)/aug-cc-pVTZ], as was done in previous studies on other systems 13,25,28,31 . From the sign and magnitudes of the uncorrected (ΔE) and basis set corrected (ΔE(BSSE)) binding energies (Table S1), it is quite apparent that the isolated H 2 O···HF dimer is less stable (<1 kcal mol −1 at all levels of theory) than in H 2 O···HF@C 70 ; this is accompanied by the O···H bond contraction. Our zeroth level symmetry adapted perturbation theory 33 calculation suggests that the interaction energy of the isolated H 2 O···HF dimer can be described in terms of the dissected energy components arising from electrostatics (−15.56 kcal mol −1 ), exchange repulsion (+14.52 kcal mol −1 ), polarization (−6.37 kcal mol −1 ) and dispersion (−2.61 kcal mol −1 ). Its encapsulation inside the C 70 markedly increases these component energies by 1.44, 1.96, 1.96 and 1.56 times, respectively, thus suggesting that the effect of exchange repulsion is increasingly dominant and is compensated by polarization and dispersive effects.
We have also estimated the ΔE of the H 2 OH + ··· − Cl and H 2 OH + ··· − Br ion-pairs using a variety of computational approaches; the results are summarized in Table S2. In doing so, the fragments such as OH 3 + and X − (X = Cl, Br) were taken as the two fragments of the complex. As expected, very large values for ΔE and ΔE(BSSE) resulted for both the complexes, confirming the formation of ion-pairs. Such large values are expected given that complex formation is predominantly electrostatically driven, with a significant transfer of charge between the interacting partners. Binding energies of ion-pairs of similar magnitudes have been reported in many recent studies 34-41 . Nature of charge rearrangement and charge transfer. The isolated H 2 O···HF, H 2 O···HCl and H 2 O···HBr dimers have stabilization energies of −8.57, −5.18 and −3.57 kcal mol −1 with CCSD(T), respectively, indicative of the presence of weak-to-medium strength hydrogen bonded interactions 34,42 . The charge transfer for these dimers is 0.056, 0.062 and 0.064 e, respectively. The result of a second-order perturbation theory analysis of the Fock matrix in natural bond orbital basis 43 suggests this charge transfer occurs from the oxygen lone-pair orbital to a σ*(H-F) anti-bonding orbital of HX2.
Significant charge rearrangement accompanying encapsulation of the guest dimers is notable in Figs S2-S5. Depending on the charge polarity of a particular atom in the guest, a specific portion of the cage is polarized with an opposite sign so as to facilitate a Coulombic attraction between them. Because of this, the charge-transfer increases significantly with cage residence, suggesting a further transfer of charge between the monomers inside the host. This is evidenced by the increase in the partial charge of the H atom by 0.0044, 0.1080 and 0.2024 e for H 2 O···HF@C 70 , H 2 OH + ··· − Cl@C 70 and H 2 OH + ··· − Br@C 70, respectively, compared to those in the isolated dimers (Fig. S2). The net charge transfers between the C 70 host and the H 2 O···HF, H 2 O···HCl and H 2 O···HBr guests is 0.107, 0.094 and 0.102 e in the encapsulated complexes; the absence of a trend in the charge transfer values is due to the integration error. This result suggests that all the three individual dimers in C 70 are Mulliken inner type complexes 34,42 .
Arrhenius theory suggests that when the acid HA reacts with the base B(H 2 O) n this causes the formation of neutral H 2 O molecules and a charge-neutral salt, n(AH) 44 . This widely applicable concept is true for the isolated H 2 O···HF, H 2 O···HCl and H 2 O···HBr dimers that are found to be electrically neutral. However, this is not so when these dimers are enclosed inside C 70 . Our result suggests that the H 2 O···HF, H 2 O···HCl and H 2 O···HBr dimers are no longer electrically neutral inside C 70 (Figs S3-S5), in agreement with previous reports on non-fullerene-based systems (for example 45 ). Although the former is a hydrogen bonded dimer and the latter two are locally ion-pairs inside C 70 46 ; similar results have been reported for non-fullerene based compounds as well 47,48 . www.nature.com/scientificreports www.nature.com/scientificreports/ Effect of electrostatic screening on the polarity of the dimers. The formation of the ion pairs inside the C 70 cage is evidence of the unusually strong charge-assisted hydrogen bond energies (Table S2). It was therefore expected that the polarity of the ion-pair would increase significantly. This is in agreement with a previous study 23 ; the dipole moment of the H 3 N···HCl system increased from 4.15 to 9.82 D upon the intermolecular proton transfer that enabled the formation of the Cl − ···HNH 2 + ion-pair salt. However, we did not notice any such enhancement in the dimer dipole moments for either H 2 O···HF@C 70 , H 2 OH + ··· − Cl@C 70 , or H 2 OH + ··· − Br@C 70 , even though this is a common characteristic on the formation of binary complexes of Mulliken inner and outer types 34 . The anomaly could thus be attributed to the effect of the cage interior that effectively screens the dipolar electric field of an entrapped species (irrespective of whether this is an ion-pair, a simple hydrogen bonded dimer, or simply an H 2 O molecule).

Effect of encapsulation on the emergence of significant intermolecular interactions: A QTAIM rationalization. A fundamental question in endohedral systems is whether the guest species is innocent
inside the cage interior, or whether it experiences significant intermolecular interactions 49 . In several studies reported over the last seven years 6,11,16 , it was shown that a guest species, such as H 2 O, residing near the center of the cage, is not involved either in hydrogen bonding or in π···H interactions with the interior carbons of the five-and six-membered rings of the C 60 and C 70 cages. We examined the bond critical point (bcp) and bond path topologies of potential bonding interactions ( Fig. 3) using Bader's quantum theory of atoms in molecules (QTAIM) 50,51 . We found that the charge density, ρ b , at the O···H bond critical point increased from 0.047 to 0.078 a.u. on passing from H 2 O···HF to H 2 O···HF@C 70 , with a concomitant change in the Laplacian of the charge density (∇ 2 ρ b ) of +0.143 to +0.179 a.u. These signatures are consistent with the contraction of the O···H bond in H 2 O···HF upon its encapsulation, suggesting that there is no H 3 O + ···F − ion-pair formed inside the C 70 cage, as claimed previously 16 .
In H 2 O···HF@C 70 , the H atoms of the H 2 O molecule are found not to be inert. Each H atom in H 2 O is engaged in a C π ···H-O hydrogen bonding interaction with an adjacent carbon of the C 70 cage (Fig. 3). The ρ b , ∇ 2 ρ b , and the total energy density (H b ) at the C π ···H-O bcps are calculated to be 0.012, +0.044 and +0.0021 a.u., respectively. In addition, the lone-pair bonding orbital of the O atom is observed to act as a donor of electron density, which is evidenced by the formation of O···C interactions for which ρ b , ∇ 2 ρ b and H b are approximately 0.0083, +0.031, and +0.0010 a.u., respectively.
There is a similar involvement of the H atoms of H 2 O with the carbon cage in H 2 OH + ··· − Cl@C 70 and H 2 OH + ··· − Br@C 70 . For these two cases, the C π ···H-O interactions are somewhat stronger and are very asymmetric. This can be seen from the values of ρ b , ∇ 2 ρ b and H b , in which, the two C π ···H-O hydrogen bonds in H 2 OH + ··· − Cl@C 70 have ρ b = 0.022 and 0.019 a.u., ∇ 2 ρ b = 0.063 and 0.0717 a.u., and H b = +0.0017 and +0.0024, respectively. The corresponding values for H 2 OH + ··· − Br@C 70 are 0.020 and 0.025, +0.067 and +0.079, +0.0023 and +0.0010 a.u., respectively. Likewise, the O···C interactions with the cage in H 2 OH + ··· − Cl@C 70 and H 2 OH + ··· − Br@C 70 are also found to be strengthened.
In each of these systems, the halogen is involved in a five-center, six bond topology, forming a cone of intermolecular interactions reminiscent of that found for cyclopentadienyl and substituted dienyl bonding to titanium 52 . Such a topology is of significance in the understanding of non-classical metal-to-saturated-carbon atom interactions. Even though the fluorine in H-F resides at a shorter distance from the interior carbon wall of the C 70 cage compared to the chlorine and bromine atoms in H-X (X = Cl, Br), it is comparatively more weakly bound to the carbons of the five-membered aromatic ring which it faces. The attractive engagement between the F in the H-F fragment and the C 5 ring of the cage can be recognized as π···σ hole (negative) interactions 53 , whereas for the other two this could be standard π···X − (X = Cl, Br) interactions. The strength of the C π ···X interactions is gauged from the values of ρ b : 0.010 a.u. for C···F; 0.014 for C···Cl; and 0.015 a.u. for C···Br. The ∇ 2 ρ b values for the corresponding interactions are +0.048, +0.0053 and +0.055 a.u., respectively, indicative of closed-shell interactions.

Nature of intermolecular interactions: A reduced density gradient and independent gradient model perspective.
Reduced Density Gradient Non-Covalent Index (RDG-NCI) 54 and Independent Gradient Model (IGM) 55,56 analyses suggest that regardless of the nature of the guest species, they are involved in an attractive engagement with the interior carbons of the C 70 cage. This is true not only for the isolated guest molecules (viz. H 2 O and HX (X = F, Cl, Br)), but also for the dimers formed by them. Figure 4 manifests this for the three endohedral systems, revealed using the sign(λ 2 ) × ρ vs. RDG plots of the RGD-NCI analysis (λ 2 is the second eigenvalue of the Hessian second derivative charge density matrix). There are several spikes in each of these plots. The green (for H 2 O···HF@C 70 ) and bluish-green spikes (for HOH + ··· − Cl@C 70 and HOH + ··· − Br@C 70 ) are prominent in the sign(λ 2 ) < 0 regions, and the red ones in the sign(λ 2 ) > 0 regions. The first represents attraction between the host and guest species and the second corresponds to repulsion in the RDG picture. Although the nature of the (green-and bluish-green) spikes are consistent with the level of attraction observed between the host and the guest species, it is largest for HOH + ··· − Br@C 70 compared to the other two (i.e., in the order H 2 O···HF@ C 70 < HOH + ··· − Cl@C 70 < HOH + ··· − Br@C 70 ). Figure 5  www.nature.com/scientificreports www.nature.com/scientificreports/ monomers are characterized by a single blue spike in the range −0.05 a.u. < sign(λ 2 ) × ρ < −0.045 a.u. for all isolated dimers, which are fingerprints of attractive interactions. The nature of the charge density and the location of the spikes reflected in these plots are marginally different from one another, consistent with the trend in the stability of the hydrogen bonds found in the dimers (H 2 O···HF > H 2 O···HCl > H 2 O···HBr).
The encapsulation of H 2 O···HF led the equilibrium O···H distance to contract, which is consistent with the enhanced electron density in the critical bonding region. Hence the peak position of the RDG spike shifted significantly to a high-density region, −0.07 a.u. < sign(λ 2 ) × ρ < −0.08 a.u. (Fig. 5a). The transfer of the proton from the HX moiety to the H 2 O molecule and subsequent formation of the hydronium ion in the equilibrium geometries of the H 2 OH + ··· − Cl@C 70 and H 2 OH + ··· − Br@C 70 systems explains why there are two RGD spikes observed in Fig. 5b and c, respectively. These provide evidence that the O-H + bond is relatively stronger than the H + ··· − X (X = Cl, Br) interactions in these two complexes. Nevertheless, the H + ··· − X interactions are not only predominantly electrostatic, but also contain a non-negligible amount of covalency 57 . Figure 6 compares the IGM results for the isolated and encaged systems. This innovative approach automates the identification of intermolecular interactions using actual and pro-molecular electron densities. A descriptor of the model, called δg inter/intra (inter-and intra-molecular contributions), reveals interactions beyond the presence of the bond critical points. As can be seen from a)-c) and d)-f), the isosurfaces (blue, green and mixed color) are centered around the bond critical point regions in the 3D maps of all the isolated and endohedral systems, signifying stabilizing interactions between interacting basins. The broad black spikes of the sign(λ 2 ) × ρ vs. δg inter/intra plots are associated with covalent bonds, appearing in the sign(λ 2 ) × ρ < 0 region and corresponding to high δg intra values. The hydrogen bonds and van der Waals interactions correspond to the spikes at low-density that are www.nature.com/scientificreports www.nature.com/scientificreports/ associated with low δg inter values. These are represented by two small and sharp spikes (red); one is due to attraction and corresponds to λ 2 < 0 and the other, due to repulsion, corresponds to λ 2 > 0 in reduced gradient density picture. The sign(λ 2 ) × ρ vs. RDG plots for H 2 O···HF@C 70 , H 2 OH + ··· − Cl@C 70 and H 2 OH + ··· − Br@C 70 illustrated in d)-f) of Fig. 6 corroborates the inferences from the RDG analysis (Figs 5 and 6). However, the isosurface plots suggest that attraction between the guest and the fullerene C 70 cage is significantly dispersed, and is in excellent agreement with the QTAIM description of bond path topologies that suggest many-fold interactions between the host and the guest in each of the endohedral systems. Figs S6-S8 provide detailed insight into the nature of these interactions evaluated using the X-ray crystal geometries of the H 2 O@C 70 and H 2 O···HF@C 70 systems without and with the presence of the octaethylporphyrinato-nickel(II) moiety.
The C 70 molecule is represented as an ellipsoid in 3D; the equilibrium position of the H 2 O···HX dimers was near at the center of the cage, oriented along the major axis. This positioning of the dimer is observed experimentally for H 2 O···HF@C 70 16 ; our calculations for this and for the H 2 O···HCl and H 2 O···HBr dimers inside the cage concurs with the experimental observation. However, our geometry optimization of the endohedral system with the dimer along the minor axis forced the dimer to rotate back to a configuration illustrated in Fig. 1, which must be due to the significant steric constraints of the cage. The results of QTAIM, RDG-NCI and IGM analyses are in good agreement with this view. It is quite unambiguous from Figs 4-6 that the dimer in its alignment along the major axis maximizes intermolecular interactions between it and interior of the host. The results of our binding energy calculation (assuming the entire guest and the C 70 as two monomers of each endohedral system) are in line with this (see Text S2). This also explains why the positional exchange of the H 2 O and HF was not detected experimentally at room temperature 16 . The 13 C NMR spectra  www.nature.com/scientificreports www.nature.com/scientificreports/ of H 2 O···HF@C 70 displayed nine signals due to its C 5v symmetry; this is in contrast with what could be expected for the averaged D 5h symmetry of both H 2 O@C 70 and HF@C 70 if dynamic motion of the encapsulated species occurred 16 .

Effect of encapsulation on the frontier orbital energies, ionization potential and electron affinity.
We calculated the Kohn-Sham HOMO-LUMO gap for C 70 to be 1.71 eV, in excellent agreement with the DFT-Generalized Gradient Approximation value of 1.714 eV reported previously 58 . From the Tauc plot, the optical energy gap of the C 70 thin film was reported to be 1.66 eV, exhibiting a semiconductor-like behavior in optical absorption properties, notwithstanding its molecular character [59][60][61][62] . The encapsulation of H 2 O···HF, H 2 O···HCl, and H 2 O···HBr inside the C 70 cage reduces the HOMO-LUMO gap, with values of 1.69, 1.65, and 1.64 eV, respectively; the result is qualitatively in agreement with that obtained using the well-known functional B3LYP (Table S3). This change is predominantly due to the energy of the LUMO level that becomes more negative 63 . This is not unexpected since the C-C bonds of the C 70 cage are affected by various H···C, O···C and C···X (X = F, Cl, Br) intermolecular interactions with the guest species. This geometric modification increases the overlap between molecular states and reduces the gap between the frontier orbitals. The decrease of the HOMO-LUMO gap due to encapsulation is consistent with previous findings 28,58,64 . It is also consistent with the marginal decrease in the vertical ionization potential and fundamental and optical gaps (lowest singlet excitation energies), and a negligible increase in the vertical electron affinity of the encaged C 70 systems compared to isolated C 70 (Table S4). Figure 7 illustrates the frontier orbitals of C 70 and H 2 O···HF@C 70 . Both the HOMO and LUMO are localized on the skeletal framework of isolated C 70 . Encapsulation of H 2 O···HF slightly shifts the energy levels of HOMO and LUMO, reducing the energy gap between them. However, this does not lead to a change in the nature of both these frontier orbitals, as evidenced in the density of states (DOS) spectra (Fig. 7a). In other words, our analysis of the atom-projected DOS spectra suggests that the encapsulated species does not contribute to the frontier orbitals (see Fig. 7c for H 2 O···HF@C 70 , as example). This is also the case for the other close-lying orbitals (HOMO-1, HOMO-2, LUMO + 1 and LUMO + 2, etc.) of the complexes; they are completely derived from orbitals states of the C 70 only, and hence the charge transfers between the cage and the guest species is negligibly small. This might explain why the experimentally observed electronic transition properties of the H 2 O···HF@C 70 and C 70 are very www.nature.com/scientificreports www.nature.com/scientificreports/ similar, with only marginal changes, compared to isolated C 70 16 . That the guest species does not contribute to the frontier orbitals is analogous with the organic-inorganic methylammonium lead triiodide (MAPbI 3 ) semiconductor. In this, the organic cation MA, entrapped inside the inorganic cage, does not contribute either to the conduction band minimum (CBM) or to the valence band maximum (VBM) that are built from significant contributions from s and p orbitals of the Pb and I atoms. Yet the cation plays a vital role in providing geometrical stability to the MAPbI 3 system, as well in reducing the bandgap of the material compared to PbI 2 (1.6 vs 2.3 eV).

Effect of encapsulation on vibrational frequencies and intensities. The elongation of the HF,
HCl and HBr bond distances upon the formation of isolated H 2 O···HF, H 2 O···HCl and H 2 O···HBr complexes ( Fig. 1) are accompanied by a significant red-shift in their fundamental vibrational frequencies and a concomitant increase in the intensities of the corresponding IR bands. The red-shift of the HX bond stretching frequency is 457, 506 and 476 cm −1 for HF, HCl and HBr, respectively. The increase in the intensity of the corresponding bands is approximately 11, 42 and 217 times larger than the corresponding values of 75, 31 and 5 km mol −1 for the isolated HX molecules, respectively.
Encapsulation of H 2 O···HCl and H 2 O···HBr inside the C 70 cages further increases the HCl and HBr bond lengths because of proton transfer (vide supra). As is seen in the simulated IR spectra shown in Fig. 8, there is no HX (X = Cl, Br) normal mode vibration in the spectra. Rather, we observe an (H 2 )O-H + stretching vibration for these two endohedral systems that occur at frequencies of 1720 and 1584 cm −1 , respectively, with the intensities of the corresponding IR bands of 134.6 and 106.8 km mol −1 , respectively; these are indeed measurable.
For H 2 O···HF@C 60 , however, the proton transfer is not as facile. Whereas the HF stretching vibration is centered at 3102 cm −1 , the intensity of the corresponding IR band is 45.3 km mol −1 . Compared to the intensities of the isolated HF molecule (75 km mol −1 ) and the isolated H 2 O···HF dimer (810 km mol −1 ), the intensity of the IR band in H 2 O···HF@C 60 is appreciably diminished; this does not explain the apparent correlation between intensity reduction and elongation of the H-F bond. The redshift of the band compared to the isolated H 2 O···HF and HF systems is 446 and 903.6 cm −1 , respectively.
The H 2 O species in these three endohedral systems offer very interesting IR characteristics. The O-H symmetric and asymmetric vibrations of the H 2 O coordinated to HX species without the cage (i.e., in the isolated dimers) occur between 3690 and 3805 cm −1 . Within the cage, the frequencies of these vibrations are red-shifted and occur at 3500 and 3720 cm −1 (Fig. 8). According to the literature 65 , the intensities of the corresponding vibrational modes should increase. However, we found them to decrease significantly. For instance, the respective intensities of these two IR bands for isolated H 2 O···HF were 11.6 and 58.5 km mol −1 , but 2.0 and 6.5 km mol −1 for H 2 O···HF@C 70. For H 2 O···HCl, they were respectively 19.3 and 62.0 km mol −1 , but 9.8 and 13.7 km mol −1 for H 2 OH + ··· − Cl@C 70 . Similarly, for H 2 O···HBr and H 2 OH + ··· − Br@C 70 , the intensities decreased from 21.0 and 58.9 km mol −1 to 11.9 and 16.5 km mol −1 .
Because of this apparent anomaly in the intensity profile, analytical second derivative calculations were performed on the endohedral geometries of all the three guest dimers by removing the C 70 cage. Positive frequencies were found in all cases, indicative of these being in meta-stable states. The simulated IR spectra are shown in Fig. S9. Although the frequency centers in the simulated spectra were somewhat shifted because of the removal of the cage, the intensities of most of the vibrational bands were found to be markedly larger than those shown in Fig. 8. For instance, the HF stretching frequency (intensity) in [H 2 O···HF] occurs at 3197 cm −1 (1273 km mol −1 ), whereas that of the O-H + band occurs at 1894 cm −1 (3175.1 km mol −1 ) and 1754 cm −1 (2992 km mol −1 ) for [H 2 OH + ··· − Cl] and [H 2 OH + ··· − Br], respectively (Fig. S9). This provides unequivocal proof that the reduction in the band intensity is a consequence of the profound electrostatic shielding effect of the C 70 cage. It explains why in many previous experimental studies the stretching and bending vibrational bands of HF and H 2 O, as well as www.nature.com/scientificreports www.nature.com/scientificreports/ those associated with the intermolecular interactions of the encaged species 6 , were not detected. Specifically, the IR bands associated with the HO-H and F-H bonds and those between H 2 O and the cage in H 2 O···HF@C 70 16 , as well as the corresponding vibrational bands of H 2 O in H 2 O@C 60 6 and in H 2 O@C 70 11 , could not be identified in experimental studies.
Common in both the isolated and complexed C 70 systems are the normal modes of the C-C stretching vibrations, which have often been identified experimentally and used to demonstrate that the encaged species residing near the center of the cage is inert 6,7,11,16 . Our calculations suggest that there are three such reasonably strong vibrations in isolated C 70 . Two of them are nearly degenerate and appear at a frequency (intensity) of 1436 cm −1 (81 km mol −1 ), and the other at 1470 cm −1 (9.8 km mol −1 ). These are in excellent agreement with those of the experimentally measured values of 1430 and 1459 cm −1 , with the first weaker in intensity 16 . The first mode for H 2 O···HF@C 70 was reasonably strong and was observed at 1427 cm −1 ; the latter was only observed as a shoulder (no value reported) 16 . In the simulated complex geometries (Fig. 7), these modes showed up at frequencies (intensities) of 1436 and 1462 cm −1 (34.0 and 8.4 km mol −1 ) for H 2 O···HF@C 70 , respectively; evidently, DFT predicts these spectra very accurately. The corresponding values for H 2 OH + ··· − Cl@C 70 were 1425 and 1447 cm −1 (27 and 6.4 km mol −1 ), and 1419 and 1440 cm −1 (26.4 and 2.0 km mol −1 ) for H 2 OH + ··· − Br@C 70 . Because of only a marginal change in the C-C vibrational spectra between the isolated and complexed C 70 , it was misleadingly concluded that the guest species is inert inside the cage for H 2 O···HF@C 70 -a conclusion that was drawn for several other such endohedral systems as well, including H 2 O@C 60 6,7,11 . Effect of encapsulation on the volume of C 70 . The centroid-to-centroid distance between the two five-membered rings of the C 70 cage placed at the opposite extremes of the major axis (Fig. 1)   www.nature.com/scientificreports www.nature.com/scientificreports/ respectively. Irrespective of the nature of the chemical interaction between the encaged species and the cage, the volume of H 2 O···HX(X = F, Cl, Br), compared to the volume of the isolated species, is reduced on encapsulation. Our calculations show that the reduction in the (0.001 a.u. isodensity envelope) volume of the guest species is 51. 32 70 16 , in which, the study did not involve any real determination of the volume of either the free or encaged C 70 ; yet it was speculated that the contraction of the O···H bond in H 2 O···HF is a result of compression of the H 2 O and HF inside the limited space of the cage environment.
Effect of encapsulation on the polarizability and dielectric properties of C 70 . Polarizability α is another landmark property of fullerene chemistry since fullerenes are highly polarizable molecules 67 . It is calculated to be 637.51 a.u. for free C 70 . Encapsulation of H 2 O···HF, H 2 O···HCl and H 2 O···HBr marginally increases α by 2.92, 6.12 and 7.31 a.u., respectively, which describes the change (Δα(guest@C 70 ) = α(guest@C 70 ) − α(C 70 )) in the overall polarizability upon the addition of the dimer. Considering the additive scheme, we found that the deviation Δα (Δα(guest@C 70 ) = α(guest@C 70 ) − (α(guest) + α(C 70 ))) to be negative: −7.32 a.u. for H 2 O···HF@ C 70 , −12.24 a.u. for H 2 OH + ··· − Cl@C 70 and −16.16 a.u. for H 2 OH + ··· − Br@C 70 . This suggests a decrease in polarizability (i.e., compression of electron clouds), analogous to what was observed for noble gas endofullerenes with C 36 ,C 50 , and C 60 68 . It is a result of the prevailing pressure of carbon cage that compresses the guest, resulting in a decrease in the total polarizability 69 .
The change α affects the dielectric constant ε of C 70 , a property that is often used to characterize the optical efficiency of semiconducting materials. Within the framework of the well-known Clausius-Mossotti model 70,71 , ε can be expressed as: ε = + −  (Table S7). This suggests a marginal decrease in the dielectric constant of the encapsulated species. Previous studies demonstrate that provided the dipole cannot reorient fast enough in composites, this causes the dielectric constant to decrease 72 . Others propose that the low dielectric constant of organic photovoltaics assists the exciton to present at larger distances 73 .

Conclusion
Using DFT and first-principles investigations we have demonstrated the reactivity of three fundamentally important complex systems, H 2 O···HX (X = F, Cl, Br), inside the fullerene C 70 cage. The geometries and vibrational characteristics of the H 2 O···HF system obtained at the PBE level of theory adequately explain the experimentally-reported features. The PBE stabilization energies for the dimers were shown to be comparable with those evaluated with CCSD(T). At the same time, we have demonstrated using geometrical, electron density and vibrational characteristics that the C 70 cage interior serves as a super-catalyst for HCl and HBr bond dissociation, enhancing the acidity of these acids by bond cleavage and assisting in complete proton transfers to H 2 O. This leads to the formation of OH 3 + X − (X = Br, Cl) ion pair complexes within the cage. This was not the case for the H 2 O···HF dimer inside C 70 , however. This is certainly the result of the significant basicity of F − compared to Cl − and Br − .
We summarize our observations as follows.
1. The C 70 cage interior is not hydrophobic and the guest species are not inert, as often contended 6,7,11,16 . A single water molecule when encaged will experience significant hydrogen bonding; this occurs regardless of whether we are dealing with an isolated water dimer 11 or an H 2 O···HX dimer. This observation is in line with others in which the effective pair interaction is not hydrophobic, yet the solvation properties are; hence fullerenes serve as an example in which hydrophobic interaction cannot be deducted from hydrophobic solvation 74 . 2. The C 70 cage interior has a tendency to donate charge density to the guest species. 3. Because the C 70 cage interior is polarizable, its effect plays a vital role that is largely responsible for the development of many anomalous features that were undetectable experimentally in previous studies (e.g., IR characteristics of intermolecular interactions). 4. Because of its nature as stated in 2 and 3, C 70 prefers to serve as a cationic host upon encapsulation of two guest molecules. 5. It has the potential to provide a terrace to the guest species for facilitating efficient proton transfer reactions between them. 6. The C 70 cage interior has the ability to screen the electrostatic (dipolar) field of the guest species, hence limiting the observability of many vibrational bands that are IR active. This explains why experiments often fail to show the IR spectra of the entrapped species inside the host. 7. Charge rearrangement, bond polarization and ion-pair formation are likely consequences of an accommodation of a host species, especially for dimer molecules.
The C 70 cage interior provides an elegant and innovative terrace to electrogenerate reactive species between H 2 O and HX. This is fully consistent with our similar investigations for other species (such as F 2 and mixed dihalogen derivatives), as well as other systems such as NH 3 and HX (X = F, Br, Cl), in the presence of H 2 O, that merit investigation to delineate the novelty of the catalytic profile of the C 70 cage interior. In fact, our preliminary results, to be reported elsewhere, show that F 2 inside the C 70 cage can be completely dissociated into two F − ions (2019) 9:10650 | https://doi.org/10.1038/s41598-019-46725-4 www.nature.com/scientificreports www.nature.com/scientificreports/ in presence of a single molecule of H 2 O. The mechanistic details involved in this and other similar systems will surely uncover the novel physical chemistry and catalytic detail of these materials.

Materials and Methods
The Gaussview 05 package 75 was used for generating the geometries of isolated molecules (HF, HCl, HBr, H 2 O and C 70 ), and their dimers (H 2 O···HX (X = F, Cl, Br) and (H 2 O···HX@C 70 (X = F, Cl, Br)). The endohedral complexes (H 2 O···HX@C 70 (X = F, Cl, Br)) were generated using the same procedure. The Cartesian coordinates of these systems were used for energy-minimizations using the Gaussian 09 26 package. For reasons discussed in the Results and Discussion section, the PBE functional implanted in Gaussian 09, together with the 6-311G(d,p) basis set, was used. The reliability of the functional in obtaining the experimental geometry and vibrational spectra of the H 2 O···HF@C 70 system is discussed in both Fig. S1 and the Results and Discussion section. The same method was used for calculations to obtain the eigenvalues of the Hessian second derivative matrix; in all cases, positive eigenvalues were obtained. This confirmed that the various monomer, dimer and endohedral geometries discussed in this study are all local minima. Tight convergence and ultrafine integration grids were used.
The potential energy surface (PES) rigid scans were performed on the isolated and C 70 complexed H 2 O···HCl (and H 2 O···HCl) geometries with MP2/aug-cc-pVTZ to show the effect of C 70 cage on the nature of proton transfer from the HCl (and HBr) species to the H 2 O molecule. The PBE energy-minimized geometries were used.
The binding energies associated with the H 2 O···HX dimers were evaluated at various levels of theory (PBE, MP2 and CCSD(T)) in conjunction with the aug-cc-pVTZ basis set; the PBE/6-311G(d,p) optimized geometries were used. To demonstrate the qualitative (and quantitative) reliability of the results of the PBE function, other functionals such as M06-2X, M06-2X-D3, ωB97XD, B97D3, PBE-D3, and PBE0-D3 (as implemented in Gaussian 09) were invoked to estimate the energies of the H 2 O···HF dimer and compared. Similar calculations for H 2 O···HF@C 70 , H 2 OH + ··· − Cl@C 70 and H 3 OH + ··· − Br@C 70 were performed with B97D3 to estimate the binding energy between the host and the guest species.
The charge density (ρ b ), the Laplacian of the charge density (∇ 2 ρ b ) and the total energy density (H b ) at bond critical points and the delocalization indices between various atom-atom pairs were evaluated within the framework of the Quantum Theory of Atoms in Molecules (QTAIM) 50 using the AIMAll package 51 . The atomic volumes v i are calculated within the framework of this theory, which were used to calculate the total volume of each of the monomers, dimers and endohedral systems using ν ν = ∑ = i N i 1 , where the sum is over all atoms in a given molecule. The mean static dipolar polarizability α of each of these systems was calculated using α α α α = + + ( ) xx yy zz 1 3 , where the α values are the three principal (diagonal) components of the 3 × 3 matrix given below. The Independent Model (IGM) 55 and Reduced Density Gradient Noncovalent Interaction (RDG-NCI) 54 analyses were performed to examine the nature of intermolecular interactions in the isolated dimers and between it and the cage interior. Depending on the evaluation of specific electron density properties and visualization, software packages such as Multiwfn 76 , AIMAll 51 and VMD 77 , as well as in-house codes, were used. The HOMO, LUMO and the Kohn-Sham gap (LUMO-HOMO) energies for the isolated and complexed C 70 systems obtained with B3LYP/6-31G* were compared with those evaluated with PBE/6-311G**; the geometries obtained for these systems at the latter level of theory was used. Further calculations were performed at these levels of theory to examine the nature of the total and atom-projected density of states spectra, vital for a fundamental understanding of the role the guest species plays in building the frontier orbitals. The Gausssum 03 and Origin 2018 78,79 packages were used for this purpose.
The geometries of the H 2 O···HF, H 2 O···HCl and H 2 O···HBr isolated dimers were also optimized in the first excited state using time-dependent density functional theory (TD-DFT) with PBE/6-311G(d,p). Similar calculations were performed for these systems by providing an additional electron. The goal for these calculations was to examine the extent of proton transfer feature in the excited and anionic states, and to compare this with that revealed when these dimers were encaged inside the C 70 in their electronic ground states.