Literature DB >> 30556625

Isolation of the Metalated Ylides [Ph3 P-C-CN]M (M=Li, Na, K): Influence of the Metal Ion on the Structure and Bonding Situation.

Christopher Schwarz1, Lennart T Scharf1, Thorsten Scherpf1, Julia Weismann2, Viktoria H Gessner1.   

Abstract

The isolation and structural characterization of the cyanido-substituted metalated ylides [Ph3 P-C-CN]M (1-M; M=Li, Na, K) are reported with lithium, sodium, and potassium as metal cations. In the solid-state, most different aggregates could be determined depending on the metal and additional Lewis bases. The crown-ether complexes of sodium (1-Na) and potassium (1-K) exhibited different structures, with sodium preferring coordination to the nitrogen end, whereas potassium binds in an unusual η2 -coordination mode to the two central carbon atoms. The formation of the yldiide was accompanied by structural changes leading to shorter C-C and longer C-N bonds. This could be attributed to the delocalization of the free electron pairs at the carbon atom into the antibonding orbitals of the CN moiety, which was confirmed by IR spectroscopy and computational studies. Detailed density functional theory calculations show that the changes in the structure and the bonding situation were most pronounced in the lithium compounds due to the higher covalency.
© 2019 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.

Entities:  

Keywords:  alkali metals; bond theory; lithium; solid-state structures; structure elucidation; ylides

Year:  2019        PMID: 30556625      PMCID: PMC6519153          DOI: 10.1002/chem.201805421

Source DB:  PubMed          Journal:  Chemistry        ISSN: 0947-6539            Impact factor:   5.236


Introduction

Bisylidic compounds such as carbodiphosphoranes (CDP) and carbodicarbenes (CDC), have gained intense research interest over the past years.1 Due to their unique electronic structures, these compounds have been used as strong donor ligands in different research areas, such as for the coordination and stabilization of small and/or reactive molecules2 or in the preparation of catalytically active transition‐metal complexes.3 The electronic structure itself has been the subject of most controversial discussions.4 Although initial debates on contributions of the double bonded ylene (R3P=CR2) and the ylidic ionic structures have ended in favor of the ylide description (R3P+−C−R2), recent computational studies argued for a predominant contribution of a dative‐bonding situation (R3P→CR2).5 This captodative bonding mode led to the conclusion that the carbon atom in ylidic compounds can serve as a “central atom” similar to a metal in a complex, which also led to the denotation of bisylides and related compounds as carbones with the carbon atom in the zero oxidation state (Figure 1).6
Figure 1

The metallic nature of carbon in ylidic compounds.

The metallic nature of carbon in ylidic compounds. Our group has reported on the anionic congeners of bisylidic compounds, that is, α‐metalated ylides (yldiides).7 These compounds also exhibit two lone pairs of electrons at the central ylidic carbon atom and are thus strong donor ligands, which, for example, allowed the stabilization of a series of borenium ions or the synthesis of highly electron‐rich phosphine ligands.8, 9 Despite their huge synthetic potential, very little is known about their electronic structures. Only four yldiides have been isolated and structurally authenticated so far.8, 10 In this regard, yldiide 1 seemed to be especially interesting because it can be viewed as an anionic carbone, in which the carbon atom is coordinated like a metal in a cyanido‐phosphine complex (Figure 1). Thus, comparable to bisylides, the yldiide may exhibit different bonding situations, which also may be influenced by the counter cation (typically M+ with M=Li, Na, or K). Therefore, cation coordination in 1‐M can either occur at the cyanido nitrogen or at the α‐carbon atom. These different coordination modes have been shown to influence the reactivity of various metalated nitriles11, 12 and are expected to affect the bonding situation. Comparable to CDPs, the electronic structure in yldiides can be described by different bonding situations. In the past, the ylidic form 1A (Figure 2) was used to account for the two lone pairs at the central carbon atom and the high nucleophilicity, which have both been experimentally proven by the use in double‐functionalization and cascade reactions or in the formation of digold complexes.13, 14 However, the dative fragmentation pattern 1B has also been proposed in analogy to transition‐metal complexes.15 Previous DFT studies have already revealed that often different (also unsymmetrical) structures must be considered to fully describe the electronic distribution in metalated ylides.16, 17 Yet, no influence of the metal ion has been considered to date. For example, yldiide 1 can potentially be described by either of the six Lewis structures shown in Figure 2, in which additionally the metal cation can coordinate to either the nitrogen or the carbon atom (hence, 12 structures in total). Thus, to obtain initial insights into the structure and bonding situation of 1, we set out to isolate and structurally characterize different metal complexes of the yldiide. We were particularly interested in answering the following questions: i) Which kind of structures are formed with different metal cations? ii) Does the bonding situation change with different metals and the different coordination modes of the metal (N‐ vs. Ccoordination)? iii) Do changes in the bonding situation correlate with experimental data (bond lengths, C−N stretching frequencies)?
Figure 2

Possible Lewis structures for yldiide 1.

Possible Lewis structures for yldiide 1.

Results and Discussion

Synthesis and NMR and IR spectroscopic studies

As a starting point of our investigations, we first addressed the isolation of yldiide 1 complexed to different alkali metal cations. The metalated ylide 1‐Na was already reported by Bestmann and co‐workers and was used in situ in a variety of organic transformations.13 The yldiide is readily accessible by a stepwise deprotonation procedure. At first, ylide 1‐H is obtained in high yields by deprotonation of the phosphonium salt 1‐H with one equiv of sodium hydride or n‐butyllithium in THF at room temperature (Scheme 1). The use of the lithium base always resulted in the formation of a lithium bromide adduct of 1‐H, in which the lithium is coordinated by the nitrile moiety (see the Supporting Information for XRD analyses of 1‐H and 1‐H), whereas in contrast, the sodium base gave way to the “free” ylide 1‐H.
Scheme 1

Preparation of the metalated ylides 1‐M (M=Li, Na, K) from phosphonium salt 1‐H.

Preparation of the metalated ylides 1‐M (M=Li, Na, K) from phosphonium salt 1‐H. Deprotonation of the ylide 1‐H was performed using alkali metal hexamethyldisilazide (HMDS) bases. Although complete conversion from 1‐H to 1‐M (with M=Li, Na or K) was accomplished with all HMDS bases according to the NMR spectroscopic data (directly obtained after mixing the reagents), only the potassium salt could be isolated in a pure form. For the lithium and sodium yldiides, incorporation of LiHMDS and NaHMDS, respectively, into the structures of 1‐M (see below) excluded their complete removal from the yldiide, and thus complicated the isolation of the pure compounds. This is particularly true for 1‐Li⋅LiHMDS, which could be isolated only in very low yield (9 %). The isolation of 1‐Na⋅NaHMDS was found to be slightly easier, but it gave strongly varying yields with the best result providing 77 % of the product. In contrast, pure 1‐K was obtained as a yellow solid in 84 % yield by simple washing of the formed solid with diethyl ether. Overall, the metalated ylides are highly reactive compounds independent of the metal cations or the use of further coordinating ligands. They tend to get easily protonated to 1‐H, particularly in THF and other ethereal solvents, which also complicated spectroscopic studies in solution.18 NMR spectroscopic studies were performed on the isolated metaled ylides, both with and without the addition of crown ethers. Only the crown‐ether complex of 1‐Li was examined after in situ preparation because its isolation proved to be particularly difficult. The successful formation of the yldiides was confirmed by a distinct high‐field shift of the signal in the 31P{1H} NMR spectrum from δ=23.2 ppm for 1‐H to δ=−1.9 to −19.9 ppm for 1‐M (see Table 1). Such a high‐field shift has also previously been observed for the formation of other yldiides.7 Most interestingly, the extent of the high‐field shift was found to be dependent on the nature of the metal cation and on the presence of additional donor bases, particularly crown ethers, thus suggesting an influence of the metal on the electronics of the yldiide. This feature can be impressively seen through the more pronounced high‐field shift in the series Licrown‐ether addition to 1‐K [e.g., δ(31P)=−10.5 ppm for 1‐K and −19.9 ppm for 1‐K⋅(18‐crown‐6)].
Table 1

NMR and IR spectroscopic properties of the metalated ylides 1‐M and their corresponding protonated congeners [1‐H and 1‐H (12c4=12‐crown‐4; 15c5=15‐crown‐5; 18c6=18‐crown‐6).

δ P [ppm] δ C1 [ppm] 1 J PC [Hz] δ C2 [ppm] 2 J PC [Hz] ν˜ (CN)(THF) [cm−1]
[1‐H2]Br 21.917.8; 54.8112.2; 9.4
1‐H 23.2−2.9; 135124.7; 7.52157
1‐Li⋅(LiHMDS)−1.9−3.4 (br)134.4; 20.11989
1‐Li⋅(12c4)−3.1−2.6; 131.8134.6; 12.61995
1‐Na⋅(NaHMDS)−5.2−5.4; 68.2142.5; 12.62008
1‐Na⋅(15c5)−10.9−0.7; 72.7141.1; 15.62023
1‐K −10.50.83; 68.0140.9; 14.42001
1‐K⋅(18c6)−19.95.4; 59.2137.3; 15.62014
NMR and IR spectroscopic properties of the metalated ylides 1‐M and their corresponding protonated congeners [1‐H and 1‐H (12c4=12‐crown‐4; 15c5=15‐crown‐5; 18c6=18‐crown‐6). It must be noted that, in contrast to all other complexes, the NaHMDScontaining sodium salt 1‐Na exhibited only a broad signal in the 31P{1H} NMR spectrum at room temperature. This signal somewhat sharpened upon heating up to 60 °C, and, most interestingly, gave rise to a splitting into four sharp signals at −30 °C. This observation indicates the formation of a complex, in which the yldiides possess four different environments. Exchange on the NMR timescale between these different environments is fast at higher temperatures, whereas it is slow at −30 °C. Four different phosphorus sites were also found in the crystal structure (see below, Figure 3), thus suggesting that the structures are the same both in the solid state as well as in solution.
Figure 3

Molecular structures of [(1‐Na)8(NaHMDS)2] and [(1‐Li)3(LiHMDS)5] in the solid state. Hydrogen atoms are omitted for clarity; displacement parameters drawn at the 30 % (for [(1‐Na)8(NaHMDS)2]) and 50 % probability level (for [(1‐Li)3(LiHMDS)5]). Selected bond lengths and angles are given in Table 2 and structural details in the Supporting Information.

In addition to the changes of the 31P{1H} NMR shift, the decrease in the 1 J PC coupling constant in 1‐M (59–72 Hz) compared to ylide 1‐H (135 Hz) was characteristic of the successful α‐metalation. This decrease can be explained by the higher p‐character of the P−C linkage due to the increase of s‐electron density at the ylidic carbon atom to stabilize the two lone pairs of electrons.19 The 13C NMR shift of the nitrile carbon atom proved to be most indicative to distinguish between N‐ versus Cmetalation in α‐metalated nitriles, because the chemical shift is highly sensitive to the local environment.20 Typically, N‐metalated nitriles show signals between δ(13C) = 140 and 157 ppm, whereas the carbon in the corresponding Cmetalated nitriles resonates at higher field between δ(13C)=115 and 138 ppm.21 Yet, comparison of the 13C NMR shifts of the different metal salts of 1 and their crown‐ether complexes shows no clear evidence for either coordination mode (Table 1). All shifts are approximately in between the two ranges, thus suggesting either a different or both types of coordination in solution. IR spectroscopy of the different yldiides 1‐M in THF solution proved to be highly informative. The C−N stretching frequencies decrease considerably upon each deprotonation step from =2250 cm−1 in 1‐H to 2137 cm−1 in 1‐H and approximately 2000 cm−1 in 1‐M. This clearly indicates a decrease in the C−N bond order and a transition from a C−N triple to a double bond.22 Considering the different possible Lewis structures of 1‐M, this observation can only be explained by structures with a covalent C=C double bond (1D or 1E, Figure 2) or by the captodative bonding mode (1B or 1F) with strong π back‐donation from the central carbon atom into the CN moiety (analogous to an electron‐rich metal in CO complexes). Two further interesting observations can be made from the IR data provided in Table 1 with respect to the impact of metal coordination on the structure of 1‐M: i) The CN stretching frequency is significantly lower for the lithium compounds 1‐Li and 1‐Li⋅(12c4) compared to their sodium and potassium analogues and ii) the crown‐ether complexes exhibit higher frequencies (approx. 10–15 cm−1) than their crown‐ether‐free analogues. These observations correlate with the increase in the electropositive character of the metal and the decrease of covalent character in their bonds to other elements. It seems that the more tightly the metal is bound to the yldiide the weaker is the C−N bond and thus, it is the weakest for 1‐Li (see below). A similar observation has also been reported for metalated and polymetalated nitriles.23

Crystallographic studies

To gain further insights into the structures of 1‐M, crystallization of the different metal salts, including their crown‐ether complexes, was attempted. Although all attempts to crystallize pure 1‐K failed, the NaHMDS/LiHMDScontaining sodium and lithium compounds could be crystallized by diffusion of n‐hexane into a benzene solution of 1‐Na and by storage of an n‐hexane solution of 1‐Li, respectively. The structure of the sodium compound was only of poor quality due to the presence of highly disordered solvent molecules in the structure, impeding a detailed discussion of the structure parameters (see the Supporting Information for details). Nonetheless, it unambiguously confirms the coordination of two NaHMDS to eight molecules of 1‐Na, which are arranged in an unusual, non‐symmetrical fashion. In [(1‐Na)8(NaHMDS)2], the sodium atoms feature contacts to the ylidic carbon atom as well as to the cyanido nitrogen atom, in accordance with the 13C NMR shift of C1 lying between the typical shifts for C‐ and N‐coordination (see above). The same holds true for the analogous lithium compound, which forms an unsymmetrical aggregate of the composition [(1‐Li)3(LiHMDS)5] in the solid state. The molecular structure is shown in Figure 3 and structural details are given in Table 2 and the Supporting Information. The aggregate is formed by two parts consisting of one yldiide 1‐Li and two LiHMDS molecules, which are connected by a further molecule of 1‐Li and LiHMDS. A possible C‐symmetry of the compound is broken solely by the orientation of the yldiide in the central unit (highlighted in blue in Figure 3). In addition to the structures of the lithium and sodium complexes of 1, a series of different LiHMDS and NaHMDS structures with varying geometries and aggregation motifs has been reported in the literature.24
Table 2

Comparison of structural properties of the metalated ylides 1‐M and their corresponding protonated congeners [1‐H and 1‐H; bond lengths are given in Å, angles in °; for 1‐K⋅(18c6), average values of two crystallographically independent molecules in the asymmetric unit are given.

P−CC−CC−NP‐C‐CC‐C‐N
[1‐H2]Br1.811(2)1.453(2)1.136(2)112.1(1)179.1(2)
1‐H 1.693(2)1.388(3)1.153(3)123.0(2)178.4(3)
1‐H⋅LiBr1.704(2)1.388(3)1.158(3)120.9(2)178.2(2)
[(1‐Li)3(LiHMDS)5]1.666(2)1.691(2)1.681(2)1.344(3)1.358(3)1.349(3)1.188(2)1.181(3)1.183(3)128.5(1)113.7(1)121.8(2)172.7(2)174.6(2)174.3(2)
[[(1‐Li)4⋅(18c6)] 1.654(3)1.642(3)1.344(3)1.342(4)1.188(3)1.171(3)121.9(2)131.2(3)173.7(3)171.9(3)
[1‐Na⋅(15c5)]1.637(3)1.370(5)1.167(4)125.3(3)174.1(3)
[1‐K⋅(18c6)]1.650(2)1.377(3)1.179(3)120.9(2)173.7(3)
Molecular structures of [(1‐Na)8(NaHMDS)2] and [(1‐Li)3(LiHMDS)5] in the solid state. Hydrogen atoms are omitted for clarity; displacement parameters drawn at the 30 % (for [(1‐Na)8(NaHMDS)2]) and 50 % probability level (for [(1‐Li)3(LiHMDS)5]). Selected bond lengths and angles are given in Table 2 and structural details in the Supporting Information. Comparison of structural properties of the metalated ylides 1‐M and their corresponding protonated congeners [1‐H and 1‐H; bond lengths are given in Å, angles in °; for 1‐K⋅(18c6), average values of two crystallographically independent molecules in the asymmetric unit are given. A comparison of structural parameters of the central PCCN moiety in the series of compounds 1‐H, 1‐H and 1‐Na⋅(HMDS) and [(1‐Li)3(LiHMDS)5] shows that the P−C and the C1−C2 bond shorten after each deprotonation step, the P−C bond goes from 1.811(2) Å to about 1.68 Å and the CC bond from 1.45 to 1.35 Å, thus in line with a transition from a single to a double bond. In contrast, the C−N bond length increases from 1.136(2) Å in 1‐H to about 1.18 Å in the lithium compound[(1‐Li)3(LiHMDS)5], which is consistent with a reduction of the C−N bond order and an C=C=N structure. The P‐C1‐C2 and C1‐C2‐N angles change only slightly compared to the ylide, however, it is noticeable that the CC‐N linkage clearly deviates from linearity. To further analyze the oligomeric structures and to probe whether the different alkali metals exhibit any preference for coordination to the nitrogen‐end or the ylidic carbon atom, crown ethers were added to 1‐M for additional complexation of the metal. Unfortunately, crystallization attempts with the small 12‐crown‐4 repeatedly failed for the lithium compound. However, single crystals were obtained with both 18‐crown‐6 (18c6) and 15‐crown‐5 (15c5) (Figure 4).25 Compound [(1‐Li)3⋅(15c5)]2 forms a unique C‐symmetric dimeric structure consisting of six yldiides and two crown ethers. The central structural motif is formed by a (Li−CC−N)2 eight‐membered ring connected to two Li2N2 four‐membered rings. Due to the large size of 15‐crown‐5, the crown ether serves as a bridging ligand for coordinating the two lithium cations. In the solid state, all yldiides exhibit contacts of the metal to both the ylidic carbon atom and the nitrogen of the nitrile moiety. In contrast to [(1‐Li)3⋅(15c5)]2, the 18‐crown‐6 compound of 1‐Li, [(1‐Li)4⋅(18c6)]∞, exhibits a polymeric structure consisting of two different dimeric subunits that contain a Li2N2 four‐membered ring. In one of these units the lithium atoms are coordinated solely by the nitrile moiety and the crown ether, whereas they also bind to the ylidic carbon atom in the second fragment. Thus, one of the yldiides again features mixed C/N‐coordination of the metal, whereas the other only binds the metal at the nitrogen end (the different coordination environments are highlighted in blue in Figure 4). This results in significant structural differences in these two yldiides. As such, the exclusively N‐coordinated yldiide features a significantly shorter C−N bond (1.171(3) versus 1.188(3) Å) and a wider P‐CC angle (131.2(3) and 121.9(2)°) compared to its C,N‐complexed analogue (see explanation below).
Figure 4

(Top) Molecular structure of [(1‐Li)3⋅(15c5)]2 and (middle) [(1‐Li)4⋅(18c6)]∞ in the solid state. Hydrogen atoms and solvent molecules are omitted for clarity; 50 % (30 % for [(1‐Li)3⋅(15c5)]2 displacement parameters. Selected bond lengths and angles are given in Table 2.25

(Top) Molecular structure of [(1‐Li)3⋅(15c5)]2 and (middle) [(1‐Li)4⋅(18c6)]∞ in the solid state. Hydrogen atoms and solvent molecules are omitted for clarity; 50 % (30 % for [(1‐Li)3⋅(15c5)]2 displacement parameters. Selected bond lengths and angles are given in Table 2.25 The crown ether complexes of 1‐Na and 1‐K were obtained by diffusion of n‐hexane into a solution of 1‐Na in toluene and n‐pentane into a THF solution of 1‐K, respectively. In contrast to 1‐Li, the 15‐crown‐5 and 18‐crown‐6 complexed sodium and potassium compounds, respectively, exhibit monomeric structures (Figure 5). Yet, they differ in the coordination of the yldiide to the respective metal cation. Although the harder sodium atom coordinates solely to the nitrogen atom of the nitrile, the softer potassium prefers the interaction with the ylidic and nitrile carbon atom. To the best of our knowledge, such a η 2‐coordination mode has never been observed for any other α‐metalated nitrile. Until now, only C‐ or N‐metalation have been observed, in which the alkali metals usually prefer N‐coordination.26 The K−C bond lengths vary between 2.899(2) and 3.293(2) Å and are in the range of known organopotassium compounds (bond lengths correspond to both independent molecules in the asymmetric unit).27 Despite the slightly bent CC‐N linkage (173.7(3)°), the K−N distance amounts to 3.749(2) Å and is thus longer than typical K−N bond lengths (2.77–3.04 Å).27c, 28 However, in contrast to transition‐metal olefin complexes, the nitrogen atom in the structure of 1‐K⋅(18c6) points towards the metal center suggesting that nonetheless some weaker K−N interaction might be present in the complex.
Figure 5

(a) Molecular structure of [1‐Na⋅(15c5)] and [1‐K⋅(18c6)] in the solid state. [1‐K⋅(18c6)] crystallizes with two crystallographically independent molecules in the asymmetric unit, only one of which is shown here. Hydrogen atoms and THF solvent molecules are omitted for clarity; 50 % displacement parameters. Selected bond lengths and angles are given in Table 2.

(a) Molecular structure of [1‐Na⋅(15c5)] and [1‐K⋅(18c6)] in the solid state. [1‐K⋅(18c6)] crystallizes with two crystallographically independent molecules in the asymmetric unit, only one of which is shown here. Hydrogen atoms and THF solvent molecules are omitted for clarity; 50 % displacement parameters. Selected bond lengths and angles are given in Table 2. With respect to the bonding situation, the η 2‐coordination is well in line with the IR data and further confirms the change in the bonding situation from a single to a double bond in 1‐K compared with the phosphonium salt. This assumption is further supported by the structural changes in the yldiides, that is, the shortening of the P−C and C1−C2 bond as well as the lengthening of the C−N bond. Hence, the crystal structures further argue for the predominance of the bonding situation C=C=N. Overall, the solid‐state structures show that the coordination modes in 1‐M are rather complex and diverse and strongly depend on the metal and additional Lewis bases. Both N‐ and Ccoordination are observed, although N‐coordination seems to be slightly more preferred for the harder cations (see 1‐Na⋅(15c5) and [(1‐Li)4⋅(18c6)]∞).

Computational studies

To gain insights into the different coordination modes of the alkali metals in yldiide 1‐M and their impact on the electronic distribution, computational studies were performed using monomeric systems (comparable to the structures of 1‐K and 1‐Na with crown ether) and dimethyl ether (Me2O) for completion of the coordination sphere of the metal. At first, the energetics of C‐ and N‐coordination of the yldiide were calculated for each metal (Li, Na, K), showing that N‐coordination in 1‐M is strongly preferred for lithium (ΔG=−28.4 kJ mol−1) and less favored for sodium (ΔG=−24.5 kJ mol−1) and particularly less favored for potassium (ΔG=−9.1 kJ mol−1). This trend is expected according to the HSAB (hard and soft acids and bases) concept. Despite the preference of N‐coordination for potassium, the calculations also reflect the preferred coordination of potassium to both carbon atoms of the α‐metalated nitrile, as found in the crystal structure (Figure 5). Energy optimizations of geometries with the metal solely bound to the α‐carbon atom always resulted in a shift of the metal towards the CN moiety. Investigation of the two highest occupied molecular orbitals (HOMO) showed that they are essentially independent of the metal coordination. This is expected because of the mostly electrostatic nature of the interactions between the yldiide and the metal cations. Thus, even for the most covalently bound lithium, the Wiberg bond indices (WBI) of the C−Li and N−Li bonds in 1‐Li amount to only 0.06 and 0.08, respectively. In addition, the energy of the orbitals is also unaffected by the alkali metal. Overall, the HOMO and HOMO‐1 are mostly undisturbed compared to the “metal‐free” yldiide. Both orbitals are distributed over the whole CC−N linkage thus confirming the delocalization of one of the lone pairs at carbon into the CN substituent (Figure 6, top). This is again well in line with the decrease in the C−N stretching frequency as well as the lengthening of the C−N bond and shortening of the CC bond upon metalation.
Figure 6

(a) Representations of the HOMO‐1 with and without coordination of lithium; (b) NPA charges and (c) WBIs depending on metal coordination.

(a) Representations of the HOMO‐1 with and without coordination of lithium; (b) NPA charges and (c) WBIs depending on metal coordination. In contrast to the orbitals, the introduction of a metal ion to the free yldiide has a significant impact on the charges in the molecule. As shown by the calculated natural population analysis (NPA) charges in Figure 6 b, the lithium atom shifts the charge density towards itself. This can clearly be seen by comparison with the “free” yldiide. Hence, a more negative charge is present at the ylidic carbon atom (−1.07 vs. −1.29) in the case of Ccoordination, and similarly the negative charge increases at the nitrogen atom (−0.57 vs. −0.83) upon N‐coordination. Interestingly, this charge transfer is most pronounced for lithium. For example, in case of N‐metalation, the charge at the nitrogen atom becomes more positive, from −0.83 to −0.78 and −0.75 when going from lithium to sodium and potassium, respectively. This is somewhat counter intuitive to the more electropositive nature of the heavier metals and their higher tendency for ionic bonding. However, the higher covalency of the bonds to lithium results in stronger orbital interactions and thus probably in a more efficient transfer of electron density from the ylidic carbon atom into the π orbitals of the C−N bond. This corroborates the decrease in the WBI of the C−N bond from potassium (WBI(C‐N) = 2.41) to lithium (WBI(C‐N) = 2.32) (Figure 6 c and the Supporting Information) and the particularly low stretching frequencies for the lithium compounds observed in the IR spectra (Table 1). As a next step, energy decomposition analysis with natural orbitals for chemical valency (EDA‐NOCV)29 calculations for all six different fragmentation patterns were performed including C‐ and N‐coordination of the metal cation (see the Supporting Information). To this end, the yldiide was fragmented once by cleavage of the P−C bond and once by cleavage of the CC bond to independently evaluate the respective bonding situation. The different types of bonding arise from different types of covalent interactions (captodative, electron‐sharing and ylidic) and thus indicate different chemical bonding situations. The importance of each bonding situation was judged based on the orbital energy value (ΔE Orb). The fragmentation with the lowest absolute ΔE Orb value was the best choice for describing the bonding situation in a molecule because the smallest alteration of the electronic charge distribution is required to yield the electronic structure of the molecule.5b,30 As shown by the results depicted in Table 3 for the P−C bond, the dative bonding is the energetically most favored bonding situation for the free yldiide and becomes even more favorable when changing from the unmetalated to the Ccoordinated form [ΔE Orb=−353 (free) vs. −296 kcal mol−1 (C‐form)]. In contrast, when the lithium is coordinated to the nitrogen atom, the dative binding mode becomes significantly less favorable [ΔE Orb=−353 (free) vs. −378 kcal mol−1 (N‐form)]. This is mostly caused by the higher electron density shift from the phosphine lone pair to the empty orbital of the carbon atom upon N‐coordination (see the Supporting Information). The changes caused by C‐ and N‐metalation are well in line with the determined NPA charges. As such, Ccoordination results in a charge accumulation at the ylidic carbon atom and thus in its reluctance to accept more electron density. This supports a dative interaction because here, less electron density is transferred than for a double bond or an ylidic interaction. Likewise, when the nitrogen atom is coordinated by the lithium atom, the whole fragment becomes more electron‐accepting, leading to a more favorable ylidic interaction relative to the dative bonding. The π backbonding is largely unaffected by the different coordination modes (see Supporting Information). Similarly, the double bonded form is unaffected by the different coordination modes.
Table 3

EDA‐NOCV results for the P−C and the C−C bond (values in kcal mol−1).

P−CC−C
dativedoubleylidicdativedoubleylidic
“free yldiide”
ΔE Pauli 536366555577424415
ΔE Elstat −266−228−476−251−211−217
ΔE Orb −353 −395 −439 −459 −457 −298
ΔE Int −83−256−359−133−244−100
N‐coordination
ΔE Pauli 563336413627462473
ΔE Elstat −272−205−324−267−229−298
ΔE Orb −378 −376 −380 −481 −462 −442
ΔE Int −86−245−291−121−229−267
C‐coordination
ΔE Pauli 416433447497397384
ΔE Elstat −205−234−327−304−284−206
ΔE Orb −296 −384 −353 −392 −421 −332
ΔE Int −86−185−233−199−308−154
EDA‐NOCV results for the P−C and the CC bond (values in kcal mol−1). In contrast to the P−C bond, the CC bond is clearly ylidic in the non‐metalated as well as in the Cmetalated form. Analogous to the P−C bond, the dative interaction between the two carbon atoms becomes more favorable when the lithium atom is coordinated by the central carbon atom due to the increased negative charge. This changes drastically in case of N‐coordination, for which all bonding situations are almost equal in energy. Most likely, this is also related to the “electron‐withdrawing” effect of the lithium atom, rendering the fragment less capable of forming an ylidic bond so that the C=C=N bonding mode becomes equal in energy. The dative interaction is largely unaffected by the N‐coordination of the metal. Thus, the EDA calculations suggest that, for C‐ and N‐coordination of the metal, the ylidic bonding situation in the CCN fragment is favorable, but that N‐coordination results in a higher preference for the double‐bonded structure. This observation is interesting because the calculations also predicted that N‐coordination is the most favored for lithium, which therefore should show the highest preference for the Ph3PC=C=N−M structure and for the weakest C−N bond. This is consistent with the IR data, which showed the lowest C−N stretching frequencies for the lithium compounds, and also with the calculated WBIs. The preferred structures according to the EDA analysis are depicted in Figure 7. Thus, no analogous metal‐like behavior of carbon analogous to bisylidic compounds is found here (see Figure 1).
Figure 7

Favored bonding situations depending on the metal coordination.

Favored bonding situations depending on the metal coordination. Altogether, our data show that the bonding situation in metalated ylides depends on the charge distribution in the molecule, which in turn depends on the nature of the metal and on the position where the alkali metal binds to the molecule. Despite the mostly electrostatic nature of the metal interactions (WBI close to zero), higher covalency in the M−X (X=C, N) bonds leads to more pronounced changes than the metal‐free system. The fact that the introduction of a cation (or simple cationic charge) has an impact on the bonding situation has also been reported for carbon monoxide, CO.31 These results suggest that in polar molecules, also other weak interactions, such as those of solvent molecules, may influence the bonding situation. This is particularly true for compounds with large charge accumulations. In fact, this is not surprising because solvent molecules have often changed the reactivity of compounds particularly in polar reagents like alkali or alkaline earth metal compounds.32 Indeed, the first systematic reactivity studies of the metalated ylides 1‐M with chlorophosphines indicated that the N‐ versus C‐attack of the phosphine strongly depends on the nature of the metal, but also on the solvent [coordinating (THF) versus weakly coordinating (toluene)]. These structure/solvent reactivity relationships are part of ongoing detailed investigations, but clearly indicate that the reactivity of the metalated ylide is complex and strongly influenced by the structure formed under the respective reaction conditions and hence on the resulting changes in the bonding situation.

Conclusions

We have reported a detailed study of the structure and bonding situation of the cyanido‐functionalized, metalated ylide [Ph3PCCN]M (1‐M, M=Li, Na, K). The compounds show a remarkable structural diversity depending on the metal and additional Lewis bases. In the crown‐ether‐ligated complexes, different monomeric structures were found for sodium and potassium, in which sodium prefers the coordination at the nitrogen of the cyanido group, whereas potassium exhibits an unusual η 2‐coordination mode binding to the CC linkage. IR spectroscopic studies demonstrate that the bonding situation in the CCN linkage changes upon metalation and is influenced by the nature of the metal. The changes are most pronounced for lithium, which exhibits the weakest C−N bond due to its more covalent interaction with the CCN fragment. Calculations on the bonding situation and electronic structure show that coordination of the metal to 1 significantly affects the NPA charges in the molecule. Compared with the metal‐free yldiide, the largest charge transfer is found for lithium, which efficiently pulls electron density towards itself. Energy decomposition analyses suggest that the bonding between the ylidic carbon atom and the cyanido group lies between an ylidic and a classical double bond. The latter becomes more important in case of the coordination of the metal at the nitrogen end rather than at the ylidic carbon atom.

Experimental Section

General methods and materials

All experiments were carried out under a dry, oxygen‐free argon atmosphere using standard Schlenk techniques. Solvents were dried using an MBraun SPS‐800 (THF, toluene, Et2O, CH2Cl2, n‐pentane, n‐hexane) or following standard procedures. All chemicals and reagents were purchased from Sigma–Aldrich, ABCR, Rockwood Lithium, or Acros Organics and used without further purification. Compound [1‐H]B was prepared in analogy to literature reports.33 Preparation of 1‐H: Phosphonium salt [1‐H]Br (29.0 g, 75.9 mmol, 1 equiv) and sodium hydride (1.91 g, 79.7 mmol, 1.05 equiv) were suspended in THF (250 mL ) and stirred for six days. The solvent was removed in vacuo and the remaining solid was extracted with CH2Cl2 (4×50 mL). The combined CH2Cl2 fractions were layered with 200 mL n‐pentane. After 3 weeks, colorless crystals formed, which were filtered and dried under vacuum, giving the product as a colorless solid (18.3 g, 60.7 mmol, 80 %). Spectroscopic data are in accordance with literature reports.13 Preparation of [(1‐Li): LiHMDS (1.44 g, 8.61 mmol, 15 equiv) was suspended in toluene (15 mL) and added through a filter cannula to 1‐H (500 mg, 1.67 mmol, 3.0 equiv). After stirring for 14 h, the solvent was removed in vacuo and n‐hexane (20 mL) was added. Upon standing for three days, the product precipitated from the solution. The supernatant solvent was removed with a syringe, the remaining solid was subsequently dried in vacuo and gave the product as a colorless solid (91 mg, 51.0 μmol, 9 %; 1 equiv). 1H NMR (400.1 MHz, [D8]THF): δ=7.85–7.60 (m, 18 H; CH ), 7.49–7.17 (m, 27 H; CH ), −0.16 ppm (s, 90 H; NSi(CH 3)3); 31P{1H} NMR (162.1 MHz, [D8]THF): δ=−1.9 ppm; 13C{1H} NMR (100.7 MHz, [D8]THF) δ=137.2 (d, 1 J CP=92.4 Hz; C , ), 134.4 (d, 2 J CP =20.1 Hz; CCN), 133.0 (d, 2 J CP=9.04 Hz; CH , ), 130.3 (s; CH , ), 128.4 (d, 3 J CP=11.0 Hz; CH , ), 6.0 (s; NSi(CH3)3), −3.4 ppm (br; PCCN); FT‐IR (KBr, in THF): =1989 cm−1. Preparation of [(1‐Na): NaHMDS (0.856 g, 4.67 mmol, 11.2 equiv) was dissolved in toluene (15 mL) and the supernatant solution was added through a filter cannula to 1‐H (1.00 g, 3.32 mmol, 8.0 equiv). After stirring the reaction mixture for 3 h, n‐hexane (20 mL) was added, resulting in the precipitation of a yellow solid from the solution. The suspension was filtered with a fritted glass funnel and dried under vacuum to give the product as yellow solid (0.945 g, 2.56 mmol, 77 %, 1 equiv). Single crystals suitable for X‐ray diffraction were obtained by slow diffusion of n‐hexane into a benzene solution. 1H NMR (250.1 MHz, [D8]THF) δ=7.9–7.6 (m, 48 H; CH ), 7.4–6.8 (m, 72 H; CH ), −0.2 ppm (s, 36 H; NSi(CH 3)3; 31P{1H} NMR (162.1 MHz, [D8]THF): δ=−5.1 ppm; 13C{1H} NMR (75.5 MHz, [D8]THF) δ=142.5 (d, 2 J CP=12.6 Hz; CCN), 137.9 (d, 1 J CP=84.5 Hz; C , ), 132.9 (d, 2 J CP=9.1 Hz; CH , ), 129.8 (d, 4 J CP=2.7 Hz; CH , ), 128.3 (d, 3 J CP=11.1 Hz; CH ,), 6.7 (s; NSi(CH3)3), −5.4 ppm (d, 1 J CP=68.2 Hz; PCCN); FT‐IR (KBr, in THF): =2008 cm−1. Preparation of 1‐K: Ylide 1‐H (4.00 g, 13.3 mmol, 1 equiv) and KHMDS (2.78 g, 13.9 mmol, 1.05 equiv) were suspended in Et2O (250 mL) and stirred for two days. The suspension was filtered with a glass frit, washed with Et2O (2×10 mL) and dried in vacuum to give the product as yellow solid (3.77 g, 11.1 mmol, 84 %; 1 equiv). 1H NMR (400.3 MHz, [D8]THF): δ=7.89–7.56 (m, CH ), 7.33–7.00 ppm (m, 9 H; CH ); 31P{1H} NMR (162.1 MHz, [D8]THF): δ=−10.5 ppm. 13C{1H} NMR (100.7 MHz, [D8]THF): δ=140.9 (d, 2 J CP=14.4 Hz; CCN), 138.6 (d, 1 J CP=83.1 Hz; C , ), 132.8 (d, 2 J CP=8.9 Hz; CH , ), 129.6 (s; CH , ), 128.3 (d, 3 J CP=10.9 Hz; CH , ), 0.8 ppm (d, 1 J CP=68.0 Hz; PCCN); FT‐IR (KBr, in THF): =2001 cm−1. Preparation of [1‐Na⋅(15c5)]: The compound [(1‐Na)8(NaHMDS)2] (400 mg, 1.08 mmol, 0.13 equiv) and 15‐crown‐5 (273 mg, 0.246 mL, 1.24 mmol, 1.2 equiv) were suspended in toluene (15 mL). After stirring for 3 h at room temperature, a red solid precipitated from the reaction mixture. The suspension was filtered with a glass frit and the obtained solid was dried under vacuum to give the product as a red solid (340 mg, 0.600 mmol, 55 %; 1 equiv). Single crystals suitable for X‐ray diffraction were obtained by slow diffusion of n‐hexane into a toluene solution. 1H NMR (400.3 MHz, [D8]THF) δ=7.94–7.61 (m, 6 H; CH ), 7.38–7.04 (m, 9 H; CH ), 3.56 ppm (s, 24 H; CH 2); 31P{1H} NMR (162.1 MHz, [D8]THF): δ=−10.9 ppm; 13C{1H} NMR (100.7 MHz, [D8]THF): δ=141.1 (d, 2 J CP=15.6 Hz; CCN), 139.1 (d, 1 J CP=83.0 Hz; C , ), 132.9 (d, 2 J CP=8.6 Hz; CH , ), 129.4 (d, 2 J CP=2.7 Hz; CH , ), 128.1 (d, 3 J CP=10.8 Hz; CH , ), 70.5 (s; CH215c5), −0.7 ppm (d, 1 J CP=72.7 Hz; PCCN); FT‐IR (KBr, in THF): =2023 cm−1. Preparation of [1‐K⋅(18c6)]: The compound 1‐K (100 mg, 29.5 μmol, 1 equiv) and 18‐crown‐6 (93.0 mg, 35.4 μmol, 1.2 equiv) were suspended in THF (2 mL). After stirring for 3 h at room temperature, all solids were completely dissolved. Subsequent addition of n‐hexane (6 mL) resulted in the precipitation of a yellow solid, which was filtered with a glass frit and dried under vacuum. The product was obtained as a yellow solid (35.1 mg, 58.0 μmol, 20 %, 1 equiv). Single crystals suitable for X‐ray diffraction were obtained by slow diffusion of n‐pentane into a THF solution. NMR spectra in the presence of excessive 18‐crown‐6: 1H NMR (400.1 MHz, [D8]THF) δ=7.85–7.65 (m, 6 H; CH ), 7.26–7.11 (m, 9 H; CH ); 31P{1H} NMR (162.1 MHz, [D8]THF): δ=−19.9 ppm; 13C{1H} NMR (100.7 MHz, [D8]THF): δ=140.6 (d, 1 J CP=77.8 Hz; C , ),), 137.3 (d, 2 J CP=15.6 Hz; CCN), 132.9 (d, 2 J CP=8.5 Hz; CH , ), 128.6 (s; CH , ), 127.5 (d, 3 J CP=10.5 Hz; CH , ), 5.4 ppm (d, 1 J CP=59.2 Hz; PCCN); FT‐IR (KBr, in THF): =2014 cm−1.

Single‐crystal X‐ray diffraction

Data collection of all compounds was conducted either with a Bruker X8‐APEX II ([1‐H]Br, 1H⋅LiBr), STOE‐IPDS 2 ([(1‐Na)8(NaHMDS)2]), Oxford SuperNova (1‐H, [(1‐Li)3(LiHMDS)5], [(1‐Li)3(15c5)]2, [1‐Na⋅(15c5)], [1‐K⋅18c6]) or Oxford Synergy diffractometer [(1‐Li)4(18c6)]∞. Crystals for investigation were coated in an inert oil (perfluoropolyalkyl ether) and mounted on a fiber loop and placed in a cold N2 stream on the diffractometer. The data were collected using the APEX 2 (Bruker), X‐Area (Stoe) or CrysAlisPro (Oxford) software and the crystal structure determinations were accomplished at 100 K using CuKα radiation (λ=1.54184 Å) or MoKα radiation (λ=0.71073 Å). The structures were solved using direct methods, refined with the SHELX software package34 and expanded using Fourier techniques. CCDC 1875234 ([(1‐Na)8(NaHMDS)2]), 1875235 (1H⋅LiBr), 1875236 (1H), 1875237 ([1‐H]Br), 1875238 ([1‐K⋅18c6]), 1875239 ([(1‐Li)3(15c5)]2), 1875240 ([(1‐Li)4(18c6)]∞), 1875241 ([(1‐Li)3(LiHMDS)5]) and 1875242 ([1‐Na⋅15c5]) contain the supplementary crystallographic data for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre. Further information including details to the structure solutions of disorderd parts and solvent molecules is provided in the Supporting Information.

NMR spectroscopy

1H, 13C{1H}, and 31P{1H} NMR spectra were recorded on Avance‐400 or DPX‐250 spectrometers at 25 °C if not stated otherwise. All values of the chemical shift are in ppm regarding the δ‐scale. All spin‐spin coupling constants (J) are printed in Hertz (Hz). To display multiplicities and signal forms correctly the following abbreviations were used: s=singlet, d=doublet, t=triplet, m=multiplet, br=broad signal. Signal assignment was supported by DEPT, APT, HSQC and HMBC experiments and by literature studies on similar compounds.8 NMR spectra of all compounds are depicted in the Supporting Information.

Vibrational spectroscopy

IR spectra were recorded on a Thermo Nicolet iS5 FT‐IR in transmission mode with a Specac “Omni‐cell” with KBr plates and a 0.1 mm spacer at 22 °C. The samples were measured according to standardized procedure: 5 mg of the substance were dissolved in 1 mL of THF in a glovebox (in case of the crown ether complexes 1 equiv of crown ether was added). The solution was added into the IR cell using a syringe. The cell was closed, taken outside the glovebox and an IR spectrum was recorded. Geometry optimizations without symmetry constraints were carried out with the Gaussian09 program package35, 36 at the B3LYP37/6‐31+G*38 level of theory. Initial structures were obtained by using GaussView 5.0.9.39 Stationary points were characterized as minima by calculating the Hessian matrix analytically at this level of theory.40 Thermodynamic corrections and Kohn–Sham orbitals were taken from these calculations. The standard state for all thermodynamic data is 298.15 K and 1 atm. Single‐point energies at the B3LYP/6‐31+G* optimized geometries were calculated at the B3LYP/6‐311++G** level. The NBO41 analyses were carried out with the internal module of Gaussian0942 at the B3LYP/6‐31+G* level of theory. Coordinates of all optimized structures are given in the Supporting Information. The EDA‐NOCV calculations were carried out with the program package ADF2017.111.43 The B3LYP/6‐31+G* geometries were used. BP8644 was chosen with uncontracted Slater type orbitals (STOs) as basis functions.45 The latter basis sets for all elements have triple‐ζ quality augmented by two sets of polarization functions (ADF basis set TZ2P). Core electrons (i.e., 1s for second‐ and [He]2s2p for third‐row atoms) were treated by the frozen‐core approximation. This level of theory is denoted BP86/TZ2P. An auxiliary set of s, p, d, f and g STOs was used to fit the molecular densities and to represent the Coulomb and exchange potentials accurately in each SCF cycle.46 Scalar relativistic effects were incorporated by applying the zeroth‐order regular approximation (ZORA) in all ADF calculations.47 The interatomic interactions were investigated by means of an energy‐decomposition analysis (EDA) developed independently by Morokuma48 and by Ziegler and Rauk49 in conjunction with the natural orbitals for chemical valence (NOCV).29,50 Further details are provided in the Supporting Information.

Conflict of interest

The authors declare no conflict of interest. As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors. Supplementary Click here for additional data file.
  60 in total

1.  A short history of SHELX.

Authors:  George M Sheldrick
Journal:  Acta Crystallogr A       Date:  2007-12-21       Impact factor: 2.290

2.  Carbenes stabilized by ylides: pushing the limits.

Authors:  Alois Fürstner; Manuel Alcarazo; Karin Radkowski; Christian W Lehmann
Journal:  Angew Chem Int Ed Engl       Date:  2008       Impact factor: 15.336

3.  Chemoselective Alkylations with N- and C-Metalated Nitriles.

Authors:  Xun Yang; Dinesh Nath; Fraser F Fleming
Journal:  Org Lett       Date:  2015-10-02       Impact factor: 6.005

4.  Metalated nitriles: organolithium, -magnesium, and -copper exchange of alpha-halonitriles.

Authors:  Fraser F Fleming; Zhiyu Zhang; Wang Liu; Paul Knochel
Journal:  J Org Chem       Date:  2005-03-18       Impact factor: 4.354

5.  Ylide Ligands as Building Blocks for Bioactive Group 11 Metal Complexes.

Authors:  Alice Johnson; Isabel Marzo; M Concepción Gimeno
Journal:  Chemistry       Date:  2018-07-13       Impact factor: 5.236

6.  β-Diketiminate calcium hydride complexes: the importance of solvent effects.

Authors:  Andrea Causero; Gerd Ballmann; Jürgen Pahl; Christian Färber; Julia Intemann; Sjoerd Harder
Journal:  Dalton Trans       Date:  2017-02-14       Impact factor: 4.390

7.  Phosphido-borane and phosphido-bis(borane) complexes of the alkali metals, a comparative study.

Authors:  Keith Izod; James M Watson; William Clegg; Ross W Harrington
Journal:  Inorg Chem       Date:  2013-01-15       Impact factor: 5.165

8.  Synergic sedation of sensitive anions: alkali-mediated zincation of cyclic ethers and ethene.

Authors:  Alan R Kennedy; Jan Klett; Robert E Mulvey; Dominic S Wright
Journal:  Science       Date:  2009-10-30       Impact factor: 47.728

9.  Coordination chemistry at carbon: the patchwork family comprising (Ph3P)2C, (Ph3P)C(C2H4), and (C2H4)2C.

Authors:  Hubert Schmidbaur; Annette Schier
Journal:  Angew Chem Int Ed Engl       Date:  2012-12-03       Impact factor: 15.336

10.  Structurally defined potassium-mediated zincation of pyridine and 4-R-substituted pyridines (R = Et, iPr, tBu, Ph, and Me2N) by using dialkyl-TMP-zincate bases.

Authors:  William Clegg; Ben Conway; David V Graham; Eva Hevia; Alan R Kennedy; Robert E Mulvey; Luca Russo; Dominic S Wright
Journal:  Chemistry       Date:  2009-07-20       Impact factor: 5.236

View more
  5 in total

1.  Ylide-Functionalization via Metalated Ylides: Synthesis and Structural Properties.

Authors:  Christopher Schwarz; Thorsten Scherpf; Ilja Rodstein; Julia Weismann; Kai-Stephan Feichtner; Viktoria H Gessner
Journal:  ChemistryOpen       Date:  2019-05-15       Impact factor: 2.911

2.  Single-Site and Cooperative Bond Activation Reactions with Ylide-Functionalized Tetrylenes: A Computational Study.

Authors:  Henning Steinert; Julian Löffler; Viktoria H Gessner
Journal:  Eur J Inorg Chem       Date:  2021-10-26       Impact factor: 2.551

3.  Isolation of the Metalated Ylides [Ph3 P-C-CN]M (M=Li, Na, K): Influence of the Metal Ion on the Structure and Bonding Situation.

Authors:  Christopher Schwarz; Lennart T Scharf; Thorsten Scherpf; Julia Weismann; Viktoria H Gessner
Journal:  Chemistry       Date:  2019-01-29       Impact factor: 5.236

4.  Synthesis, Crystal and Electronic Structures of a Thiophosphinoyl- and Amino-Substituted Metallated Ylide.

Authors:  Mike Jörges; Alexander Kroll; Leif Kelling; Richard Gauld; Bert Mallick; Stefan M Huber; Viktoria H Gessner
Journal:  ChemistryOpen       Date:  2021-09-27       Impact factor: 2.630

Review 5.  Phosphorus-ylides: powerful substituents for the stabilization of reactive main group compounds.

Authors:  Abir Sarbajna; V S V S N Swamy; Viktoria H Gessner
Journal:  Chem Sci       Date:  2020-07-17       Impact factor: 9.825

  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.