Literature DB >> 26751824

The Synthesis of Organic Molecules of Intrinsic Microporosity Designed to Frustrate Efficient Molecular Packing.

Rupert G D Taylor1, C Grazia Bezzu2, Mariolino Carta2, Kadhum J Msayib2, Jonathan Walker1, Rhys Short1, Benson M Kariuki1, Neil B McKeown3.   

Abstract

Efficient reactions between fluorine-functionalised biphenyl and terphenyl derivatives with catechol-functionalised terminal groups provide a route to large, discrete organic molecules of intrinsic microporosity (OMIMs) that provide porous solids solely by their inefficient packing. By altering the size and substituent bulk of the terminal groups, a number of soluble compounds with apparent BET surface areas in excess of 600 m(2)  g(-1) are produced. The efficiency of OMIM structural units for generating microporosity is in the order: propellane>triptycene>hexaphenylbenzene>spirobifluorene>naphthyl=phenyl. The introduction of bulky hydrocarbon substituents significantly enhances microporosity by further reducing packing efficiency. These results are consistent with findings from previously reported packing simulation studies. The introduction of methyl groups at the bridgehead position of triptycene units reduces intrinsic microporosity. This is presumably due to their internal position within the OMIM structure so that they occupy space, but unlike peripheral substituents they do not contribute to the generation of free volume by inefficient packing.
© 2015 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.

Entities:  

Keywords:  adsorption; aromatic nucleophilic substitution; microporous materials; molecular packing; triptycene

Year:  2016        PMID: 26751824      PMCID: PMC4755154          DOI: 10.1002/chem.201504212

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


Introduction

There is increasing interest in the study of microporous materials made using organic components,1 driven by potential applications in gas storage,2 selective gas separation membranes,3 heterogeneous catalysis,4 and nanoparticle encapsulation.5 Broadly, microporous materials can be subdivided into two structural classes: ordered materials, the porosity of which typically arises from their crystalline structure (for example, zeolites,6 metal–organic frameworks (MOFs),4a, 7 covalent organic frameworks (COFs),8 and porous molecular crystals9), and amorphous materials, the porosity of which arises from a disordered framework (for example, activated carbons,10 hypercrosslinked polymers,11 and other porous polymer networks12). The majority of both ordered and amorphous microporous materials are based on network structures, which are inherently insoluble. However, demand for solution‐processable porous materials has led to interest in non‐network polymeric and molecular materials. For example, microporosity can be generated from the inability of soluble polymer chains to pack together efficiently in the solid state, as demonstrated by the polymers of intrinsic microporosity (PIMs).3, 12, 13 Ordered porous molecular materials may be obtained from the self‐assembly of the component molecules into a crystalline porous packing arrangement,9, 14 or from the crystallisation of molecular cages, which act as prefabricated pores. For example, the groups of Cooper and Mastalerz have demonstrated that molecular cages, produced from multiple Schiff base condensations between simple amine and aldehyde precursors, can produce soluble materials with apparent BET surface areas (SABET) approaching 1400 m2 g−1.15 Recently, Mastalerz et al. have created highly porous cages from triptycene‐based precursors that were suitably functionalised to undergo Schiff base condensations16 or boronic ester formation.17 Using the latter, a crystalline molecular material with a SABET of 3758 m2 g−1 was produced. The prefabricated pore structure of macrocycles18 and cages15c,15h, 19 can also generate microporosity from amorphous packing in the solid state. However, amorphous microporous materials derived from discrete organic molecules that are neither cages nor macrocycles are still relatively rare.20 We are engaged in a joint modelling21 and synthesis20b programme to investigate organic molecules of intrinsic microporosity (OMIMs), which we define as discrete molecular compounds that are designed to generate microporous materials solely from their inability to pack efficiently in the solid state. Theory suggests that for both two‐22 and three‐dimensional23 objects, the most inefficient packing is produced when the constituent shapes possess highly concave faces. Accordingly, our design strategy is to combine rigid, aromatic precursors (Figure 1) to form large molecules with multiple concavities (Figure 2). Precursors are composed of core units (1, 2) that possess ortho‐difluorine functionality and terminal units (3‐9) with catechol functionality, allowing for efficient combination by a double nucleophilic aromatic substitution to give dioxan fused units, a reaction previously utilised successfully in the synthesis of PIMs.13, 24 We selected 4,4′‐dicyano‐2,2′,3,3′,5,5′,6,6′‐octafluorobiphenyl (1)25 as an OMIM core owing to the two electron‐withdrawing nitrile groups, encouraging efficient dioxan formation. By serendipity, 2,2′,2′′,3,3′′,5,5′,5′′,6,6′,6′′‐undecafluoro‐[1,1′:3′,1′′‐terphenyl]‐4,4′,4′′‐tricarbonitrile (2) was obtained as a by‐product during the synthesis of 1, and was investigated owing to its higher reactive functionality, allowing access to highly substituted materials. Catechol (3)‐ and naphthalene (4)‐based arms were selected as base‐line controls to examine the effect of increasing arm length and adding substituents. Triptycene (5 and 6)26, 27 and bulkier arm groups, based on spirobifluorene (7),28 propellane (8),29 and hexaphenylbenzene (9),30 were selected owing to their high internal free volume26b and proven porosity‐enhancing properties in PIMs and other microporous polymers.
Figure 1

Precursors 1–9 (top) and general reaction scheme (box) for reaction of catechol terminal units with fluorinated cores 1 (left) and 2 (right).

Figure 2

Structures of compounds 10–39, including OMIM‐1 to OMIM‐18.

Precursors 1–9 (top) and general reaction scheme (box) for reaction of catechol terminal units with fluorinated cores 1 (left) and 2 (right). Structures of compounds 10–39, including OMIM‐1 to OMIM‐18.

Results and Discussion

Synthesis

The biphenyl fluorinated core precursor 1 was prepared in moderate yield (34 %) by the hexaethylphosphorous triamide mediated coupling of pentafluorobenzonitrile by adapting a procedure described in a patent.25 The terphenyl core precursor 2 was isolated as an unexpected by‐product in low yield (6 %) but in sufficient quantity to facilitate its use as a core for adduct synthesis. Each adduct was prepared by mixing the core precursor 1 or 2 with an excess of an arm unit chosen from 3–9 (Figure 1) in anhydrous DMF at 65 °C for 48 h in the presence of potassium carbonate. Each adduct (Figure 2) was purified either by column chromatography or by repeated trituration to remove more soluble impurities. Owing to the lack of symmetry of many of the substituted precursors, several of the reported OMIMs/adducts were isolated as mixtures of inseparable regioisomers (Table 1). The molecular mass, homogeneity and discrete nature of each adduct was confirmed by matrix‐assisted laser desorption ionisation mass spectrometry (MALDI‐MS) and gel permeation chromatography (GPC). For each adduct, the polydispersity index obtained from GPC was less than 1.1.
Table 1

OMIM/adduct data for compounds 10–39.

Compound numberOMIM number[a] PrecursorsApparent BET surface area [m2 g−1][b] Pore volume [cm3 g−1]Yield [%]Isomerically pure (Y/N)CCDC number
10 [c] 1+3 a 70.0374Y1406071
11 [c] 1+3 b 410.1580N
12 [c] 1+3 c 670.2370N
13 1+3 d 510.1390Y
14 2+3 a 130.0587Y1406073
15 2+3 b 70.0164N
16 2+3 c 1020.3788N
17 [c] 1+4 a 250.0544Y1406072
18 [c] 1+4 b 2600.3364N
19 [c] 1+4 c 1320.2539N
20 2+4 a 290.1113Y
21 2+4 b 2590.3543N
22 OMIM‐7 2+4 c 3470.4197N
23 [c,d,e] OMIM‐1 1+5 a 4850.4073Y955894
24 OMIM‐8 1+5 b 4620.3351Y1406070
25 [d,e] OMIM‐4 1+5 c 6540.5491N
26 OMIM‐9 1+5 d 5990.4247N
27 [e] OMIM‐5 1+5 e 7020.6083Y973327
28 OMIM‐10 2+5 a 4230.4430Y
29 OMIM‐11 2+5 b 3510.3079Y
30 OMIM‐12 2+5 c 7260.7281N
31 OMIM‐13 2+5 d 6510.4765N
32 OMIM‐14 2+5 e 6980.4477Y
33 [e] OMIM‐6 1+6 6220.6486N
34 [d] OMIM‐2 1+7 3330.2890N
35 [d] OMIM‐3 1+8 5950.4556N
36 OMIM‐15 1+9 4070.4673N
37 OMIM‐16 2+6 5910.4978N
38 OMIM‐17 2+7 4710.3688N
39 OMIM‐18 2+8 6120.4984N

[a] We only classify those materials with an apparent BET surface area of more than 300 m2 g−1 as OMIMs. [b] Measured experimentally at 77 K following degassing at 135 °C under vacuum for 15 h. [c] Structure previously reported for packing simulation.21b [d] Structure previously reported for packing simulation.21a, 21c [e] Structure and experimental data previously reported.20b

OMIM/adduct data for compounds 10–39. [a] We only classify those materials with an apparent BET surface area of more than 300 m2 g−1 as OMIMs. [b] Measured experimentally at 77 K following degassing at 135 °C under vacuum for 15 h. [c] Structure previously reported for packing simulation.21b [d] Structure previously reported for packing simulation.21a, 21c [e] Structure and experimental data previously reported.20b The potential of 1 and 2 as OMIM cores was first tested by their reaction with excess 3 a to ascertain their reactivity towards aromatic nucleophilic substitution. After purification, crystals of the resultant adducts (10 and 14) were achieved by slow diffusion of methanol into chloroform solutions. X‐ray crystallography coupled with 19F NMR spectroscopy (Figure 3) revealed these adducts to be the tetra‐substituted adduct (10), with a near orthogonal relationship between the two long struts of the molecule, and the penta‐substituted adduct (14), with one residual fluorine atom on the central aryl ring. This substitution pattern was found consistently with all catechol adducts of 1 and 2 in this study. Bridged products, in which a single arm bridges two cores (Supporting Information, Figure S1) were often found as trace impurities; however, these were readily removed by column chromatography. Detailed synthetic procedures and spectroscopic data for all novel precursors and adducts are given in the Supporting Information.
Figure 3

Molecular structures of 10 (top) and 14 (bottom) obtained from XRD analysis with their respective 19F NMR spectra, indicating total fluorine substitution in 10 and the retention of a single fluorine atom on the middle phenyl ring of the terphenyl unit in 14.

Molecular structures of 10 (top) and 14 (bottom) obtained from XRD analysis with their respective 19F NMR spectra, indicating total fluorine substitution in 10 and the retention of a single fluorine atom on the middle phenyl ring of the terphenyl unit in 14.

Discussion

The simplest adduct (10) was found to be effectively non‐microporous, with an apparent BET surface area (SABET) of 7 m2 g−1 as measured by nitrogen sorption at 77 K and it was highly insoluble, suggesting that the constituent molecules are able to pack together in an efficient manner, despite its awkward displaced cruciform structure. As predicted by molecular modelling,21b the introduction of bulky tert‐butyl groups to the periphery of the adducts (11, 12) improved apparent microporosity (SABET=41 and 67 m2 g−1, respectively). Both adducts were isolated as mixtures of regioisomers and were found to be highly soluble in common solvents (for example dichloromethane and THF), which combined with the increased porosity, suggests a disruption of the cohesive interactions between the constituent molecules. To further study this effect, 2,3‐dihydroxy‐5,5,8,8‐tetramethyl‐6,6,7,7‐tetrahydronapthalene (3 d) was prepared and combined with 1 to give 13, which demonstrates comparable properties to 12. Simple adducts of core 2 were also prepared using three catechol precursors 3 a–c. Much like its biphenyl analogue, terphenyl‐based adduct 14 was poorly soluble in common solvents and possesses a negligible apparent SABET of 13 m2 g−1. However, the addition of two tert‐butyl groups per arm gave a highly soluble material (16) with a greater apparent SABET of 102 m2 g−1 over its biphenyl analogue (12, 67 m2 g−1), suggesting that the higher functionality of the terphenyl core generates a more porous material. Adducts 17 and 20, derived from the reaction between unsubstituted naphthalene‐2,3‐diol precursor (4 a) and cores 1 and 2, respectively, gave poorly soluble materials with no apparent microporosity (Table 1), similar to adducts 10 and 14 prepared using unsubstituted catechol (3 a). X‐ray crystallographic analysis of 10, 14, and 17 (Figure 4) demonstrates that for each of these adducts there are extensive π–π interactions between the aromatic arms causing the molecules to pack efficiently, resulting in poor solubility and negligible porosity. However, the use of tert‐butyl substituted naphthalene‐2,3‐diol (4 b) as a precursor gave biphenyl 18 and terphenyl 21 with significant apparent SABET of 260 and 259 m2 g−1, respectively. This enhanced porosity is attributed to the tert‐butyl substituents prohibiting close interactions between aromatic arms of adducts resulting in less efficient packing.
Figure 4

X‐ray crystal structures of 10 (top), 14 (middle), and 17 (bottom), revealing π–π induced staking of the constituent molecules. Protons and solvent molecules are removed for clarity.

X‐ray crystal structures of 10 (top), 14 (middle), and 17 (bottom), revealing π–π induced staking of the constituent molecules. Protons and solvent molecules are removed for clarity. Given its success as a component for making organic microporous materials,26 triptycene derived arms were also investigated for making OMIMs. Hence, cores 1 and 2 were combined with triptycene‐2,3‐diol (5 a)27 to give 23 (OMIM‐1)20b and 28 (OMIM‐10), respectively. Both OMIM‐1 and OMIM‐10 proved to be highly soluble and gave amorphous solids with apparent SABET of 485 and 423 m2 g−1, respectively. The addition of bulky terminal substituents to give triptycene precursors 5 c and 5 e enhances porosity in the resulting OMIMs, irrespective of whether they are attached to a biphenyl (OMIM‐4 and OMIM‐5) or terphenyl core (OMIM‐12 and OMIM‐14), with each demonstrating a high apparent SABET in the range 654–726 m2 g−1. OMIM‐12 possesses intrinsic microporosity greater than many PIMs and comparable to that of the archetypal PIM‐1 (ca. 780 m2 g−1).13a OMIM‐12 is of a similar porosity to a π‐extended triptycene recently reported by the Mastalerz group (SABET=754 m2 g−1),20c which were similarly designed to pack inefficiently. Notably, triptycene peripheral units substituted with methyl groups at their bridgehead positions (5 b, d) gave OMIMs with lower apparent values of SABET relative to their non‐methyl containing counterparts. For example, OMIM‐11 has an apparent SABET of 351 m2 g−1, as compared to 423 m2 g−1 for OMIM‐10. The values of SABET for the tert‐butyl substituted analogues (OMIM‐12 and OMIM‐13) also differ by a similar amount (75 m2 g−1). It appears that the space adjacent to the bridgehead in triptycene terminated OMIMs is directly contributing to the porosity of the material, hence, filling this space with a methyl group reduces the amount of intrinsic microporosity that can be generated during the amorphous packing of the molecules. Indirect evidence for this feature of amorphous packing comes from the single‐crystal XRD analysis of OMIM‐1 (Figure 5), which shows solvent‐filled channels defined partially by the triptycene bridgehead positions. Similar local ordering may occur within the amorphous packing of triptycene‐based OMIMs. Further XRD analysis of OMIM‐1, OMIM‐5, and OMIM‐8 is presented in the Supporting Information, Figure S2–S4.
Figure 5

Space filling view (top) and expanded view (bottom) of channels running through OMIM‐1.

Space filling view (top) and expanded view (bottom) of channels running through OMIM‐1. An investigation to determine the suitability of various rigid structural building units as peripheral arms for OMIM construction was performed by combining catechol derivatives of benzenopentacene (6)20b spirobifluorene (7),28 propellane (8),29 and hexaphenylbenzene (9)30 with cores 1 and 2. Benzenopentacene containing OMIM‐6 and OMIM‐16, were found to possess apparent SABET in the range 10–15 % lower than their more symmetrical isomeric triptycene‐based counterparts OMIM‐5 and OMIM‐14. It is probable that this is due to greater intermolecular interactions between the long struts of the benzenopentacene arms as compared to those of triptycene. From an OMIM comparison, it can be concluded that the efficiency of unsubstituted structural units for generating microporosity is in the following order: propellane>triptycene>hexaphenylbenzene>spirobifluorene>naphthyl=phenyl. This is consistent with findings from previously reported packing simulation studies.21

Conclusions

Molecular adducts, designed to possess well‐defined concavities, were synthesised using the reaction between fluorinated biphenyl and terphenyl cores and peripheral arms chosen from a diverse range of rigid structural units. The intrinsic microporosity generated by the inefficient packing of these adducts was evaluated using nitrogen sorption at 77 K. The use of small, planar, non‐substituted arms gave insoluble materials with negligible surface areas (<30 m2 g−1). However, by using arms composed of bulky rigid structural units, OMIMs were prepared with apparent SABET in the range 333–612 m2 g−1. The efficiency of unsubstituted structural units for generating microporosity is in the following order: propellane>triptycene>hexaphenylbenzene>spirobifluorene>naphthyl=phenyl. Substitution of these arms with bulky groups further enhanced microporosity (up to 726 m2 g−1), which is presumably due to reducing intermolecular cohesive interactions. In contrast, the introduction of methyl groups at the bridgehead position of triptycene units reduced intrinsic microporosity. In this case, the internal position of the methyl groups within the OMIM structure means that they occupy space but, unlike peripheral substituents, cannot contribute to the generation of free volume by frustrating packing.

Experimental Section

Example synthesis of adduct 10

4,4′‐Dicyano‐2,2′,3,3′,5,5′,6,6′‐octafluorobiphenyl (1; 0.230 g, 0.661 mmol) and catechol (3 a) were added to an oven‐dried flask and purged with nitrogen. Anhydrous DMF (10 mL) was then added via syringe and the reaction mixture heated to form a solution, at which point, oven dried potassium carbonate (0.306 g, 2.78 mmol) was quickly added, the reaction sealed under nitrogen flow, heated to 65 °C and left to stir for 48 h. After cooling to room temperature, the reaction mixture was poured into water (200 mL), acidified with 2 m HCl, and allowed to stir as a suspension for 2 h. The crude product was then collected by filtration, washed with water (200 mL) and methanol (200 mL), then dried under suction. Purification of the crude material was achieved by column chromatography (dichloromethane/hexane, 1:1, R f=0.4) to give 10 (0.307 g, 74 %) as a yellow powder (m.p.>300 °C); IR (CH2Cl2 film) 2234, 1495, 1440, 1309, 1272, 1253 cm−1; 1H NMR (500 MHz, CDCl3): δ=7.00 (d, J=8.0 Hz, 4 H, ArH), 6.96–6.93 (m, 4 H, ArH), 6.88–6.85 (m, 4 H, ArH), 6.66 ppm (d, J=8.0 Hz, 4 H, ArH); 13C NMR (125 MHz, CDCl3): δ=141.1, 140.8, 140.0, 136.2, 125.3, 125.2, 117.0, 114.4, 110.6, 91.7 ppm (one carbon missing); HRMS (EI+, m/z) calc. for C38H16N2O8: 628.0907 [M +], found 628.0903; GPC analysis (CHCl3) M=548, M w=578 g mol−1 relative to polystyrene, M w/M n=1.055; BET surface area=7 m2 g−1; total pore volume=0.03 cm3 g−1 at p/p o=0.98. Crystallography data (CHCl3/MeOH): triclinic, space group=P , a=11.0410(6) Å, b=11.3324(7) Å, c=14.2766(7) Å, α=113.045(5), β=90.659(4), γ=106.663(5), V=1559.14 Å3, Z=2, R=4.01. Full experimental details and spectroscopic data for all of the new compounds are given in the Supporting Information. CCDC 955894, 973327, 1406070, 1406071, 1406072, 1406073, and 1406074 contain the supplementary crystallographic data for this paper. These data are provided free of charge by The Cambridge Crystallographic Data Centre. 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.
  52 in total

1.  Optimal packings of superdisks and the role of symmetry.

Authors:  Y Jiao; F H Stillinger; S Torquato
Journal:  Phys Rev Lett       Date:  2008-06-20       Impact factor: 9.161

2.  Porous organic cages.

Authors:  Tomokazu Tozawa; James T A Jones; Shashikala I Swamy; Shan Jiang; Dave J Adams; Stephen Shakespeare; Rob Clowes; Darren Bradshaw; Tom Hasell; Samantha Y Chong; Chiu Tang; Stephen Thompson; Julia Parker; Abbie Trewin; John Bacsa; Alexandra M Z Slawin; Alexander Steiner; Andrew I Cooper
Journal:  Nat Mater       Date:  2009-10-25       Impact factor: 43.841

3.  A permanent mesoporous organic cage with an exceptionally high surface area.

Authors:  Gang Zhang; Oliver Presly; Fraser White; Iris M Oppel; Michael Mastalerz
Journal:  Angew Chem Int Ed Engl       Date:  2014-01-08       Impact factor: 15.336

4.  Energy-efficient hydrogen separation by AB-type ladder-polymer molecular sieves.

Authors:  Bader S Ghanem; Raja Swaidan; Xiaohua Ma; Eric Litwiller; Ingo Pinnau
Journal:  Adv Mater       Date:  2014-07-19       Impact factor: 30.849

5.  Ultra-microporous triptycene-based polyimide membranes for high-performance gas separation.

Authors:  Bader S Ghanem; Raja Swaidan; Eric Litwiller; Ingo Pinnau
Journal:  Adv Mater       Date:  2014-03-11       Impact factor: 30.849

6.  Microporous polycarbazole with high specific surface area for gas storage and separation.

Authors:  Qi Chen; Min Luo; Peter Hammershøj; Ding Zhou; Ying Han; Bo Wegge Laursen; Chao-Guo Yan; Bao-Hang Han
Journal:  J Am Chem Soc       Date:  2012-03-30       Impact factor: 15.419

7.  Hexaphenylbenzene-based polymers of intrinsic microporosity.

Authors:  Rhys Short; Mariolino Carta; C Grazia Bezzu; Detlev Fritsch; Benson M Kariuki; Neil B McKeown
Journal:  Chem Commun (Camb)       Date:  2011-05-17       Impact factor: 6.222

8.  Metal-organic framework materials with ultrahigh surface areas: is the sky the limit?

Authors:  Omar K Farha; Ibrahim Eryazici; Nak Cheon Jeong; Brad G Hauser; Christopher E Wilmer; Amy A Sarjeant; Randall Q Snurr; SonBinh T Nguyen; A Özgür Yazaydın; Joseph T Hupp
Journal:  J Am Chem Soc       Date:  2012-08-31       Impact factor: 15.419

9.  Polymers of intrinsic microporosity (PIMs): robust, solution-processable, organic nanoporous materials.

Authors:  Peter M Budd; Bader S Ghanem; Saad Makhseed; Neil B McKeown; Kadhum J Msayib; Carin E Tattershall
Journal:  Chem Commun (Camb)       Date:  2003-12-05       Impact factor: 6.222

10.  Hexabenzo[4.4.4]propellane: a helical molecular platform for the construction of electroactive materials.

Authors:  Paromita Debroy; Sergey V Lindeman; Rajendra Rathore
Journal:  Org Lett       Date:  2007-09-15       Impact factor: 6.005

View more
  3 in total

1.  Polymer ultrapermeability from the inefficient packing of 2D chains.

Authors:  Ian Rose; C Grazia Bezzu; Mariolino Carta; Bibiana Comesaña-Gándara; Elsa Lasseuguette; M Chiara Ferrari; Paola Bernardo; Gabriele Clarizia; Alessio Fuoco; Johannes C Jansen; Kyle E Hart; Thilanga P Liyana-Arachchi; Coray M Colina; Neil B McKeown
Journal:  Nat Mater       Date:  2017-07-31       Impact factor: 43.841

2.  The Synthesis of Organic Molecules of Intrinsic Microporosity Designed to Frustrate Efficient Molecular Packing.

Authors:  Rupert G D Taylor; C Grazia Bezzu; Mariolino Carta; Kadhum J Msayib; Jonathan Walker; Rhys Short; Benson M Kariuki; Neil B McKeown
Journal:  Chemistry       Date:  2016-01-11       Impact factor: 5.236

3.  Arene Selectivity by a Flexible Coordination Polymer Host.

Authors:  James S Wright; Iñigo J Vitórica-Yrezábal; Stephen P Thompson; Lee Brammer
Journal:  Chemistry       Date:  2016-08-02       Impact factor: 5.236

  3 in total

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