Heterobimetallic complexes composed of bismuth and lithium carboxylates as polyurethane catalysts – alternatives to organotin compounds †

Organotin compounds are important catalysts for the synthesis of polyurethanes consisting of aliphatic isocyanates and polyols. Due to their toxicity, however, it has been a long-standing goal to develop more environmentally benign catalysts. Bismuth octoates and neodecanoates are such well-known alternatives, but their catalytic activity is insu ﬃ cient for many applications. Herein we show that the catalytic activity of bismuth carboxylates can be enhanced signi ﬁ cantly by the addition of lithium carboxylates. Structural and spectroscopic results reveal the spontaneous formation of heterobimetallic complexes consisting of two bismuth and four lithium carboxylates that show dynamic behaviour in solution. The mechanism of the bismuth-catalyzed urethane reaction was elucidated in detail using quantum chemical calculations, paving the way for a rational catalyst design.


Introduction
Polyurethanes (PUs) represent one of the most important classes of polymer materials with numerous technical applications, including elastomers, rigid foams, soft foams, coatings and adhesives. 1 These materials account for 7.9 wt% of the world polymer 2 production due to their versatile applications at low cost. A key advantage of PU technology is the simple construction of material composites in a single processing step via a polyaddition reaction of diisocyanates and diols. 1 Depending on the desired end product, the polymerization is controlled by task-specific catalysts. 3,4 Tertiary amines and Lewis-acidic metal complexes are commonly employed to mediate this polyaddition reaction. Recent developments in organocatalytic PU formation include the use of amidines/guanidines, 5 NHCs 6-9 and organic acids. [10][11][12] For coating applications involving aliphatic polyisocyanate components, 13,14 catalysts with high catalytic activity at very low concentrations are required that accelerate the gelling reaction between isocya-nate and alcohol, but do not promote the isocyanate-water reaction as this would lead to blistering. Further important catalyst properties are sufficient stability against light and air as well as tolerance against influences from other ingredients of the coating formulation. Organotin compounds, in particular dibutyltin dilaurate (DBTL), remain the most important catalysts in aliphatic urethane synthesis due to their remarkable catalytic activity and advantageous properties. 4,15 However, owing to the toxicity of organotin compounds and the difficulties to remove them from the final products, 16,17 it has been a long-standing goal to develop more environmentally benign alternative catalysts. 18 In efforts to replace organotin catalysts, numerous Lewis-acidic metal complexes were identified as effective catalysts for the isocyanate-hydroxyl reaction, [19][20][21][22][23][24][25][26][27] including some inorganic tin complexes with amino/alkoxide ligands that exhibit thermolatent behaviour. [28][29][30][31][32][33][34][35] However, apart from DBTL, only a few metal complexes based on bismuth, aluminium, and zirconium are among the promising candidates for clear coat applications, as discoloration of the coating formulation is highly undesirable. 21,26,27 While the practical use of Zr catalyst is hampered by their sensitivity towards hydrolysis and low tolerance against functional groups, 27 bismuth carboxylates are recognized as non-toxic alternatives to DBTL catalysts. 15,[36][37][38][39] In addition, there is evidence that combinations of bismuth and zinc carboxylates or bismuth and lithium carboxylates can accelerate the drying time of polyurethane coating formulations. 40 3 . [42][43][44] The synthesis of bismuth carboxylates via salt metathesis reactions often results in salt contaminations due to the variable and high coordination number of bismuth(III) complexes. 45 We therefore prepared bismuth pivalate Bi( piv) 3 (1), bismuth 2,2-dimethylbutanoate Bi(dmb) 3 (2) and bismuth 2-phenylisobutyrate Bi ( pib) 3 (3) according to the procedure reported by Andrews. 44 Heating a mixture of BiPh 3 and stochiometric amounts of carboxylic acid at 110°C for 16 hours afforded bismuth carboxylates 1-3 in quantitative yield (Scheme 1). The progress of the reaction was monitored by 1 H NMR spectroscopy based on the formation of benzene and the decrease in BiPh 3 signals. Complexes 1-3 were obtained as white solids that are soluble in common polar and nonpolar solvents. In the solid-state structures of bismuth carboxylates, the bismuth centres are typically nine coordinate and comprise the carboxylates in chelating and bridging (µ 2 and µ 3 ) binding modes. 46 51,52 The almost identical steric bulk of the 2,2-dimethylbutanoate ligands and the observation that complex 2 is soluble in nonpolar solvents, suggests that 2 is also a tetramer in the solid state. The introduction of a phenyl group in 2-phenylisobutyrate further increases the steric bulk of the carboxylate ligand. An X-ray diffraction study of a single crystal grown from a saturated THF solution revealed that 3 is dimeric in the solid state ( Fig. S57 †). Each Bi atom is coordinated by three chelating pib ligands and one THF molecule. These monomers are connected by two pib ligands in tridentate bridging coordination mode forming a Bi 2 O 2 four-membered ring with long intermolecular Bi-O bonds of 2.87 Å. Next, we explored the influence of lithium carboxylates on the coordination behaviour of bismuth complexes 1-3. Heating a suspension of complex 1-3 and two equivalents of the corresponding lithium carboxylate in acetonitrile gave a clear solution, which indicates the formation of a heterobimetallic compound as Li( piv) is poorly soluble in acetonitrile. The heterobimetallic complexes 7-9 were readily obtained as white crystalline solid in excellent yield after allowing the reaction mixture to cool to room temperature (Scheme 1). Note that complexes 7-9 also crystallize from an equimolar mixture of bismuth and lithium carboxylates, which demonstrates the propensity for the formation of heterobimetallic complexes. The asymmetric and symmetric CO 2 stretching frequencies in the IR spectra of complexes 1-3 and 7-9 appear at very similar frequencies ( Table 1). The separation Δν of these bands can provide information about the coordination mode of the carboxylate ligand. 53 For complexes 1-3 and 7-9 Δν values between 160 and 186 cm −1 are obtained, which are in the typical range of chelating and bridging carboxylates (150-200 cm −1 ). 53 Single crystals X-ray diffraction studies revealed the formation of isostructural complexes 7, 7′, 8, 9 and 9′ with the general formula 2 ], consisting of two bismuth carboxylates, four lithium carboxylates and two acetonitrile ligands ( Fig. S58-S62 †). The solid-state structure of the pivalate complex is depicted in Fig. 1. Compound 7′ and 9′ are two polymorphs containing one acetonitrile solvate molecules in the asymmetric unit (see the ESI † for details). Their geometrical parameters are very similar to those of 7 and 9, respectively, and thus are not included in the discussion. Selected bond lengths for the heterobimetallic complexes 7, 8 and 9 are summarized in Table 2. In the hexanuclear complexes 7-9, two Bi(O 2 CR) 3 units are linked by a Li 4 (O 2 CR) 4 (MeCN) 2 ladder-type structure with a crystallographic inversion centre located in the Li 2 (O 2 CR) 2 six-membered ring. Each Bi is eight coordinate in a disordered square antiprismatic environment (Fig. 2)

NMR-study of heterobimetallic carboxylate complexes
To study the behavior of the Bi carboxylate complexes 1-3 and 7-9 in solution we performed 1 H, 13 C, 7 Li and 1 H DOSY NMR experiments. The 1 H and 13 C NMR spectra of Bi complexes 1-3 in CD 2 Cl 2 show only one set of signals for the carboxylate ligands. The characteristic 13 C signal of the carboxylate groups appears in the range of 190.0-187.6 ppm. By means of 1 H DOSY NMR spectroscopy the hydrodynamic radius r H of 1 was determined, which is in good agreement with the radius of the tetramer in the solid state (Table 3). However, three sets of signals would be expected for the carboxylate ligands of the tetramer according to the solid-state structure of 1. Thus, the NMR data clearly indicates that in solution complexes 1-3 form dynamic aggregates. Dynamic behavior was also observed for the heterobimetallic complexes 7-9 in dichloromethane solutions, as indicated by the appearance of only one set of signals for the carboxylate ligands in the 1 H or 13 C NMR spectra. The carboxylate 13 C NMR resonances of 7-9 (189.7-187.5 ppm) are significantly broadened and appear in the same range of those of 1-3. In the 7 Li NMR spectra two signals were observed for complexes 8   Table 3 Hydrodynamic radius r H from 1 H-DOSY-NMR and averaged radius from the molecular structure of Bi( piv) 3

(1) and heterobimetallic complex 7 in [Å]
Entry Bi(piv) 3 (1) Bi 2 Li 4 (piv) 10 (9) 3.181 (1) (0.01 ppm, −0.5 ppm) and 9 (−0.1 ppm, −1.5 ppm) while 7 shows only one signal at −0.1 ppm. The hydrodynamic radius r H of 7 is significantly larger than that of 1 and agrees with the radius for the presence of dynamic aggregates in solutions of 7-9 was provided by low temperature 1 H NMR experiments. Upon cooling a solution of 7 in CD 2 Cl 2 to 193 K, the singlet for the methyl groups split into more than eight resonances, indicating the presence of several different aggregates rather than the discrete complex 7 in solution (Fig. S30 †). Several signals for the methyl protons were also observed in the low temperature 1 H NMR spectra of complexes 8 and 9 ( Fig. S31 and S32 †), which, however, can also be attributed to the presence of several rotamers. It should be noted that Li( piv) is very poorly soluble in CD 3 CN, but dissolves quickly upon addition of Bi( piv) 3 . Mixtures of Bi( piv) 3 and Li( piv) with 1 : 1 and 1 : 2 stoichiometry in CD 3 CN both show a singlet at 1.13 ppm in the 1 H NMR spectra, which appears downfield compared to that of Bi ( piv) (1.06) and downfield relative to that of Li( piv) (1.16 ppm). Collectively, the NMR studies show that the addition of Li carboxylates to bismuth carboxylates leads to the spontaneous formation of heterobimetallic complexes in solution. However, these species are not well defined, but rather constitute a highly dynamic mixture of aggregates with an average size comparable to that of 7-9.

Catalysis of the isocyanate/alcohol reaction
The catalytic performance of complexes 1-3 and 7-9 in the urethane reaction was examined using ethanol and the aliphatic isocyanate 2-ethylhexyl-(6-isocyanatohexyl)-carbamates (Desmodur LD). The latter was chosen as model for 1,6-hexamethylene diisocyanate (HDI) owing to its lower volatility and toxicity. The urethane reactions were performed at ambient temperature with 0.05 mol% catalyst loading using xylene as solvent (Scheme 2). DBTL was included in our study as a reference. Using ATR-IR spectroscopy the reaction progress was monitored by measuring the decreasing intensity of the isocyanate band at 2260 cm −1 . The selectivity of the urethane reaction was verified by the simultaneous increase of the characteristic ν(CvO) (1692 cm −1 ) and ν(CHN) (1519 cm −1 ) frequencies of the alkyl urethane group. 57 From these experiments, summarized in Fig. 3, the Bi carboxylates 1-3 show moderate catalytic activity and accelerate the urethane reaction significantly compared to the non-cata-lyzed reaction. Among the bismuth compounds the catalytic activity was found to increase with the steric bulk of the carboxylate ligands pib > dmb > piv. Strikingly, the heterobimetallic complexes 7-9 clearly emerge as superior catalysts compared to their lithium free counterparts and compound 8 was found to be almost as active as DBTL, leading to 90% conversion after only 15 minutes.
We next wanted to explore the influence of lithium carboxylates on the catalytic performance of the bismuth complexes in more detail. Catalyst mixtures containing different amounts of bismuth and lithium carboxylates were prepared by adding the respective carboxylates to the xylene/ethanol mixture immediately before catalysis, one after the other (Fig. 4). Owing to the sufficient solubility of Li( pib) in xylene, the pib based catalyst system was chosen for this study. As expected, Li ( pib) does not accelerate the isocyanate conversion at room Scheme 2 Catalysis conditions for isocyanate alcohol reaction. Catalyst concentration is based on the amount of bismuth.  temperature. However, the catalytic activity of Bi( pib) 3 is slightly increased upon adding 0.5 equiv. of Li( pib), while the addition of two or more equiv. of Li( pib) resulted in reaction rates that are identical to that of complex 9, which suggests that complex 9 is the major species in these reaction mixtures and heterobimetallic compounds with higher Li content are not formed. Surprisingly, the best catalytic performance was observed using a 1 : 1 mixture of Bi( pib) 3 and Li( pib). Collectively, these results reveal the superior performance of the 1 : 1 mixture while the stability of complex 9 might be responsible for the lower reaction rate with higher lithium content.

Mechanistic investigation and density functional theory calculations
In an attempt to evaluate the reactivity of the urethane catalysts with alcohol and isocyanate functions, complexes 1 and 7 were mixed with stoichiometric amounts of ethanol or 1,6hexamethylene diisocyanate (HDI) and the resulting 1 H NMR spectra were compared to those of the individual compounds ( Fig. S33-S36 †). The 1 H NMR spectrum of the mixtures of 1 or 7 with HDI show the resonances of the individual components and thus do not indicate a reaction. The same is observed for the mixture of 1 and ethanol. However, the CH 2 resonance of ethanol in the mixture with 7 is significantly broadened and is not shifted to lower/higher frequency, suggesting that ethanol is involved in the dynamic formation of coordination compounds in solution, although the formation of free pivalic acid was not detected. 58,59 Preliminary kinetic studies on the bismuth pivalate catalysed urethane reaction reveal a 0.4 order kinetics with respect to ethanol (see the ESI †).
So far, mechanistic studies on metal complex catalyzed PU reaction have been performed mainly for organotin catalysts, 58,60-65 leaving critical questions still unanswered. It remains for example unclear what factors determine the selectivity and speed of the isocyanate-hydroxyl reaction when catalyzed by tin, zirconium, bismuth, or zinc catalysts. We investigated the Bi-catalyzed urethane formation by means of quantum chemical computations, providing the basis for further mechanistic investigations. These computations confirm that bismuth carboxylates successfully catalyze urethane formation for both, aliphatic and aromatic isocyanates, speeding up the reaction by a factor of 10 4 to 10 5 . Complex kinetics govern the bismuth pivalate catalysis, albeit partially following pseudo 0 th -order kinetics (cf. Fig. 5). The turnover-limiting step in the catalytic cycle is the migratory insertion of the isocyanate into the Bi-O bond of the alkoxy ligand, finally forming the carbamate that remains coordinated to the metal center in a κ 1 -or κ 2 -fashion. The urethane product is released from the catalyst by protonation at the carbamate nitrogen-or carbonyl oxygen atom. Bismuth pivalate serves as model for gaining conceptional mechanistic insight into the Bi-catalysed urethane formation. Propyl isocyanate and ethyl alcohol serve as model substrates. Further information on the condensed phase representation and computational details in general can be found in the ESI. † After briefly addressing the uncatalyzed isocyanate-hydroxyl reaction we will take a closer look at its mechanistic spectrum in the context of bismuth catalysis.
The uncatalyzed reaction between aliphatic alcohol and aliphatic or aromatic isocyanate requires a proton shuttle in the 6-membered ring transition-state (TS) structure to give urethane product in reasonable time (Scheme 3). In the absence of a suitable proton shuttle in the reaction mixture, the alcohol itself acts as such. This stoichiometrytwo alcohols and one isocyanateleads to pseudo second-order kinetics with a strong dependence on the concentration of the alcohol. The computed reaction barriers for our model system are 110 and 100 kJ mol −1 for aliphatic and aromatic isocyanates, respectively, in the infinite dilution scenario. The reaction rates are consistent with the experimental reaction times, which are in the range of days and hours, respectively, at ambient temperature. Note that the effective rate of the urethan bond forming reaction can be influenced by various phenomena such as catalysis, autocatalysis, solvent effects and molecular interactions. 63,[66][67][68][69] Computations reveal that cyclic, tetrameric bismuth tricarboxylate clusters form thermodynamically driven by a free energy difference of −40 and −18 kJ mol −1 at 25°C and 60°C, respectively (Fig. 6). Contrast this with the endergonic formation of di-and trimeric bismuth-carboxylate complexes, their free energy differences amounting to roughly +35 kJ mol −1 , starting from the monomer at 60°C. Because of the huge computational effort of using the tetramer as catalyst, we Scheme 3 Uncatalyzed model reaction between propyl isocyanate and ethyl alcohol. A urethane species containing also an ester functionality represents the reaction medium. On the right side the 6-membered ring transition-state structure is shown.  3 . This curve initially follows pseudo first-order kinetics, deviating though after a while from the simulated first-order kinetics profile and appearing similar to 0 th -order kinetics. We assume complex kinetics because of this change of kinetic profile. limit our discussion to the monomer. Consider that the concentration of Bi( piv) 3 in solution is low (0.05 mol%) and that entropy (and diffusion) hamper the formation of clusters.
Considering the multiple coordination forms of Bi, we first examined the alcohol-isocyanate coupling in more detail. In general, the coupling reaction is accelerated by facilitating the nucleophilic attack of the alcohol at the carbon atom of the isocyanate thereby lowering the activation barriers. Coordination of alcohol to bismuth lowers the alcohol's pK a , facilitating its deprotonation and enhancing its nucleophilic character. Several pathways of urethane formation have been considered with the insertion reaction either proceeding concertedly with protonated or unprotonated alcohol (coordinated to the catalyst's metal center) or following an outer-sphere type of mechanism (cf. Fig. 7). Clearly, the higher nucleophilicity of coordinated alkoxides in concert with the higher electrophilicity of the isocyanate carbon atom by interaction of the isocyanate nitrogen atom with a Lewis acid is most effective in lowering the reaction barrier. This reaction, by taking place on a metal center, is reminiscent of a 1,2-migratory insertion. Migratory insertion rates are strongly dependent on the solvent or reaction matrix. 70 Consequently, the polarity of ligands (forming the local matrix around the Bi center) affects the rate of 1,2-migratory insertion.
The initial ligand exchange most likely follows an associative mechanism because of the sterically accessible bismuth center in Bi( piv) 3 . There is no significant electronic stabilization of the pivalate's negative charge formed upon initial departure from the bismuth center, which additionally supports this mechanism. Conversely, conjugate base stabilization hints at a dissociative ligand exchange mechanism. Computations additionally confirm that the associative mechanism is dominant. Hence, the ligand exchange rate law is first order in the entering alcohol ligand. The expected large negative entropy of activation is the reason for an additional contribution to the solvent dependence of the hydroxyl-isocyanate reaction. It turns out that Bi( piv) 3 prefers exchanging pivalate for alcohol only once (see Fig. 8 and 9). The concerted 1,2migratory insertion of the alkoxide ligand is identified by quantum chemical computations as predominant reaction mechanism (compare the reaction barriers in Fig. 7).
The efficiency of a catalytic cycle depends on the barrier of the turnover-limiting step but most importantly though on the Fig. 6 Computed Bi( piv) 3 -dimer, -trimer, and -tetramer structures. The highly symmetric Bi( piv) 3 -tetramer structure consists of Bi atoms in tetrahedral a geometry and agrees with the solid-state structure. 42 Fig. 7 Potential modes of reaction between isocyanate and alcohol on the bismuth center. These modes are shown in increasing order of the effective barrier height. Concerted protonation of the forming carbamate can be ruled out. Protonated alcohol also is unlikely to react with isocyanate in a Bi complex. Fig. 8 Energy diagram of the Bi( piv) 3 catalytic cycle. Migratory insertion is turnover limiting. The resting state consists of Bi( piv) 2 (OEt)(PrNCO) (cat-II). The effective equilibria are shifted towards the resting state. Fig. 9 Catalytic cycle starting with Bi( piv) 3 as catalyst precursor. EtO-Bi ( piv) 2 is the active form of the catalyst, its being regenerated by the alcohol component. Highest efficiency is achieved by following this catalytic cycle (heavy arrows).
absence of easily accessible off-cycle species with siphoning character. Favorable equilibria towards off-cycle species retard the catalytic cycle and open the door to the possible formation of side products such as allophanates (Fig. 9). Clearly, the Bi ( piv) 3 catalytic cycle consists of easily accessible off-cycle equilibria. These efficiency-lowering equilibria have been computationally identified (Fig. 8 and 9). Their presence additionally is supported by the overall kinetic profile, obtained by experiments, hinting at complex kinetics (cf. Fig. 5).
The kinetics of Bi( piv) 3 -catalysed urethane formation is mainly determined by which catalyst intermediate is regenerated after the last step of the catalytic cycle, that is the step releasing the product from the catalyst. The intermittent 0 thorder kinetic character hints at a resting state located at Bi ( piv)(OEt)(PrNCO) (cat-II). Up to this point we have not considered the computed energy diagram in Fig. 8. The relative energy levels of Bi( piv) 3 (cat), Bi( piv) 2 (OEt) (cat-I), and cat-II seem to contradict this hypothesis. However, taking into account the uncertainty of quantum chemical computations in condensed phase, the energy level of cat-II might be lower than the one of cat-I. Seeing this, we should not be surprised to find cat-I being regenerated in the catalytic cycle. Computations show, however, that both the alcohol and free pivalic acid ligand can release the urethane product from the catalyst via proton transfer. Whether the alcohol or the carboxylic acid coordinates to the catalyst depends on their instantaneous concentration. Off-cycle equilibria become important in case of the regeneration of cat, retarding the efficiency of the catalytic cycle (Fig. 9). The most efficient catalytical path, having the least probability of off-cycle intermediates, follows the regeneration of cat-I. As the reaction proceeds, urethane increasingly forms and may immobilize the catalyst ligand pivalic acid, as outlined in Fig. 10. This immobilization redirects the catalytic cycle towards from cat being regenerated to the more efficient pathway regenerating only cat-I.
Computations reveal that protonation of carbamate is thermoneutral and, in addition, does not discriminate between carboxylic acid and alcohol as protonating agents in terms of kinetics (cf. Fig. 8). The corresponding barriers amount to 30 and 33 kJ mol −1 , respectively. The alcohol and carboxylic acid concentrations determinenext to steric demandwhich species favorably coordinates prior to protonating the carbamate to release the product. Let us have a closer look at the two scenarios. Firstly, when the alcohol protonates the carbamate (via associative ligand exchange), cat-I is regenerated, leading to the catalytic cycle consisting of less off-cycle equilibria (cf. Fig. 9). Secondly, when free carboxylic acid takes over protonation, catalyst precursor cat is regenerated, which appears to be the non-active form of the catalyst. In this latter scenario, ligand exchange with alcohol must take place, competing with off-cycle equilibria already at the beginning of the catalytic cycle or perhaps before forming the active catalyst species cat-I. Pre-equilibria including alcohol corroborate the experimental finding of the non-integer partial reaction order of the alcohol reagent in the rate law.
The selectivity of this catalyzed isocyanate-hydroxyl reaction depends on the height of the barrier of the turnover-limiting step compared to the ones that lead to side products. Formation of allophanate requires conquering a barrier of as much as 46 kJ mol −1 , compared to 41 kJ mol −1 for the turnover-limiting migratory insertion. The extent of allophanate formation mainly depends on the relative height of the reaction barrier forming allophanates compared to that of carbamate protonation. As outlined in Fig. 8, isocyanate may coordinate to cat-P, competing with the coordination of alcohol (or pivalic acid). Latter coordination is clearly favored by 10-15 kJ mol −1 over that of isocyanate. Moreover, the barrier to form allophanate amounts to as much as 46 kJ mol −1 and, thus, lies higher than the ones leading to product release. These results corroborate the absence of allophanates in experiments.
Because the 1,2-migratory insertion is turnover limiting the resting state of the Bi( piv) 3 catalyzed reaction is cat-II. Consequently, the barrier of the liberation of the product from the catalyst is not turnover limiting. Generally spoken, if computations result in this scenario, the overall kinetics most probably exhibits pseudo 0 th -order. The case in which product release is turnover limiting also leads to pseudo 0 th -order kinetics. But now the resting state consists of the catalyst-product (cat-P) intermediate, raising the probability of side-product formation. Let us assume now that the barrier of ligand exchange with alcohol is higher compared to the one of the insertion step. In this scenario, the reaction will follow pseudo 1 st -order kinetics and the resting state would be the catalyst precursor species (cat).
Often humidity cannot be excluded during the isocyanatehydroxyl reaction, depending on the application. Therefore, water may be present during synthesis, raising the question of the hydrolytic stability of Bi( piv) 3 . Water coordinates favorably to cat and cat-I (Fig. 8). However, the formation of bismuth hydroxides is unlikely because the reaction pathway is consistently uphill in energy. The computations further show that the barrier for the reaction of water with isocyanate is higher than that for the reaction with alcohol. This high barrier explains why Bi( piv) 3 does not promote the isocyanate-water reaction. Apparently, water blocks the catalyst and removes it from the catalytic cycle by siphoning. Unfortunately, this bismuth-aqua complex is no longer available for catalyzing the hydroxyl-isocyanate reaction.
The mechanistic concepts can be applied (in a first approximation) to multi-nuclear bismuth complexes (e.g. to Bi-Li clusters). Preliminary quantum chemical computations reveal that in contrast to mononuclear bismuth catalysis, isocyanate ligand association is thermoneutral or even exergonic when lithium is involved in the catalyst complex. The substrate iso- cyanate then preferably coordinates to lithium, whereas the deprotonated alcohol substrate is bridging between lithium and bismuth. Presumably, the higher Lewis acidity of lithium compared to that of the bismuth species leads to a higher protonating potential of alcohols and free carboxylic acids, resulting in a more efficient release of product from the dinuclear complex.

Conclusions
In conclusion, a series of Bi carboxylates 1-3 and heterobimetallic complexes 7-9 have been synthesized and characterized. The isostructural clusters 7-9 consist of two bismuth carboxylates and four lithium carboxylates and are spontaneously formed upon mixing bismuth and lithium carboxylates in solution. 1 H DOSY NMR experiments and variable temperature NMR studies reveal that solutions of 1-3 and 7-9 contain a dynamic mixture of aggregates with an average size comparable to that of their solid-state structures.
Studies on the catalytic performance of the Bi compounds in the urethane forming reaction show that the Bi-catalyzed reaction can be significantly accelerated by using heterobimetallic mixtures of Bi and Li carboxylates, with 8 almost reaching the activity of DBTL without promoting the isocyanate-water reaction. Theoretical studies show that urethane formation utilizing Bi(piv) 3 follows complex kinetics that exhibit intermittent pseudo 0 th -order kinetics. The turnover-limiting step in the catalytic cycle is the migratory insertion of the isocyanate into the Bi-O bond of the alkoxy ligand, eventually forming the carbamate. The urethane product is released from the catalyst by proton transfer from either the alcohol substrate or the free pivalic acid ligand. Collectively, heterobimetallic complexes composed of bismuth and lithium carboxylates are a promising alternative to the toxic organotin-containing catalysts in PU synthesis. Successive studies on the synthesis of even more active bismuth catalysts are under investigation in our laboratories.

Conflicts of interest
There are no conflicts to declare.