# Impact of POM's coordination mode and Mo-hybrid constituents on the binding, stability, and catalytic properties of hybrid (pre)catalysts

Jana Pisk, Višnja Vrdoljak, Mirna Mandarić, Tomica Hrenar, Dominique Agustin, Mirta Rubčić

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Jana Pisk, Višnja Vrdoljak, Mirna Mandarić, Tomica Hrenar, Dominique Agustin, et al., Impact of POM's coordination mode and Mo-hybrid constituents on the binding, stability, and catalytic properties of hybrid (pre)catalysts. RSC Advances, 2024, 14 (27), pp.19029 - 19040. 10.1039/d4ra03563a. hal-04632363

HAL Id: hal-04632363

https://hal.science/hal-04632363v1

Submitted on 2 Jul 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Impact of POM's coordination mode and Mo-hybrid constituents on the binding, stability, and catalytic properties of hybrid (pre)catalysts†

Jana Pisk, ID a Višnja Vrdoljak, ID *a Mirna Mandarić, a Tomica Hrenar, ID a Dominique Agustin ID bc and Mirta Rubčić ID a

The assembly of MoVI O₂²⁺ and methoxy-substituted salicylaldehyde nicotinoyl hydrazone ligands afforded two classes of hybrid polyoxometalates (POMs). In the Class I architectures, [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉]·xCH₃COCH₃ (D = CH₃COCH₃ or H₂O, x = 0 or 2, and L¹⁻³ = ligands bearing the OMe group at position 3, 4 and 5, respectively), the main driving force for self-assembly is the electrostatic interaction between the components. Class II architectures are composed of a POM anion covalently linked to two Mo-complex units through the terminal O⁺ or bridging µ₂-OₚOM oxygen atoms, as found in Lindqvist-based hybrids [{MoO₂(HL¹⁻³)}₂Mo₆O₁₉]·xCH₃CN (x = 0 or 2) and the asymmetrical β-octamolybdate-based hybrid [{Mo₂O₄(HL²)(H₂L)}{MoO₂(HL²)}₂Mo₈O₂₆]·CH₃CN·H₂O. Quantum chemical calculations were applied to evaluate the impact of the POM hybrid constituents on the hybrid-type stability, showing that it strongly depends on the ligand substituent position and ancillary ligand nature. Hybrids were tested as catalysts for cyclooctene epoxidation using tert-butyl hydroperoxide (TBHP in water or decane) and with or without the addition of acetonitrile (CH₃CN) as an organic solvent. The catalytic results provided by the use of TBHP in decane are the best ones and classify all the prepared catalysts as very active, with the conversion of cyclooctene >90%, and high selectivity towards epoxide, >80%. We also examined the influence of the ligand structure, POM's hybrid type, and coordination mode on the Mo-hybrid activity and selectivity.

Received 15th May 2024

Accepted 6th June 2024

DOI: 10.1039/d4ra03563a

rsc.li/rsc-advances

## Introduction

Polyoxometalates (POMs) have attracted attention in many fields of science due to their structural diversity, compositional adjustability, and numerous chemical features.¹⁻⁵ They have particularly received attention in catalysis,⁶,⁷ contributing to the atomically efficient and energetically sustainable objectives. Namely, POMs are highly effective catalysts in various relevant processes, such as in oxidation reactions of olefins,⁸⁻¹⁰ alcohols,¹¹,¹² aldehydes,¹³ or amines.¹⁴ The development of new POM

functional materials with unique properties is therefore essential for the continuous progress in that research field. In this context, their derivatization by organic components¹⁵⁻¹⁹ or transition metal complexes²⁰⁻²⁵ has proven to be a convenient method.

Among various catalysts, MoVI complexes have been particularly useful for the epoxidation of cyclic olefins.²⁶⁻³⁰ The metallosupramolecular species with long Mo-Nligand bonds were particularly effective.³¹⁻³⁴ It was, therefore, reasonable to combine MoVI coordination compounds with POMs into organic-inorganic hybrid architectures. In principle, POMs should be a suitable substitute for any ligand, but due to steric hindrance effects, one might expect their coordination to Mo with a long Mo-O bond. However, their controlled assembly in organic media is less straightforward. To date, only two examples of crystalline species [{MoO₂(HL)}₂{Mo₆O₁₉}] have been structurally characterized and in both of them, POM links two Mo-complex fragments through POM terminal oxygen atoms.³⁵,³⁶

The coordination of Lindqvist anion to Mo-complex through the bridging µ₂-OₚOM is not yet known, although the bridging oxygen atoms are more nucleophilic and basic than the O⁺ terminal sites.³⁷ This drawback may be due to steric hindrance imposed by POM and organic ligands. Only one example of the

<figure><img src="image_4.png" /></figure>

<footnote>ªUniversity of Zagreb, Faculty of Science, Department of Chemistry, Horvatovac 102a, 10000 Zagreb, Croatia. E-mail: visnja.vrdoljak@chem.pmf.hr</footnote>

<footnote>ᵇLCC-CNRS (Laboratoire de Chimie de Coordination), 205 Route de Narbonne, BP44099, CEDEX 4, 31077 Toulouse, France</footnote>

<footnote>Department of Chemistry, IUT Paul Sabatier, Université Paul Sabatier, University of Toulouse, Av. G. Pompidou, CS20258, 81104 Castres, France</footnote>

<footnote>† Electronic supplementary information (ESI) available: (1) Powder diffraction patterns, (2) ATR-IR spectra, (3) TGA curves, (3) additional figures for compounds, (4) tables of selected bond distances and angles and of hydrogen bonds parameters. Crystallographic data sets for the structures 1H⁺·2CH₃COCH₃, 1H⁺, 1H, 2H⁺·2CH₃COCH₃, 2H, 3H⁺·2CH₃COCH₃, 2O·CH₃CN·H₂O, 3H, 3H·2CH₃CN. CCDC 2335829-2335837. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4ra03563a</footnote>
<figure><img src="image_5.png" /></figure>

Lindqvist µ₂·OₚOM coordination to Cu-complex has been reported so far.³⁸ Additionally, a few examples of such α-Keggin-type hybrid structures with transition metal complex moieties are known.³⁹⁻⁴² Furthermore, according to the Cambridge Structural Database,⁴³ the formation of asymmetrically anchored Mo hybrid systems has never been reported.

So far, more results have been achieved in the synthesis of various MoVI POM hybrid salts. POM anions tend to associate with cationic Mo-complexes *via* multiple non-covalent interactions and give ionic hybrids (Class I).⁴⁴⁻⁴⁶ However, investigation of their potential was limited by their low stability. Coordination of POM to metal centres offers a way to overcome this issue and allows the discovery of new POM coordination modes.

Herein, we report the successful self-assembly of individual building units into ionic hybrids [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉]·xsolv (Class I) consisting of coordination cations and POM anions, and their conversion into the corresponding covalently linked architectures [{MoO₂(HL¹⁻³)}₂Mo₆O₁₉]·xsolv (Class II). As organic ligands, we selected methoxy-substituted salicylaldehyde nicotinoyl hydrazones H₂L¹⁻³, (Scheme 1), considering the feasibilities for the fine-tuning of structural features and their tendency towards reversible protonation and deprotonation.

All hybrids were characterized using the single-crystal and powder X-ray diffraction (XRD) methods, infrared (IR-ATR) spectroscopy, thermogravimetric analysis (TGA), and elemental analysis. In the case of Class II assemblies, both coordination modes through the terminal Mo = Ot and bridging µ₂-O oxygen atoms to two peripheral Mo-coordination subunits were successfully obtained. The β-octamolybdate hybrid [{Mo₂O₄(HL²)(H₂L²)}{MoO₂(HL²)}]Mo₈O₂₆]·2CH₃CN·H₂O corresponds to the first example of an asymmetrical arrangement of Mo-fragments around the POM anion.

Quantum chemical calculations were employed to support experimental results and estimate the influence of the ancillary ligand, hydrazone substituent, and POM coordination type on the hybrid-class formation.

We envisaged those POM-hybrid architectures, combining favourable properties of coordination and inorganic moieties, as promising catalysts for the epoxidation of olefins. The covalently anchored ones would be particularly appealing because the activity of Mo-complexes could be directly related to the POM bonding nature. In our research, the catalytic activity of hybrids in the selective oxidation of cyclooctene was evaluated using tert-butyl hydroperoxide (TBHP) as the oxidant (in water or decane) and with or without the addition of the organic solvent (CH₃CN).

<figure><img src="image_6.png" /><figcaption>**Scheme 1** Nicotinoyl hydrazonato ligands in neutral zwitterionic and singly-deprotonated (HL⁻) forms bearing OMe group at positions 3 (H₂L¹), 4 (H₂L²), and 5 (H₂L³).</figcaption></figure>

# Results and discussion

## Synthesis and characterization of hybrids

Ligands H₂L¹⁻³ (Scheme 1) were prepared by the reaction of isonicotinyl hydrazine with an appropriate aldehyde according to the published procedures.²⁶,⁴⁷ The synthetic pathways for POM-based hybrids and a list of all hybrids are provided in Scheme 2. Almost all assemblies accommodate lattice solvent molecules. Their formulae are consistent with analytical, spectroscopic and X-ray diffraction data. A corresponding designation of compounds can be found in Scheme 3.

Reactions of [MoO₂(acac)₂] (acac = acetylacetonate) and hydrazone ligands H₂L¹⁻³ in a molar ratio of 3.3 : 1 in acetone resulted initially in orange suspensions, indicating the MoVI hydrazone complex formation. A significant quantity of the orange solids was separated. A comparison of the PXRD patterns identified them to be a metallosupramolecular species [MoO₂(L¹⁻³)]ₙ, known from our previous investigation.²⁶,³¹ Upon standing the filtrate at room temperature, the condensation of the remaining metal oxo-units led to the formation of the Class I assemblies [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉]·2CH₃COCH₃, where D = acetone in 1Hᵃ·2CH₃COCH₃, or water in 2HW·2CH₃COCH₃ and 3HW·2CH₃COCH₃. The red crystals deposited after about one to two weeks.

The reactions of [MoO₂(acac)₂] with H₂L¹ and H₂L³ in a molar ratio of 10 : 1 in acetonitrile proceeded similarly. Despite the presence of a higher amount of [MoO₂(acac)₂], which is more favourable for the formation of octamolybdates,⁴⁵ the reactions with H₂L¹ and H₂L³ always resulted in the formation of hexamolybdate hybrids [MoO₂(HL¹)(H₂O)]₂[Mo₆O₁₉] (1HW) and [(MoO₂(HL³))₂Mo₆O₁₉]·2CH₃CN (3H·2CH₃CN). This could be due to the better stability of the hexamolybdate anion in acetonitrile in the presence of Mo-complexes with these two ligands. The resulting hybrids 1HW and 3H·2CH₃CN, belong to Class I and Class II, respectively. However, in the case of H₂L², the condensation always resulted in the asymmetrical Class II hybrid [{Mo₂O₄(HL²)(H₂L²)}{MoO₂(HL²)}]Mo₈O₂₆·CH₃CN·H₂O (2O·CH₃CN·H₂O). To our knowledge, this is the first example of a POM-based hybrid with octamolybdate covalently expanded by three Mo-hydrazone complex units.

Traces of moisture in the solvents are crucial for the polyoxometalate anion formation. The addition of acetic acid during the reactions increased the purity of the compounds, but there was no difference in final hybrid product formation. The powder X-ray diffraction (PXRD, Fig. S1, see ESI†) patterns of the bulk products are well-matched with the calculated ones. However, the yields depend on the quantity of the formed coordinated polymers. It should be mentioned that the reaction in acetone in a molar ratio of 10 : 1 yielded each time a gelatinous reaction mixture, whereas the synthesis in acetonitrile in a molar ratio of 3.3 : 1 gave a mixture of products.

Solvothermal treatment of the corresponding hybrid salts 1Hᵃ·2CH₃COCH₃, 1HW, 2HW·2CH₃COCH₃, and 3HW·2CH₃COCH₃ in dry acetonitrile resulted in the exclusion of coordinated solvent molecules (acetone or water) and POM anion immobilisation. This approach resulted in [{MoO₂(HL¹⁻³)}₂Mo₆O₁₉]·xCH₃CN, x = 0 or 2
<figure><img src="image_7.png" /><figcaption>**Scheme 2** Summary of POM-based hybrids (Class I and Class II). Black arrows depict the syntheses carried out at room temperature in acetone or acetonitrile. Orange and blue arrows stand for the solvothermal transformations carried out in acetone and acetonitrile, respectively.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Scheme 3** Compounds designation – N = 1, 2, and 3 correspond to hybrids with L¹, L², and L³ within the [MoO₂(HL)(D)]⁺ cation. Hexamolybdates and octamolybdate are abbreviated as H and O, respectively, whereas a and w represent ancillary ligand D, acetone and water.</figcaption></figure>

(1H, 2H, and 3H·2CH₃CN), Class II hybrids. The condensation under solvothermal conditions in dry acetone was successful in the transformation of 1Hᵃ·2CH₃COCH₃ and 3Hʷ·2CH₃COCH₃ into 1H and 3H, respectively. In the case of 1Hʷ and 2Hʷ·2CH₃COCH₃, such conversion failed to proceed.

Attenuated Total Reflectance Infrared (ATR-IR) spectra for all compounds are consistent with the formation of hybrids (Fig. S2–S4, see ESI†). Whereas, the vibrations around 950 cm⁻¹ could be assigned to the v(Mo=O) in the POM anion and the bands in the region of 942–886 cm⁻¹ corresponds to the terminal v(Mo=O) of the Mo-hydrazonato complex. The bands in the range of 789 cm⁻¹ to 712 cm⁻¹ may be assigned to the v(Mo–O–Mo) vibrations of the hexamolybdate Mo-oxo groups, while those in the range of 689 cm⁻¹ to 642 cm⁻¹ are related to v(Mo–O–Mo) vibrations of the octamolybdate cluster.⁴⁸,⁴⁹

Crystals of solvates 2Hʷ·2CH₃COCH₃, 3Hʷ·2CH₃COCH₃, and 3H·2CH₃CN exhibit excellent stability and retain their structural integrity at room temperature for long periods without decomposing. According to the thermal analysis (Fig. S2 and S3, see ESI†), they lose the coordinated and/or non-coordinated solvent molecules in one mass loss step in the range 122–161 °C (2Hʷ·2CH₃COCH₃), 158–187 °C (3Hʷ·2CH₃COCH₃), and 167–222 °C (3H·2CH₃CN). In the second step, mass loss started around 290 °C.

However, in the case of 1Hᵃ·2CH₃COCH₃ and 2O·CH₃CN·H₂O, a complete loss of lattice solvent occurred at room temperature after a few minutes. The analysis of residuals indicates that they are of the composition 1Hᵃ and 2O. Thermogravimetric curves of 1Hᵃ and 1Hʷ exhibit two separate mass loss steps (Fig. S2 and S3†). The first mass loss between 130 °C and 157 °C and between 182 °C and 213 °C corresponds to the loss of coordinated acetone and water molecules, respectively.

Solvent-free hybrids 1, 2, 2O, and 3 are very stable. As shown in Fig. S3 and S4†, the mass loss starts at 275 °C, 243 °C, 295 °C, and 317 °C, respectively, with the decomposition of the metal-oxo cluster. In all the cases, the final step mass losses agree with the calculated values, considering the product of the decomposition to be the MoO₃.

## Description of molecular and crystal structures

Within the structures of the ionic type one can distinguish those where the complex cation is of the [MoO₂(HL)(D)]⁺ kind (Fig. 1), where HL⁻ denotes a singly deprotonated form of the ligand and D represents ancillary ligand, acetone or water molecule, fulfilling the deformed octahedral coordination environment of the Mo atom (1Hᵃ·2CH₃COCH₃, 3Hʷ·2CH₃COCH₃, 2Hʷ·2CH₃COCH₃ and 1Hʷ). In all structures of ionic type of hybrids, polyoxometalate anions are classical Lindqvist ones, Mo₆O₁₉²⁻.

Within the class of covalently linked hybrids, one finds those in which the Lindqvist anion Mo₆O₁₉²⁻ acts a bridging unit between two [MoO₂(HL)]⁺ fragments (1H, 2H, 3H, and 3H·2CH₃CN), and the one where β-octamolybdate bridges [Mo₂O₄(HL)(H₂L)]³⁺ and [MoO₂(HL)]⁺ fragments, 2O·2CH₃CN·H₂O.

It should be highlighted that in all structurally characterized compounds, except the 2O·2CH₃CN·H₂O one, the hydrazone ligand coordinates in the monoanionic form *via* O1, N1 and O2 atoms, which is evident from the relevant bond distances (Scheme 1, Tables S1 and S2 in the ESI†), and also consistent with previously reported examples.⁴⁵ It is also interesting to note that coordinated ligands assume in most cases an essentially

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>Fig. 1 Molecular structures of 1H<sup>a</sup>·2CH<sub>3</sub>COCH<sub>3</sub> with the atom-labeling scheme as illustrative example of ionic structures described in this work. Ellipsoids are drawn at 30% probability level. Hydrogen atoms are presented as spheres of arbitrary small radii. The O<sub>11</sub> oxygen atom is situated on an inversion center, and for clarity, only one of the two equally occupied positions for the non-coordinated acetone molecule is shown. Dihedral angles between the planes of the pyridyl and the phenyl ring of the aldehyde residue amount to 5°.</figcaption></figure>

planar conformation (Fig. 1, 2 and S5, see ESI†), while the largest deviations from planarity are observed for the ligands in 3H·2CH<sub>3</sub>CN and 2O·2CH<sub>3</sub>CN·H<sub>2</sub>O.

Molecular structures of ionic hybrids 1H<sup>a</sup>·2CH<sub>3</sub>COCH<sub>3</sub>, 3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub>, 2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> and 1H<sup>w</sup> contain Lindqvist anions, Mo<sub>6</sub>O<sub>19</sub><sup>2-</sup>, and complex cations of the type [MoO<sub>2</sub>(HL)(D)]<sup>+</sup> (D = CH<sub>3</sub>COCH<sub>3</sub> or H<sub>2</sub>O), in which the Mo atoms assume distorted octahedral environments (Fig. 1 and S5b-d†). As expected for octahedral coordination geometry, the longest distances within the coordination sphere of the Mo atom of the [MoO<sub>2</sub>(HL)(D)]<sup>+</sup> units are those positioned *trans* to Mo=O functionalities, namely Mo1–O6 and Mo1–N1 bonds. The presence of additional donor molecules at the sixth coordination site of molybdenum in [MoO<sub>2</sub>(HL)(D)]<sup>+</sup> cations along with the presence of solvent molecules in 1H<sup>a</sup>·2CH<sub>3</sub>COCH<sub>3</sub>, 3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub>, 2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> lead to diverse packing scenarios (Fig. S6–S9, see ESI†). The crystal structure of 1H<sup>a</sup>·2CH<sub>3</sub>COCH<sub>3</sub> unveils strong hydrogen bonding between the protonated pyridine moiety and non-coordinated acetone molecule in a discrete D<sub>1</sub><sup>1</sup>(2) graph-set motif. In the absence of further strong hydrogen bond donors, hexamolybdate anions remain involved only in C–H···O interactions with the [MoO<sub>2</sub>(HL)(D)]<sup>+</sup> units and the solvent molecules within the crystal structure (Fig. S6†). The same, evidently robust D<sub>1</sub><sup>1</sup>(2), hydrogen bond motif, established between the protonated pyridine moiety and the acetone molecule is also found in the structures of 2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> and 3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> (Fig. S7 and S8†).

In contrast to 1H<sup>a</sup>·2CH<sub>3</sub>COCH<sub>3</sub>, the presence of additional hydrogen bond donors, i.e. coordinated water molecules in [MoO<sub>2</sub>(HL)(H<sub>2</sub>O)]<sup>+</sup>, affords a complex hydrogen bonded networks in both cases. Namely, in the 2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> the coordinated water molecule forms hydrogen bonds with the Lindqvist anion, and with the acetone molecule yielding thus complex supramolecular chains along the crystallographic *a*-axis (Fig. S8, see ESI†). On the other hand, in 3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub> one observes supramolecular dimers of [MoO<sub>2</sub>(HL)(H<sub>2</sub>O)]<sup>+</sup> units,

<figure><img src="image_11.png" /><figcaption>Fig. 2 Molecular structures of: (a) 3H·2CH<sub>3</sub>CN and (b) 2O·2CH<sub>3</sub>CN·H<sub>2</sub>O, and (c) 3H with the atom-labeling scheme. Ellipsoids are drawn at 30% probability level, whereas hydrogen atoms are shown as spheres of arbitrary small radii. In (a) O11 and in (c) O9 atoms are situated on an inversion center. For the clarity of presentation, in (b) solvent molecules are omitted. Dihedral angles between the planes of the pyridyl and the phenyl ring of the aldehyde residue for 3H·2CH<sub>3</sub>CN and 3H, respectively amount to ca. 33.3° and ca. 2.2° (for 2O·2CH<sub>3</sub>CN·H<sub>2</sub>O are discussed below).</figcaption></figure>

formed via O6–H6A···N2 hydrogen bonds, which are bridged by Mo<sub>6</sub>O<sub>19</sub><sup>2-</sup> anions. Crystal structures in both cases are stabilized by plethora of C–H···O interactions (Fig. S7, S8 and Table S3, see ESI†). Supramolecular architecture of 1H<sup>w</sup> differs from the previously described ones, as there are no solvent molecules present in the structure. Here the [MoO<sub>2</sub>(HL)(H<sub>2</sub>O)]<sup>+</sup> units associate *via* O–H···O and N–H···O interactions into complex assemblies which are joined by hydrogen-bond anchored Lindqvist anions (Fig. S9†). It is interesting to note that in the absence of the solvent molecules in 1H<sup>w</sup>, i.e. acetone, the protonated pyridine moiety establishes
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>Fig. 3 Supramolecular assemblies observed in the crystal structures of: (a) 1H, and (b) 3H·2CH₃CN. In (a) N3–H3···O4 hydrogen bonds yield a R₂²(16) motif, and N3–H3···O15 a C₁⁺(12) motif. In (b) N3–H3···O4 and N3–H3···N4 hydrogen bonds yield R₂²(16) and D₁⁺(2) motifs.</figcaption></figure>

hydrogen bond with one of the oxygen atoms of the {MoO₂}²⁺ core of the [MoO₂(HL)(H₂O)]⁺ unit.

In the structures of neutral complexes polyoxomolybdate anions, Mo₆O₁₉²⁻ in 1H and 3H·2CH₃CN and Mo₈O₂₄⁴⁻ in 2O·2CH₃CN·H₂O, act as bridging entities (Fig. 2). Polyoxoanions in these structures coordinate via their terminal oxygen atoms to molybdenum atoms of the {MoO₂(HL)} fragments, thus providing the Mo atoms of the later octahedral environment.

Structure of 2O·2CH₃CN·H₂O is unique as the Mo₈O₂₆⁴⁻ unit bridges one dinuclear [Mo₂O₄(HL)(H₂L)]³⁺ and one mononuclear [MoO₂(HL)]⁺ entity (Fig. 2b). In the dinuclear unit one finds the coordinated ligand both in singly-deprotonated HL⁻ as well as in the neutral zwitterionic form H₂L (Scheme 1), whereas in mononuclear unit [MoO₂(HL)]⁺ ligand is coordinated in its singly-deprotonated form (Scheme 1). In this structure one can observe substantial differences in the conformation of zwitterionic form of H₂L with respect to singly-deprotonated ones, HL⁻. Namely, the angle between the pyridyl and aldehyde rings of the zwitterionic form of H₂L amounts to ca. 2.5°, in the case of HL⁻ in the [Mo₂O₄(HL)(H₂L)]³⁺ unit it is ca. 35.4° and for the [MoO₂(HL)]⁺ is ≈21.9°. It is worthwhile to note that singly-deprotonated forms of the ligand are found both with syn- and anti- oriented pyridyl moieties with respect to the coordinating ONO compartment of HL⁻. The crystal structure reveals a rich assortment of hydrogen bonds. Namely, hydrogen atoms of the protonated pyridyl moieties of all three coordinated ligands afford bifurcated hydrogen bonds with the octamolybdate anions. Additional hydrogen bonding involves crystal water molecules and the octamolybdate anions, as well as the hydrazinic N–H group (zwitterionic form H₂L) and the oxygen atoms of the {MoO₂}²⁺ core, forming thus a complex supramolecular network (Fig. S10†). It should be noted that acetonitrile molecules are not involved in strong intermolecular interactions, which explains the observed positional disorder, related to acetonitrile molecules, in the structure. So far, only Class I octamolybdate hybrids with coordinated solvent molecules on the {MoO₂}²⁺ core have been studied, offering valuable insight into their structural properties.⁴⁵,⁴⁶

When comparing the structures of 1H and 3H·2CH₃CN, one observes the difference in the engagement of the protonated pyridine moiety in the hydrogen bonding due to the presence of solvent molecules in 3H·2CH₃CN (Fig. 3). Namely, in 1H

hydrogen atom of the protonated pyridine moiety participates in bifurcated hydrogen bond, N3–H3···O4 and N3–H3···O15, forming thus supramolecular layers in the ab-plane (Fig. S11, see ESI†).

Molecular and crystal structures of the hybrids 2H and 3H are similar. Namely, in both cases molecular structures consist of the [MoO₂(HL)]⁺ type of cations which are bridged by Lindqvist anions, Mo₆O₁₉²⁻, to afford overall neutral species (Fig. 2c and S5a†). It is interesting to observe that in both 2H and 3H bridging oxygen atoms of the {Mo₆O₁₉} units coordinate to Mo1 atoms of the complex unit {MoO₂(HL)}, with the corresponding distances being 2.723(5) for 2H and 2.705(6) for 3H (in both cases Mo1–O6 distances), Table S4, see ESI†.

Similarity in molecular structures of 2H and 3H is also reflected in their packing arrangements. In both cases, such complexes associate via N–H···O hydrogen bonds involving both {[MoO₂(HL)]} and {Mo₆O₁₉} fragments (Fig. 4), which are further associated through C–H···O interactions (Fig. S12, S13 and Table S4, see ESI†).

Finally, it is interesting to observe that in the Lindqvist containing hybrids, 2H and 3H have their pyridyl moieties-oriented anti with respect to the coordinating ONO compartment of the HL⁻, whereas in the structures of (1Hᵃ·2CH₃COCH₃, 3Hʷ·2CH₃COCH₃, 2Hʷ·2CH₃COCH₃, 1Hʷ, 1H and 3H·2CH₃CN) one observes both syn and anti-arrangements.

## Quantum chemical calculations

In order to explain possible paths of crystallization, structures of [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉] and [{MoO₂(HL¹⁻³)}₂Mo₆O₁₉] were built starting from the crystallographically determined structures, and their geometries were optimized. For each optimized geometry, calculations of harmonic frequencies were performed and thermodynamic properties were computed. Relative differences in standard Gibbs energies of binding Δ<sub>b</sub>G'ₜₚel (Tables 1 and 2) were calculated using the calculated standard Gibbs energies of formation.

Analysis of calculated structures for [MoO₂(HL¹)(D)]₂[Mo₆O₁₉] revealed that the most stable complex should be the one where the metal ion is coordinated with MeOH. This complex has the lowest relative standard Gibbs energy of binding Δ<sub>b</sub>G'ₜₚel. However, for
<figure><img src="image_14.png" /><figcaption>Fig. 4 Supramolecular assembly observed in the crystal structure of 2H. {[MoO₂(HL)]} units of the 2H structure is engaged in N3–H3···O4 hydrogen bonds forming C₁⁺(8) chains (along the b-axis), while the {Mo₆O₁₉} fragments participate in N3–H3···O11 hydrogen bonds shaping C₁⁺(10) motif.</figcaption></figure>

other complexes [MoO₂(HL²⁻³)(D)]₂[Mo₆O₁₉] the complexes with the lowest ΔbG°rel were the ones with coordinated H₂O molecule.

Depending on the orientation of the pyridyl moieties of the HL⁻ ligand, being oriented *anti* or *syn* (Fig. 5), the hybrids have different Gibbs binding energies. By considering the Gibbs binding energy for the structures having *syn*-orientation, Fig. 5a, in the case of Mo complex containing the ligand HL¹, coordination of POM *via* the terminal oxygen atom is preferred over the bridging one, while in the case of ligands HL² and HL³, the coordination *via* the bridging oxygen atom is more stable.

Those differences in Gibbs energies amount to 31.1 kJ mol⁻¹. Even so, hybrid 3H·2CH₃CN (coordinated through Ot) crystalized depending on the reaction conditions. This could probably be due to the stabilization through the hydrogen bonding implied by the presence of solvent molecules in 3H·2CH₃CN and sterical hindrance directly related to the size of the POM. Nevertheless, in the case of the nicotinic hydrazone moiety *anti*-orientation, shown in Fig. 5b, the POM coordination *via* the bridging oxygen atom is more favourable. The results are in line with the experimental results for 1H, 2H, and 3H. The strong

<table><caption>**Table 1** Calculated relative differences in standard Gibbs energies of binding at 298.15 K and 1 atm for [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉] relative to the complex of minimal energy (B3LYP-D3/Def2TZVP level of the theory). *Syn*- and *anti*-orientation of pyridyl moieties relative to coordinating ONO compartment of HL⁻ in the complex cation</caption><thead><tr><th rowspan="2"></th><th colspan="3">ΔbG°rel/kJ mol⁻¹</th><th colspan="3"></th></tr><tr><th colspan="3">[MoO₂(HL¹⁻³)(D)]⁺ syn-orientation</th><th colspan="3">[MoO₂(HL¹⁻³)(D)]⁺ anti-orientation</th></tr></thead><tbody><tr><td>D</td><td>HL¹</td><td>HL²</td><td>HL³</td><td>HL¹</td><td>HL²</td><td>HL³</td></tr><tr><td>CH₃CN</td><td>9.22</td><td>3.31</td><td>4.24</td><td>3.46</td><td>4.17</td><td>3.75</td></tr><tr><td>H₂O</td><td>4.86</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td><td>0.00</td></tr><tr><td>CH₃OH</td><td>0.00</td><td>5.05</td><td>5.15</td><td>8.76</td><td>2.77</td><td>4.90</td></tr><tr><td>CH₃COCH₃</td><td>14.28</td><td>6.80</td><td>7.58</td><td>12.26</td><td>7.46</td><td>6.99</td></tr><tr><td>Experimental</td><td><b>1H<sup>W</sup></b></td><td><b>2H<sup>W</sup>·2CH₃COCH₃</b></td><td></td><td><b>1H<sup>a</sup>·2CH₃COCH₃</b></td><td></td><td><b>3H<sup>W</sup>·2CH₃COCH₃</b></td></tr></tbody></table>

<table><caption>**Table 2** Calculated relative differences in standard Gibbs energies of binding at 298.15 K and 1 atom for [MoO₂(HL¹⁻³)Mo₆O₁₉] relative to the complex of minimal energy (B3LYP-D3/Def2TZVP level of the theory). *Syn* and *anti*-orientation of pyridyl moieties relative to coordinating ONO compartment of HL⁻ in the complex cation</caption><thead><tr><th rowspan="2"></th><th colspan="3">ΔbG°rel/kJ mol⁻¹</th><th colspan="3"></th></tr><tr><th colspan="3">[MoO₂(HL¹⁻³)]₂(Mo₆O₁₉)] syn-orientation</th><th colspan="3">[MoO₂(HL¹⁻³)(Mo₆O₁₉)] anti-orientation</th></tr></thead><tbody><tr><td>POM donor atom</td><td>HL¹</td><td>HL²</td><td>HL³</td><td>HL¹</td><td>HL²</td><td>HL³</td></tr><tr><td>μ₂-O<sub>POM</sub></td><td>555.21</td><td>507.65</td><td>496.88</td><td>521.87</td><td>508.11</td><td>514.77</td></tr><tr><td>O<sup>t</sup>-POM (terminal)</td><td>535.94</td><td>530.90</td><td>527.98</td><td>536.42</td><td>527.50</td><td>530.88</td></tr><tr><td>Experimental</td><td><b>1H</b></td><td></td><td><b>3H·2CH₃CN</b></td><td></td><td><b>2H</b></td><td><b>3H</b></td></tr></tbody></table>
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>Fig. 5 Representation of (a) syn- and (b) anti-orientation of pyridyl moieties relative to coordinating ONO compartment of HL⁻ in the complex cation.</figcaption></figure>

hydrogen bonds identified in those structures involve terminal oxygen atoms.

## Catalytic studies

Those hybrid species are interesting since they bear two active parts, one with the cationic [MoO₂(HL)]⁺ complex and the other one with the anionic [Mo₆O₁₉]²⁻. The different interaction between both entities could potentially influence their reactivity. To test the catalytic activity of the prepared catalysts in the cyclooctene epoxidation reactions, the following parameters were studied: (i) influence of the oxidant (TBHP in water vs. TBHP in decane) and (ii) addition of the organic solvent (CH₃CN) to the reaction media. The purpose of parameter variations was to evaluate the impact of green chemistry principles, based on the use of ecologically more appropriate oxidant (TBHP in water) and no addition of organic solvent to the reaction media, on the catalytic activity of the tested catalysts.

Relevant catalytic parameters are summarized in Table 3. As expected, the catalytic results provided by the use of TBHP in decane are the best ones and classify all the prepared catalysts as very active, conversion of cyclooctene >93%, and highly selective towards epoxide, >82% (Table 3 and Fig. S5, see ESI†). This is following the fact that solubility in the organic phase of the MoO₂(HL) moiety is favoured and the expected reactivity goes mainly through the transition metal complex. It is worth noting that epoxide can undergo further transformation into the corresponding diol, which is predominantly found in the water phase but also partially in the organic phase. Our analysis *via* GC primarily targets the organic layer of the reaction mixture, where the substrate, oxidizing agent, and desired product naturally reside. Given that TBHP is utilized in the water phase, the diol would be distributed across both phases, potentially affecting the accuracy of our determinations. Hence, considering the various oxidizing agents tested, our research was tailored specifically towards the epoxide.

The activity of the complex 3H·2CH₃CN should be highlighted, Fig. 6 and 7. TOF₂₀min is extremely high (value of 6546), implying very fast and abrupt activation of the complex towards active species. The observation correlates with the results obtained for the theoretical calculations above, where 3H·2CH₃CN is less stable than 3H, faster transforming into catalytically active species. For all the catalysts the first 60 min of the reaction were carefully followed and different trends can be observed. For the reactions with the complex 2Hw·2CH₃COCH₃ induction time could be distinguished, implying slower activation than in the case of other studied catalysts. This could be due to the presence of the coordinated water

molecule to molybdenum. The same type of phenomenon was seen in the case of 1Hw, with the slowest activity in the first 20 min of the reaction (Table 3).

In general, the use of TBHP in water, as environmentally friendly media, resulted in very good activity and high

<table><caption>Table 3 Relevant catalytic results for the cyclooctene epoxidation. Reaction conditions: time, 6 h; temperature, 80 °C, catalyst/cyclooctene/TBHP molar ratio: 0.025/100/200 for all compounds</caption><thead><tr><th>Catalyst</th><th>Conversion<sup>a</sup>/%</th><th>Selectivity<sup>b</sup> (%)</th><th>TOF<sup>c</sup></th><th>TON<sup>d</sup></th></tr></thead><tbody><tr><td colspan="5"><strong>Oxidant: TBHP in decane, no organic solvent added</strong></td></tr><tr><td>1H<sup>a</sup></td><td>93</td><td>89</td><td>952</td><td>1233</td></tr><tr><td>1H<sup>w</sup></td><td>97</td><td>87</td><td>1090</td><td>1292</td></tr><tr><td>1H</td><td>99</td><td>82</td><td>1605</td><td>1332</td></tr><tr><td>2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>99</td><td>85</td><td>876</td><td>1325</td></tr><tr><td>2H</td><td>99</td><td>83</td><td>1256</td><td>1329</td></tr><tr><td>2O</td><td>98</td><td>86</td><td>850</td><td>1300</td></tr><tr><td>3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>99</td><td>84</td><td>1441</td><td>1313</td></tr><tr><td>3H·2CH<sub>3</sub>CN</td><td>99</td><td>84</td><td>6546</td><td>1329</td></tr><tr><td>3H</td><td>99</td><td>65</td><td>1520</td><td>1321</td></tr><tr><td colspan="5"><strong>Oxidant: TBHP in water, no organic solvent added</strong></td></tr><tr><td>1H<sup>a</sup>.</td><td>45</td><td>54</td><td>413</td><td>582</td></tr><tr><td>1H<sup>w</sup></td><td>78</td><td>95</td><td>245</td><td>1032</td></tr><tr><td>1H</td><td>78</td><td>68</td><td>800</td><td>1034</td></tr><tr><td>2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>87</td><td>95</td><td>570</td><td>1164</td></tr><tr><td>2H</td><td>85</td><td>75</td><td>570</td><td>1121</td></tr><tr><td>2O</td><td>61</td><td>87</td><td>408</td><td>706</td></tr><tr><td>3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>87</td><td>84</td><td>492</td><td>1157</td></tr><tr><td>3H·2CH<sub>3</sub>CN</td><td>98</td><td>96</td><td>587</td><td>1309</td></tr><tr><td>3H</td><td>97</td><td>80</td><td>494</td><td>1268</td></tr><tr><td colspan="5"><strong>Oxidant: TBHP in water, CH<sub>3</sub>CN added</strong></td></tr><tr><td>1H<sup>a</sup>.</td><td>41</td><td>44</td><td>410</td><td>548</td></tr><tr><td>1H<sup>w</sup></td><td>69</td><td>71</td><td>646</td><td>926</td></tr><tr><td>1H</td><td>80</td><td>71</td><td>264</td><td>1067</td></tr><tr><td>2H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>63</td><td>97</td><td>966</td><td>836</td></tr><tr><td>2H</td><td>53</td><td>87</td><td>568</td><td>710</td></tr><tr><td>2O</td><td>52</td><td>59</td><td>503</td><td>705</td></tr><tr><td>3H<sup>w</sup>·2CH<sub>3</sub>COCH<sub>3</sub></td><td>65</td><td>95</td><td>500</td><td>879</td></tr><tr><td>3H·2CH<sub>3</sub>CN</td><td>84</td><td>84</td><td>560</td><td>1119</td></tr><tr><td>3H</td><td>96</td><td>56</td><td>512</td><td>1098</td></tr></tbody></table>

<footnote>a Conversion of cyclooctene = n(cyclooctene converted)/n(cyclooctene engaged) calculated after 6 h. b n(Cyclooctene oxide formed)/n(cyclooctene converted) after 6 h. c TOF = n(cyclooctene transformed)/3 × n(catalyst)/time at 20 minutes. All potential “active species” were considered, i.e. three potential active “sites” per species. The concept of three potential sites was adopted from our earlier published research, which considers two active sites originating from two mononuclear coordination complexes, along with one active site arising from POM.⁵¹ d TON is defined at the end of the reaction n(cyclooctene transformed)/3 × n(catalyst).</footnote>
<figure><img src="image_17.png" /><figcaption>Fig. 6 Kinetic profile for cyclooctene conversion when using TBHP in decane as oxidant with the catalysts 3H·2CH₃CN (red), 3Hw (blue) and 3H (black).</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>Fig. 7 Kinetic profile for cyclooctene conversion when using TBHP in water as oxidant with the catalysts 3H·2CH₃CN (red), 3Hw (blue) and 3H (black).</figcaption></figure>

selectivity, for most of the tested catalysts, justifying the green chemistry concept pursued in this research (Fig. 7). The complex 3H, obtained from H₂L³ ligand, operates the best under eco-friendly conditions. It can be assumed that the water is promoting epoxide ring opening and the formation of the diol that was not quantified through GC measurements since only the organic layer was analysed. On the other hand, the addition of CH₃CN surprisingly inhibited the catalytic activity, but the selectivity towards epoxide did not seem strongly affected. The reverse observation was noted for the tungsten metal complexes coordinated by a similar class of ligands.⁵⁰

Having in mind that the self-assembly of organic–inorganic POM-based salts in the presence of hydrazone ligands is rare, we compared the catalytic performance of the investigated compounds with the one previously reported, [MoO₂(HL)]₂Mo₆O₁₉·2CH₃CN, where HL presents monodeprotonated ligand obtained by the condensation of pyridoxal and nicotinic acid hydrazide.⁵¹ The compound [MoO₂(HL)]₂Mo₆O₁₉·2CH₃CN provides the cyclooctene conversion of 72% with selectivity towards epoxide of 58% when using TBHP in water, without adding organic solvent. For the presented

investigation, no clear distinction and correlation could be concluded for the nature of ligands (position of OMe group) influencing the catalytic performance of the tested catalysts. Further, all the tested complexes in this research show much better activity and selectivity than the previously tested “referent” system (Bu₄N)₂Mo₆O₁₉,⁵¹ justifying the role of [MoO₂L] species (conversion parameter 30%, selectivity towards epoxide 53%, TOF 42, TON 123, under the same reaction conditions).

Further comparison was made considering Mo tetranuclear and polynuclear coordination compounds obtained from H₂L¹⁻³ ligands (Table S7, see ESI†).²⁶,³¹ With TBHP in decane conversion was 74% and selectivity towards epoxide 85%, while with TBHP in water, the investigated parameters were much lower (27% and 56%, respectively). In this case, it seems that the presence of polyoxomolybdate anion positively influences the catalytic performance of the tested catalyst, supporting the fact that the Lindqvist-type anion is one active site of the catalyst.

# Experimental section

## Materials and instruments

The commercially available acetylacetone, nicotinic acid hydrazide, 2-hydroxy-3-methoxybenzaldehyde, 2-hydroxy-4-methoxybenzaldehyde, 2-hydroxy-5-methoxybenzaldehyde, and organic solvents were purchased from Alfa Aesar or Aldrich. Unless otherwise mentioned, all solvents were used without further purification. [MoO₂(acac)₂] (acac = acetylacetonate),⁵² and ligands²⁶,⁴⁷ were synthesised following the previously published procedures. Dry solvents were used for the solvothermal reactions. Acetonitrile was dried by the use of phosphorus pentoxide and then distilled. Acetone was dried with anhydrous potassium carbonate for at least 48 hours and then distilled prior to use. Light was excluded from all reactions using foil to avoid Mo reduction.

Attenuated Total Reflectance Infrared (ATR-IR) spectra were recorded at room temperature on a PerkinElmer Spectrum One spectrometer fitted with an ATR reflectance attachment. All data were processed in OMNIC software. NMR data for organic compounds was recorded on a Bruker Avance 400 MHz at room temperature.

Thermal stability was examined by conducting thermogravimetric analysis (TGA) on a Mettler-Toledo TGA/SDTA851 at a heating rate or 5 °C min⁻¹. The samples were heated from 25 °C to 600 °C in an oxygen atmosphere (flow rate: 20 mL min⁻¹) using Al₂O₃ crucibles. The results were developed by applying the Mettler STARᵉ evaluation software. Elemental analysis (C, H, and N) was determined by the Analytical Services Laboratory of the Ruder Bošković Institute, Zagreb.

## X-ray crystallography

The powder X-ray diffraction data were collected using a Malvern Panalytical Aeris powder diffractometer in the Bragg-Brentano geometry with a PIXcel¹D detector, using CuKα radiation (λ = 1.5418 Å). Samples were contained on a Si sample holder. Powder patterns were collected at room temperature (5 <

<figure><img src="image_19.png" /></figure>
<figure><img src="image_20.png" /></figure>

$2\Theta < 40^{\circ}$). The data were collected and visualized utilizing the Malvern Panalytical HighScore Software Suite.⁵³

High-quality single crystals of 1H<sup>a</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 2H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 3H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 1H<sup>w</sup>, 1H, 2H, 2O.2CH<sub>3</sub>CN·H<sub>2</sub>O, 3H, 3H.2CH<sub>3</sub>CN were grown by crystallization from the corresponding solutions. Diffracted intensities were collected on Rigaku XtaLAB Synergy diffractometer equipped with Dualflex source using CuK<sub>α</sub> radiation (λ = 1.54184 Å) for 1H<sup>a</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 3H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 2H, 2O.2CH<sub>3</sub>CN·H<sub>2</sub>O, 3H and 3H.2CH<sub>3</sub>CN, as well as MoK<sub>α</sub> (λ = 0.71073) for 2H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 1H<sup>w</sup>, 1H and HyPix detector via ω-scans. The crystals were kept at 150 K or 170 K during data collection (see Table S1 in the ESI†). Data were processed using the CrysAlis program package. A summary of general and crystal data, intensity data collection and final refinement parameters are presented in Table S5 and S6 (see ESI†). The structures were solved with dual space methods using SHELXT. The refinement procedure by full-matrix least-squares methods based on *F*<sup>2</sup> 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 the use of the riding model with *U*<sub>iso</sub> = 1.2*U*<sub>eq</sub> of the connected carbon atom or as ideal CH<sub>3</sub> groups with *U*<sub>iso</sub> = 1.5*U*<sub>eq</sub>. Acetonitrile molecules in 2O.2CH<sub>3</sub>CN·H<sub>2</sub>O are heavily disordered and for one of them hydrogen atoms of the methyl group could not be satisfactory modelled and were thus omitted from the model. In the structure of 1H<sup>a</sup>.2CH<sub>3</sub>COCH<sub>3</sub> acetone molecule in the crystal lattice was disordered over two positions with equal occupancies 0.5. Hydrogen atoms attached to heteroatoms were located in the difference Fourier maps at the final stages of the refinement procedure except for one of the coordinated ligands in 2O.2CH<sub>3</sub>CN·H<sub>2</sub>O, where hydrogen atom on the N<sub>2</sub> was ascribed based on geometrical considerations and the charge, while considering also the most reasonable hydrogen bond geometry. All refinements were performed using SHELXL, which was operated within the Olex2 suite. Geometrical calculations and molecular graphics were done with Mercury and Platon.⁵⁴,⁵⁵

## Catalytic studies

The Agilent 7820A chromatograph, outfitted with a Flame Ionization Detector (FID) and an HP5-MS capillary column (dimensions: 30 m × 0.32 mm × 0.25 µm), was utilized for chromatogram production. Parameters for the gas chromatography (GC) analysis were established through the examination of genuine reactant and product samples. The analysis focused on quantifying the transformation of olefins into their respective epoxides, employing calibration curves for precise measurement. Acetophenone served as the internal standard throughout this process, ensuring accuracy in quantification.

## Quantum chemical calculations

Optimizations of geometries for all computed complexes were performed using the hybrid functional B3LYP⁵⁶ with the D3 version of Grimme's dispersion⁵⁷ in combination with the Def2TZVP⁵⁸,⁵⁹ basis set. Initial geometries were taken from crystallographically determined structures. To confirm that the

obtained geometries were local minima, harmonic frequency calculations were performed.⁶⁰,⁶¹ The standard Gibbs energies of formation were calculated at *T* = 298.15 K and *p* = 101.325 Pa. Strengths of the intramolecular interactions were estimated by the calculation of relative differences in the standard Gibbs energies of binding Δ<sub>b</sub>G<sub>rel</sub> according to the formula:

$$ \Delta_b G^\circ_{\text{rel}}(\text{MoLS}_1) - \Delta_b G^\circ_{\text{rel}}(\text{MoLS}_2) = \Delta_f G^\circ(\text{MoLS}_1) - \Delta_f G^\circ(\text{S}_1) \\ -\Delta_f G^\circ(\text{MoLS}_2) + \Delta_f G^\circ(\text{S}_2) $$

All quantum chemical calculations were carried out using the Gaussian 16 program package.⁶²

**Procedure I – synthesis of 1H<sup>a</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, 2H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, and 3H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>.** Hydrazone H<sub>2</sub>L<sup>1–3</sup> (0.05 g, 1.84 × 10<sup>−4</sup> mol) was added to a solution of [MoO<sub>2</sub>(acac)<sub>2</sub>] (0.2 g, 6.13 × 10<sup>−4</sup> mol) in acetone (80 mL), and the reaction mixture was left to stir at room temperature for three hours and left overnight. The deposited solid was separated, the filtrate evaporated under vacuum to 1/3 of its volume and left at room temperature for about one week. Slow crystallization yielded crystals of hybrids which were separated by filtration.

[MoO<sub>2</sub>(HL<sup>1</sup>)(CH<sub>3</sub>COCH<sub>3</sub>)]<sub>2</sub>[Mo<sub>6</sub>O<sub>19</sub>].2CH<sub>3</sub>COCH<sub>3</sub> (1H<sup>a</sup>.2CH<sub>3</sub>COCH<sub>3</sub>). Orange product was obtained. Isolated 59 mg, 35.8% yield based on H<sub>2</sub>L<sup>1</sup>. Calc. for 1H<sup>a</sup>, C<sub>34</sub>H<sub>36</sub>Mo<sub>8</sub>N<sub>6</sub>O<sub>31</sub> (1792.19): C, 22.79; H, 2.02; N, 4.69%. Found: C, 22.54; H, 1.78; N, 4.36%. TG: Calc. CH<sub>3</sub>COCH<sub>3</sub>, 6.48; MoO<sub>3</sub>, 64.25%. Found: CH<sub>3</sub>COCH<sub>3</sub>, 6.66; MoO<sub>3</sub>, 63.66%. Selected IR data (cm<sup>−1</sup>): 1632 (CO), 1614 (CN), 1348 (CO), 1256 (CO), 957 (MoO<sub>4</sub>)<sub>POM</sub>, 942, 920, 908 (MoO<sub>2</sub>), 777, 755 (MoO<sub>3</sub>), 600 (δ(O<sub>b</sub>MoO<sub>4</sub>)).

[MoO<sub>2</sub>(HL<sup>2</sup>)(H<sub>2</sub>O)]<sub>2</sub>[Mo<sub>6</sub>O<sub>19</sub>].2CH<sub>3</sub>COCH<sub>3</sub> (2H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>). Dark red product was obtained. Isolated 71 mg, 42.3% yield based on H<sub>2</sub>L<sup>2</sup>. Calc. for 2H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, C<sub>34</sub>H<sub>40</sub>Mo<sub>8</sub>N<sub>6</sub>O<sub>33</sub> (1828.22): C, 22.34; H, 2.21; N, 4.60%. Found: C, 22.12; H, 1.98; N, 4.43%. TG: Calc. CH<sub>3</sub>COCH<sub>3</sub>, 6.35; Calc. H<sub>2</sub>O, 1.97; MoO<sub>3</sub>, 62.99%. Found: CH<sub>3</sub>COCH<sub>3</sub> + H<sub>2</sub>O, 8.24; H<sub>2</sub>O; MoO<sub>3</sub>, 62.99%. Selected IR data (cm<sup>−1</sup>): 1669 (CO), 1607 (CN), 1357 (CO), 1225 (CO), 959 (MoO<sub>4</sub>)<sub>POM</sub>, 939, 926, 908 (MoO<sub>2</sub>), 789, 761 (MoO<sub>3</sub>), 591 (δ(O<sub>b</sub>MoO<sub>4</sub>)).

[MoO<sub>2</sub>(HL<sup>3</sup>)(H<sub>2</sub>O)]<sub>2</sub>[Mo<sub>6</sub>O<sub>19</sub>].2CH<sub>3</sub>COCH<sub>3</sub> (3H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>). Dark red product was obtained. Isolated 59 mg, 35.1% yield based on H<sub>2</sub>L<sup>3</sup>. Calc. for 3H<sup>w</sup>.2CH<sub>3</sub>COCH<sub>3</sub>, C<sub>34</sub>H<sub>40</sub>Mo<sub>8</sub>N<sub>6</sub>O<sub>33</sub> (1828.22): C, 22.34; H, 2.21; N, 4.60%. Found: C, 22.10; H, 2.12; N, 4.45%. TG: Calc. CH<sub>3</sub>COCH<sub>3</sub>, 6.35; Calc. H<sub>2</sub>O, 1.97; MoO<sub>3</sub>, 62.99%. Found: CH<sub>3</sub>COCH<sub>3</sub> + H<sub>2</sub>O, 8.03; MoO<sub>3</sub>, 62.87%. Selected IR data (cm<sup>−1</sup>): 1666 (CO), 1605 (CN), 1338 (CO), 1244 (CO), 957 (MoO<sub>4</sub>)<sub>POM</sub>, 927, 903, 895 (MoO<sub>2</sub>), 776, 712 (MoO<sub>3</sub>), 584 (δ(O<sub>b</sub>MoO<sub>4</sub>)).

**Procedure II – synthesis of 1H<sup>w</sup>, 2O·CH<sub>3</sub>CN·H<sub>2</sub>O, and 3H·2CH<sub>3</sub>CN.** Hydrazone H<sub>2</sub>L<sup>1–3</sup> (0.054 g; 0.2 mmol) was added to a solution of [MoO<sub>2</sub>(acac)<sub>2</sub>] (0.65 g; 2.0 mmol) in a mixture of 80 mL acetonitrile and 200 µL of acetic acid. The reaction mixture was left to stir at room temperature for 3 hours and left overnight. The deposited solid was separated by filtration after 2 days, and the clear solution left at room temperature. Slow crystallization yielded crystals within two weeks.
[**MoO₂(HL¹)(H₂O)**]₂[**Mo₆O₁₉**] (**1H**<sup>W</sup>). Dark red product was obtained. Isolated 32 mg, 18.1% yield based on H₂L¹. Calc. for **1H**<sup>W</sup>, C₂₈H₂₈Mo₈N₆O₃₁ (1712.06): C, 19.64; H, 1.65; N, 4.91%. Found: C, 19.54; H, 1.32; N, 4.77%. TG: Calc. H₂O, 2.10; MoO₃, 67.25%. Found: H₂O, 2.42; MoO₃, 66.89%. Selected IR data (cm⁻¹): 1607 (CN), 1349 (CO), 1261 (CO), 950 (MoO₇ₜPOM, 938, 922, 909, 886 (MoO₂), 776 759, 730 (MoO₆b), 600 (δ(O₆bMoO₇t)).

[**{MoO₂(HL²)}₃Mo₈O₂₆**]·**CH₃CN·H₂O** (2O·**CH₃CN·H₂O**). Brownish-orange product was obtained. Isolated 25 mg, 15.8% yield based on H₂L². Calc. for 2O, C₄₂H₃₇Mo₁₁N₉O₄₁ (2379.11): C, 21.20; H, 1.57; N, 5.30%. Found: C, 21.03; H, 1.24; N, 5.13%. TG: Calc. MoO₃, 66.55%. Found: MoO₃, 66.43%. Selected IR data (cm⁻¹): 1616 (CN), 1354 (CO), 1228 (CO), 948 (MoO₇ₜPOM, 934, 907 (MoO₂), 856, 838 (O=Mo-O), 734, 689, 664, 649, 642 (MoO₆b), 555 (δ(O₆bMoO₇t)).

[**{MoO₂(HL³)}₂Mo₆O₁₉**]·**2CH₃CN** (**3H·2CH₃CN**). Dark red product was obtained. Isolated 68.5 mg, 40.7% yield based on H₂L³. Calc. for **3H·2CH₃CN**, C₃₂H₃₀Mo₈N₈O₂₉ (1757.87): C, 21.86; H, 1.72; N, 6.37%. Found: C, 21.68; H, 1.53; N, 4.72%. TG: Calc. **CH₃CN**, 4. 88; Calc. MoO₃, 65.50%. Found: **CH₃CN**, 4. 84; MoO₃, 65.28%. Selected IR data (cm⁻¹): 1607 (CN), 1343 (CO), 1271 (CO), 959, 951 (MoO₇ₜPOM, 916, 907, 891 (MoO₂), 777 (MoO₆b), 575 (δ(O₆bMoO₇t)).

**Procedure III – synthesis of 1H, 2H, and 3H.** Procedure I was repeated. The obtained crystals of **1H**·2**CH₃COCH₃**, **2H**<sup>W</sup>·2**CH₃COCH₃**, and **3H**<sup>W</sup>·2**CH₃COCH₃** were filtered in glove box and transferred in 30 mL Teflon-lined autoclave together with 25 mL of dry acetonitrile or acetone and acetic acid (100 µL) and heated at 125 °C for 24 h. Crystals were obtained overnight.

[**{MoO₂(HL¹)}₂Mo₆O₁₉**] (**1H**). Red product was obtained from acetonitrile. Isolated 39 mg, 25.2% yield based on H₂L¹. Calc. for **1H**, C₂₈H₂₄Mo₈N₆O₂₉ (1676.03): C, 20.07; H, 1.44; N, 5.01%. Found: C, 19.79; H, 1.15; N, 4.87%. TG: Calc. MoO₃, 68.70%. Found: MoO₃, 68.50%. Selected IR data (cm⁻¹): 1604 (CN), 1350 (CO), 1265 (CO), 952 (MoO₇ₜPOM, 918, 901 (MoO₂), 869 (O=Mo-O), 782 754, 731 (MoO₆b), 577 (δ(O₆bMoO₇t)).

[**{MoO₂(HL²)}₂Mo₆O₁₉**] (**2H**). Red product was obtained from acetonitrile. Isolated 37 mg, 23.9% yield based on H₂L². Calc. for **2H**, C₂₈H₂₄Mo₈N₆O₂₉ (1676.03): C, 20.07; H, 1.44; N, 5.01%. Found: C, 19.90; H, 1.16; N, 4.92%. TG: Calc. MoO₃, 68.70%. Found: MoO₃, 68.40%. Selected IR data (cm⁻¹): 1606 (CN), 1354 (CO), 1292 (CO), 955, 951 (MoO₇ₜPOM, 931, 908, 891 (MoO₂), 862 (O=Mo-O), 761, 729 (MoO₆b), 558 (δ(O₆bMoO₇t)).

[**{MoO₂(HL³)}₂Mo₆O₁₉**] (**3H**). Red product was obtained from acetone. Isolated 32 mg, 20.0% yield based on H₂L³. Calc. for **3H**, C₂₈H₂₄Mo₈N₆O₂₉ (1676.03): C, 20.07; H, 1.44; N, 5.01%. Found: C, 19.87; H, 1.22; N, 4.91%. TG: Calc. MoO₃, 68.70%. Found: MoO₃, 68.42%. Selected IR data (cm⁻¹): 1609 (CN), 1357 (CO), 1277 (CO), 953 (MoO₇ₜPOM, 921, 910, 886 (MoO₂), 866 (O=Mo-O), 764, 722 (MoO₆b), 573 (δ(O₆bMoO₇t)).

General procedure for the epoxidation of cyclooctene using
TBHP (in water or decane)

In a typical reaction, cyclooctene (2.76 mL, 20 mmol), aqueous TBHP (70% w/w, 40 mmol) or TBHP in decane (5.5 mol dm⁻³, 40 mmol) and acetophenone (0.1 mL) and molybdenum (pre)catalyst,

0.005 mmol, were used, and the reaction mixture was stirred for 5 h at 80 °C. The reaction was monitored by withdrawing small aliquots at certain time intervals and analyzing them by GC.

General procedure for the epoxidation of cyclooctene using
TBHP in water with the addition of CH₃CN

The procedure was the same as without CH₃CN, except the 2 mL of CH₃CN was added to the reaction mixture before the catalyst was added.

Conclusions

Mo-POM organic-inorganic hybrids were successfully prepared by conventional and solvothermal methods. In the Class I hybrids [MoO₂(HL¹⁻³)(D)]₂[Mo₆O₁₉], POM anions tend to associate with cationic Mo-complexes *via* multiple non-covalent interactions. However, in Class II architectures [{MoO₂(HL¹⁻³)}₂Mo₆O₁₉], POM serves as a sterically hindered ligand and completes the coordination sphere of two Mo atoms from peripheral coordination subunits either through the terminal or bridging oxygen atoms. **2H** and **3H** are the first examples of polyoxometalates that incorporate bridging POM coordinated architectures. The preference of Mo-complex cations for coordination with solvent and terminal or bridging oxygen groups on the POM moiety was estimated by calculation of complexes' relative standard Gibbs energies of binding. Density functional theory was utilized for the calculation of standard Gibbs energies of formation, which were used for the estimation of binding energies. In most of the complexes, there is a clear preference for coordination with water molecules, whereas the bridging oxygen group of the POM is preferred in coordination with POM moiety.

The assembly process is additionally sensitive to the ligand's methoxy group position. Thus, the condensation with H₂L² resulted in the unsymmetrical octamolybdate-based hybrid [{MoO₂(H₂L²)}{MoO₂(HL²)}]Mo₈O₂₆{MoO₂(HL²)}], covalently expanded by three Mo-hydrazone complex units. All Class II hybrids have long bonds Mo-O-POM due to sterical hindrance.

In this study, the catalytic activity of all compounds was investigated and compared. With TBHP in decane as an oxidant, all the prepared catalysts showed excellent catalytic activity towards the cyclooctene oxidation, where >90% conversion and >80% selectivity could be obtained. Notably, **3H·2CH₃CN** catalyst operates the best under eco-friendly conditions. The addition of CH₃CN surprisingly inhibited the catalytic activity when TBHP in water is used as oxidant, but the selectivity towards epoxide remained almost the same.

Data availability

The data supporting this article have been included as part of the ESI.† Crystallographic data for **1H**·2**CH₃COCH₃**, **1H**<sup>W</sup>, **2H**<sup>W</sup>·2**CH₃COCH₃**, **2H**, **3H**<sup>W</sup>·2**CH₃COCH₃**, **2O·CH₃CN·H₂O**, **3H**, **3H·2CH₃CN** has been deposited at the Cambridge Structural Data base under with deposition numbers CCDC 2335829–2335837.
## Author contributions

Jana Pisk: investigation, supervision, formal analysis, writing – original draft preparation. Mirna Mandarić: investigation and formal analysis. Tomica Hrenar: quantum chemical calculations, data curation, and writing – original draft, visualization. Dominique Agustin: supervision, review and editing. Mirta Rubčić: investigation, writing – original draft preparation, review and editing. Višnja Vrdoljak: conceptualization, supervision, investigation, writing – original draft preparation, visualization, review and editing.

## Conflicts of interest

There are no conflicts to declare.

## Acknowledgements

This work has been fully supported by the Croatian Science Foundation under the project (IP-2016-06-4221). We acknowledge the support of project CIuK co-financed by the Croatian Government and the European Union through the European Regional Development Fund–Competitiveness and Cohesion Operational Programme (Grant KK.01.1.1.02.0016). LCC CNRS and Chemistry Dept of IUT Paul Sabatier are acknowledged for providing the equipment for catalysis experiments. We thank Matea Pajski and Silvija Mrkonja for their technical help with catalytic studies.

## Notes and references

1. J. J. Walsh, A. M. Bond, R. J. Forster and T. E. Keyes, Coord. Chem. Rev., 2016, **306**, 217.

2. Y. Dong, Y. Yang, L. Qiu, G. Dong, D. Xia, X. Liu, M. Li and R. Fan, ACS Appl. Energy Mater., 2019, **2**, 4224.

3. A. Bijelic, M. Aureliano and A. Rompel, Angew. Chem., Int. Ed., 2019, **58**, 2980.

4. L. S. van Rompuy and T. N. Parac-Vogt, Curr. Opin. Biotechnol., 2019, **58**, 92.

5. E. Coronado, C. Gimenez-Saiz, C. J. Gomez-Garcia and S. C. Capelli, Angew. Chem., Int. Ed., 2004, **43**, 3022.

6. S.-S. Wang and G.-Y. Yang, Chem. Rev., 2019, **115**, 4893.

7. S. Roy, V. Vemuri, S. Maiti, K. S. Manoj, U. Subbarao and S. C. Peter, Inorg. Chem., 2018, **57**, 12078.

8. K. Kamata, K. Yonehara, Y. Sumida, K. Yamaguchi, S. Hikichi and N. Mizuno, Science, 2003, **300**, 964.

9. N. Mizuno, K. Yamaguchi and K. Kamata, Coord. Chem. Rev., 2005, **249**, 1944.

10. L. Salles, J.-Y. Piquemal, R. Thouvenot, C. Minot and J.-M. Brégeault, J. Mol. Catal. A: Chem., 1997, **117**, 375.

11. L. Jing, J. Shi, F. Zhang, Y. Zhong and W. Zhu, Ind. Eng. Chem. Res., 2013, **52**, 10095.

12. D. Sloboda-Rozner, P. L. Alsters and R. Neumann, J. Am. Chem. Soc., 2003, **125**, 5280.

13. N. Mizuno and M. Misono, Chem. Rev., 1998, **98**, 199.

14. Y. M. A. Yamada, H. Tabata, M. Ichinohe, H. Takahashi and S. Ikegami, Tetrahedron, 2004, **60**, 4087.

15. D. Schaming, C. Allain, R. Farha, M. Goldmann, S. Lobstein, A. Giraudeau, B. Hasenknopf and L. Ruhlmann, Langmuir, 2010, **26**, 5101.

16. Y. Wang, X. Liu, W. Xu, Y. Yue, B. Li and L. Wu, Inorg. Chem., 2017, **56**, 7019.

17. H. Karoui and C. Ritchie, New J. Chem., 2018, **42**, 25.

18. D. Schaming, C. Costa-Coquelard, I. Lampre, S. Sorgues, M. Erard, X. Liu, J. Liu, L. Sun, J. Canny, R. Thouvenot and L. Ruhlmann, Inorg. Chim. Acta, 2010, **363**, 2185.

19. H. Kwen, V. G. Young Junior and E. A. Maatta, Angew. Chem., Int. Ed., 1999, **38**, 1145.

20. S. Schönweiz, M. Heiland, M. Anjass, T. Jacob, S. Rau and C. Streb, Eur. J. Chem., 2017, **23**, 15370.

21. N. K. Mishra, D. Bansal and S. Supriya, ACS Omega, 2022, **7**(35), 31403.

22. C. Yao, L.-K. Yan, W. Guan, C.-G. Liu, P. Song and Z.-M. Su, Dalton Trans., 2010, **39**, 7645.

23. R. M. Smith, I. Colliard, M. Amiri, M. A. Galindo and M. Nyman, Cryst. Growth Des., 2022, **22**, 2294.

24. A. Yokoyama, T. Kojima, K. Ohkubo, M. Shiro and S. Fukuzumi, J. Phys. Chem. A, 2011, **115**, 986.

25. J.-H. Li, X.-L. Wang, G. Song, H.-Y. Lin, X. Wang and G.-C. Liu, Dalton Trans., 2020, **49**, 1265.

26. M. Mandarić, E. Topić, D. Agustin, J. Pisk and V. Vrdoljak, Int. J. Mol. Sci., 2024, **25**(3), 1503.

27. J. Pisk, D. Agustin, V. Vrdoljak and R. Poli, Adv. Synth. Catal., 2011, **353**, 2910.

28. R. Bikas, V. Lippolis, N. Noshiranzadeh, H. Farzaneh-Bonab, A. J. Blake, M. Siczek, H. Hosseini-Monfared and T. Lis, Eur. J. Inorg. Chem., 2017, **6**, 999.

29. D. Biswal, N. R. Pramanik, S. Chakrabati, M. G. B. Drew, B. Sarkar, M. R. Maurya, S. K. Mukherjee and P. Chowdhury, New J. Chem., 2017, **41**, 4116.

30. M. R. Maurya, L. Rana and F. Avecilla, Inorganica Chim. Acta, 2015, **429**, 138.

31. V. Vrdoljak, M. Mandarić, T. Hrenar, I. Đilović, J. Pisk, G. Pavlović, M. Cindrić and D. Agustin, Cryst. Growth Des., 2019, **19**, 3000.

32. D. Cvijanović, J. Pisk, G. Pavlović, D. Šišak-Jung, D. Matković-Čalogović, M. Cindrić, D. Agustin and V. Vrdoljak, New J. Chem., 2019, **43**, 1791.

33. J. Pisk, M. Rubčić, D. Kuzman, M. Cindrić, D. Agustin and V. Vrdoljak, New J. Chem., 2019, **43**, 5531.

34. J. Pisk, D. Agustin and V. Vrdoljak, Catal. Commun., 2020, **142**, 106027.

35. V. Vrdoljak, B. Prugovečki, D. Matković-Čalogović and J. Pisk, CrystEngComm, 2011, **13**, 4382.

36. V. Vrdoljak, B. Prugovečki, D. Matković-Čalogović, R. Dreos, P. Siega and C. Tavagnacco, Cryst. Growth Des., 2010, **10**, 1373.

37. A. Dolbecq, A. Guirauden, M. Fourmigué, K. Boubekeur, P. Batail, M.-M. Rohmer, M. Bénard, C. Coulon, M. Sallé and P. Blanchard, J. Chem. Soc., Dalton Trans., 1999, 1241.

38. B. Gai, H. He, Y. Zhao, Z. Mao and H. Fu, Chem. Res. Chin. Univ., 2016, **32**, 527.

39. S. Wang, Z. Sun, C. Zhang, L. Ni, C. Wang, Y. Gao, L. Lv, J. Chang and W. Hao, Inorg. Chem. Commun., 2014, **41**, 47.

<figure><img src="image_21.png" /></figure>
<figure><img src="image_22.png" /></figure>

40 Y. Wang, Z. Zhang, E. Wang, Y. Qi and S. Chang, *Aust. J. Chem.*, 2008, **61**, 874.

41 M. Yuan, Y. Li, E. Wang, C. Tian, L. Wang, C. Hu, N. Hu and H. Jia, *Inorg. Chem.*, 2003, **42**, 3670.

42 H. Jin, Y. Qi, E. Wang, Y. Li, X. Wang, C. Qin and S. Chang, *Cryst. Growth Des.*, 2006, **6**, 2693.

43 C. R. Groom, I. J. Bruno, M. P. Lightfoot and S. C. Ward, The Cambridge Structural Database, *Acta Crystallogr.*, 2016, **B72**, 171.

44 E. Topić, I. Landripet, M. Duguin, J. Pisk, I. Đilović, V. Vrdoljak and M. Rubčić, *New J. Chem.*, 2020, **44**, 13357.

45 N. Bebić, E. Topić, M. Mandarić, T. Hrenar and V. Vrdoljak, *CrystEngComm*, 2021, **23**, 6349.

46 M. R. Maurya, S. Dhaka and F. Avecilla, *Polyhedron*, 2014, **67**, 145-215.

47 J. Pisk, T. Hrenar, M. Rubčić, G. Pavlović, V. Damjanović, J. Lovrić, M. Cindrić and V. Vrdoljak, *CrystEngComm*, 2018, **20**, 1804.

48 A. J. Bridgeman and G. Cavigliasso, *Chem. Phys.*, 2002, **279**, 143.

49 B. Modec, J. V. Brenčič and J. Zubieta, *J. Chem. Soc., Dalton Trans.*, 2002, 1500.

50 V. Vrdoljak, J. Pisk, B. Prugovečki, D. Agustin, P. Novak and D. Matković-Čalogović, *RSC Adv.*, 2016, **6**, 36384.

51 J. Pisk, B. Prugovečki, D. Matković-Čalogović, T. Jednačak, P. Novak, D. Agustin and V. Vrdoljak, *RSC Adv.*, 2014, **4**, 39000.

52 J.-J. Chen, J. W. McDonald and W. E. Newton, *Inorg. Chem.*, 1976, **15**, 2612.

53 T. Degen, M. Sadki, E. Bron, U. König and G. Nénert, The HighScore Suite, *Powder Diffr.*, 2014, **29**, S13.

54 C. F. Macrae, I. Sovago, S. J. Cottrell, P. T. A. Galek, P. McCabe, E. Pidcock, M. Platings, G. P. Shields,

J. S. Stevens, M. Towler and P. A. Wood, *J. Appl. Crystallogr.*, 2020, **53**, 226.

55 A. L. Spek, *Acta Crystallogr.*, 2009, **D65**, 148.

56 (a) C. Lee, W. Yang and R. G. Parr, *Phys. Rev. B: Condens. Matter Mater. Phys.*, 1998, **37**, 785; (b) R. Peverati and D. G. Truhlar, *Phys. Chem. Chem. Phys.*, 2021, **14**, 16187.

57 S. Grimme, S. Ehrlich and L. Goerigk, *J. Comput. Chem.*, 2011, **32**, 1456.

58 F. Weigend and R. Ahlrichs, *Phys. Chem. Chem. Phys.*, 2005, **7**, 3297.

59 F. Weigend, *Phys. Chem. Chem. Phys.*, 2006, **8**, 1057.

60 T. Hrenar, I. Primožič, D. Fijan and M. Majerić Elenkov, *Phys. Chem. Chem. Phys.*, 2017, **19**, 31706.

61 I. Primožič, T. Hrenar, K. Baumann, L. Krišto, I. Križić and S. Tomić, *Croat. Chem. Acta*, 2014, **87**, 153.

62 M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery Jr, J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman and D. J. Fox, *Gaussian 16, Revision A.03*, Gaussian, Inc., Wallingford CT, 2016.