Exploring Chimeric Calix[4]tetrolarene Molecular Scaffolds

Apr 4, 2018 - The enhanced solubility of the modified calix[4]tetrolarene macrocycles has been further explained from diminutive cooperative hydrogen ...
0 downloads 4 Views 2MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

A: Molecular Structure, Quantum Chemistry, and General Theory

Exploring Chimeric Calix[4]tetrolarene Molecular Scaffolds: Theoretical Investigations Dipali N. Lande, and Shridhar P. Gejji J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.8b01686 • Publication Date (Web): 04 Apr 2018 Downloaded from http://pubs.acs.org on April 4, 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 35 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

Hybrid structures of calixarene (lower rim) and pyrogallolarene (upper rim) 194x193mm (96 x 96 DPI)

ACS Paragon Plus Environment

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

Exploring Chimeric Calix[4]tetrolarene Molecular Scaffolds: Theoretical Investigations Dipali N. Lande and Shridhar P. Gejji* Department of Chemistry, Savitribai Phule Pune University, 411 007, India

Corresponding author: [email protected] Fax No.: +91-20-225691728, Telephone No.: +91-20 2560124

1 ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35 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

Abstract Structure and spectral characteristics of chimeric mixture

of calixarene and

pyrogallolarene (usually referred as calix[4]tetrolarene) and its derivatives are studied employing the M06-2X based density functional theory. Different conformers viz., cone, partial cone 1,2-alternate and 1-3-alternate, were identified as the stationary point structures on their potential energy surfaces. Amongst these, the symmetric C4V cone conformer is found to be energetically favorable which can be attributed to cyclic array of hydrogen bonding network in the calix[4]tetrolarene or its thia analogue. The substitution of methoxy groups at the upper rims of calix[4]tetrol- and thicalix[4]tetrol -arenes influences significantly the cooperative hydrogen bonding network and conformational behavior of these hosts. The methoxy substituted macrocycles showed lowering in symmetry from C4v to C2v engendering the pinched cone conformer as the lowest energy structure. The enhanced solubility of the modified calix[4]tetrolarene macrocycles have further been explained from diminutive cooperative hydrogen bonding in its top rim than the pyrogallolarene, which is evidenced from the quantum theory of atoms in molecule and non-covalent interaction reduce density gradient method. Discernibly the underlying cooperative hydrogen bonding emerge with signature in the characteristic vibrational patterns of the calixarene based molecular scaffolds. The chemical shift parameters of their 1H

NMR spectra have further been characterized.

2 ACS Paragon Plus Environment

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 35

1. Introduction Arene based molecular scaffolds such as calix[n]arene,1 resorcin[n]arenes, 2pyrogallol[n]arenes,3

pillar[n]arenes,4

biphen[n]arenes,5

corona[n]arenes,6

asar[n]arenes,7 oxatub[4]arenes8-10 are play a pivotal role in the domain of host-guest chemistry. The structural attributes for example their shape, size, and electronic properties facilitate their use in countless scientific field encompassing: molecular machine,11 supramolecular polymer,12 catalysis,13 and drug delivery14 and makes them fascinating. Amongst these one of the more explored cavitand is calix[n]arenes which has chalice-like structure and derived by condensation of para-substituted phenol derivatives in the presence of formaldehyde. A family of this host usually is endowed with three dimensional hydrophobic cavity that can accommodate a variety of organic and cationic guests.15-18 Modification of these hosts by nucleophilc and electrophilic aromatic substitutions at the endo- or exo- positions of its central annulus makes them largely flexible and paves a way to develop a strategy to devise the newer cavitands which exhibit efficient molecular recognition.19-22 Besides the methylene linkage of calix[n]arenes has led to newer families of molecular containers categorized as hetero- and hetera-calixarenes.23-26 A particular branch of calixarenes designed by varying the nature of starting phenol by 1,2,3-trihydroxybenzene is Pyrogallol[4]arene. To characterized the pyrogallolarenes possessing a concave surface and 12 peripheral hydroxyl groups, the solid state analytical methods particularly, X-ray crystallography, are utilized.

To this direction Mattay,27

Atwood,28 and Rissanen29 inaugurated structures of the “cone” or “rccc” conformers for the isolated scaffolds. Pursuance to this Dakanale30 and coworkers elucidated the conformation and configuration of these cavitands using the NMR experiments. 3 ACS Paragon Plus Environment

The

Page 5 of 35 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

preparation, physico-chemical characterization, and mesomorphic properties of such new class of pyrogallol[4]arene-based columnar liquid crystals have also been illustrated in the literature.30-32 The pioneering works of Cram et al.33-38 considered pyrogallol[4]arenes as the building block for the synthesis of hemicarcerands in which the covalently bonded dimeric capsules provide impetus for further investigations accompanying the host–guest complexation. With this view the strategies to modulate the binding behavior of pyrogallolarene have been developed in the literature, for instance, Reinhoudt et al.39-43 combined the calix[4]arene and pyrogallol[4]arene subunits via incorporation of linkers which yielded rigid container endowed with nanometer-size cavity, usually referred by ‘carceplexes’. These authors further studied immobilization of pyrogallol[4]arene based carceplexes on a gold surface and subsequently analyzed their complexes with the neutral as well as cationic guests. Discernibly blooming of these synthetic receptors has further opened up plethora of applications in interdisciplinary areas and accelerated expansion of the field of supramolecular chemistry.

In this light the chimeric calixtetrolarenes

synthesized by Cohen and coworkers44 which may as well be considered as a fruit of the perfect combination of calixarene (lower rim) and pyrogalloarene (upper rim) molecular receptors. This chimeric cavitand was obtained by simple and efficient one step reaction between partially methylated 1,2,3,5-benzenetetrol with paraformaldehyde using trifluroacetic acid as a catalyst using dichloroethane as the solvent. These authors combined the NMR, HRMS analyses and X-ray diffraction experiments and elucidated its structure. A detailed molecular level understanding of the electronic structure, conformational behavior, cavity size (dimensions) and charge distributions in such synthetic macrocycles is expected to shed light on the rationale on its complexation 4 ACS Paragon Plus Environment

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

behavior. These hosts offer certain advantages for design of efficient experimental systems. Inspired by this in the present endeavor we seek to explore the impact of chimeric mixture of calixarene and pyrogallolarene. The following questions have been addressed How substitution at the bridging positions or at the upper rim of the calix[4]tetrolarene influence the conformational energies and the electronic structure of these macrocycles? How topography of scalar fields such as molecular electrostatic potential characterize the charge distribution within their cavity and shed light on the complexation apriorily? How underlying cooperative hydrogen bonding interactions reflect in the vibrational and 1H NMR spectra? Schematic representation of these macrocycles structures are displayed in scheme 1. The computational strategy adopted is outlined below.

Scheme 1. Atomic labeling and abbreviations used herein.

2. Computational Strategy The present work analyze the energy rank order of calix[4]tetrolarene conformers and their derivatives. The density functional theory (DFT) incorporating the global hybrid meta GGA (M06-2X) functional, which simulates well the hydrogen bonding and adequately typifying other weak interactions, have been applied in conjunction with the 65 ACS Paragon Plus Environment

Page 6 of 35

Page 7 of 35 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

311++G(d,p) basis set (the polarized and diffuse functions being added to all atoms).45-50 Stationary point structures exhibiting qualitatively different hydrogen bonded patterns were optimized using the Gaussian 09 suite of programs.51 Normal vibrational frequencies in various conformers were computed for the stationary point structures which were confirmed to be the local minima on multi-dimensional potential energy surfaces. The normal vibrations were assigned by visualization of the displacements of atoms around their equilibrium (mean) positions using the GAUSSVIEW-5 program.52 The electrophilic and nucleophilic reactive sites were predicted from the MESP comprising of the bare nuclear and electronic contributions which bring forth effective electron-rich regions in the host.53-56 Subsequent MESP topography was mapped by examining eigenvalues of the Hessian matrix where the gradient of V(r) vanishes. Quantum Theory of Atoms in Molecules (QTAIM), proposed by Bader57 was employed to study the molecular electron density topography and accordingly the bond critical points (bcp) were characterized with the AIMAll software.58

The nonlocalized dispersion as well as repulsive nonbonded

interactions are certain from the NCI-RDG method.59 The RDG defined through s = 

|∇ |

( )/ /

within the QTAIM approach and obtained using the Multiwfn program60 and

visualized through the visual molecular dynamics51 software. The chemical shift parameters (δH) in NMR spectra were obtained for the lowest energy conformers by subtracting the nuclear magnetic shielding tensors of protons in the host from those in the tetramethylsilane (reference) within the framework of the gauge-independent atomic orbital (GIAO) method.62 Furthermore the effect of solvent (chloroform) on electronic structure and 1H NMR chemical shifts were simulated using the self-consistent reaction field (SCRF) theory that incorporates the polarization continuum model (PCM).63 6 ACS Paragon Plus Environment

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 35

3. Results and Discussion 3.1 Structural characteristics Calix[4]tetrolarene and its thia derivatives comprising of benzenetetrol (1,2,3,5tetrahydroxybenzene) monomer units can be envisaged as the chimeric mixture of calixarene and pyrogallolarene hosts. Like calix[4]arenes such artificial receptors possess four stereo-centers at the bridging position and different isomeric arrangements led to various conformers showing different arrangements for the hydroxyl groups on the lower (narrow) rim: all four hydroxyls pointing along the same direction (cone structure); one hydroxyl pointing in the opposite direction (partial cone); two neighboring monomers directed in different directions (1,2 conformer) and lastly, two radially opposing monomers pointing in different directions (1,3 conformer), have been optimized. The structures thus derived were subjected to earlier screening using the M06-2X theory with the lower 6-311G basis set. The minimum energy and next low lying structures of these hosts are shown in Figure S1 to S4 of the supporting information. The hierarchy of relative stabilization energies of molecular receptors (I to V) follows the order: cone > partial cone > 1,2-alterrnate > 1,3-alternate.

Subsequently the minimal energy conformers of

calix[4]tetrolarene and its modified analogue subjected to optimizations with the extended 6-311++G(d,p) basis set. The optimized structures are displayed in Figure 1 along with the atomic labelling scheme. As can be inferred

C4v symmetrical cone structure of the

calix[4]tetrolarene (I) possess intramolecular hydrogen bonding interactions between the OH groups (oriented in clockwise direction) at its narrower lower rim

which are

effectively locked in a cyclic network of cooperative hydrogen bonding, contribute to its larger stability over the remaining conformers. Additionally the upper rim (wider) consists 7 ACS Paragon Plus Environment

Page 9 of 35 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

of hydroxyl groups oriented in anticlockwise direction which bring forth interactions between neighboring hydroxyl functionalities. The bridging sulphur atoms (-S-) in II brings about expansion of the host cavity subsequently weakening the cooperative hydrogen bonding network, engenders the shortening of C-O and O-H bonds ( ~0.014 Å and ~0.005 Å respectively) compared to the cone conformer of calix[4]tetrolarene. Substitution of meta methoxy group on the upper rim of parent calix[4]tetrolarene causes larger steric repulsions between the meta methoxy substituents in the symmetrical C4v cone conformer. Thus, III and IV reveal lowering of symmetry concomitantly minimizing steric hindrance, which in turn, transforms to the ‘pinched cone’ structure, wherein the two opposite aromatic rings orient nearly parallel with methoxy substituents directing inside (endo) the cavity whereas the remaining methoxy substituent protruding outside the host cavity; the other two aromatic rings adopt flattened positions with the substituents showing the opposite arrangement compared to parallel ones (cf. Figure 1). The flattened aromatic rings in III revealed a separation of 10.19 Å, which compares well with 10.07 Å observed in the X-ray structure reported by Zafrani and Cohen.44 3.2 Molecular electrostatic potential (MESP) In an effort to predict reactive sites in molecular scaffolds studied herein the MESP topography was used as a puissant descriptor. Figure 2 represents the maps of the electrostatic potential overlaid on 0.001 au isodensity surface of calix[4]tetrolarene and its derivatives. The electrostatic potential intensifications are displayed on the color scale as: red (more negative) < orange < yellow < green< blue (positive). A prevalence of green shows the potential midway between the two extreme ends. As may readily be discernible the negative potential has largely been localized near hydroxyl oxygens of the cavitand 8 ACS Paragon Plus Environment

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

which partly can be attributed to lone pair of oxygen and the π-electron clouds of aromatic rings.

An inner surface within its cavity is less electron rich than the calix[4]arene

receptor. Maximal brunt of a positive charge (dark blue) near the hydroxyl protons can be noticed. The bridging hetero atoms in II pulls the π-electron density of aromatic moiety concomitantly weakening the negative π potential inside its cavity, as evidenced from figure. The methoxylation of upper rim hydroxyl functionalities renders more electro-rich character to the cavity (reddish yellow) and the oxygen centers (dark red) which can be attributed to the disruption hydrogen bond network. Moreover, the lower rim oxygen atoms of the parallel rings render large negative potential. It may, therefore be concluded that MESP topography provides a concise tool for clarifying the substituent effect. Furthermore electron-rich regions around the bridging sulphur atoms in II and IV macrocycles render multiple coordinating sites for effective host-guest binding. In summary, the wide variety of guest molecules can coordinate to interior of bowl (I and II) and funnel (in III and IV) either by hydrophobic effect or various electrostatic interactions. 3.3 Quantum theory of atoms in molecules (QTAIM) An in-depth analysis of hydrogen bonding enables one to rationalize molecular recognition ability and should prove useful for design of task-specific host–guest complexes. In this sense the Bader QTAIM methodology has been proven advantageous to segregate the underlying molecular interactions which tie-up the various chemical concepts related to structure, bonding and reactivity to topological parameters of the electron density distribution function. The negative value of Laplacian of electron density ( ∇ ) at bcp indicates concentration of electronic charge in the internuclear region suggesting sharing of electronic charge between two nuclei ascertaining the covalent bond 9 ACS Paragon Plus Environment

Page 10 of 35

Page 11 of 35 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

(polar). The closed shell type interactions for example, hydrogen bond, ionic bond or van der Waal interactions are recognized from depletion of the electronic charge in the internuclear region. A demarcation of hydrogen bonding from van der Waal interactions can be carried out using the criteria set by Koch and Popelier.64 Accordingly the ρbcp parameters fall in the range 0.002− 0.035 au, whereas the corresponding Laplacian (∇2ρ) parameter lie between 0.014−0.139 au. The correlation between the ∇2ρ and energetic parameters at the bcp can be established with the use of the Virial theorem 0.25 ∇2ρ = 2 Gbcp + Vbcp, where Gbcp, denotes the kinetic energy (conventionally >0) component which is directly proportional to mobility of the electron density, whereas the Vbcp term refers to the potential energy (< 0) of the electrons at bcp and provides a measure of pressure exerted on the electrons at the bcp by other electrons. To comprise such energetic parameters in a system of classification of interatomic contacts, Cremer and Kraka65 suggested a criterion on the basis of the sign of ∇2ρ and total electron energy density Hbcp of the charge distribution which can be expressed as Hbcp = Gbcp + Vbcp . As pointed out earlier, the ∇2ρ < 0 suggests covalent and greater than zero indicates closed-shell interactions. However, sometimes the negative Vbcp value greater than the positive value of Gbcp results in negative Hbcp , but the ∇2ρ is still positive implying relatively stronger hydrogen bonding rendered with partial covalent character. Weak and medium hydrogen bonds show both positive ∇2ρ and Hbcp values. QTAIM topological parameters and the corresponding bond distances are enumerated in Table 1. The stronger cooperative hydrogen bonding (O-H···O) interactions between the hydroxyl groups at lower rim are evident from the presence of bcp with larger electron density of 0.04314 au with the corresponding Laplacian being positive, however the total electron density (Hbcp) turns out be a negative (-0.0032 au) in case of the 10 ACS Paragon Plus Environment

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 35

calix[4]tetrolarene (I) host. It may, therefore, be conjectured that these hydrogen bonding interactions have partially covalent character.

The intramolecular hydrogen bonding

arising from –OH group at the C1 and C3 positions can be ascertained from the corresponding QTAIM topological parameters as opposed to the hydroxyl group at C2 position which do not participate in the hydrogen bonding. The substitution of –S- at linking position engenders larger separation between the hydroxyl groups with smaller ρ (0.0248 au) values implying the weaker hydrogen bonding. The minimum energy pinched cone conformer of host III and IV possess varying hydrogen bond strengths at their lower rims; the parallel ring hydroxyl group reveal stronger interaction than flatten ones. As shown in Table 1, for II to IV both ∇2ρ and H parameter turn out to be > 0 which imply closed shell type interactions. Additionally the correlation of hydrogen bond separation and electron density can also be deciphered from the plot of ln(ρbcp) as a function of hydrogen bond distance d(O-H···O) portrayed in Figure 3. A linear correlation (correlation coefficient 0.98) shows the reciprocal relation between ρbcp and the interatomic distances. Further the variation of G and V as a function of interatomic separation (A···B) is portrayed in Figure 4. Beside this the strength of hydrogen bonding at the bcp was determined by using the Espinosa-Molins-Lecomte formulae (EHB =0.5 x Vbcp). Calculated strength of OH···O interactions thus were estimated to be 56 kJ mol-1 and 41 kJ mol-1 corresponding to the lower and upper rims, respectively, for the calix[4]tetrolarene receptor. The rest of macrocycles (II to IV) hydrogen bond strength falls in the range 15.3 - 26.7 kJ mol-1. These inferences are corroborated through the ellipticity parameters ( =

 

− 1) with strong

hydrogen bonding displaying relatively lower elliptical value (0.07) whilst the weaker hydrogen bonding correspond to the larger  values. 11 ACS Paragon Plus Environment

Page 13 of 35 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

3.4 Noncovalent interactions reduced density gradient (NCI-RDG) The weak intramolecular interactions should not necessarily indicate the presence of the bcp (critical point where the density gradient vanish) within the regime of QTAIM theory. To scrutinize these interactions particularly, the O1H1 and O2H2 hydroxyl group in addition to nonlocalized dispersion and repulsive non-bonded interactions, we apply the NCI-RDG method which has emerged as a powerful tool to categorize and envisaged weak interactions: stabilized (hydrogen bonding), destabilized (steric repulsion) and delocalized weak (van der Waal). The 2D plots of the reduced density gradient (s) versus the electron density multiplied by the sign of the second Hessian eigenvalue (λ2) and 3D RDG isosurfaces in real space by color mapping providing the fingerprint region of close contact within the molecular system. A color coding scheme used is as follows: strong repulsive non-bonded overlap shown in red, attractive interactions in blue and the green regions suggest electrostatic interactions. The 3D spatial visualization of NCI isosurfaces and 2D stiletto-heel shape scatter plots of host I to IV are portrayed in Figure 5. The appearance of a small, flat, pilled shaped dark blue isosurfaces between the O-H···O interactions at the lower rim of host at high critical density region in 2D plot (-0.045au) in calix[4]tetrolarene corroborates the earlier inferences on the stronger hydrogen bonding. The bluish green isosurfaces are observed between the O3-H3···O1 interactions with trough at -0.035 au. The bicolored (green and red) isosurfaces are observed between the O1H1, O2H2 and O3H3 groups of all monomer. Here the attractive component of the surface between the oxygen and hydrogen atoms, is counterbalanced by destabilization originating from the induced five-membered ring closure. The bright red cigar-shaped surfaces which are elongated along the direction of decreasing density appear between the central core of 12 ACS Paragon Plus Environment

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

aromatic rings indicating destabilized (steric crowding) interactions and reflect as trough for high ρ critical values (ρ = 0.020 au) in the 2D scatter plot. Similar inferences are drawn for the rest of molecular scaffolds. 3.5 Infrared and 1H NMR spectra Understanding the molecular mechanism of cooperative self-assembly is a key component for modeling self-assembled supramolecular architectures and analyzed through the lens of the vibrational spectroscopy. Underlying cooperative interactions reflect in the characteristic vibrations of their infrared spectra. Normal vibration frequencies of calix[4]tetrolarene and its derivatives (I to IV) were calculated within the regime of the M06-2x/6-311++G(d,p) theory. The molar absorption coefficient (or, molar absorptivity in units of 0.1 m2mol−1) versus the frequency (in cm−1) showing the -OH stretchings in I to V are illustrated in Figure 6. All these molecular receptors possess four hydroxyl groups at the lower rim. The stretching vibration corresponding to all OH groups in the parent calix[4]tetrolarene coincide; the highly symmetrical C4v cone conformer renders all benzenetetrol units to be equivalent. Theoretical calculations demonstrated that stronger dynamic interactions within all four oscillators stem from circular network of cooperative intramolecular hydrogen bonding. Consequently their vibrations are not precisely vibrations of individual hydroxyl groups and were assigned to stretching vibrations of cyclic (O–H···)4 system with four distinct frequency bands. Out of these ʋ1 near 3513 cm-1 has zero intensity by symmetry rules, the ʋ2 and ʋ3 bands are doubly degenerate intense bands those appear at 3470 cm-1, and the symmetrical stretching ʋ4 shows up near 3378 cm-1 again with the negligible intensity. Similar inferences are drawn for the upper rim hydroxyl groups. For II molecular receptor, pattern of group vibrations of cyclic 13 ACS Paragon Plus Environment

Page 14 of 35

Page 15 of 35 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

tetramer system still preserve, despite the weakening of cooperativity in hydrogen bonding network. The 3716 cm-1 band was assigned to hydroxyl vibrations of parallel aromatic rings and 3594 cm-1 arise from those of flattened rings in the pinched cone conformer (C2V) in III. The ʋ(OH) stretching vibration shows a upshift (blue shift) on modification via substitution at the linkage as in IV due to transfer of electron density from sulphur atom to aromatic rings of the host. Natural bond orbital analysis revealed diminutive electron density in its anti-bonding O-H natural orbital in consonance with the earlier corollaries. In other words the number and intensity patterns of OH vibrations in the infrared spectra serves as an indicator of cooperativity effects in the intramolecular hydrogen bonding. Moreover, the vibration near 3132 cm-1 of the spectra of calix[4]tetrolarene stem from methylene linkage and corresponds to lower wavenumber bands near the 3987 cm-1 in III. Alternatively the cooperativity of hydrogen bonding can be envisaged from the frontier molecular constituting the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). Subsequently the charge distribution within macrocycle are characterized through frontier orbitals the form of electron density (isosurface ±0.04 a. u. ) depicted in Figure 7. It may readily be discernible that HOMO largely resides over the aromatic rings and near hydroxyl oxygen in host I to IV invoking cooperative hydrogen bonding interactions arising from the hydroxyl groups. On the other hand the LUMO is delocalized over the π cloud of the aromatic rings and also near the sulphur atoms. The proton NMR chemical shifts are computed from the GIAO method which serves as criterion for probing the intramolecular interactions. To take into account bulk solvent effects the NMR calculations were performed employing the SCRF-PCM theory with 14 ACS Paragon Plus Environment

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

chloroform as the solvent. Calculated δH values of calix[4]tetrolarene and its derivatives from the M06-2X/6-311++G(d,p) theory are summarized in Table 2. The protons of receptor I to IV derivatives can be classified as: hydroxyl protons at upper rim as H1, H2 and H3 whereas at lower rim as H5; the H6 protons from methyl linkage and H7 refers to those from the methoxy group. The hydroxyl protons at the lower rim (H5) of parent calix[4]tetrolarene facilitate stronger cooperative intramolecular hydrogen bonding interactions engendering largely deshielded (10.83 ppm) signals. The hydrogen bonded H3 protons at the upper rim emerge with δH values near 8.28 ppm whilst the H2 and H1 protons show signals at the 6.11 ppm and 5.61 ppm, respectively. The distinct δH signals of methylene linkage protons, differing by ~ 0.24 ppm, were further noticed. For thia analog of host I the concerned, δH parameters follow the order H5 > H3 > H1 > H2. The H5 protons reveal shielding by 2.66 ppm which reaffirm the weakening of the cooperative hydrogen bonding. The lowest energy pinched cone conformer of host III display the δH signal for the H5 hydroxyl protons near 8.10 ppm (parallel ring); 8.98 ppm (flatten ring); the corresponding δH in the measured spectra was observed at 8.90 ppm. The methylene bridge signals emerge with doublets at 5.33 ppm and 3.44 ppm typifying respectively, the axial and equatorial protons. The methoxy protons chemical shift from the present theory turns out to be 3.75 ppm, which agrees well with its experimental counterpart (near 3.76 ppm). As compared to III, the hydroxyl protons at the lower rim (H5) in IV revealed upfield signals. It may be speculated that conformational behavior, substitution effect and underlying interactions reflect in their δH parameters.

15 ACS Paragon Plus Environment

Page 16 of 35

Page 17 of 35 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

Conclusions In nut shell, the electronic structure, conformational energetics, charge distribution, vibrational and 1H NMR spectra of the Calix[4]tetrolarene and its derivatives, which represent chimers of calixarene (lower rim) and pyrogallolarene (upper rim) have been analyzed in the present endeavor. A systematic conformational search on the molecular scaffolds yields the cone, partial cone, 1,2-alternate and 1,3-alternate conformers. The cooperativity in hydrogen bonding at the lower rim engenders the symmetrical cone conformer to be of the lowest energy for I and II while in case of receptors III and IV the pinched-cone is favored. A disruptions in cooperative hydrogen bonding network of the upper rim of these chimers than the individual pyrogallolarene host, may be expected that enhancing the solubility. Furthermore the methoxy substitution renders largely electronrich cavity in III and IV than the calix[4]tetrolarene or its thia derivative suggesting its favorable binding for the electron deficient or neutral guests. The hydroxyl protons at the lower rim facilitating circular array of hydrogen bonding network reflect as largely deshielded signals in the calculated 1H NMR spectra. Lastly the characteristic patterns in the hydrogen bonded OH vibrations are noticed from the infrared spectra. Supporting Information Different conformers of calix[4]tetrolarene and its derivatives from M06-2X/6-311G(d,p) theory. Acknowledgements SPG acknowledges support from the Research Project (37(2)/14/11/2015-BRNS) from the Board of Research in Nuclear Sciences (BRNS), India. DNL is thankful to UPE Phase-II (UPE 262 (A) 3) project from the University Grand Commission (UGC). Authors thank the Centre 16 ACS Paragon Plus Environment

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 35

for Development of Advanced Computing (CDAC), Pune for providing the National Param Supercomputing Facility.

References 1. Gutsche, C. D. Calixarenes, Monographs. In Supramolecular Chemistry; Stoddart, J. F., Ed.; The Royal Society of Chemistry: Cambridge, 1989. 2. Oshovsky, G. V.; Reinhoudt D. N.; Verboom, W. Triple-ion Interactions for the Construction of Supramolecular Capsules. J. Am. Chem. Soc., 2006, 128, 5270-5278. 3. Haworth, R. D.; Lamberton, A. H. Some Derivatives of Catechol and Pyrogallol. Part II J. Chem. Soc., 1946, 1003–1005. 4. Ogoshi, T.; Kanai, S.; Fujinami, S.; Yamagishi T.; Nakamoto, Y. Para-Bridged Symmetrical Pillar[5]arenes: their Lewis Acid Catalyzed Synthesis and Host–guest Property. J. Am. Chem. Soc., 2008, 130, 5022-5023. 5. Chen, H.; Fan, J.; Hu, X.; Ma, J.; Wang, S.; Li, J.; Yu, Y.; Jia, X.; Li, C. Biphen[n]arenes. Chem. Sci., 2015, 6 (1), 197–202. 6. Ren, W.-S.; Zhao, L.; Wang, M.-X. Functionalized O

6

–Corona[6]arenes: Synthesis,

Structure, and Fullerene Complexation Property. Org. Lett., 2016, 18 , 3126–3129. 7. Schneebeli, S. T.; Cheng, C.; Hartlieb, K. J.; Strutt, N. L.; Sarjeant, A. A.; Stern, C. L.; Stoddart, J. F. Asararenes-A Family of Large Aromatic Macrocycles. Chem. Eur. J., 2013, 19, 3860–3868. 8. Jia, F.; He, Z.; Yang, L.-P.; Pan, Z.-S.; Yi, M.; Jiang, R.-W.; Jiang, W. Oxatub[4]arene: A Smart Macrocyclic Receptor with Multiple Interconvertible Cavities. Chem. Sci. 2015, 6, 6731– 6738. 17 ACS Paragon Plus Environment

Page 19 of 35 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

9. Jia, F.; Wang, H.-Y.; Li, D.-H.; Yang, L.-P.; Jiang, W. Oxatub[4]arene: A Molecular “transformer” Capable of Hosting a Wide Range of Organic Cations. Chem. Commun. 2016, 52, 5666-5669. 10. Jia, F.; Li, D.-H.; Yang, T.-L.; Yang, L.-P.; Dang, L.; Jiang, W. Oxatub[5,6]arene: Synthesis, Conformational Analysis, and the Recognition of C60 and C70. Chem. Commun., 2017, 53, 336–339. 11. Mandl, C.; König, B. Chemistry in Motion—Unidirectional Rotating Molecular Motors. Angew. Chem. Int. Ed., 2004, 43, 1622–1624. 12. Li, C. Pillararene-Based Supramolecular Polymers: from Molecular Recognition to Polymeric Aggregates. Chem. Commun., 2014, 50, 12420–12433. 13. Jiang, J.; Meng, Y.; Zhang, L.; Liu, M. H. Self-Assembled Single-Walled Metal-Helical Nanotube (M-HN): Creation of Efficient Supramolecular Catalysts for Asymmetric Reaction. J. Am. Chem. Soc., 2016, 138, 15629–15635. 14. Sun, T.; Zhang, Y. S.; Pang, B.; Hyun, D. C.; Yang, M.; Xia, Y. Engineered Nanoparticles for Drug Delivery in Cancer Therapy. Angew. Chem. Int. Ed., 2014, 53, 12320-12364. 15. Boulet, B.; Joubert, L.; Cote, G.; Bouvier-Capely, C.; Cossonnet, C.; Adamo, C. A Combined Experimental and Theoretical Study on the Conformational Behavior of a Calix[6]arene. J. Phys. Chem. A, 2006, 110, 5782-5791. 16. Bourdelande, J. L.; Font, J.; Gonzalez-Moreno, R.; Nonell, S. Inclusion Complex of Calix[8] Arene-C60: Photophysical Properties and Its Behaviour as Singlet Molecular Oxygen Sensitiser in the Solid State. J. Photochem. Photobiol. A, 1998, 115, 69. 17. Furer, V. L.; Borisoglebskaya, E. I.; Zverev, V. V.; Kovalenko, V. I. Optical Spectroscopy and Theoretical Studies in Calixarene Chemistry. Spectrochim. Acta, 2005, 62, 483. 18 ACS Paragon Plus Environment

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

18. Arena, G.; Contino, A.; Gulino, F. G.; Magrı̀, A.; Sciotto, D.; Ungaro, R. Complexation of small Neutral Organic Molecules by Water Soluble Calix[4]arenes. Tetrahedron Lett., 2000, 41, 9327-9330. 19. Creaven, B. S.; Donlon, D. F.; McGinley, J.; Coordination Chemistry of Calix[4]arene Derivatives with Lower rim Functionalization and their Applications. Coord. Chem. Rev., 2009, 253, 893-962. 20. Van Loon, J. D.; Arduini, A.; Coppi, L.; Verboom, W.; Pochini, A.; Ungaro, R.; Harkema, S; Reinhoudt, D. N.; Selective Functionalization of calix[4]arenes at the Upper Rim. J. Org. Chem., 1990, 55, 5639-5646. 21. Ozturk, G.; Akkaya, E. U. Differential and substrate-selective reactivity of Calix[4] arene Derivatives with Cyclenyl-Zn (II) Modifications at the Upper Rim. Org. Lett., 2004, 6, 241-243. 22. Joseph, R.; Rao, C. P. Ion and Molecular Recognition by Lower Rim 1, 3-di-conjugates of Calix[4]arene as Receptors. Chem. Rev., 2011, 111, 4658-4702. 23. Sessler, J. L.; Anzenbacher Jr., P.; Shriver, J. A.; Jursíková, K.; Lynch, V. M.; Marquez, M. Direct Synthesis of Expanded Fluorinated Calix[n]pyrroles:  Decafluorocalix[5]pyrrole and Hexadecafluorocalix[8]pyrrole. J. Am. Chem. Soc., 2000, 122, 12061–12062. 24. Nagarajan, A.; Ka, J-W.; Lee, C-H. Synthesis Of Expanded Calix[n]Pyrroles and their Furan or Thiophene Analogues. Tetrahedron, 2001, 57, 7323–7330. 25. Kumar, S.; Paul, D. Singh, H. Syntheses, Structures and Interactions of Heterocalixarenes. Adv. Heterocycl. Chem., 2005, 89, 65–124. 26. Mislin, G.; Graf, E.; Hosseini, M. W.; De Cian, A.; Fischer, J. Synthesis and Structural analysis of Mercaptothiacalix[4]arene, Tetrahedron Lett. 1999, 40, 1129. 19 ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35 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

27. Gerkensmeier, T.; Agena, C.; Iwanek, W.; Frohlich, R.; Kotilad, S.; Nather, C.; Mattay, J. Z. Synthesis and Structural Studies of 5, 11, 17, 23-tetrahydroxyresorc [4]arenes. Naturforsch. 2001, 56, 1063–1073. 28. Atwood, J. L.; Barbour, L. J.; Jerga, A. Hydrogen-bonded Molecular Capsules are Stable in Polar Media. Chem. Commun., 2001, 22, 2376–2377. 29. Mansikkamaki, H.; Nissinen, M.; Rissanen, K. Noncovalent π⋅⋅⋅ π-Stacked Exo-Functional Nanotubes: Subtle Control of Resorcinarene Self-Assembly. Angew. Chem., Int.Ed. 2004, 43, 1243–1246. 30. Bonsignore, S.; DuVosel, A.; Guglielmetti, G.; Dalcanale, E.; Ugozzoli, F. Influence of Steric Interactions and Random Side Chain Variations on the Mesomorphic Properties of Bowlic Mesogens. Liq. Cryst. 1993, 13, 471–482. 31. Ricco, M.; Dalcanale, E. Molecular Conformation and Magnetic Behavior of Macrocyclic Columnar Liquid Crystals. J. Phys. Chem. 1994, 98, 9002–9008. 32. Lippmann, T.; Dalcanale, E.; Mann, G. Influence of Macrocyclic Core Configuration on the Mesomorphic Properties of Liquid Crystalline Calixarenes. Gazz. Chim. Ital. 1995, 125, 595–599. 33. Cram, D. J.; Jaeger, R.; Deshayes, K. Host-guest complexation. 65. Hemicarcerands that Encapsulate Hydrocarbons with Molecular Weights Greater than Two Hundred. J. Am. Chem. Soc. 1993, 115, 10111–10116. 34. Robbins, T. A.; Knobler, C. B.; Bellew, D. R.; Cram, D. J. Robbins, T. A., Knobler, C. B., Bellew, D. R., Cram, D. J. Host-guest complexation. 67 A highly adaptive and strongly binding hemicarcerand. J. Am. Chem. Soc., 1994, 116, 111–120.

20 ACS Paragon Plus Environment

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 35

35. Yoon, J.; Sheu, C.; Houk, K. N.; Knobler, C. B.; Cram, D. Syntheses, Binding Properties, and Structures of Seven New Hemicarcerands Each Composed of Two Bowls Bridged by Three Tetramethylenedioxy Groups and a Fourth Unique Linkage1, 2. J. Org. Chem., 1996, 61, 9323–9339 36. Helgeson, R. C.; Paek, K.; Knobler, C. B.; Maverick, E. F.; Cram, D. J. Guest-Assisted and Guest-Inhibited

Shell

Closures

Provide

Differently

Shaped

Carceplexes

and

Hemicarceplexes1, 2. J. Am. Chem. Soc., 1996, 118, 5590–5604 37. Helgeson, R. C.; Knobler, C. B.; Cram, D. J. Correlations of Structure with Binding Ability Involving Nine Hemicarcerand Hosts and Twenty-Four Guests1. J. Am. Chem. Soc. 1997, 119, 3229–3244; 38. Yoon, J.; Cram, D. J. The First Water-soluble Hermicarceplexes. Chem. Commun. 1997, 497–498. 39. Timmerman, P.; Verboom, W.; van Veggel, F. C. J. M.; van Duynhoven, J. P.; Reinhoudt, D. N. An Organic Molecule with a Rigid Cavity of Nanosize Dimensions. Angew. Chem., Int. Ed. Engl. 1994, 33, 1292–1295. 40. Timmerman, P.; van Veggel, F. C. J. M.; van Duynhoven, J. P.; Reinhoudt, D. N. A Novel Type of Stereoisomerism in Calix[4]arene-Based Carceplexes. Angew. Chem., Int. Ed. Engl., 1994, 33, 2345–2348; 41. Timmerman, P.; Nierop, K. G. A.; Brinks, E. A.; Verboom, W.; van Veggel, F. C. J. M.; van Hoorn, W. P.; Reinhoudt, D. N. Hydrophobic Concave Surfaces and cavities by Combination of Calix[4]Arenes and Resorcin[4]Arenes. Chem. d Eur. J., 1995, 1, 132– 143

21 ACS Paragon Plus Environment

Page 23 of 35 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

42. Higler, I.; Verboom, W.; van Veggel, F. C. J. M.; De Jong, F.; Reinhoudt, D. N.; Liebigs Synthesis and Application of Iso(thio)cyanate-Functionalized Calix[4]arene, Eur. J. of Org. Chem., 1997, 1577–1586. 43. Huisman, B.; van Veggel, F. C. J. M.; Reinhoudt, D. N. Pure Appl. Chem., Supramolecular Chemistry at Interfaces, 1998, 70, 1985–1992. 44. Zafrani, Y.; and Cohen, Y. Calix [4,5] tetrolarenes: A New Family of Macrocycles. Org. let., 2017, 19, 3719-3722. 45. Peverati, R.; Truhlar, D. G. Exchange–Correlation Functional with Good Accuracy for both Structural and Energetic Properties While Depending Only on the Density and Its Gradient. J. Chem. Theory Comput., 2012, 8, 2310–2319. 46. Zhao, Y.; Truhlar, D. G. Applications and Validations of the Minnesota Density Functionals. Chem. Phys. Lett., 2011, 502 , 1–13. 47. Lande D. N.; Rao S. S.; Gejji S. P.; Deciphering Noncovalent Interactions Accompanying 7, 7, 8, 8-Tetracyanoquinodimethane Encapsulation within Biphene [n] arenes: NucleusIndependent Chemical Shifts Approach. ChemPhysChem, 2016, 17, 1 – 14. 48. Goerigk L.; Grimme S.; A thorough Benchmark of Density Functional Methods for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. Phys. Chem. Chem. Phys., 2011, 13, 6670-6688. 49. Walker M.; Harvey A. J.; Sen A.; Dessent C. E.; Performance of M06, M06-2X, and M06-HF Density Functionals for Conformationally Flexible Anionic Clusters: M06 Functionals Perform Better than B3LYP for a Model System with dispersion and Ionic HydrogenBonding Interactions. J. Phys. Chem. A, 2013, 117, 12590-12600.

22 ACS Paragon Plus Environment

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

50. Yuan K.; Guo Y. J.; Yang T.; Dang J. S.; Zhao P.; Li Q. Z.; Zhao X.; Theoretical Insights into the Host–Guest Interactions between [6] Cycloparaphenyleneacetylene and its Anthracene-Containing derivative and Fullerene C70. J. Phys. Org. Chem., 2014, 27, 772782. 51. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian, Inc., Wallingford CT,2009. 52. R. Dennington, T. Keith and J. Milliam, Semichem. Inc., Shawnee Mission, KS, 2009. 53. Murray, J. S.; Shields, Z. P.; Seybold, P. G.; Politzer, P. Intuitive and Counterintuitive Noncovalent Interactions of Aromatic π Regions with the Hydrogen and the Nitrogen of HCN. J. Comp. Sci., 2015, 10, 209–216. 54. Politzer, P.; Murray, J. S.; Clark, T. Mathematical Modeling and Physical Reality in Noncovalent Interactions. J. Mol. Model., 2015, 21, 52. 55. Murry, J. S.; Seminario, J. M.; Politzer. P. A Computational Study of the Structures and Electrostatic Potentials of Some Azines and Nitroazines. J. Mol. Struct., (TheoChem) 1989, 87, 95-108. 56. Lande D. N.; Gejji S. P.; Cooperative Hydrogen Bonding, Molecular Electrostatic Potentials, and Spectral Characteristics of Partial Thia-Substituted Calix[4]arene Macrocycles. J. Phys. Chem. A, 2016, 120, 7385-7397. 57. Bader, R. F. W., Atoms in Molecule Oxford Science Publication: Oxford U. K. 1990. 58. AIMAll (Version 14.11.23), Todd A. Keith, TK Gristmill Software, Overland Park KS, USA, 2014 (aim.tkgristmill.com).

23 ACS Paragon Plus Environment

Page 24 of 35

Page 25 of 35 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

59. Johnson, E. R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A. J.; Yang, W. Revealing Noncovalent Interactions. J. Am. Chem. Soc., 2010, 132, 6498–6506. 60. Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem., 2012, 33 (5), 580–592. 61. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphics, 1996, 14, 33−38. 62. Wolinski, K.; Hinton, J. F.; Pulay, P. Efficient Implementation of the Gauge-Independent Atomic Orbital Method for NMR Chemical Shift Calculations. J. Am. Chem. Soc., 1990, 112, 8251–8260. 63. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999−3093. 64. Koch, U.; Popelier, P. L. A. Characterization of CHO Hydrogen Bonds on the Basis of the Charge Density. J. Phys. Chem., 1995, 99, 9747–9754. 65. Cremer, D., Kraka, E.: Chemical bonds without Bonding Electron Density-does the Difference Electron-Density Analysis Suffice for a Description of the Chemical Bond?, Angew. Chem., Int. Ed. Engl., 1984, 23, 627–628.

24 ACS Paragon Plus Environment

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 35

Figures: Lower Rim

Host Cavity

I

II

III

IV

Figure 1. Optimized geometry of the calix[4]tetrolarene and its derivatives.

25 ACS Paragon Plus Environment

Page 27 of 35 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

Host cavity

Lower rim

I

II

III

IV

Figure 2. Electron density isosurfaces (0.001 au) overlaid with MESP (±0.03 au)

26 ACS Paragon Plus Environment

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

Figure 3. The ρbcp parameters for O-H···O hydrogen bonds for host I to IV along with the fitted logarithmic relation is given. See text for details.

27 ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35 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 4. A plot of kinetic and potential electron density as a function of intramolecular hydrogen bond distances. See text for details

28 ACS Paragon Plus Environment

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

I

II

29 ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35 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

III

IV

Figure 5. Color-filled RDG isosurfaces depicting intramolecular interaction regions in the host I to IV.

30 ACS Paragon Plus Environment

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

Figure 6. OH vibration frequencies in the calculated infrared spectra of I-IV receptors.

31 ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35 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

HOMO

LUMO

I

II

III

IV

Figure 7. Frontier orbitals (± 0.04 a.u) in calix[4]tetrolarene and modified macrocycles.

32 ACS Paragon Plus Environment

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 34 of 35

Table 1. QTAIM parameters for the nonbonding interactions in I to IV macrocycles. Interactions

r

ρ

∇2 ρ

ε

V

G

H

O5-H5---O5

1.707

0.043

0.145019

0.067704

-0.04255

0.039401

-0.00315

O3-H3---O1

1.794

0.034

0.13001

0.081922

-0.03105

0.031776

0.000727

II

O5-H5---O5

1.931

0.025

0.096959

0.098211

-0.02036

0.022301

0.00194

O3-H3---O1

1.979

0.022

0.086768

0.110237

-0.01699

0.019341

0.00235

III

O5-H5a---O5

1.797

0.035

0.125284

0.052716

-0.03169

0.031506

-0.00018

O5-H5b--O5

1.877

0.029

0.115807

0.112546

-0.02595

0.027451

0.001502

O5-H5a---O5

2.073

0.018

0.069758

0.074231

-0.01333

0.015386

0.002053

2.257 0.014 0.052457 0.240044 and b refers to parallel and flattened position rings respectively.

-0.01027

0.011693

0.001421

Hosts I

IV

O5-H5b--O5

a

33 ACS Paragon Plus Environment

Page 35 of 35 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 2. 1H NMR chemical shifts (ppm) in calix[4]tetrolarene and their derivatives. δH H1 H2 H3 H5 H6 H7

I 6.11 5.61 8.28 10.83 3.80,3.56

II 6.20 5.45 8.47 8.17

34 ACS Paragon Plus Environment

III

IV

8.54 5.33, 3.44 3.75

7.50 3.73