THE STABILIZING ROLE OF HALIDE IONS IN ENDOHEDRAL [20

chlorines play a more relevant role into the guest encapsulation than those attached to ... for the combination of one-dimensional transport and silic...
0 downloads 3 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

C: Physical Processes in Nanomaterials and Nanostructures

THE STABILIZING ROLE OF HALIDE IONS IN ENDOHEDRAL [20]SILAFULLERANES: INSIGHTS FROM DFT CALCULATIONS TOWARDS SILICON NANOCAGES Miguel Ponce-Vargas, and Alvaro Muñoz-Castro J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 18 May 2018 Downloaded from http://pubs.acs.org on May 18, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

THE STABILIZING ROLE OF HALIDE IONS IN ENDOHEDRAL [20]SILAFULLERANES: INSIGHTS FROM DFT CALCULATIONS TOWARDS SILICON NANOCAGES

Miguel Ponce-Vargasa* and Alvaro Muñoz-Castrob* a

Université Paris-Est, Laboratoire Modélisation et Simulation, Multi Echelle, MSME UMR

8208 CNRS, 5 Bd Descartes, 77454 Marne-la-Vallée, France. E-mail: [email protected] b

Grupo de Química Inorgánica y Materiales Moleculares, Facultad de Ingeniería, Universidad

Autónoma de Chile, El Llano Subercaseaux 2801, Santiago, Chile. E-mail: [email protected]

ABSTRACT

As miniaturization of electronic devices is rapidly approaching the nanoscale, a deeper understanding of the electronic and structural properties of silicon nanoclusters, called to be the

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 26

next-generation materials for circuit design, becomes of paramount importance. Herein, a detailed DFT study of the binding forces between [20]Silafullerenes frameworks and their central halide ions is conducted, prompted by the recent synthesis of the first discrete Si20 dodecahedra stabilized by an endohedral chloride and valence saturation, [Cl@Si32Cl44]-, as well as the fabrication of the first electron transistor device based on a single silicon cluster. Although more intense stabilizing forces are obtained in the chloride-containing system, a small energetic difference with respect to bromide-centered one is found (4.76 kcal mol-1) suggesting the synthetically accessibility of the latter. An energy decomposition analysis is conducted revealing that in all cases the electrostatic term is the major contributor to the binding forces, representing about of 71%. Additionally, the higher order electrostatic terms become more relevant as the halide volume increases, and this effect is quantified through a local multipole analysis. This methodology also enables us to state that those silicon atoms directly linked to peripheral chlorines play a more relevant role into the guest encapsulation than those attached to trichlorosilyl. It is also evidenced, that the presence of such peripheral groups deeply influences the charge of the inner cluster cavity making it more positive and suitable for the encapsulation of the halide anions. We expect that this research will allow a better understanding of the driven forces of these novel structures, also contributing to experimental teams searching for novel building blocks for nanoscale transistors.

ACS Paragon Plus Environment

2

Page 3 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

INTRODUCTION

The continuous search for miniaturization in microelectronics industry has led to the development of new silicon structures able to be integrated in smaller and more power-efficient devices.1-3 In this context, silicon clusters have been extensively studied by multiple experimental and theoretical teams searching for stable frameworks at the nanoscale size.4-7 Recently, Song and coworkers have reported the fabrication of a single electron transistor device based on a silicon cluster connected to a gold break junction with a nanometer scale separation.8 Also, a field effect transistor based on undoped silicon nanowires have been successfully fabricated opening up new possibilities for the combination of one-dimensional transport and silicon technology,9 thus, positioning silicon nanostructures as promising molecular devices materials for practical applications. Unlike carbon, sp2 hybridization is unfavorable for silicon10, leading to asymmetric and reactive clusters11, however, there exists a wide consensus with respect to the main structural modifications that could contribute to the stability of the silicon assemblies, i.e., silica coating12, hydrogen capping13, phosphorus14, carbon15, and metal doping16-19, as well as the incorporation of a central metal20-22 or halide ion.23 In the case of the silicon counterpart of graphene, silicene, the functionalization with OH radicals produces the largest structural modifications due to the presence of intralayer hydrogen bonds leading to further stability.24 The recent synthesis of the first [20]Silafullerane, by Wagner and coworkers25, represents a major milestone in the field, where the key features leading to a stable cluster are the appended trichlorosilyl and chlorine groups, the endohedral chloride ion, and the overall Si-sp3 backbone. However, the nature of the related forces ensuring the encapsulation of the endohedral anion and the silicon cage has not been studied in detail.26

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 26

Herein we provide a systematic theoretical study of the nature of the cluster-anion interactions through reliable density functional theory (DFT) tools that will allow a better description of this recently synthesized [Cl@Si32Cl44]- structure. In addition, we extend our scope to similar systems containing bromide and iodide guests, also evaluating the effect of the external trichlorosilyl groups and chlorine atoms in the inner cavity charge. First, we conduct an NCI analysis27,28 in order to visualize the regions where more intense noncovalent forces arises, next we carry out an Energy Decomposition Analysis29-30 to clearly established the nature of the forces between the central ions and the silicon atoms. In parallel we carried out a Multipole Analysis31 to quantify the electronic cloud distortion experienced by the silicon atoms with the incorporation of the halide, also plotting the dipole vectors and quadrupole tensors to have a complete picture of the electronic effects caused by the endohedral halide. Finally, the 29Si-NMR properties were evaluated, in order to provide valuable and reliable information to experimental teams devoted to the synthesis of these fascinating structures.

ACS Paragon Plus Environment

4

Page 5 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

COMPUTATIONAL DETAILS

Relativistic density functional theory32 calculations were carried out by using the ADF code33, via the scalar ZORA Hamiltonian. An all electron triple-ξ Slater basis set including a polarization function (STO-TZP) was employed within the generalized gradient approximation (GGA) of Becke and Perdew (BP86).34-35 The choice of this functional is dictated by its good performance in the study of vibrational spectra of neutral36 and metal-doped silicon clusters.7,37,38 Moreover, the good performance of BP86 has also been tested for the calculation of the vertical dettachment energy of Gallium clusters39 and the optical spectra of icosahedral silver-nickel nanoclusters.40 Geometry optimizations were conducted via the analytical energy gradient method implemented by Verluis and Ziegler.41 To properly account for weak London forces, the D3 Grimme dispersion correction was incorporated.42 In D3, the dispersion coefficients are geometry-dependent as they are adjusted on the basis of local geometry. Calculations of the nuclear shielding tensors were conducted within the GIAO formalism, employing the OPBE functional along with an all-electron STO-TZ2P basis set, which were referenced to trimethylsilane (TMS) in order to account for the 29Si-NMR chemical shifts.43,44 A Noncovalent Interactions Analysis (NCI), as proposed by Yang and coworkers, was conducted by using the NCIPLOT-3.027 and NCImilano28 codes, the later using the electron density from ADF calculations. Furthermore, an Energy Decomposition Analysis was carried out following the Morokuma-Ziegler scheme.29,30 Also, the ion-multipole contribution in the studied complexes has been estimated according to a methodology based on a comparison with the hypothetical noble gas centered counterparts45,46 which has been employed in anion-host interactions studies.47,48 Finally, the Atomic Multipole Moment Analysis proposed by Swart and coworkers31, was used to study in deep the local electronic cloud variation at the silicon atoms.

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

For all isosurfaces and figures the software packages Chemcraft49 and Visual Molecular Dynamics (VMD)50 were employed.

RESULTS AND DISCUSSION

The synthesis of [Cl@Si32Cl44]-, has been recently reported by Wagner and coworkers25 and it involves a chloride-induced disproportionation reaction of hexachlorodisilane, where the chloride ion plays a structure-determining role. Some relevant structural parameters concerning [Cl@Si32Cl44]- and the related hypothetical frameworks [Br@Si32Cl44]- and [I@Si32Cl44]- are presented in Table 1 and Figure 1. The backbone silicon atoms directly linked to chloride are labelled as Si1, those bonded to trichlorosilyl as Si2, and the trichlorosilyl ones as Si3 (see Figure 1). Remarkably, a good agreement between crystallographic data and theoretical results is found for [Cl@Si32Cl44]-. Making a comparison between the hollow silicon cluster and the halide-encapsulating cages, it seems that the incorporation of the halide ion does not significantly alter the silicon cluster, since [Br@Si32Cl44]- and [I@Si32Cl44]- geometries exhibit similar bond distances, whereas in [Cl@Si32Cl44]- all Si-Si bonds are slightly reduced indicating a marginal cluster contraction. This denotes that the Si32Cl44 cage is rigid enough to overcome any further deformation upon halide encapsulation. For all structures we carry out a NCI Analysis27,28 (see Figure 2) in order to determine which atoms contributing the most in keeping the anion in the center of the cavity, This analysis is based on a 2D plot of the reduced density gradient s, and the electron density ρ, where s can be expressed as

ACS Paragon Plus Environment

6

Page 7 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

=

 1 ⁄ 

23 ⁄

When non-covalent interactions arise, there exists a noticeable variation in s, producing density critical points that generate an isosurface. To determine the type of interaction, the second eigenvalue of the electronic density Hessian (λ2) is invoked, where intense stabilizing interactions, such as hydrogen bonds, are characterized by λ2 < 0, non-bonded interactions such as steric repulsion by λ2 > 0, and weak interactions by λ2 ≈ 0. These λ2 values are denoted over the isosurface by using a color scale. In the [X@Si32Cl44]- series, the NCI analysis evidences the rise of non-covalent forces between the silicon atoms and the central anion, irrespective of the silicon environment. As it is expected, these noncovalent forces become stronger for those silicon centers in close proximity to the halide, which is denoted by blue regions in the isosurface. The nature of the cluster∙∙∙anion interactions can be explored in deep through the Energy Decomposition Analysis (EDA) proposed by Morokuma and implemented by Ziegler29,30 which separates the total interaction energy in different terms such as electrostatic (∆Eelec), Pauli repulsion (∆EPauli), and orbital overlapping (∆Eorb). In this analysis, the interaction between the frozen charge densities is accounted by ∆Eelec, which is calculated by subtracting the Coulomb integral of the fragments from that corresponding to the overall system. The product of density fragments, which is normalized but violates the Pauli principle, is antisymmetrized and renormalized to give an intermediate state. The energy difference between the states before and after the antisymmetrization-renormalization is called Pauli repulsion, ∆EPauli, and can be related to destabilizing forces between the fragments. Finally, the fragment orbitals are relaxed to yield the final state corresponding to the inclusion system. The energy decrease caused from such orbital mixing is identified as the orbital interaction, ∆Eorb. Adding ∆Eelec, ∆EPauli, and ∆Eorb we obtain the total interaction energy ∆Eint, of the fragments. To overcome the basis set

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 26

superposition error (BSSE), the counterpoise method proposed by Boys and Bernardi51 is employed denoting forces of about 3 kcal mol-1, which are negligible with respect to the intensity of the studied energies. The dispersion interaction is included according to the Grimme pairwise dispersion correction42 where an attractive energy term summed over all atomic pairs is incorporated, taking into account the weak forces associated to instantaneous fluctuations of fragment electronic densities. Finally, the energy difference between the isolated fragments geometries and the overall system is denoted as preparation energy, ∆Eprep. The Energy Decomposition Analysis (Table 2) evidences the major stability of the [Cl@Si32Cl44]- structure (-110.35 kcal mol-1), where the main contributing term is ∆ , representing about 71% of the stabilizing terms, denoting the mainly electrostatic character of the cavity∙∙∙anion interactions. Remarkably, ∆ significantly increases in line with the volume of the encapsulated ion, i.e. from -159.77 kcal mol-1 in [Cl@Si32Cl44]- to -278.21 kcal mol-1 in [I@Si32Cl44]-. However, this effect is counterbalanced by the increment of the Pauli repulsion which varies from -113.41 kcal mol-1 in [Cl@Si32Cl44]-, to -300.70 kcal mol-1 in [I@Si32Cl44]-. The calculated preparation energy, ∆Eprep, varies in agreement to the suggested stiffness from geometrical parameters (see above), showing a small ∆Eprep amounting to 1.77 for [Cl@Si32Cl44], 1.58 for [Br@Si32Cl44]- and 1.66 kcal mol-1 for [I@Si32Cl44]-. The effect of replacing the peripheral chlorine atoms by hydrogens can be evaluated through the [X@Si32H44]- series (see Table 3), where it is clear that a significant loss of stability occurs, with total binding forces now ranging from -68.53 kcal mol-1 in [Cl@Si32H44]- to -47.64 kcal mol-1 in [I@Si32H44]-. The ∆ term is the most affected with the absence of the peripheral trichlorosilyl and chlorines, unraveling that an inductive effect promoted by these groups, that make the inner cavity more electron-deficient occurs, also enhancing the electrostatic interactions towards the negative-charged central ion, and hence ensuring the latter encapsulation observed ACS Paragon Plus Environment

8

Page 9 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

for chloride anion. Furthermore, the dispersion forces contribute the least to total interaction energy in both the [X@Si32Cl44]- and [X@Si32H44]- series, representing about 2.5%. Interestingly, the obtained preparation energy, ∆Eprep, in the [X@Si32H44]- series, denotes a lower stiffness in comparison to the [X@Si32Cl44]- series, showing a larger ∆Eprep amounting to 3.24 for [Cl@Si32H44]-, 4.82 for [Br@Si32H44]- and 8.43 kcal mol-1 for [I@Si32H44]-, suggesting that the electron withdrawing -SiCl3 groups are able to induce more rigidity in comparison to -SiH3, in addition to a more favorable encapsulation energy, supporting the statement of their relevant role from the initial report on [Cl@Si32Cl44]- .25 Moreover, a small difference in the favorable encapsulation energies is found between [Cl@Si32Cl44]- and [Br@Si32Cl44]- with interaction energies of -110.35 and -105.59 kcal mol-1, respectively, suggesting that a bromine-containing silicon cluster could be stable enough to be obtained experimentally, and also [I@Si32Cl44]- despite of the lower encapsulation energy (-93.54 kcal mol-1). This exposes that the hypothetical analogous structures are suitable for further experimentation. In addition, the case of an endohedral fluorine atom is studied (see Table S2), revealing a ∆Eint of -132.65 kcal mol-1, where the ∆Eelec term contributes in 65.82%. Thus, the total interaction energy of [F@Si32Cl44]- is comparable to that observed for its Cl, Br and I counterparts. When the peripheral chlorine atoms are replaced by H, in [F@Si32H44]-, the stabilization is dramatically decreased to -3.37 kcal mol-1, highlighting the relevant effect of the peripheral chlorines in making the cluster cavity more positive, and then suitable for the incorporation of a negative-charged species. In order to acquire a deeper understanding of the inclusion forces, we decompose the main contributing term to the stabilizing forces, i.e. the ∆ term, into two groups, one encompassing ion-multipole interactions (i.e. ion-dipole, ion-quadrupole, etc.) and other covering

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

higher-order interactions (i.e. dipole-dipole, dipole-quadrupole, quadrupole-quadrupole, etc.) following a protocol recently applied to inclusion compounds displaying metallic centers.45,46,52 Then, the central halide ions (Cl-, Br-, I-) are replaced by their isoelectronic noble gas counterparts (Ar, Kr, Xe), thus enabling only the rise of the higher-order forces. The comparison   of the electrostatic terms for each series (1 − ∆ /∆ )% allows us to quantify the

contribution of ion-multipole electrostatic forces into ∆ . The results of the EDA for these hypothetical silicon clusters containing noble gases are presented in Table S1. The ion-multipole nature of the electrostatic forces is revealed, representing 83.14%, 77.24% and 69.67% of ∆ for [Cl@Si32Cl44]-, [Br@Si32Cl44]- and [I@Si32Cl44]-, respectively. Additionally, the increasing relevance of the higher order terms from the chloride system to the iodine one can be explained by the major electronic cloud distortion experienced by the anion as its volume increases. In turn, the orbital term (∆ ) contribution also increases with the halide volume, suggesting a more effective orbital overlap, which in turn allows a higher electron charge transfer from the cluster to the anion. These results are in good agreement with the atomic charges obtained from the multipole analysis (see Table 4), where a -0.44 |e-| charge amount is transferred from the cluster to the iodide ion, quite far from the +0.22 |e-| and -0.03 |e-| obtained for the chloride and bromide clusters, respectively. The influence of the peripheral halide atoms is evaluated along the [Cl@Si32X44]- (X=F, Cl, Br, I) series (Table S3), where it is clear that as more electronegative X, a stronger encapsulation energy is found, going from -202.71 kcal mol-1 in [Cl@Si32F44]- to -73.51 kcal mol1

in [Cl@Si32I44]-. In all cases it is the ∆Eelec which contributes the most to total interaction

energy, but this effect is more pronounced in [Cl@Si32F44]- and [Cl@Si32Cl44]- (~71%) in

ACS Paragon Plus Environment

10

Page 11 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

comparison to [Cl@Si32Br44]- and [Cl@Si32I44]- (~57%). These results support further explorative synthesis efforts for obtaining the related F exo-bonded derivative. To quantify the extent of the electronic cloud distortion experienced by the silicon atoms with the inclusion of the halide ion, we calculate their local dipole vectors and quadrupole tensors according to the methodology proposed by Swart and coworkers31 where the local electrostatic potential is mainly determined by the charge distribution around the respective atom, or within the atomic multipole expansion, by the atomic multipoles near to that point. This charge analysis offers an accurate description of the electrostatic potential from the charge distribution in molecules by writing the total density as a sum of atomic densities expressed in terms of atomic functions, which in turn are used to define a set of atomic multipoles. This methodology enables a simple description of charge distributions in endohedral complexes, which represents a major challenge.52,53 The results of the multipole analysis are presented in Tables 5 and 6, and the corresponding graphical representations in Figures 3 and 4. The vector analysis reveals a more pronounced electronic cloud distortion in terms of dipole moments for those silicon atoms bonded to chlorine (Si1), which is expected given the influence caused by the near electronegative halide. Then, the Si1 dipole moment magnitude varies from 0.963 D in Si32Cl44. to 0.511 D in [Cl@Si32Cl44]-, 0.485 D in [Br@Si32Cl44]-, and 0.395 D in [I@Si32Cl44]-. On the other hand, the vectors associated to those silicon atoms linked to the trichlorosilyl groups (Si2) vary to a considerable lesser extent from 0.310 D in Si32Cl44 to 0.215 D in [Cl@Si32Cl44]-, 0.235 D in [Br@Si32Cl44]-, and 0.205 in [I@Si32Cl44]-. In the case of the peripheral silicon atoms (Si3), their dipole magnitude reduction can be associated with the charge transfer toward the central halide. The quadrupole anisotropy values show a significant electronic deviation in Si1 as the size of the central halide increases, especially in the case of [I@Si32Cl44]-, where a different

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

electronic cloud distribution is found, with the quadrupole positive lobe pointing toward the center of the cavity (see Figure 4). In the case of Si2, although the presence of the chloride ion clearly modifies the electronic density, it is the incorporation of Br- and I-, which causes a more pronounced distortion, reversing the lobes charges, thus favoring the presence of the central anion. Notably, the small quadrupole tensors for the endohedral halide anions in all evaluated systems suggest that there is not a significant electronic cloud deviation for them, even in the iodide case. Finally, the 29Si-NMR properties were evaluated owing to the valuable data they offer for further experimental guidance. For [X@Si32Cl44]-, the three chemically different Si atoms were identified, namely, Si1, Si2 and Si3 (Table 7 and Figure 1) in a 8:12:12 ratio, which exhibits an almost spherical nuclear shielding tensor due to the sp3-hybridization of the silicon belonging to the structural backbone. For the experimentally characterized [Cl@Si32Cl44]- cluster, the calculated values for the isotropic chemical shift (δiso) amounts to 31.6, -59.9 and 13.3 ppm, for Si1, Si2 and Si3, respectively, in good agreement with the experimental

29

Si-NMR spectrum

values (31.1, -60.3 and 10.3 ppm).25 For the hypothetical hollow Si32Cl44 cluster, the respective values are 53.9, -63.1 and 7.9 ppm, denoting a high shielding for Si1 (~22 ppm), in contrast to the deshielding in a lesser extent for Si2 and Si3, of about 3 and 5 ppm. Such observation is in line with the more predominant role of the silicon atoms bonded to peripheral chlorines (Si1) into the guest encapsulation, in comparison to those bearing trichlorosilyl ligands, and appears as more sensitive to the chlorine encapsulation in the 29Si-NMR spectrum. Similarly, for the Br- and I- clusters, Si2 and Si3 remain in a narrow range in comparison to [Cl@Si32Cl44]-, unraveling the progressive shielding of the chemical shift for Si1 from 53.9 (Si32Cl44) < 31.6 ([Cl@Si32Cl44]-) < 25.7 ([Br@Si32Cl44]-) < 15.5 ppm ([I@Si32Cl44]-), further

ACS Paragon Plus Environment

12

Page 13 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

confirming the respective characteristic pattern in the

29

Si-NMR spectrum. We also notice that

into [I@Si32Cl44]- the peaks corresponding to Si1 and Si2 are almost overlapped.

CONCLUSIONS

A series of [20]Silafullerenes exhibiting an endohedral halide ion have been studied through a detailed DFT protocol, which reveals the preference of the Si20 framework for the chloride guest, as well as the strong electrostatic character (∆Eelec of about 71%) of the noncovalent forces in all systems. Looking in detail for the ∆Eelec, the increasing relevance of the higher order terms going from [Cl@Si32Cl44]- to [I@Si32Cl44]- is manifested, which can be explained by a major electronic cloud distortion for the silicon atoms as the guest volume increases. The small difference between the total interaction energies corresponding to the recently reported [Cl@Si32Cl44]- (-110.35 kcal mol-1) and its hypothetical counterpart [Br@Si32Cl44]- (105.59 kcal mol-1) suggests that experimental efforts could be oriented to the synthesis of the latter. Conversely, the lower stability of [I@Si32Cl44]- (-93.54 kcal mol-1) can be rationalized by the sharp rise of the Pauli repulsion term (300.70 kcal mol-1), more than twice of that corresponding to [Cl@Si32Cl44]- (113.41 kcal mol-1). A comparative study between the [X@Si32Cl44]- clusters and the hypothetical silyl counterparts with -SiH3 at the position of -SiCl3, given by [X@Si32H44]- analogues, evidences that peripheral chloride atoms and trichlorosilyl groups clearly promotes the guest incorporation as they make the inner cavity more positive, and hence more suitable for accommodating the halide guest. It supports the research of new peripheral groups able to influence the cage charge through inductive effects. Thus, the role of -SiCl3 is crucial in the initial encapsulation of chlorine

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 26

in the experimentally observed [Cl@Si32Cl44]- species, making [Br@Si32Cl44]- as a potentially synthetic target, among related [X@Si32Cl44]- structures. The electronic cloud distortion experienced by the silicon atoms with the guest inclusion can be rationalized in terms of local dipoles and quadrupole moments. This analysis reveals that silicon atoms linked to the peripheral chlorines play a more relevant role in encapsulating the endohedral guest variation in comparison to those linked to trichlorosilyl, as it is evidenced by their more pronounced multipole variation with the guest incorporation. We envisage that the stability of [Cl@Si32Cl44]- as well as the apparent synthetic accessibility to its analogues will stimulate the exploration of similar materials as an interesting approach to fabricate novel singleelectron transistors.

ACKNOWLEDGEMENTS The authors thank the financial support from FONDECYT 1180683.

ACS Paragon Plus Environment

14

Page 15 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TABLES AND FIGURES

Table 1: Selected structural parameters of the studied systems (Å and degrees). [Cl@Si32Cl44]-

Si32Cl44

Exp.a

Calc.

[Br@Si32Cl44]- [I@Si32Cl44]-

Si1-Si2

2.405

2.355

2.398

2.403

2.412

Si2-Si3

2.367

2.320

2.367

2.370

2.373

Si2-Si2

2.388

2.352

2.386

2.394

2.403

Si1···Si1

6.664

6.547

6.639

6.659

6.691

Si2···Si2

6.779

6.628

6.761

6.777

6.795

Si2-Si1-Si2

109.3

108.9

109.3

109.3

109.1

Si1-Si2-Si1

106.2

106.7

106.1

106.2

106.4

Si1-Si2-Si3

107.9

107.5

108.2

108.1

108.1

Si1···center

3.318

3.319

3.330

3.330

3.346

Si2···center

3.386

3.382

3.388

3.388

3.397

a. Experimental data from reference 22

Table 2: Energy Decomposition Analysis (EDA) (kcal mol-1) for the halide-substituted silicon clusters. [Cl@Si32Cl44]-

[Br@Si32Cl44]1.58

[I@Si32Cl44]-

∆Eprep

1.77

1.66

∆Eorb

-57.35

25.63%

-73.45

25.80%

-108.82

27.60%

∆Eelec

-159.77

71.40%

-203.94

71.64%

-278.21

70.57%

∆Edisp

-6.64

2.97%

-7.28

2.56%

-7.21

1.83%

∆EPauli

113.41

179.08

300.70

∆Eint

-110.35

-105.59

-93.54

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

Table 3: Energy Decomposition Analysis (EDA) (kcal mol-1) for the hydrogen-substituted silicon clusters. [Cl@Si32H44]-

[Br@Si32H44]4.28

[I@Si32H44]-

∆Eprep

3.24

8.43

∆Eorb

-62.30

35.55%

-76.33

32.76%

-105.89

31.97%

∆Eelec

-106.95

61.03%

-150.37

64.54%

-219.78

66.35%

∆Edisp

-5.99

3.42%

-6.30

2.70%

-5.59

1.69%

∆EPauli

106.71

170.58

283.62

∆Eint

-68.53

-62.42

-47.64

Table 4. Charge variation for selected atoms in the studied systems. [Cl@Si32Cl44]-

[Br@Si32Cl44]-

[I@Si32Cl44]-

Si1

-0.09

-0.08

-0.04

Si2

0.08

0.11

0.14

Si3

-1.11

-1.12

-1.14

X

0.22

-0.03

-0.44

Table 5: Electronic dipole vectors for silicon and halide ions (Debyes). Si32Cl44

[Cl@Si32Cl44]- [Br@Si32Cl44]- [I@Si32Cl44]-

Si1

0.963

0.511

0.485

0.395

Si2

0.310

0.215

0.235

0.205

Si3

1.161

0.688

0.708

0.750

ACS Paragon Plus Environment

16

Page 17 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 6. Anisotropy of the electronic quadrupole tensors for the representative silicon atoms and halide ions (Buckinghams). Si32Cl44

[Cl@Si32Cl44]- [Br@Si32Cl44]-

[I@Si32Cl44]-

Si1

-3.7305

-2.5918

-1.6536

0.3398

Si2

0.6075

0.5878

-0.9071

-1.1787

Si3

2.1675

1.3262

1.1482

0.7379

X

-

-0.0005

0.0006

0.0042

Table 7. Calculated 29Si-NMR chemical shifts (δ) for the studied series (ppm).

Si32Cl44 [Cl@Si32Cl44][Br@Si32Cl44][I@Si32Cl44]a

Si1 53.9 31.6 25.7 15.5

Si2 -63.1 -59.9 -58.8 -58.0

Si3 7.9 13.3 13.4 14.6

Absolute shielding referenced to TMS (σref=366.3 ppm)

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

Figure 1. Optimized structures and some relevant distances for the inclusion compounds (distances in Å)

ACS Paragon Plus Environment

18

Page 19 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. NCI Analysis of the studied systems.

Figure 3. Graphical representation of the electronic dipole vectors of the studied systems where the blue arrowhead denotes the negative region (Debyes).

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 26

Figure 4. Graphical representation of the electronic quadrupole tensors of the studied systems (orange lobes denote negative regions and the blue ones, positive regions).

Figure 5. Calculated 29Si-NMR spectrum for the studied systems.

ACS Paragon Plus Environment

20

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

BIBLIOGRAPHIC REFERENCES [1] Morton, J. J. L.; McCamey, D. R.; Eriksson, M. A.; Lyon, S. A. Embracing the Quantum Limit in Silicon Computing. Nature 2011, 479, 345-353. [2] Bai, J.; Cui, L.-F.; Wang, J.; Yoo, S.; Li, X.; Jellinek, J.; Koehler, C.; Frauenheim, T.; Wang, L.-S.; Zeng, X. C. Structural Evolution of Anionic Silicon Clusters. J. Phys. Chem. A 2006, 110, 908-912. [3] Prieto, J.; Torres, G. A.; Ghoneim, M. T.; Inayat, S. B.; Ahmed, S. M.; Hussain, A. M.; Hussain, M. M. Transformational Silicon Electronics. ACS Nano 2014, 8, 1468-1474. [4] Ross, M. W.; Castleman Jr. A. W. Influence of Group 10 Metals on the Growth and Subsequent Coulomb Explosion of Small Silicon Clusters under Strong Light Pulses. ChemPhysChem 2013, 14, 771-776. [5] Geitner, F. S.; Fässler, T. F. Low Oxidation State Silicon Clusters-Synthesis and Structure of [NHCDippCu(η4-Si9)]3−. Chem. Commun. 2017, 53, 12974-12977. [6] Grubisic, A.; Ko, Y. J.; Wang, H.; Bowen, K. H. Photoelectron Spectroscopy of Lanthanide−Silicon Cluster Anions LnSin- (3 ≤ n ≤ 13; Ln = Ho, Gd, Pr, Sm, Eu, Yb): Prospect for Magnetic Silicon-Based Clusters. J. Am. Chem. Soc. 2009, 131, 10783-10790. [7] Li, X.; Claes, P.; Haertelt, M.; Lievens, P.; Janssens, E.; Fielicke, A. Structural Determination of Niobium-doped Silicon Clusters by Far-infrared Spectroscopy and Theory. Phys. Chem. Chem. Phys. 2016, 18, 6291-6300. [8] Bai, Z.; Liu, X.; Lian, Z.; Zhang, K.; Wang, G.; Shi, S.-F.; Pi, X.; Song, F. A Silicon Cluster Based Single Electron Transistor with Potential Room-Temperature Switching. Chin. Phys. Lett. 2018, 35, 037301. [9] Weber, W. M.; Geelhaar, L.; Graham, A. P.; Unger, E.; Duesberg, G. S.; Liebau, M.; Pamler, W.; Chèze, C.; Riechert, H.; Lugli, P.; Kreupl, F. Silicon-Nanowire Transistors with Intruded Nickel-Silicide Contacts. Nano Lett. 2006, 6, 2660-2666. [10] Barman, S.; Sen, P.; Das, G. P. Ti-Decorated Doped Silicon Fullerene: A Possible Hydrogen-Storage Material. J. Phys. Chem. C 2008, 112, 19963-19968. [11] Li, Y.; Lyon, J. T.; Woodham, A. P.; Fielicke, A.; Janssens, E. The Geometric Structure of Silver Doped Silicon Clusters. ChemPhysChem 2014, 15, 328-336.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 26

[12] Zhang, D.; Guo, G.; Liu, C. New Route for Stabilizing Silicon Fullerenes. J. Phys. Chem. B 2006, 110, 14619-14622. [13] Qin, W.; Lu, W.-C.; Xia, L.-H.; Zhao, L.-Z.; Zang, Q.-J.; Wang, C. Z.; Ho, K. M. Theoretical Study on the Structures and Optical Absorption of Si172 Nanoclusters. Nanoscale 2015, 7, 14444-14451. [14] Pham, H. T.; Tam, N. M.; Jeilani, Y. A.; Nguyen, M. T. Structural Evolution and Bonding of Phosphorus-doped Silicon SinPm-/0/+ with n = 1-10, m = 1,2. Comput. Theor. Chem. 2017, 1107, 115-126. [15] Chu, Q.-Y.; Li, B. X.; Yu, J. Stability of Carbon-doped Silicon Clusters: FP-LMTO Molecular Dynamics Calculations. J. Mol. Struct. THEOCHEM 2007, 806, 67-76. [16] Grubisic, A.; Ko, Y. J.; Wang, H.; Bowen, K. H. Photoelectron Spectroscopy of LanthanideSilicon Cluster Anions LnSin- (3 ≤ n ≤ 13; Ln = Ho, Gd, Pr, Sm, Eu, Yb): Prospect for Magnetic Silicon-Based Clusters. J. Am. Chem. Soc. 2009, 131, 10783-10790. [17] Reveles, J. U.; Khanna, S. N. Electronic Counting Rules for the Stability of Metal-Silicon Clusters. Phys. Rev. B 2006, 74, 035435. [18] Claes, P.; Ngan, V. T.; Haertlet, M.; Lyon, J. T.; Fielicke, A.; Nguyen, M. T.; Lievens, P.; Janssens, E. The Structures of neutral Transition metal Doped Silicon Clusters, SinX (n = 6-9; X = V, Mn). J. Chem. Phys. 2013, 138, 194301. [19] Ziella, D. H.; aputo, M. C.; Provasi, P. F. Study of Geometries and Electronic Properties of AgSin Clusters using DFT/TB. Int. J. Quantum Chem. 2011, 111, 1680-1693. [20] Jin, X.; Arcisauskaite, V.; McGrady, J. E. The Structural Landscape in 14-Vertex Clusters of Silicon, M@Si14: When Two Bonding Paradigms Collide. Dalton Trans. 2017, 46, 11636-11644. [21] Jaeger, J. B.; Jaeger, T. D.; Duncan, M. A. Photodissociation of Metal-Silicon Clusters: Encapsulated versus Surface-Bound Metal. J. Phys. Chem. A 2006, 110, 9310-9314. [22] Ngan, V. T.; Pierloot, K.; Nguyen, M. T. Mn@Si14+: A Singlet Fullerene-like Endohedrally Doped Silicon Cluster. Phys. Chem. Chem. Phys. 2013, 15, 5493-5498. [23] Pichierri, F.; Kumar, V.; Kawazoe, Y. Encapsulation of Halide Anions in Perhydrogenated Silicon Fullerene: X-@Si20H20 (X = F, Cl, Br, I). Chem. Phys. Lett. 2005, 406, 341-344. [24] Denis, P. A. Stacked Functionalized Silicene: A Powerful System to adjust the Electronic Structure of Silicene. Phys. Chem. Chem. Phys. 2015, 17, 5393-5402.

ACS Paragon Plus Environment

22

Page 23 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[25] Tillmann, J.; Wender, J. H.; Bahr, U.; Bolte, M.; Lerner, H.-W.; Holthausen, M. C.; Wagner, M. One-Step Synthesis of a [20]Silafullerane with an Endohedral Chloride Ion. Angew. Chem. Int. Ed. 2015, 54, 5429-5433. [26] Pei, Y.; Gao, Y.; Zeng, X. C. Exohedral Silicon Fullerenes: SiNPtN/2 (20 ≤ N ≤ 60). J. Chem. Phys. 2007, 127, 044704. [27] Contreras-García, J.; Johnson, E. R.; Keinan, S.; Chaudret, R.; Piquemal, J.; Beratan, D. N.; Yang, W. NCIPLOT: A Program for Plotting Noncovalent Interaction Regions. J. Chem. Theory Comput. 2011, 7, 625- 632. [28] Saleh, G., Lo Presti, L., Gatti, C.; Ceresoli, D. NCImilano: An Electron-Density-Based Code for the Study of Noncovalent Interactions. J. Appl. Cryst. 2013, 46, 1513–1517. [29] Morokuma, K. Molecular Orbital Studies of Hydrogen Bonds. J. Chem. Phys., 1971, 55, 1236-1244. [30] Su, P.; Li, H. Energy decomposition analysis of covalent bonds and intermolecular interactions. J. Chem. Phys. 2009, 131, 014102/1-15. [31] Swart, M.; Van Duijnen, P.; Snidjers, J. G. A Charge Analysis Derived from an Atomic Multipole Expansion. J. Comput. Chem. 2001, 22, 79- 88. [32] Parr, R.; Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: London, U.K., 1989. [33] Baerends, E. J.; Ziegler, T.; Autschbach, J.; Bashford, D.; Bérces, A.; Bickelhaupt, F. M.; Bo, C.; Boerrigter, P. M.; Cavallo, L.; Chong, D. P. et al. ADF2014; Vrije Universiteit: Amsterdam, The Netherlands, 2014. [34] Becke, A. D. Density-functional Exchange Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098-3100. [35] Perdew, J. P. Density-functional Approximation for the Correlation Energy of the Inhomogeneous Electron Gas. Phys. Rev. B 1986, 33, 8822-8824. [36] Haertelt, M.; Lyon, J. T.; Claes, P.; de Haeck, J.; Lievens, P.; Fielicke, A. Gas-phase Structures of Neutral Silicon Clusters. J. Chem. Phys. 2012, 136, 064301. [37] Ngan, V. T.; Gruene, P.; Claes, P.; Janssens, E.; Fielicke, A.; Nguyen, M. T.; Lievens, P. Disparate Effects of Cu and V on Structures of Exohedral Transition Metal-Doped Silicon Clusters: A Combined Far-Infrared Spectroscopic and Computational Study. J. Am. Chem. Soc. 2010, 132, 15589-15602.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 26

[38] Gruene, P.; Fielicke, A.; Meijer, G.; Janssens, E.; Ngan, V. T.; Nguyen, M. T.; Lievens, P. Tuning the Geometric Structure by Doping Silicon Clusters. ChemPhysChem 2008, 9, 703-706. [39] Zhao, Y.; Xu, W.; Li, Q.; Xie, Y.; Schaefer III, H. F. Gallium Clusters Gan (n = 1-6): Structures, Thermochemistry, and Electron Affinities. J. Phys. Chem. A. 2004, 108, 7448-7459. [40] Harb, M.; Rabilloud, F.; Simon, D. Structural, Electronic, Magnetic and Optical Properties of Icosahedral Silver-Nickel Nanoclusters. Phys. Chem. Chem. Phys. 2010, 12, 4246-4254. [41] Verluis, L.; Ziegler, T. The Determination of Molecular Structures by Density Functional Theory. The Evaluation of Analytical Energy Gradients by Numerical Integration. J. Chem. Phys. 1988, 88, 322-328. [42] Grimme, S. Density Functional Theory with London Dispersion Corrections. WIREs Comput. Mol. Sci. 2011, 1, 211-228. [43] Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868 [44] Handy, N. C.; Cohen, A. J. Left-Right Correlation Energy. Mol. Phys. 2001, 99, 403–412. [45] Ponce-Vargas, M.; Muñoz-Castro, A. A Study on the Versatility of Metallacycles in Host– guest chemistry: Interactions in Halide-centered Hexanuclear Copper (II) Pyrazolate Complexes. Phys. Chem. Chem. Phys. 2014, 16, 13103-13111. [46] Ponce-Vargas, M.; Muñoz-Castro, A. Metal containing Cryptands as Hosts for Anions: Evaluation of Cu (I)⋯ X and π⋯ X interactions in Halide–tricopper (I) Complexes through Relativistic DFT Calculations. Phys. Chem. Chem. Phys. 2015, 17, 18677-18683. [47] Ortolan, A.; Garamori, G. F.; Bickelhaupt, F. M.; Parreira, R. L. T.; Muñoz-Castro, A.; Kar, T. How the Electron-deficient Cavity of Heterocalixarenes recognizes Anions: Insights from Computation. Phys. Chem. Chem. Phys. 2017, 19, 24696-24705. [48] Ortolan, A. O.; Østrøm, I.; Caramori, G. F.; Parreira, R. L. T.; Da Silva, E. H.; Bickelhaupt, M. Tuning Heterocalixarenes to Improve Their Anion Recognition: A Computational Approach. J. Phys. Chem. A 2018, 122, 3328-3336. [49] Zhurko, G. A.; Zhurko, D. A. Chemcraft. Version 1.7. [50] Humphrey, W.; Dalke, A.; Schulten, K. VMD – Visual Molecular Dynamics. J. Molec. Graph. 1996, 14.1, 33-38.

ACS Paragon Plus Environment

24

Page 25 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

[51] Boys, S. F.; Bernardi, F. The Calculation of Small Molecular Interactions by the Differences of Separate Total Energies. Some Procedures with Reduced Errors. Mol. Phys. 1970, 19, 553566. [52] Ponce-Vargas, M.; Muñoz-Castro, A. Tiara-like Complexes acting as Iodine Encapsulating Agents: The Role of M···I Interactions in [M(µ-SCH2CO2Me)2]8⊂I2 (M = Ni, Pd, Pt) Inclusion Compounds. J. Phys. Chem. C 2016, 120, 23441-23448. [53] Denis, P. A. Chemical Reactivity of Lithium-doped Fullerenes. J. Phys. Org. Chem. 2012, 25, 322-326.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 26

TOC GRAPHIC

ACS Paragon Plus Environment

26