Bhaskar
Mondal
a,
Robin D.
Wilkes
b,
Jonathan M.
Percy
*a,
Tell
Tuttle
*a,
Richard J. G.
Black
b and
Christopher
North
b
aWestCHEM, Department of Pure and Applied Chemistry, University of Strathclyde, Glasgow, G1 1XL, UK. E-mail: tell.tuttle@strath.ac.uk; jonathan.percy@strath.ac.uk; Tel: +44 (0)141 548 2290
bPhosphonicS Ltd, 44c Western Avenue Milton Park, Abingdon, Oxfordshire OX14 4RU, UK
First published on 8th October 2013
Dispersion corrected density functional theory (DFT-D) has been applied to understand the performance of several palladium metal scavengers. Nine different sulfur-based ligands and three different palladium metal sets have been investigated in detail. Based on a thorough analysis of the thermodynamic binding parameters ΔH, ΔG and ΔS, we have identified the best binding modes for all scavenger ligands. Bis-monodentate coordination is favoured over chelation in ΔH and ΔG values for most of the scavenger ligands. Special attention has been paid to the ligand strain energies, which account for the structural changes of the ligands upon complexation indicating that small (5-membered) chelates are considerably less favourable than expected. Some ligands can use their longest chain (>7-atoms) to yield trans chelates, which ligands with shorter chains (≤6-atoms) are unable to form. A secondary amino nitrogen (RR′NH) is found to be the best donor with highest binding enthalpy for Pd(II) metal systems. In terms of the strength of the initial binding interactions, –SMe > –SH; capping thiols (–SH) as thioethers (–SMe) is therefore suggested to be an effective strategy in scavenger design. These observations mark the beginning of a knowledge base of the full range of possible interactions, leading to the construction of a sulfur ligand database for the design of scavenger systems.
A number of processes have been developed for the removal of metal impurities from pharmaceutical intermediates and products.4,5 Standard methods of metal removal, such as distillation, use of ion exchange resins, precipitation, extraction, and crystallisation usually have significant disadvantages and residues can persist even beyond rigorous purification steps of this type.6 Functionalised materials which can scavenge residual metals from APIs or waste streams represent a strategic solution to these problems. One of the key advantages of using these materials is that the functionality can be designed to have a very high affinity towards specific metal residues. The molecular backbones of APIs often bear a number of functional groups and stereogenic centres situated in close proximity; these functional arrays are created to ensure binding to molecular targets. However, strong binding of precious metal atoms can be an unwanted outcome; the design of an effective metal scavenger therefore requires the incorporation of combinations of ligands or functional groups that have an overall higher affinity for the metal than that of the API.
A wide range of multifunctional and complex metal scavengers are required when the APIs and related intermediates are structurally diverse. PhosphonicS has developed and applied a range of successful multifunctionalised materials based on mesoporous silica, which possess siloxane (Si–O–Si) bonds linking organic functionality to the silica framework, for the removal of palladium residues from product streams;7 some of these have been deployed at scale.8 Materials based on polystyrene matrices are also available commercially. Sulfur-based organic functional groups which offer potential strong metal–ligand interactions, and which can be exploited in metal sequestration are privileged in the design of these materials.
The final product of a palladium metal-catalysed reaction may be contaminated by impurities in which the metal is present in different oxidation states.9 For our study, Pd(0) and Pd(II) have been selected as the initial metal centres for investigation because of their ubiquity in catalysis, and the availability of detailed knowledge about their catalytic mechanisms.10
Scavenger materials are represented by generic structures 1 and 2. Linkers may be purely hydrocarbon, or contain a thioether linkage, depending on the method of synthesis.11 Species 3 (Silicycle SiliaMetS Thiol), 4 (PhosphonicS SPM32) and 5 (Biotage MP-TMT) represent the types of materials available commercially and in current use.
In this contribution, we have studied detailed thermodynamic aspects of the basic chemistry exploited by scavenger materials during the scavenging process. This approach provides a quantitative understanding of the performance of several scavengers containing sulfur-based functionality based on electronic structure calculations of binding parameters for metal–scavenger binding. We have used three different metal sets that represent palladium species ubiquitous in coupling reaction chemistry. We have replaced the linkage to the silica matrix (Fig. 1) for each scavenger material with a methyl group to simplify the calculations. This represents an approximation but it is unlikely that a longer tether will significantly affect the Lewis basicity of the sulfide. A key goal of the study is the exposure of the effectiveness of the internal thioether ligand, a key design feature of the materials.
The binding between the metal sets and the ligands was investigated using the reactions shown in Scheme 1. Reaction (1) is used for the calculation of binding parameters (ΔH, ΔG and ΔS) for bis-monodentate coordination of the functional groups to the metal centres, and reaction (2) is used for chelation between the metal centres and bidentate ligands.
![]() | ||
Scheme 1 Model metal–ligand binding reactions; L is the arbitrary ligand associated with the metal sets, L1 and L2 represent donor sites of a multidentate ligand. |
The aim of this work is to model a range of established palladium scavengers which contain a thioether sulfur as part of a linkage to a solid support, to study the conformational preferences of the ligand arrays they present, and to assess the importance of the thioether moiety as a participating centre in a chelate array. Based on these observations, important initial principles and findings are presented which begin the construction of a sulfur ligand knowledge base for the design of the next generation of metal scavenger systems.12
In the absence of fortuitous error cancellation, the larger the basis set, the more accurate it is. In this case, the largest basis set BS3 can therefore be taken as the benchmark among the three. However, the computational time associated with BS3 will be significantly higher than BS2 and BS1; the use of BS3 is therefore not practical for the investigation of the full range of ligands and binding motifs described in Fig. 4. We therefore investigated the use of the smaller basis sets for this task. The methods were compared based on consistency of structural parameters with the X-ray crystal structure, consistency in relative enthalpies with the highest level of theory tested, and computational cost (CPU time). In this study the benchmark level of theory is the B2PLYP functional in conjunction with the BS3 basis set.
The polarising effect of solvent on the binding energy of the systems is expected to be minimal given the neutral nature of the ligands. Nonetheless, we have carried out exploratory calculations to confirm this for ligands 10 and 12, which indicate the effect to be a slight reduction in the binding energy by ∼2 kcal mol−1. The implicit solvent effect of methanol (selected because it is the simplest alcohol solvent, and alcohol solvents are often present in reactions and streams for scavenging) has been calculated using the polarizable continuum model (PCM) as implemented in Gaussian 09.22 Some aspects of the electronic structures have been analysed using Natural Bond Orbital (NBO) calculations at the same level of theory using NBO 3.1 program23 which is implemented in Gaussian 09. All computations are performed with the Gaussian 09 program.24
The results are presented fully in Tables S1 and S2 respectively in the ESI.† Based on the binding enthalpy/energy values obtained, the B97-D functional in conjunction with basis BS2 provides a reliable level of theory for current systems at an acceptable balance between accuracy and computational cost.
There were no significant structural differences between the complex geometries obtained with the three different functionals.
The accuracy of the structures generated was confirmed with the help of known crystal structures (complex S1, S2 and S3 in the ESI†); we found that both M06-L and B97-D functionals in conjunction with BS2 could reproduce bond lengths to within 0.05 Å of values from crystallographic data. The small deviations may be due to the fact that the model approach is not capable of accounting for crystal-packing effects in this case. Because the B97-D functional in conjunction with BS2 was able to reproduce binding enthalpies (ΔH) relative to the benchmark data, and structural features relative to the crystallographic data accurately, it was selected for further investigation of a wider range of complexes. A complete discussion of the method evaluation is presented in the ESI.†
The adoption of the tetrahedral geometry for the complexes of 8 is due to the metal's occupancy of its zero oxidation state, instead of the +2 oxidation state in 6 and 7. These two different geometrical arrangements can be explained with standard ligand field theory as shown in Fig. 3. The Pd(II) system has d8 electronic configuration which favours a square planar geometry with an empty d orbital. This empty d orbital (conventionally designated the dx2−y2) offers a path (along x and y axes) for each of four ligands to approach the metal and avoid electron–electron repulsions. In contrast, Pd(0) prefers to adopt a tetrahedral geometry to limit electron–electron repulsion between the ligand and the 10 d electrons of the Pd as far as possible. The tetrahedral geometry has the smallest difference between the highest and lowest energy orbitals; with 10 electrons, all orbitals are filled and the tetrahedral geometry offers the lowest energy arrangement possible. Tetrakis(triphenylphosphino)-palladium(0) is one of the best examples in the literature of a tetrahedrally coordinated Pd(0) metal centre.25
![]() | ||
Fig. 3 Electronic distribution for Pd(0)/tetrahedral complexes (a) and Pd(II)/square planar complexes (b) formed from 8 and 6 respectively, with ligand 9. |
Monodentate ligand | H rel |
---|---|
9 | −9.2 |
10a | −8.7 |
11a | −8.3 |
12a | −7.2 |
13a | −2.0 |
14a | −2.2 |
15a | −2.1 |
16a | −7.9 |
Table 1 shows that all of the ligands prefer to form a trans complex with 6. The preference for the trans arrangement is much greater for ligands 9, 10a to 12a, and 16a (about 7–9 kcal mol−1), but still significant for 13a to 15a (about 2 kcal mol−1). The values presented in Table 1 effectively separate all species into two groups, which could reflect the presence of the amide moiety in ligands 13a to 15a. In the cis arrangements, the ligands that contain an amide moiety are able to form two hydrogen bonds between the amide N–H and the Cl ligands attached to the metal centre. In contrast, occupancy of the trans orientation allows only one of these hydrogen bonds to form and is therefore destabilised relative to the cis complexes.
Complexes | ΔH | ΔG | TΔS |
---|---|---|---|
6a/9 | −73.9 | −50.8 | −23.1 |
Type 10 | |||
6/10a | −82.7 | −59.4 | −23.3 |
6/10b | −71.7 | −51.5 | −20.3 |
6/10c | −60.1 | −46.1 | −13.9 |
6/10d | −65.8 | −51.5 | −14.3 |
Type 12 | |||
6/12a | −82.6 | −56.8 | −25.8 |
6/12b | −74.6 | −50.9 | −23.7 |
6/12c | −60.2 | −44.4 | −15.8 |
Type 15 | |||
6/15a | −78.7 | −54.1 | −24.6 |
6/15b | −44.9 | −32.9 | −12.0 |
6/15c | −53.8 | −39.3 | −14.5 |
6/15d-cis | −72.2 | −54.9 | −17.3 |
6/15d-trans | −77.8 | −60.1 | −17.7 |
6/15e-cis | −67.9 | −49.9 | −18.1 |
6/15e-trans | −87.8 | −69.5 | −18.2 |
Type 17 | |||
6/17a | −40.8 | −28.1 | −12.6 |
6/17b | −36.6 | −23.9 | −12.6 |
We have computed the binding parameters (enthalpies (ΔH), free energies (ΔG) and entropies (ΔS)) for all of the ligands in their different binding modes with metal system 6. The mono coordinated complexes were all trans, while cis and trans chelated complexes were investigated where accessible.
A number of interesting features were revealed by the study of the first ligand type 10; trans-monodentate coordination (10a, 10b) via thiol (–SH) or thioether (–SMe) is enthalpically favoured over chelation. Thioether coordination is over 10 kcal mol−1 stronger than thiol coordination enthalpically, a significant and unexpected difference; the same behaviour was observed for ligand type 12. The 10c/10d comparison shows that the advantage of thioethers is also realised within chelates. This finding supports the design of commercial scavenger materials which privilege the internal thioether linkage.
Both ligand types form 5-membered chelates; the most favourable chelate (with 10d) is still less stable than the trans-10a complex by 8 kcal mol−1. It would appear that 5-membered chelates are not favourable for these simple bis-thiol or bis-sulfide ligands. The correct evaluation of the free energy of ligand binding requires that entropy changes are accounted for correctly. Taking the average value of −21.8 kcal mol−1 for TΔS for the formation of trans-10a and trans-10b yields a value of 73 eu for ΔS (at 298 K). This corresponds to a ΔS of ca. 37 eu per ligand, very close to typical values for bimolecular reactions.26 The corresponding value of TΔS for chelate formation is 14 kcal mol−1, or ΔS ≈ 47 eu. The difference of 10 eu represents the cost of immobilising rotors in the ligand following initial monocoordination; as the cost of immobilising one σ-bond rotor is ca. 4 eu, the quantitative treatment of binding entropy would appear to be accurate, because three rotors must be frozen when a 5-membered chelate forms (Scheme 2).
The difference in binding entropies between chelates from 10c/10d, and 12c, is consistent with the presence of one additional C–C bond around which free rotation can occur in 12 (at 298 K, a difference of ∼1 kcal mol−1 corresponds to ca. 4 eu associated with the internal rotation).
While monocoordination appears strongly competitive in these model studies, the local concentration of thiol on the scavenger materials is likely to be quite low because they are based on rigid matrices which are loaded relatively lightly with ligands, and this relatively high local dilution of ligands will erode any advantage of monocoordination significantly. Materials which bear chelating ligands which offer multiple coordination sites avoid this problem of very low local ligand concentration.27
The binding enthalpy (ΔH) order runs parallel to the binding free energy (ΔG) as displayed in Fig. S7 (ESI†) where the variation of binding parameters (ΔH and ΔG) are plotted for all metal–ligand complexes. Subsequent discussion is therefore restricted to the relative enthalpies in most cases. The effects of the presence of different types of chalcogenide donor centres were tested with ligand type 13 (see Table S10 in the ESI†). No particularly favourable modes of binding were discovered but sulfur centres displayed a distinct advantage. For example, a thioether/thiocarbonyl donor pair (13c) was 17 kcal mol−1 more effective than the thioether/carbonyl sites in 13b. Increasing the chain length from n = 1 (13c) to n = 2 (13d) had an effect of about 2 kcal mol−1 on chelation.
Ligand 15 presents a rich array of potential complexes of different types and geometries. trans Bis-monocoordination of the thioether (15a) was significantly more favourable than the formation of the 5-membered chelate from either 15b or 15c. The latter, which features amide NH as the second donor centre, was unexpectedly 9 kcal mol−1 more stable than the complex between 6 and 15b which involves a CO⋯Pd interaction. Ligand 15d forms an 8-membered chelate which can be cis- or trans configured; to our surprise, the trans-chelate was 6 kcal mol−1 more favourable than the cis-stereoisomer and only 1 kcal mol−1 less favourable than trans bis-monocoordination of the thioether in 15a, so for the first time, a chelated mode is directly competitive with monodentate binding. The formation of the trans 11-membered chelate in which ligand 15e was complexed to 6 was very strongly favoured, by 20 kcal mol−1 over the cis-stereoisomer, and 9 kcal mol−1 over trans bis-monocoordination. The more flexible longer chain ligand therefore shows for the first time, a decisive and unambiguous advantage over the binding of two separate ligands. Despite the relatively high cost of the loss of internal rotational modes, the advantage over trans bis-monocoordination is preserved in a free energy difference of 15 kcal mol−1 over 6/15a. Fig. 5 displays the chelate complexes of 14d and 15e with 6, along with the relative enthalpies for cis and trans isomers in kcal mol−1; the behaviour of 14d also represents (at lower conformational complexity) the complexes available from 15d in cis- and trans-8-membered chelates.
![]() | ||
Fig. 5 cis and trans Chelates of 6 with 14 and 15: (a) cis-6/14d; (b) trans-6/14d; (c) cis6/15e; (d) trans-6/15e. The relative enthalpies (Hrel in kcal mol−1) are shown in bold. |
A similar set of calculations was carried out for metal set 7; Table 3 presents the results. A complete table (Table S11†) containing the binding parameters for 7 with the full set of ligands is presented in the ESI.† Species 7 represents a catalytically active monocationic palladium complex obtained when an acetate ligand dissociates from palladium acetate. The binding enthalpies and free energies shown in Table 3 reflect the considerably increased Lewis acidity of the palladium centre in 7 relative to that in 6. The values obtained for TΔS are very similar for 6 and 7, which is to be expected as the entropy of binding is a function of molecularity and ligand flexibility, and not the metal centre. As discussed previously, cis binding is forced in the complexes of 7 because of the nature of the acetoxy ligand so all the chelates represented in Table 3 are the cis-stereoisomers. Nevertheless, the 8-membered cis-chelate from 15d is competitive with monocoordination, and the 11-membered cis-chelate from 15e is strongly favoured over monocoordination, showing the presence of the internal thioether linkage to advantage.
Complexes | ΔH | ΔG | TΔS |
---|---|---|---|
Type 10 | |||
7/10a | −121.9 | −98.9 | −22.9 |
7/10b | −110.5 | −88.7 | −21.8 |
7/10c | −102.7 | −88.4 | −14.3 |
Type 12 | |||
7/12a | −121.7 | −97.3 | −24.3 |
7/12b | −111.7 | −87.8 | −23.8 |
7/12c | −104.3 | −88.9 | −15.3 |
Type 15 | |||
7/15a | −127.3 | −103.0 | −24.3 |
7/15b | −97.3 | −84.2 | −13.2 |
7/15c | −94.8 | −80.3 | −14.5 |
7/15d | −121.3 | −103.9 | −17.5 |
7/15e | −128.5 | −110.9 | −17.6 |
Type 17 | |||
7/17a | −85.5 | −72.1 | −13.4 |
7/17b | −81.8 | −68.6 | −13.2 |
Complexes | ΔELstrain | ΔEcomplex |
---|---|---|
6/15b | 11.5 | −46.3 |
6/15c | 11.3 | −55.5 |
6/15d-cis | 6.3 | −75.2 |
6/15d-trans | 5.6 | −80.3 |
6/15e-cis | 8.8 | −70.7 |
6/15e-trans | 4.6 | −90.5 |
From Table 4, it is clear that the longest chain species 15e favours trans complex formation with the lowest contribution from the strain energy of the ligand (4.6 kcal mol−1). This effect is also observed in the 8-membered chelate complexes with 15d, although the difference between the strain energy in the cis and trans configurations is reduced. The highest strain energies for 15 results when the 5-membered ring chelates of 6/15b and 6/15c are formed (∼11 kcal mol−1, Table 4).
Table S12† shows that strain energy from the ligand can contribute 2–12 kcal mol−1 to the complexation energy depending upon the type of binding. With a constant chain length (S/O or S/N coordination), the strain energies are about 10 kcal mol−1 for 13 to 15 irrespective of binding modes. For 16, S–S chelation incurs almost double the strain energy which arises from S–N chelation. Strain energies for ligands 10 to 12 are in the range of 2–5 kcal mol−1 with 3 kcal mol−1 variation due to the difference in chain length between 10 and 12.
Complexes | ΔH | ΔG | TΔS |
---|---|---|---|
Type 10 | |||
8/10a | −28.4 | −7.1 | −21.2 |
8/10b | −22.0 | −3.6 | −18.4 |
8/10c | −19.0 | −7.1 | −11.9 |
8/10d | −21.2 | −8.8 | −12.3 |
Type 12 | |||
8/12a | −26.0 | −2.4 | −23.5 |
8/12b | −21.1 | 0.2 | −21.3 |
Type 15 | |||
8/15a | −26.5 | −3.3 | −23.2 |
8/15e | −26.2 | −10.5 | −15.7 |
Type 17 | |||
8/17 | −20.4 | −9.6 | −10.8 |
For complexes 8/15, monocoordination of two molecules through –SMe has the highest binding enthalpy (−26.5 kcal mol−1), which is comparable to the strength of the 11-membered chelate from 15e (−26.2 kcal mol−1). Structures for the alternative metal–ligand chelates S/O, S/N and S/Nmid were investigated but none of the attempts yielded a tetrahedral Pd(0) chelate complex. At every attempt, the tetrahedral chelate complex minimised to a tricoordinate complex in which the ligand is bound to the metal centre via thioether-S, indicating a more favourable interaction for the thioether-S binding relative to the alternative coordination through either the carbonyl O or secondary amine. Similar observations were also made with tetrahedral chelation of 12c through a 6-membered ring. For 10a, the monocoordination through thioether sulfur is about 7 kcal mol−1 stronger than chelation through two thioether sulfurs in 10d and the monocoordination through thiol in 10b is about 3 kcal mol−1 better than chelation through thiol and thioether sulfur in 10c. For both 10 and 12, thioether sulfur coordination to the metal centre is more favourable than thiol coordination. This observation is consistent with the trends observed for square planar Pd(II) metal complexes.
Both attempts to find a tetrahedral complex between 8 and 17via chelation involving nitrogen and sulfur (17a or 17b) resulted in a tricoordinate complex, in which the ligand is coordinated through the nitrogen alone (17c). Fig. S8 in the ESI† shows two different chelated starting geometries that minimise to a structure in which only N-coordination is observed.
All the binding enthalpies for the Pd(0) complexes from 8 are lower than those for Pd(II) species 6 and 7 suggesting very strongly that palladium removal with scavenger materials should be considerably easier when the material is present in the higher oxidation state.28
A natural bond orbital (NBO) analysis23 has been performed on the 8/17c complex to explore the extent of coordination of the different atom centres to the metal and to rationalise the observation of the N-coordinated complex alone.
Second order perturbation energies from NBO analysis measure the strengths of the interactions between two different moieties. This analysis shows that the dominant interaction is between N(6) of the ligand and Pd(16) centre. In addition, two of the S-atoms (S(7) and S(8)) from the ligand are also interacting with the metal centre which further stabilises the tricoordinate complex (Fig. 6).
![]() | ||
Fig. 6 Selected second order perturbation energies E(2) (kcal mol−1) and selected interatomic distances (Å) from NBO analysis for 8/17c complex. Atom numbers are shown in red. |
Interacting lone pair orbitals of S(7) and S(8) from the ligand are presented together in Fig. 7. These lone pairs are predominantly (67%) 3s in character, with a small contribution from the hybrid sp2-orbital which forms the C–S bond.
Complexes | ΔH | ΔG | TΔS |
---|---|---|---|
Type 15 using thioether-S | |||
6/15a | −78.7 | −54.1 | −24.6 |
6/15b | −44.9 | −32.9 | −12.0 |
6/15c | −53.8 | −39.3 | −14.5 |
6/15d-cis | −72.2 | −54.9 | −17.3 |
6/15d-trans | −77.8 | −60.1 | −17.7 |
6/15e-cis | −68.0 | −49.9 | −18.1 |
6/15e-trans | −87.8 | −69.5 | −18.2 |
Type 15 not using thioether-S | |||
6/15f | −66.4 | −40.8 | −25.6 |
6/15g | — | — | — |
6/15h | −95.1 | −70.3 | −24.8 |
6/15i | −69.4 | −37.2 | −32.2 |
6/15j | −76.4 | −59.1 | −17.3 |
6/15k | −103.4 | −76.6 | −26.8 |
Of the six different possible ways of binding to the metal centre which do not involve thioether-S which we explored, monocoordination through the secondary amino group in 15k results in the highest binding enthalpy at −103.4 kcal mol−1. Binding through the primary amino group in 15h results in a comparable binding enthalpy at −95.1 kcal mol−1. These two nitrogen sites offer the strongest binding of all the ligand atoms investigated. The most stable trans chelate from 15e lies 16 kcal mol−1 above mono coordinated complex 15k.
We could not locate the 8-membered chelate complex which involves the amide nitrogen and the terminal primary amino group (15g). Every attempt led to the formation of the most stable trans chelate 15e, which is presented in Fig. S9 in the ESI.† However, the inability to form the 8-membered chelate is a reflection of the preference for these systems to form large (≥10) ring chelates (vide supra) rather than the binding ability of the amino groups.
Finally, our treatment of potentially acidic binders (RSH, RNH2, RR′NH) requires comment. We have treated these species exclusively as neutral ligands; we have not evaluated the strengths of the interactions of their conjugate bases (RS−, RNH−, RR′N−) with metal centres. Mechanistic models of thioether formation proposed by Hartwig29 and explored experimentally by Jutand30 involve initial coordination of the neutral thiol to the metal, followed by deprotonation; in the absence of amine bases, the deprotonation is slow. Jutand and co-workers could not rule out a classical mechanism involving proton transfer between thiol and amine to afford a thiolate/ammonium ion pair entirely, but presented strong evidence for metalation/deprotonation. The strength of the initial binding of the metal to scavenger materials is therefore evaluated most accurately by considering neutral species, rather than their conjugate bases.
Thioether (–SMe) ligand sites were generally more effective than the corresponding thiols while a secondary amino nitrogen (RR′NH) is the best donor type, with the highest binding enthalpy for Pd(II) metal systems (6 and 7). While bis monocoordination appears strongly competitive in these model studies, the local concentration of any supported ligand on the scavenger materials is likely to be quite low, and this relatively high local dilution of ligands will erode any advantage of bis monocoordination significantly. The design of materials which bear chelating ligands based on multiple coordination sites therefore represents an astute and effective solution to this problem, with the internal thioether providing a significant contribution to a strong interaction.
Footnote |
† Electronic supplementary information (ESI) available: Details of computational method evaluation, binding parameters (ΔH, ΔG and TΔS) for the full set of ligands, decomposition of complexation energies, structural changes during optimisation, optimised Cartesian coordinates at B97-D/BS2 level of theory. See DOI: 10.1039/c3dt52282b |
This journal is © The Royal Society of Chemistry 2014 |