Synthesis of Cycloveratrylene Macrocycles and Benzyl Oligomers Catalysed by Bentonite under Microwave/Infrared and Solvent-Free Conditions

Tonsil Actisil FF, which is a commercial bentonitic clay, promotes the formation of cycloveratrylene macrocycles and benzyl oligomers from the corresponding benzyl alcohols in good yields under microwave heating and infrared irradiation in the absence of solvent in both cases. The catalytic reaction is sensitive to the type of substituent on the aromatic ring. Thus, when benzyl alcohol was substituted with a methylenedioxy, two methoxy or three methoxy groups, a cyclooligomerisation process was induced. Unsubstituted, methyl and methoxy benzyl alcohols yielded linear oligomers. In addition, computational chemistry calculations were performed to establish a validated mechanistic pathway to explain the growth of the obtained linear oligomers.


Introduction
The development of general, catalytic and selective procedures for benzyl alcohol C-O bond activation is highly desirable and would constitute a broadly applicable set of transformations in organic synthesis [1]. Because a benzyl alcohol C-O bond is more reactive than that of an aliphatic alcohol, benzyl alcohol C-O activation is favoured in acid-catalysed reactions because of the greater stability of the intermediate formed [2,3]. This activation reaction has been used as a strategy in the synthesis of cyclotriveratrylenes [4] and oligotoluenes [5,6]-two groups of molecules with attractive chemical properties. Cyclotriveratrylenes (CTVs) [7][8][9][10] are cyclic molecular hosts that have been extensively employed in host-guest chemistry as a supramolecular scaffold [11][12][13][14][15]. CTVs can be prepared in three different ways: (a) by the acid-catalysed condensation of 1,2-disubstituted benzenes with formaldehyde; (b) through the condensation of diphenylmethane with 1,2-disubstituted benzenes and (c) most commonly through the use of dimethoxy-substituted benzyl alcohols under strongly acidic conditions (H 2 SO 4 /CH 3 COOH at 90 °C, H 3 PO 4 at 80 °C or Sc(OTf) 3 at 110 °C) [16][17][18][19][20][21]. Oligotoluenes are of industrial interest due to their application as insulating oils in high-voltage electrical devices [22], as reagents in the production of termiticide emulsions that demonstrate good penetration and emulsion stability [23] and as reagents in the preparation of corrosion-protection products [24]. The synthesis of these compounds has been developed following two different strategies: (a) from a mixture of toluene and benzyl chloride under ultrasound (20 MHz, 6 h) and thermal energy conditions (105 °C, 2 h) [5,6] or (b) from benzyl alcohol with carbon disulfide [2] or dichloromethane [3] used as solvents. Both processes were performed in the presence of bentonitic clay.
In recent years, a growth in the use of clays as reaction media for organic transformations has been observed. Natural and modified clays are safer, easy to handle and environmentally attractive compared with mineral acid solutions or metallic catalysts [25][26][27][28][29][30]. Bentonitic clays have been an excellent and versatile acidic Lewis-Brönsted [31] catalyst in the synthesis of Biginelli and Hantzsch esters [32], phenylmethanes [33,34], triphenylbenzenes and triphenylpyrylium salts [35], linear and cyclic ethers [36,37], thiocetals [38], anthracene and its radical cation [39] and in the generation of H 2 O 2 and hydroxyl radicals [40]. Presently, we desire to extend that knowledge to the development of a process that is more efficient, eliminates the use of solvent and halogenated starting materials, reduces reaction times and avoids tedious workup protocols. With these ideas in mind, we synthesised cycloveratrylene macrocycles and benzyl oligomers from the combination of Tonsil Actisil FF (TAFF), a commercial bentonite clay previously characterized by our group [5,6], with microwave-assisted organic synthesis (MAOS) [41][42][43][44][45] or infrared irradiation [46][47][48][49][50][51][52] in the absence of a solvent. In addition, computational chemistry calculations were performed to establish a validated mechanistic pathway to explain the growth of the obtained linear oligomers.

Results and Discussion
The initial experiments were performed with alcohols 1a-3a and TAFF under microwave and infrared irradiation; the results are summarised in Table 1. Thus, when compounds 1a-3a were treated with TAFF/MW at 85 °C (100 W) until the disappearance of the starting material after 1.50 min, the corresponding cyclic trimers, namely cyclotripiperonylene (1b, CTP, 85%), cyclotriveratrylene (2b, CTV, 90%) and 1,2,3,6,7,8,11,12,13-nonamethoxy-10,15-dihydro-5H-tribenzo[a,d,g]cyclononene (3b, NDTC, 80%), were obtained (Table 1, entries 1 and 3). These molecules are the condensation products of the corresponding benzylic cations. The cyclic trimers should proceed stepwise via the mono-, diand trimeric cation species and should be followed by a closure step to give the compounds 1b-3b [3]. The reactions with alcohols 1a and 2a invariably produced a variety of other minor products, which could be a disproportionation product (Table 1, compound 1c) or are a typical example of an acid-catalysed reaction of a primary alcohol to produce ethers (Table 1, compounds 1d, 1e, 2c). The formation of 1f might be explained by a subsequent electrophilic aromatic substitution reaction between the ether 1e and the respective benzyl cation. In the case of alcohol 3a, 1,2,3,6,7,8,11,12,13nonamethoxy-10,15-dihydro-5H-tribenzo[a,d,g]cyclononene (3b) was obtained as the major product, whereas the compounds 3c and 3d were the minor products. These products were formed (5% and 15%) as a consequence of the reaction conditions because, when the reaction was performed in carbon disulphide at reflux in the presence of TAFF for 7 h, 3b was obtained in lesser yield (7%) [36,37].
As shown in Table 1, when the reactions of alcohols 1a-3a were performed under infrared irradiation at 95 °C (375 W), similar results were obtained. The principle difference was the reaction time, which changed for each alcohol. Notwithstanding, a short irradiation time (3-7 min) was required for the complete disappearance of the starting material; this short time was ideal for our goal of developing a time-efficient process.
Next, we considered the possibility of using benzyl alcohols and testing the potential of TAFF/MW and TAFF/IR conditions for inducing the formation of benzyl oligomers. Thus, benzyl alcohol 4a yielded the oligomers 4c and 4d under MW conditions (Entry 1, Table 2). When the reaction was performed using infrared irradiation, oligomers 4c-e were formed (Entry 1, Table 2). In the same context, we also investigated the reactivity of 2-methylbenzyl alcohol (5a), which generated benzyl oligomers 5c-f when the reaction was performed under MW conditions and benzyl ethers oligomers 5g-i when the reaction was performed under infrared irradiation (Entry 2, Table 2). Finally, 2-methoxybenzyl alcohol (6a) yielded the benzyl oligomers 6c and 6d with both MW heating and infrared irradiation (Entry 3, Table 2). In all cases, short irradiation times were necessary for complete reaction of the starting materials (4-10 min).  Given that it is not practical to purify the reaction mixtures by conventional chromatographic methods (CC or HPLC), the presence of these benzyl oligomers in the reaction mixtures was determined by the analysis of GC-EIMS and HRMS spectra following the protocol previously described by our group [5,6]. Thus, when the reaction was analysed by GC-EIMS, each spectrum showed a set of fragments that were assigned unequivocally to the molecular ion of each group of isomers. In addition, each molecular ion was confirmed using high-resolution experiments; consequently, the corresponding elemental composition was obtained (Table 3). Notably, the respective elemental composition change was seven units for carbon and six units for hydrogen; this fact is in agreement with the difference of a benzyl moiety between the groups. We next considered a theoretical analysis to rationalise the observed products and some of the key reaction steps in the oligomerisation reaction of benzyl alcohol 4a. We performed a detailed calculation of the molecular structure and electronic properties of the ortho-, meta-and para-trimers, tetramers and pentamers (Table 4) and predicted the relative reactivities and regioselectivities using density functional theory (or Conceptual DFT). We used DFT to analyse the chemical potential, electronegativity, hardness and condensed Fukui functions for an electrophilic aromatic substitution.   4c-S1o 4c-S1m 4c-S1p 4d-S1oo 4d-S1om 4d-S1op 4d-S1mm 4d-S1mp 4d-S1pp 4d-S2om 4d-S2op 4d-S2mm 4e-S1ooo 4e-S1oom 4e-S1oop Table 4. Cont.
4e-S1omo 4e-S1omm 4e-S1omp 4e-S1opo 4e-S1opm 4e-S1opp 4e-S1mom 4e-S1mop 4e-S1mmm 4e-S1mmp 4e-S1mpm 4e-S1mpp 4e-S1pop 4e-S1pmp 4e-S1ppp 4e-S2ooo Thus, the geometry of the benzyl alcohol (4a) (monomer) was determined using B3LYP/6-311G(d,p). The calculated results of the most stable conformation reproduced the bond distances, C-C Ar , within 1.393-1.398 Å; this result is in agreement with the values previously reported for the same compound. Traetteberg and co-workers [53] obtained a d(C-C Ar ) value of 1.394 Å using the gas electron diffraction method. This molecule has a conformational gauche form with the OH group oriented toward the phenyl plane (not stable, Figure 1(a)) and a trans form with the OH group oriented away from the phenyl ring (not stable, Figure 1(b)). When the hydrogen atoms on CH 2 are in anti or syn orientations to the hydrogen atom of the OH group, the energy difference between the two states is 0.36 kcal/mol. Additionally, the benzyl alcohol has both the gauche (Figure 1(c)) and planar (Figure 1(d)) orientations of the OH group, but with hydrogen atoms in semiperpendicular orientations toward the phenyl plane. The energy difference between the two states is 1.64 kcal/mol. Both conformational gauche forms are stable, but the state with hydrogen atoms in a semiperpendicular orientation toward the phenyl plane is more stable, by 2.09 kcal/mol. The benzyl alcohol monomer in its stable form has dihedral angles CCCO = −32.2° and CCCOH = −58.0°. The minimum energy molecules of the trimer, tetramer and pentamer can be substituted at the ortho-, meta-and para-positions (Table 4). These systems should exist as different isomers, with each isomer in its respective conformation. The substituted isomers display either a regular or irregular symmetry distribution. All of these structures are stable. The calculated energies of these isomeric compounds are shown in Table 5. Of the three possible isomers of the trimer, the meta-isomer appears to be the most stable, followed by the para-isomers (the difference between the smallest and greatest is 1.19 kcal/mol). For the nine isomers of the tetramer, meta-isomers are more stable than para-isomers (the difference between the smallest and greatest is 3.56 kcal/mol). In the case of 39 isomers of the pentamer, the meta substitution is more stable than para-isomers (the difference between the smallest and greatest is 5.72 kcal/mol). The differences in the relative energies of the isomers of the trimers and tetramers are small and suggest that these isomers are almost identically stable. Additionally, a larger but still relatively small (given that they are conformers) difference in the energy is observed for the pentamers-isomer.  Atomic charge values of the dimer, trimer, tetramer and pentamer are depicted in Figure 2. For a particular atom type, differences in its atomic charge can be roughly correlated to differences in its nucleophilic power. The general pattern displayed by NPA [54][55][56][57][58][59][60][61][62][63][64][65][66][67] atomic charges collected in Figure 2 gives the carbon atom in the ortho isomer a negative value of approximately −0.208 e, although in this position, steric effects will be present. Table 5 shows the HOMO and LUMO energies and the interfrontier molecular orbital energy gaps, ∆E = E LUMO − E HOMO , for the considered isomers. The band gap energy values also provide an indication of the stability of a system. For large band gaps, a greater amount of excitation energy is needed to remove an electron from the valence band. From the difference between the smallest and greatest data to the band gap energy values, it is possible to see that the trimer isomers 4c shows a gap of 0.12 eV, the tetramer isomers 4d has a value of 0.13 eV, and the pentamer isomers 4e has a value of 0.2 eV. As evident from the results in Table 5, the ortho isomer has the highest HOMO/LUMO gap because of the reduced torsional potential and interaction between the substituted groups. When the substitutions are in the meta and para positions, the difference in the band gap energy increases with the number of rings; the gap energy is 0.06 eV for the trimer, 0.10 eV for the tetramer and 0.12 eV for the pentamer. Therefore, for an aromatic compound, the HOMO-LUMO gaps must be sufficiently large to prevent electron localisation. This observation is a theoretically correct answer according to the rules of aromaticity-the smaller the HOMO-LUMO band gap, the less aromatic the system. Thus, a linear molecule is a non-aromatic, non-ring structure. If the substitutions are in the ortho-and parapositions, the difference in the band gap values between them is 0.12 eV for the trimers and tetramers and 0.16 eV for pentamers. The B3LYP/6-311G(d,p)-calculated HOMO, LUMO and gap energies of the studied oligomers are shown in Figure 3. The 4e-S3omp pentamer has the lowest band-gap energy among this group at 5.75 eV, which means that this molecule is the most reactive of all the structures (Table 5). Figure 4 shows the HOMO and LUMO contour molecular orbitals of the more stable oligomers.

Benzyl alcohol 4a
Dimer 4b Trimer 4c-S1m Tetramer 4d-S2mm Pentamer 4e-S2mmm The results of the calculations of the reactivity indices, ionisation potential (IP), electron affinity (EA), electronegativity (χ) (as the negative of the chemical potential, µ), chemical potential and hardness (η) for all molecules investigated in this work, as obtained with the B3LYP/6-311G(d,p) model, are presented in Table 5. The ionisation potential of a stable molecule is always positive. An inspection of Table 5 reveals that all molecules present a positive value of the ionisation potential (trimer 7.76-8.01 eV, tetramer 7.49-7.78 eV and pentamer 7.34-7.67 eV). Thus, this result is an indication of the stability of the set of molecules. All of the structures present a negative value for the electron affinity, which indicates that energy is released when a neutral species becomes an anion. The difference between the smallest and greatest energies are −0.84 to 0.87 eV for the trimer, −0.66 to −0.74 eV for the tetramer and −0.53 to −0.62 eV for the pentamer. The released energy decreased with an increase in the number of rings. The monomer has the greatest value of electron affinity, which indicates that the anion is more stable relative to the neutral system. The electronegativity shows how the electrons will flow from regions of high electronic density in a molecule to other sites of lower electronic density. As such, this value is an important index of reactivity for a given system. As evident from the results in Table 5, the electronegativities are given in the following order: 3.6 eV (monomer) > 3.51 eV (dimer) ≈ 3.58-3.44 eV (trimer) ≈ 3.55-3.47 eV (tetramer) ≈ 3.53-3.48 (pentamer). These results indicate that the monomer is more prone to attract electrons during the interaction with another chemical compound. However, we note that the chemical potential of the oligomers derivatives is localised in the following intervals: trimers −3.44 to −3.58 eV, tetramers −3.40 to −3.86 eV and pentamers −3.37 to −3.56 eV. The smallest chemical potential corresponds to substitutions in either ortho or meta positions; this result clarifies that the flow of charge transfer is from the para position to the electrophile during an electrophilic aromatic substitution process. The global hardness of the oligomers, given as the difference between the smallest and largest values, is 4.31-4.43 eV for the trimers, 3.87-4.23 eV for the tetramers and 3.96-4.11 eV for the pentamers.
According to the principle of maximum hardness, more reactive systems will show low hardness values, and less reactive systems will show high hardness values [58]. Molecules arrange themselves to maximise hardness. A high value of chemical hardness indicates high kinetic stability and low reactivity, and, thus, this parameter was found to be a cardinal index for molecular structure, bonding and reactivity. When we consider that the stability of the species is directly related to its hardness, then the stability of the pentasubstituted isomers is lower compared to those of the tetrasubstituted and trisubstituted isomers. Among the trisubstituted, tetrasubstituted and pentasubstituted isomers, the ortho isomer was generally found to be the hardest, and the meta-and para-isomers were the softest. Thus, the hardness measure indicates that the meta-and para-isomers are the most reactive positional isomers.
To complement the charge analysis and find the active sites of the molecules, we calculated the Fukui functions. The results of the calculations of the condensed Fukui functions, f -, for the benzylic alcohol monomers and for some oligomers are presented in Table 6. The benzylic alcohol is more active towards an electrophilic attack through atom Cp (carbon atom in para-position). The same atom is inactive for a radical attack or a nucleophilic attack. The dimer has the same activity for either of the two attacks at its para carbon atom. The oligomer increases the number of benzylic alcohol units, and its Fukui index decreases. For the other carbon atoms, f − decreased. The charge and HOMO-LUMO molecular orbitals allow us to establish that the growth of the oligomer could occur through the meta-position, whereas the Fukui indices indicate that the reactivity is at the para-position. According to Figure 5 and Table6, para-carbon in oligomers is the most nucleophilic.

Benzyl alcohol 4a Dimer 4b
Trimer 4c-S1m Tetramer 4d-S2mm Pentamer 4e-S2mmm Scheme 1. Proposed reaction mechanism for the oligomerisation process. Finally, from a mechanistic point of view, the catalytic action of TAFF should enhance the electrophilic character of the benzyl alcohol and facilitate the electrophilic aromatic substitution reaction that yields the dimer 4b, which is the key intermediate to understand the formation of oligomers 4c-e. The interaction between TAFF and benzyl alcohol might be due to the protonated and unprotonated active sites that correspond to the acidic Lewis character of the clay (Scheme 1).

Materials
Tonsil Actisil FF, is cheap (US $1. 7. X-ray thermodiffractograms show that the laminar structure is unstable above 150 °C. Quartz and cristobalite are also important components in the clay composition, as observed by X-ray diffraction. The corresponding BET surface area is 198.718 m 2 g −1 , and the pore volume and average pore diameter are 32.04 × 10 −2 cm 3 g −1 and 77.8 Å, respectively. It is worth mentioning that a detailed characterization of the clay ( 29 Si and 27 Al MAS-NMR, SEM, IR-Py, DTA, TG, and Ho) has already been performed and reported by our research group [5,6]. The particle size is 325 mesh. The benzyl alcohols were purchased from Aldrich and were used without further purification. The solvents were also acquired from Aldrich and were purified by standard methods prior to use. All 1 H-and 13 C-NMR spectra were recorded on a Varian Gemini (300 MHz) spectrometer using CDCl 3 as solvent and TMS as an internal reference. EIMS (70 eV) spectra, HRMS data and GC-MS analysis were obtained using a JEOL JMS AX505HA mass spectrometer. The product distribution was determined by a Varian gas chromatograph Stard 3400 equipped with flame ionization detector and a 30m × 0.53 mm column packed with polyethylene glycol, the relative proportion of the products were calculated assuming that the detector gave equal response for each compound. Thin-layer chromatographic analyses were performed using Merck silica gel 60 F 254 (0.25 mm) pre-coated plates, while products were purified on flash chromatographic columns of silica gel 60 (70−230 mesh). Microwave reactions were performed in a CEM Discover LabMate instrument. The middle infrared irradiation was performed using a Phillips IR lamp (375 W/220 V) integrated to an infrared reactor designed by our research group [46][47][48], and validated by a wide number of applications.

Typical Procedure for the Catalytic Reaction under Microwave Conditions
A mixture of the benzyl alcohols 1a-6a (2 mmol) and Tonsil Actisil FF (20 mg) were reacted in a CEM Discover LabMate in a sealed vessel at 85 °C (100 W and 1 bar of pressure) until disappearance of the starting material. The reaction was conveniently monitored by TLC and the experiments were repeated in various times. Then, the clay was removed by filtration through Celite and washed with ethyl acetate (3 × 5 mL). The combined filtrates were dried on anhydrous Na 2 SO 4 , and the solvent was eliminated under reduced pressure. To the benzyl alcohols 1a-3a, the residue was subjected to chromatography on a silica-gel column using n-C 6 H 14 /EtOAc as the eluent, affording the compounds 1b-f, 2b-c and 3b-d. To the benzyl alcohols 4a-6a, the presence of benzyl oligomers in the reaction mixtures was determinate by the analysis of GC-EIMS and HRMS spectra, followed the protocol previously described for our group [5,6]. Table 3 shows the elemental composition and high resolution data of oligotoluenes.

Typical Procedure for the Catalytic Reaction under Infrared Conditions
A mixture of the benzyl alcohols 1a-6a (2 mmol) and Tonsil Actisil FF (20 mg) were thoroughly mixed in a round-bottomed flask (10 mL). The mixtures were exposed to infrared irradiation with an infrared lamp at 95 °C (375 W) until disappearance of the starting material. The reaction was conveniently monitored by TLC. Then, the clay was removed by filtration through Celite and washed with ethyl acetate (3 × 5 mL). The combined filtrates were dried on anhydrous Na 2 SO 4 , and the solvent was eliminated under reduced pressure. To the benzyl alcohols 1a-3a, the residue was subjected to chromatography on a silica-gel column using n-C 6 H 14 /EtOAc as the eluent, affording the compounds 1b-f, 2b-c and 3b-3d. To the benzyl alcohols 4a-6a, the presence of benzyl oligomers in the reaction mixtures was determinate by the analysis of GC-EIMS and HRMS spectra (Table 1), followed the protocol previously described for our group. Table 3 shows the elemental composition and high resolution data of oligotoluenes.

Synthesis of Cyclotriveratrylene (2b) [4]
Following the general procedures, the reaction was carried out starting from (3,4- In this way, compound 2c was also obtained. Following the general procedures, the reaction was carried out starting from (3,4,5-trimethoxyphenyl) methanol (3a, 396 mg, 2 mmol) and Tonsil Actisil FF (20 mg) under solvent-free conditions. When the reaction was finished (MW = 5 min, IR = 7 min) the reaction crude was purified by flash column chromatography on silica gel using n-C 6

Synthesis of Oligotoluenes 4c-e
Following the general procedures, the reaction was carried out starting from phenylmethanol (4a, 216 mg, 2 mmol) and Tonsil Actisil FF (20 mg) under solvent-free conditions. When the reaction was finished (MW = 10 min, IR = 10 min) the presence of benzyl oligomers 4c-e in the reaction mixture was determined by GC-EIMS analysis and HRMS spectroscopy.
Synthesis of Oligotoluenes 5c-i: Following the general procedures, the reaction was carried out starting from (2-methylphenyl)methanol (5a, 244 mg, 2 mmol) and Tonsil Actisil FF (20 mg) under solvent-free conditions. When the reaction was finished (MW = 5 min, IR = 7.5 min) the presence of benzyl oligomers 5c-i in the reaction mixture was determined by GC-EIMS analysis and HRMS spectroscopy.
Synthesis of Oligotoluenes 6c-d: Following the general procedures, the reaction was carried out starting from (2-methoxylphenyl)methanol (6a, 276 mg, 2 mmol) and Tonsil Actisil FF (20 mg) under solvent-free conditions. When the reaction was finished (MW = 4 min, IR = 5 min) the presence of benzyl oligomers 6c-d in the reaction mixture was determinate by GC-EIMS analysis and HRMS spectroscopy (Table 3).

Theoretical Calculations
Calculations were performed using the Gaussian 09 system [61]. The geometries of the isolated benzyl alcohol molecule and its dimer, trimers, tetramers and pentamers complexes were fully optimized at the B3LYP [62,63] level using the 6-311G(d,p) basis set [64,65]. Harmonic vibrational frequency calculations have been performed at the same computational level to confirm that the structures are energetic minima. In the present work we analyze the theoretical structure and energetic of the ortho, meta, and para trimmers, tetramers and pentamers (Table 4). Natural population analyses (NPA) were also performed based on the optimized geometry [66]. We have also examined HOMO and LUMO levels; the energy gap is evaluated as the difference between the HOMO and LUMO energies. The quantitative definitions for chemical potential (µ) [67] and chemical hardness (η) [68,69] for an N-electron system with total energy E can respectively be given as: and: where IP and EA are ionization potential and electron affinity of the system, respectively. Using the ∆SCF finite difference approach, we can calculate the IP and EA for the N-electron system as follows: kui functions are common descriptors of site reactivity [70][71][72]. They are defined as the derivative of the electron density with respect to the total number of electrons N in the system, at the constant external potential. : The condensed Fukui functions can also be employed to determine the reactivity of each atom in the molecule. The corresponding condensed functions are given by: (for nucleophilic attack) (for electrophilic attack), and (for radical attack), where is the gross charge of atom in the molecule. A high value of implies a high reactivity of that site . Besides, the type of condensed Fukui function whose value is highest at a particular site, predicts the type of attack that predominates at that site.

Conclusions
In summary, we have developed an efficient and solvent-free protocol for the synthesis of cycloveratrylene macrocycles and benzyl oligomers under mild conditions in the presence of Tonsil Actisil FF (TAFF), a commercial bentonite clay, under microwave heating and infrared irradiation. Compared with the traditional method, this method possesses at least four advantages: (i) the energy supply (MW and IR) is harmless; (ii) the reaction was performed in the absence of carbon disulphide or dichloromethane, which means that no volatile organic compounds (VOCs) were used in the present protocol; (iii) the reaction times are short; and (iv) the use of any halogenated starting materials, metal compounds or corrosive acids that may cause pollution are avoided. The only catalyst used in this reaction is bentonite clay, which can be easily removed from the reaction mixture. Notably, the catalytic reaction is sensitive to the type of substituent on the aromatic ring. Thus, when benzyl alcohol was substituted with a methylenedioxy, two methoxy or three methoxy groups, a cyclooligomerisation process was induced. Unsubstituted, methyl-and methoxybenzyl alcohols yielded linear oligomers.
In addition, to predict the relative reactivities and regioselectivities of the oligomerisation processes of benzyl alcohol, we calculated the energies and optimised the structures of trimer, tetramer and pentamer isomers. The meta-isomer appears to be more stable than the para-isomers. After the reactivity descriptors calculated for the oligomer derivatives were analysed to identify the system with the lowest band gap and the largest ionisation potential, we calculated these values for the penta-oligomers. The penta-oligomers also showed a lower electron affinity. The chemical hardness indicates that as the number of benzylic alcohol units is increased, the oligomer becomes less stable. These molecules are more reactive. Additionally, the position of substitution is important. The electronegativity of the oligomer derivatives is almost the same. Finally, the charge indicates that reactivity is highest at the ortho position but also that steric effects are present. The HOMO-LUMO energy gap allows us to establish that the growth of the oligomer could occur at the meta-or paraposition. The Fukui indices, however, indicate that the highest reactivity is at the para-position.