Efﬁcient Molybdenum Hydrazonato Epoxidation Catalysts Operating under Green Chemistry Conditions: Water vs. Decane Competition

: Molybdenum compounds containing benzaldehyde-based hydrazones were obtained. The reaction in MeOH resulted with monomeric Mo complexes, [MoO 2 (L)(MeOH)], while the reaction in dichloromethane (DCM) provided oligomeric complexes, [MoO 2 (L)] n . The solid-state structures of the obtained compounds were investigated through Infrared Spectroscopy - Attenuated Total Reﬂec-tion (IR-ATR), Thermogravimetric analysis (TGA), and via X-ray diffraction. The prepared molybdenum species were employed as cyclooctene epoxidation catalysts. TBHP ( tert -butylhydroperoxide) in water and TBHP in decane were employed and compared as oxidants, with 0.25 mol% [Mo]. The catalyst activity and selectivity towards epoxide is >90% for all the reactions. The results have been linked to theoretical calculations, showing the importance of the ﬁrst step, i.e., the transformation of [MoO 2 (L)(MeOH)] into the pentacoordinate [MoO 2 (L)]. found 6.8%, calc. for MoO 3 32.4%, found 31.5%. Selected IR data (cm − 1 ): 1623 (C=N), 1378 (C–Ophenolate), 1345 (C–Oenolate), 942, 927 (MoO 2 ), 894 (O=Mo–O MeOH ).


Introduction
Epoxides are known to be crucial raw materials widely employed in everyday life. One of the main routes for epoxide production is alkene epoxidation, where alkene is converted into epoxide when subjected to a peroxy acid of general formula RCOOOH. Peroxy acids usually applied are peroxy trifluoroacetic acid, peroxyacetic acid, and metachloroperoxybenzoic acid [1]. From an environmental point of view, the introduced pathway bears severe drawbacks (i.e., long and tedious post-treatment and high generation of waste). Nowadays, TBHP (tert-butylhydroperoxide) makes a significant contribution, especially at the industrial scale, as it produces less waste and the byproduct formed upon the reaction can be reused [2,3]. For metal-catalyzed epoxidations with alkyl hydroperoxides (ROOH) used as oxidants, highly active catalysts show Lewis acidity and low oxidation potentials in their highest oxidation states. Their relative activity is usually Mo > W > Ti~V. The commercially available molybdenum compounds [MoO 2 (acac) 2 ] (acac = C 5 H 7 O 2 ) and Mo(CO) 6 have been tested as very active and selective epoxidation catalysts [4,5]. Within the development of the research in this area, the systems containing {MoO 2 } 2+ core and coordinated Schiff base ligand showed a great epoxidation potential [6,7]. In most of the reactions dealing with such Mo complexes, 1,2-dichloroethane (DCE) or other chlorinated solvents were the best solvents for the reaction. Focus is also done on the minimization or complete elimination of the solvents used as reaction media, following the green chemistry policies and solvent-less chemistry. For instance, the Montreal Protocol on substances depleting the ozone layer is a global agreement to protect the Earth's ozone layer by phasing green chemistry policies and solvent-less chemistry. For instance, the Montreal Protocol on substances depleting the ozone layer is a global agreement to protect the Earth's ozone layer by phasing out the chemicals that diminish it. Regarding it carbon tetrachloride was substituted in 1989 by the halogenated solvents chloroform and dichloromethane (DCM) [8]. Since DCM today is recognized to have cancerogenic properties and to deplete the ozone layer, the idea is to replace DCM with other, less harmful solvents [9]. Due to it, modern solvent selection tools, driven by Pfizer and Sanofi, classify water and alcohol as preferred and recommended solvents [10].
In the proposed quest, our research follow-up was to obtain efficient molybdenum epoxidation catalysts containing benzaldehyde-based hydrazones (Scheme 1). The choice of ligands lies in the fact that the similar molybdenum complexes showed good activity and were classified as efficient catalytic systems [11][12][13]. In the presented research, the stress was on the use of 2-hydroxy-benzhydrazide derivatives, that could provide catalytically active molybdenum complexes. In addition, the presence/position and position of the methoxy group on the ligand has been a studied factor, to see its influence on the catalytic activity. The catalytic results have been compared to theoretical calculations, linking the activity to a mechanistic study. Scheme 1. Ligands H2L H, 3OMe , and 4OMe employed for the complex and catalyst preparation.

Catalysts Preparation and Solid-State Characterization
All the ligands were prepared by the conventional procedure, in methanol or ethanol. With the intention of synthetic improvement, a solvent-less technique was applied as well. The stoichiometric amounts of appropriate aldehyde (salicylaldehyde, 3-methoxy-or 4 methoxy-salicylaldehyde, and 2-hydroxy-benzhydrazide were employed, and liquid assisted grinding (LAG) was applied. The existence of the same ligands, prepared by different routes, was confirmed by IR-ATR and PXRD.

Catalysts Preparation and Solid-State Characterization
All the ligands were prepared by the conventional procedure, in methanol or ethanol. With the intention of synthetic improvement, a solvent-less technique was applied as well. The stoichiometric amounts of appropriate aldehyde (salicylaldehyde, 3-methoxyor 4 methoxy-salicylaldehyde, and 2-hydroxy-benzhydrazide were employed, and liquid assisted grinding (LAG) was applied. The existence of the same ligands, prepared by different routes, was confirmed by IR-ATR and PXRD.  [14,15]. The appearance of the stretching frequencies characteristic of ν(MoO 2 ) (found at 940 cm −1 , 917 cm −1 and 900 cm −1 ) characterized unambiguously the monomeric compounds (Supplementary Materials Figure S1). In addition, the bands typical for bands characterizing the N-H and C=O groups of the ligand in the neutral form, at around 3190 cm −1 and 1662 cm −1 , respectively) are absent in the spectra of molybdenum complexes, implying keto-enol tautomerization and ligand coordination to the molybdenum center through deprotonated oxygen atoms. Addition support for this hypothesis is a band at 1345 cm −1 . The vibration band belonging to C=N imine of the neutral ligands at 1650 cm −1 was shifted to 1625 cm −1 in the spectra of the molybdenum complexes.
Oligomeric complexes [MoO2L]n were easily distinguished from monomeric by the simple comparison of the IR-ATR spectra. The wide absorption band at around 810 cm -1 , attributed to Mo=O … Mo polymerization mode, vanished in the spectra of the corresponding monomeric complex [MoO2L(MeOH)] [14,15]. The appearance of the stretching frequencies characteristic of ν(MoO2) (found at 940 cm -1 , 917 cm -1 and 900 cm -1 ) characterized unambiguously the monomeric compounds (Supplementary Materials FigureS1). In addition, the bands typical for bands characterizing the N-H and C=O groups of the ligand in the neutral form, at around 3190 cm -1 and 1662 cm -1 , respectively) are absent in the spectra of molybdenum complexes, implying keto-enol tautomerization and ligand coordination to the molybdenum center through deprotonated oxygen atoms. Addition support for this hypothesis is a band at 1345 cm -1 . The vibration band belonging to C=Nimine of the neutral ligands at 1650 cm -1 was shifted to 1625 cm -1 in the spectra of the molybdenum complexes.
The  [18] were already reported in the literature. The studied  [17] and [MoO 2 (L 4OMe )(MeOH)]·MeOH [18] were already reported in the literature. The studied complexes exhibit robust coordinative and geometrical features which were preserved throughout the whole class of hydrazonato dioxomolybdenum(VI) species. The cis-dioxomolybdenum {MoO 2 } 2+ core was found to be coordinated by doubly deprotonated ligand in its hydrazidato form if considering hydrazone functional group (Supplementary Materials Table S2). The coordination is achieved through a geometrically convenient ONO coordination pocket built from aryl 2-oxo group, imine nitrogen atom, and hydrazone group oxygen atom (Figure 1a,b). The sixth coordination place of {MoO 2 } 2+ core was coordinated by methanol hydroxy group, even when competing for nitrogen auxiliary is introduced in case of [MoO 2 (L 4OMe )(MeOH)]·4,4-bpy. Molecules of complexes were connected via hydrogen bonding into supramolecular dimers (Figure 1c,d), which were furthermore densely packed, forming crystal structures with no solvent-accessible voids found (Supplementary Materials Figures S4 and S5). All aforementioned observations were in line with previously reported structures of similar compounds. However, it was beneficial to compare the geometrical features of complexes and try to connect them with observed performance, to build knowledge for the development of materials with desirable properties (Supplementary Materials Figure S6). If one considers octahedral distortion of {MoO 2 } 2+ core and planarity of the coordinated ligand as relevant geometrical parameters, it can be seen (Table 1) that these parameters significantly deviate depending on the ligand substituent and solvent coordinated at the sixth coordination place, even though coordination geometry remains relatively similar (Supplementary Materials Table S1). Indeed, this sample of crystal structures was too small for definitive structure-property correlation, but gives additional motivation for further investigation of this type of complexes. thermore densely packed, forming crystal structures with no solvent-accessible voids found (Supplementary Materials Figures S4 and S5). All aforementioned observations were in line with previously reported structures of similar compounds. However, it was beneficial to compare the geometrical features of complexes and try to connect them with observed performance, to build knowledge for the development of materials with desirable properties (Supplementary Materials FigureS6). If one considers octahedral distortion of {MoO2} 2+ core and planarity of the coordinated ligand as relevant geometrical parameters, it can be seen (Table 1) that these parameters significantly deviate depending on the ligand substituent and solvent coordinated at the sixth coordination place, even though coordination geometry remains relatively similar (Supplementary Materials Table S1). Indeed, this sample of crystal structures was too small for definitive structure-property correlation, but gives additional motivation for further investigation of this type of complexes.

Catalytic Results
All obtained, monomeric and oligomeric complexes, were used as epoxidation catalysts, using cyclooctene as a substrate, and tested oxidants were TBHP in decane or TBHP in water. No other organic solvents were added to the reaction mixture, justifying the green chemistry protocols. Kinetic profiles for all the complexes are presented in Figure 2 and the obtained results are compiled in Table 2.

Catalytic Results
All obtained, monomeric and oligomeric complexes, were used as epoxidation catalysts, using cyclooctene as a substrate, and tested oxidants were TBHP in decane or TBHP in water. No other organic solvents were added to the reaction mixture, justifying the green chemistry protocols. Kinetic profiles for all the complexes are presented in Figure 2 and the obtained results are compiled in Table 2.  As previously observed with other pentacoordinate [MoO2L] complexes bearing ligands with very close geometries than those studied here, TBHP in decane provided better results in terms of conversion and selectivity than TBHP in water. However, for all the complexes used, the cyclooctene conversion after 6 h was almost complete, and selectivity towards the epoxide being higher than 90%. At the difference with previously published results, it was surprising to observe that the process using TBHP in water as an oxidant gave herein excellent results with the studied complexes. Indeed, all complexes showed conversion of cyclooctene in the range 90-99%, the monomeric complex [MoO2(L 4OMe )(MeOH)] giving the highest conversion and the oligomeric complex  As previously observed with other pentacoordinate [MoO 2 L] complexes bearing ligands with very close geometries than those studied here, TBHP in decane provided better results in terms of conversion and selectivity than TBHP in water. However, for all the complexes used, the cyclooctene conversion after 6 h was almost complete, and selectivity towards the epoxide being higher than 90%. At the difference with previously published results, it was surprising to observe that the process using TBHP in water as an oxidant gave herein excellent results with the studied complexes. Indeed, all complexes showed conver- . Thus, to compare the activity from both complexes, cyclooctene conversion and epoxide production were followed in the first hour of the reaction and are presented in Figure 3. As seen, even if induction periods can be noticed for both catalysts in the first 20 min of the reaction, the monomeric catalyst [MoO 2 (L 3OMe )(MeOH)] became more active in a shorter time. In relationship with the postulated mechanism, experimental data confirmed that the energy needed to decoordinate the solvent seems lower than for the oligomeric compound. This was observed already with similar species.
relative to the oxidant carrier (water vs. decane). Highest TOFs observed for [MoO2(L 3OMe )(MeOH)] and [MoO2(L 3OMe )]n show that the catalytic activity is ligand-dependent. Another factor to discuss comes from the complete proposed postulated mechanism, where the active species was a pentacoordinate species [MoO2L] and its formation had to be reached from the complexes added into the reaction media. According to it, [MoO2(L 3OMe )(MeOH)] and [MoO2(L 3OMe )]n might have the same active species, i.e., [MoO2(L 3OMe )]. Thus, to compare the activity from both complexes, cyclooctene conversion and epoxide production were followed in the first hour of the reaction and are presented in Figure 3. As seen, even if induction periods can be noticed for both catalysts in the first 20 min of the reaction, the monomeric catalyst [MoO2(L 3OMe )(MeOH)] became more active in a shorter time. In relationship with the postulated mechanism, experimental data confirmed that the energy needed to decoordinate the solvent seems lower than for the oligomeric compound. This was observed already with similar species. The reported research with the monomeric complexes [MoO2(L)(EtOH)], where L represents the coordinated ligand obtained from 4-hydroxybenzhydrazide, showed great activity of the catalyst, but moderate selectivity towards epoxide, when using TBHP in water [11]. This implied a positive effect of the 2-hydroxybenzhydrazide, applied in this research, and easier formation of the pentacoordinate active species, [MoO2(L)], that participated in the catalytic cycle. TOF20min values for the reported complexes varied between 350 and 469, while TON parameter was in the range 354 and 368, similar to the results reported herein. Furthermore, similar investigations employing polymeric catalysts, [MoO2(L)]n, obtained from the ligands containing nicotinic acid hydrazide moiety, showed much lower activity and selectivity towards epoxide with TBHP in water, while by applying TBHP in decane the catalytic results improved drastically.

Theoretical Considerations
The experimental results show that the 3-methoxy group on the ligand seems to enhance considerably the activity of the molybdenum complexes, with the best TOFs observed among the three ligands studied (for monomeric and oligomeric complexes). In addition, decane was proved earlier to be a superior TBHP carrier than water and to favor the reactivity of the Mo complexes. The presence of water (from aqueous TBHP) can potentially inhibit the catalytic system, but the observed conversion difference was minimal in the present investigation. Since activity was quite high and conversion important, such behavior was visible here only within the TOF where the discrepancy of values between ligands was higher for TBHP in water than for TBHP in decane. As a reminder, the postulated mechanism and a corresponding DFT study gave interesting insight. Previously, The reported research with the monomeric complexes [MoO 2 (L)(EtOH)], where L represents the coordinated ligand obtained from 4-hydroxybenzhydrazide, showed great activity of the catalyst, but moderate selectivity towards epoxide, when using TBHP in water [11]. This implied a positive effect of the 2-hydroxybenzhydrazide, applied in this research, and easier formation of the pentacoordinate active species, [MoO 2 (L)], that participated in the catalytic cycle. TOF 20min values for the reported complexes varied between 350 and 469, while TON parameter was in the range 354 and 368, similar to the results reported herein. Furthermore, similar investigations employing polymeric catalysts, [MoO 2 (L)] n , obtained from the ligands containing nicotinic acid hydrazide moiety, showed much lower activity and selectivity towards epoxide with TBHP in water, while by applying TBHP in decane the catalytic results improved drastically.

Theoretical Considerations
The experimental results show that the 3-methoxy group on the ligand seems to enhance considerably the activity of the molybdenum complexes, with the best TOFs observed among the three ligands studied (for monomeric and oligomeric complexes). In addition, decane was proved earlier to be a superior TBHP carrier than water and to favor the reactivity of the Mo complexes. The presence of water (from aqueous TBHP) can potentially inhibit the catalytic system, but the observed conversion difference was minimal in the present investigation. Since activity was quite high and conversion important, such behavior was visible here only within the TOF where the discrepancy of values between ligands was higher for TBHP in water than for TBHP in decane. As a reminder, the postulated mechanism and a corresponding DFT study gave interesting insight. Previously, several Mo complexes have been examined through a systematic DFT study. As an example, on the [MoO 2 (SAP)] complexes and different substitutions on ligand showed that the presence of OH group (on a specific position on the aldehyde part) [19] or NO 2 substituent (on aminophenol part) [20] on the ligand provoked higher activity. The presence of 3-and 4-methoxy groups on the ligand was studied several times with other Mo complexes [21]. The trend observed experimentally among similar complexes could be confirmed theoretically in all the cases, but it was not possible to establish a general rule linking the activity of the complexes and the methoxy position/presence. Thus, it was interesting to study the presented complexes based on the previous works to see the relationship between experiment and theory among the three pentacoordinate species. Since those complexes were very efficient under the same conditions applied as in previous studies, the research interest was to correlate the activity and thermodynamic values, for ex-ample, the TS values. To be consistent, calculations were limited to the part responsible for the oxygen transfer, i.e., starting from the monomeric species [MoO 2 (L)(MeOH)]. Thus, the different studied steps include the complex desolvation to reach the pentacoordinate species (step 0), the formation of the [MoO 2 (L)(TBHP)] adduct (step I), and the approach of the olefin to the O β (step II) that was implicated in the epoxide formation. The calculations were limited to the conditions in the gas phase since it was demonstrated earlier that trends obtained through PCM corrections and gas-phase were similar. Those data are gathered in Table 3. The adducts and the TS were indicated in Scheme 3. All optimized geometries are collected in Supplementary Materials Tables S4-S6. several Mo complexes have been examined through a systematic DFT study. As an example, on the [MoO2(SAP)] complexes and different substitutions on ligand showed that the presence of OH group (on a specific position on the aldehyde part) [19] or NO2 substituent (on aminophenol part) [20] on the ligand provoked higher activity. The presence of 3-and 4-methoxy groups on the ligand was studied several times with other Mo complexes [21]. The trend observed experimentally among similar complexes could be confirmed theoretically in all the cases, but it was not possible to establish a general rule linking the activity of the complexes and the methoxy position/presence. Thus, it was interesting to study the presented complexes based on the previous works to see the relationship between experiment and theory among the three pentacoordinate species. Since those complexes were very efficient under the same conditions applied as in previous studies, the research interest was to correlate the activity and thermodynamic values, for example, the TS values.
To be consistent, calculations were limited to the part responsible for the oxygen transfer, i.e., starting from the monomeric species [MoO2(L)(MeOH)]. Thus, the different studied steps include the complex desolvation to reach the pentacoordinate species (step 0), the formation of the [MoO2(L)(TBHP)] adduct (step I), and the approach of the olefin to the Oβ (step II) that was implicated in the epoxide formation. The calculations were limited to the conditions in the gas phase since it was demonstrated earlier that trends obtained through PCM corrections and gas-phase were similar. Those data are gathered in Table 3. The adducts and the TS were indicated in Scheme 3. All optimized geometries are collected in Supplementary Materials Tables S4-S6.  From the absolute values obtained by calculations for each step, the decoordination of the methanol (step 0) is more favored with [MoO 2 (L 4OMe )(MeOH)] complex, which can be as well seen experimentally through TG analysis. The stabilization of the pentacoordinate species [MoO 2 (L)] with TBHP (step I) was stronger in the case of [MoO 2 (L 3OMe )] (8.9 kcal mol −1 ). The TS energy (step II) was smaller for [MoO 2 (L H )] with 23.0 kcal mol −1 , i.e., around 0.2-0.5 kcal mol −1 more favorable than the two other complexes. Steps III and IV were not important for the discussion but are indicated in the table. The experimental TOF values showed that complexes bearing [MoO 2 (L 3OMe )] are the fastest. This observation could be done if the energy necessary to form the pentacoordinate species (step 0) is added, i.e., the energy needed for the methanol decoordination. Thus, the TS value was the lowest (in relative value) for [MoO 2 (L 3OMe )]. The difference is very low, and it has to be seen that the important discussion is about the TOF and not the better activity since it has been shown that there is not a high difference of activity between the three ligands. Those data can be linked to the geometrical data.
The Table 4 gathered relevant distances and one dihedral angle (related to ligand planarity) around the molybdenum atom, with numbering according to Scheme 4. Experimental results showed that the [MoO 2 (L 3OMe )] species seemed to be the best catalyst. Some geometrical parameters were discussed. Thus, as seen in energy diagram, Mo-O β interatomic distance parameter influenced the TS path mainly (step II). The shorter this distance, the smaller the TS. This factor favored the path with [MoO 2 (L H )]. This was certainly linked to electronic and steric effects. Interesting facts appeared with the two nitrogen atoms, i.e., N 1 (linked to Mo) and N 2 (linked to the hydrogen of the pending OH moiety). In all the steps, although if a slight difference was observed, the N 2 . . . H bond was shorter for [MoO 2 (L 3OMe )] catalyst. At the difference, the impact on the N 1 . . . Mo interatomic distance was more visible with the longest distance with the same catalyst. This seems to give an additional clue of the electronic effect and subsequently to the mechanism, showing that a stronger Mo stabilization through N 1 atom does not favor the reactivity, and this stabilization is influenced by the N 2 . . . H stabilization. Of note, the MoO 2 fragment was mostly identical with the three ligands. Table 4. Interatomic distances (Å) and dihedral angle ( • ) for the calculated structures according to Scheme 4.

Synthesis of Ligands
Ligand preparation by solution-based synthesis: Ligands H2L H , H2L 3OMe , H2L 4OMe and were prepared by condensation of aldehyde and hydrazide, according to the reported procedure [24][25][26]. Ligand preparation by mechanochemical synthesis: Liquid assisted grinding (LAG) experiments were performed using a Retsch MM200 ball mill. 0.1044 g (0.7 mmol) of benzaldehyde, 3-or 4-methoxysalicylaldehyde, 0.1041 g (0.7 mmol) 2-hydroxy-benzydrazide, and 50 μL of methanol were placed with one 10 mm milling ball in a 10 mL Teflon jar. For the preparation of H2L H no solvent was added to the mixture of salicylaldehyde, 0.085 g (0.7 mmol) and hydrazide (0.7 mmol). A sample of the powder

Synthesis of Ligands
Ligand preparation by solution-based synthesis: Ligands H 2 L H , H 2 L 3OMe , H 2 L 4OMe and were prepared by condensation of aldehyde and hydrazide, according to the reported procedure [24][25][26]. Ligand preparation by mechanochemical synthesis: Liquid assisted grinding (LAG) experiments were performed using a Retsch MM200 ball mill. 0.1044 g (0.7 mmol) of benzaldehyde, 3-or 4-methoxysalicylaldehyde, 0.1041 g (0.7 mmol) 2-hydroxy-benzydrazide, and 50 µL of methanol were placed with one 10 mm milling ball in a 10 mL Teflon jar. For the preparation of H 2 L H no solvent was added to the mixture of salicylaldehyde, 0.085 g (0.7 mmol) and hydrazide (0.7 mmol). A sample of the powder obtained by milling at 25 Hz, after 45 min, was immediately analyzed by PXRD and compared to the compound obtained by solution-based synthesis. No purification of the sample was performed.

Synthesis of Molybdenum(VI) Complexes
A mixture of [MoO 2 (acac) 2 ] (0.64 g, 2 mmol) and H 2 L H, 3OMe or 4OMe (2 mmol) in methanol (120 mL) or DCM (140 mL) was refluxed for four hours. The solution was left at room temperature and the resulting crystalline substance was filtered and dried in a desiccator up to a constant weight.

General Procedure for the Epoxidation Of Cyclooctene by TBHP
A mixture of cyclooctene (2.76 mL, 20 mmol), acetophenone as an internal standard, and molybdenum catalyst (0.05 mmol) was stirred and heated to 80 • C before the addition of TBHP in water (70% w/w, 40 mmol) or in decane (5.5 M, 40 mmol). The reaction was monitored for 6 h with withdrawal and analysis of organic phase aliquots (0.1 mL) at required times. Each withdrawn sample was mixed with 1 mL of Et 2 O and analyzed by GC.

Theoretical Calculations
The geometries of all species under investigation were optimized without any symmetry constraint with the Gaussian 09 rev. D01 program suite [27], with the DFT approach using the B3LYP three-parameter functional [28][29][30] in conjunction with the 6-31G* basis set [31][32][33][34] for the light atoms (O, N, C, H) and the CEP-31G set for the Mo atom [35,36]. The geometries of all complexes and intermediates were optimized from starting geometries determined or inspired by X-ray diffraction without any symmetry constraint. All coordinates have been listed in Supplementary Materials. Frequency analysis confirmed that the optimized geometries of all the stable compounds and intermediates were local minima. The transition states were optimized using a preliminary scan of a relevant internal coordinate, followed by full optimization of the TS guided by the knowledge of such coordinates. All optimized geometries were confirmed to be stationary points and local minima (for stable molecules or reaction intermediates) or first order saddle points (for the TSs) by frequency analyses. For all the TSs, analysis of the imaginary frequency confirmed the expected motion along the reaction coordinate. Those values and the relative schemes have been added in Supplementary Materials (Table S4-S6)) The calculated frequencies were also used to derive the thermochemical parameters at 298 K according to the standard approximations (ideal gas, rigid rotor, and harmonic oscillator)

Physical Methods
Elemental analyses were provided by the Analytical Services Laboratory of the Ruder Bošković Institute, Zagreb. Thermogravimetric (TG) analysis was performed using a Mettler TG 50 thermobalance in aluminum crucibles. All experiments were performed in oxygen atmosphere with a flow rate of 200 cm 3 min −1 , with heating rates of 10 K min −1 . Differential scanning calorimetry (DSC) measurements were carried out using a Mettler-Toledo DSC823e calorimeter and analyzed by Mettler STARe 9.01. software. IR-ATR spectra were recorded on a Perkin Elmer Spectrum One spectrometer, in the spectral range between 4500-450 cm −1 . The catalytic reactions were followed by gas chromatography on an Agilent 6890A chromatograph equipped with an FID detector and a DB5-MS capillary column (30 mm × 0.32 mm × 0.25 mm). The GC parameters were quantified using authentic samples of the reactants and products. The conversion of cis-cyclooctene and the formation of cyclooctene oxide were calculated from calibration curves (r 2 = 0.999) relative to acetophenone as an internal standard.
High  Table S3. The structures were solved with dual space methods using SHELXT [38]. The refinement procedure by full-matrix least squares methods based on F 2 values against all reflections included anisotropic displacement parameters for all non-H atoms. Hydrogen atoms bound to carbon atoms were placed in geometrically idealized positions and refined by use of the riding model with U iso = 1.2U eq of the connected carbon atom or as ideal CH 3 groups with U iso = 1.5U eq . Hydrogen atoms attached to heteroatoms were located in the difference Fourier maps at the final stages of the refinement procedure. Their coordinates were refined freely but with restrained O-H distances of 0.82(2) Å. All refinements were performed using SHELXL-2013 [39]. The SHELX programs operated within the Olex2 [40] suite. Geometrical calculations and molecular graphics were done with Mercury [41].

Conclusions
Oligomeric [MoO 2 L] n and monomeric [MoO 2 L(MeOH)] molybdenum complexes were prepared by the conventional synthetic procedures, employing DCM or MeOH, as solvents. Formation of oligomeric complexes from monomeric ones could be caused directly through heating. This provided a more convenient synthetic route, eliminating the use of the environmentally harmful DCM solvent. All prepared molybdenum complexes were tested as epoxidation catalysts, with the assistance of TBHP in water or TBHP in decane. For the first time, monomeric and oligomeric complexes showed both excellent activity and selectivity towards epoxide. Similar results were obtained, no matter the oxidant used, giving priority on the use of TBHP in water, in terms of the green process. DFT calculations have extracted trends with this family of complexes rationalizing the experimental values in regard to the ligand structure. Further research will be oriented to the mechanochemical synthesis and solvent-less methods for the complex preparation to respect green chemistry principles. In terms of the catalytic process, improvements will be done towards the use of different model substrates.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10 .3390/catal11070756/s1, Figure S1: Comparison of IR-ATR spectra of the complex [MoO2(L 3OMe ) (MeOH)] (pink) and [MoO2(L 3OMe )]n (blue). Figure S2: Enumeration of atoms in the crystal structure of [MoO 2 (L 3OMe )(MeOH)], Figure S3: Enumeration of atoms in the crystal structure of [MoO 2 (L 4OMe ) (MeOH)]·4,4-bpy, Figure S4: Crystal packing in the crystal structure of [MoO 2 (L 3OMe )(MeOH)] shown parallel to crystal a-axis, Figure S5: Crystal packing in the crystal structure of [MoO 2 (L 4OMe ) (MeOH)]·4,4-bpy shown parallel to crystal b-axis, Figure S6 Table S1: Geometrical parameters of the octahedral environment of molybdenum atom in the title complexes, Table S2: Bond lengths (in Å) and angles (in • ) for the prepared complexes. Atom enumeration is shown in Figures S1 and S2, Table S3: Summary of general and crystallographic data for the prepared complexes, Table S4: Simplified images of all the calculated Mo-containing species, Table S5: Imaginary frequency (cm −1 ) of the calculated transition states shown in Table S1, Table S6: DFT Coordinates of all optimized species. CCDC 2085847-2085848 contains the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via Available online: www.ccdc.cam.ac.uk/structures.