掺杂对铒、钕蒽基配合物的光二聚反应及光磁性能调控的影响

2024-01-20 03:56麻秀芳秘鹏飞鲍松松郑丽敏
无机化学学报 2024年1期
关键词:化工学院南京大学鹏飞

麻秀芳 秘鹏飞 鲍松松 郑丽敏

(南京大学化学化工学院,配位化学国家重点实验室,南京 210023)

0 Introduction

Light-responsive single molecule magnets (SMMs)have been a major research hotspot in the field of molecular magnets in the past few decades, which have broad research and application prospects[1-6].Lanthanide-based SMMs are appealing candidates for highperformance SMMs owing to their significant single-ion magnetic anisotropy arising from strong spin-orbit coupling and crystal-field effects[7-9].However, due to the difference in the number of 4fnelectrons and thus the magnetic anisotropies, different lanthanide complexes exhibit distinct SMM behaviours in the same ligand field[10-12].Not only that, the characteristic luminescence of lanthanide complexes is also highly dependent on specific lanthanide ions[13-14].Generally, direct excitation of LnⅢions is not efficient as a result of the parity-forbiddenf-ftransitions and low molar absorptivity.Therefore, it is necessary to sensitize the lanthanide ions through the antenna ligands, and the sensitization efficiency is very much related to the energy gap between the ligand3Tstate and the emissive excited state of LnⅢ[15].

We have been interested in dysprosium-anthracene complexes showing light-responsive magnetic and luminescent properties due to the photocycloaddition reaction of the face-to-faceπ-πinteracted anthracene pairs[16-17].A typical example is [DyⅢ(SCN)2(NO3)(depma)2(4-hpy)2] (1Dy)[18], where depma stands for 9-diethylphosphono methylanthracene and 4-hpy is 4 - hydroxypyridine, which has the highest effective energy barrier (Ueff=277 K) among the known photoresponsive Ln-SMMs.1Dy shows a rapid photo-induced structural transformation to form compound [DyⅢ(SCN)2(NO3)(depma2)(4-hpy)2]n(2Dy) with chain structure,accompanied by changes in SMM and luminescent properties.Since the photocycloaddition of anthracene involves the formation of an excimer in which one anthracene in the excited singlet state isπ-πinteracted with another in the ground state, we envision that when energy transfer between the excimer and the lanthanide ion is effective, the process is interfered with by ions of different lanthanide elements.Another question then arises:is it possible to modulate with light the SMM and luminescent properties of lanthanideanthracene complexes that are unable to undergo a photocycloaddition reaction?

In this paper, we report three new complexes,[LnⅢ(SCN)2(NO3)(depma)2(4-hpy)2] [Ln=Er (1Er), Nd(2Nd), Y (3Y)] that are isomorphic to 1Dy.NdⅢor ErⅢions were chosen because their complexes may exhibit interesting SMM behaviour depending on the specific coordination environment[19-22].Meanwhile, it has been found that energy transfer between anthracene and NdⅢor ErⅢions is efficient, resulting in the near-infrared(NIR) emission from metal ions[23-24].Indeed, 1Er and 2Nd exhibit field-induced SMM behavior and NIR luminescence but fail to undergo photocycloaddition reaction under 395 nm UV light irradiation.In contrast, the diamagnetic complex 3Y can undergo photocycloaddition (Scheme 1).Interestingly, by doping ErⅢor NdⅢions in the 3Y lattice, we obtained diluted samples Ln0.1Y0.9(SCN)2(NO3)(depma)2(4-hpy)2(Ln=Er(1Er@Y), Nd (2Nd@Y)), which not only can undergo photocycloaddition but also are accompanied by changes in magnetic and optical properties.

Scheme 1 Photocycloaddition of anthracene groups occurs in 3Y,1Er@Y and 2Nd@Y,but not in pure 1Er and 2Nd

1 Experimental

1.1 Materials and methods

The 9-diethyl-phosphonomethylanthracene (depma) was synthesized according to the literature[25].All other starting materials were obtained from commercial sources without any purification.PE 240C analyzer was used to do the elemental analyses for C, H and N.The metal content in the doped samples was determined by inductively coupled plasma atomic emission spectrometer (ICP-AES,Avio500).The infrared spectra were collected on a Bruker Tensor 27 spectrometer in 4 000-400 cm-1region in pressed KBr pellets.Thermogravimetric analysis was recorded on a Mettler-Toledo TGA STARe thermal analyzer in the range of 30 -600 ℃under a nitrogen flow at a heating rate of 5 ℃·min-1.The instrument to collect powder X-ray diffraction (PXRD) data was Bruker D8 advance diffractometer with CuKαradiation(λ=0.154 06 nm,U=40 mV,I=40 mA) in a range of 5°-50° at room temperature.The1H NMR spectra were recorded on a BRUKER AVANCEⅢ400 MHz spectrometer with samples dissolved in DMSO-d6.The steady fluorescence spectra were attained at Bruker Spectrofluorimeter LS55.Timeresolved fluorescent decays were carried out on a Fluorolog TCSPC spectrofluorometer (Horiba Scientific)equipped with laser exciting at 365 nm.The dc (direct current) and ac (alternating current) magnetic susceptibility data were obtained using polycrystalline samples by a Quantum Design MPMS SQUID VSM magnetometer.

1.2 Synthesis of Er(SCN)2(NO3)(depma)2(4-hpy)2(1Er)

46.1 mg (0.1 mmol) Er(NO3)3·6H2O and 19.4 mg(0.2 mmol) KSCN were added to 2 mL CH3CN solution and stirred for 30 min at room temperature.The white precipitate was removed by filtration.Then 19.0 mg(0.2 mmol)4-hpy and 3 mL CH3CN solution of 65.6 mg(0.2 mmol) depma were added and stirred for about 15 min, and the yellow precipitate was filtered.The clear yellow filtrate was left to evaporate slowly, and block yellow crystals were precipitated after 2 d.Yield: 44.4%.Elemental Anal.Calcd.for C50H52ErN5O11P2S2(% ): C,50.37; H, 4.40; N, 5.87.Found (%): C, 50.54; H, 4.58;N,5.79.

1.3 Synthesis of Nd(SCN)2(NO3)(depma)2(4-hpy)2(2Nd)

Compound 2Nd was prepared by adopting an analogous method to 1Er except for the replacement of Er(NO3)3·6H2O with Nd(NO3)3·6H2O (34.8 mg, 0.1 mmol).Yield: 42.5%.Elemental Anal.Calcd.for C50H52NdN5O11P2S2(%): C, 51.36; H, 4.48; N, 5.99.Found(%):C,51.69;H,4.37;N,6.07.

1.4 Synthesis of Y(SCN)2(NO3)(depma)2(4-hpy)2(3Y)

Compound 3Y was prepared by adopting an analogous method to 1Er except for the replacement of Nd(NO3)3·6H2O with Y(NO3)3·6H2O (38.3 mg, 0.1 mmol).Yield: 51.4%.Elemental Anal.Calcd.for C50H52YN5O11P2S2(%): C, 53.91; H, 4.71; N, 6.29.Found(%):C,53.43;H,4.54;N,6.40.

1.5 Synthesis of Er0.1Y0.9(SCN)2(NO3)(depma)2(4-hpy)2(1Er@Y)

The synthesis of the doped 1Er@Y complex was performed following a similar procedure to 1Er except that a mixture of Er(NO3)3·6H2O (0.01 mmol) and Y(NO3)3·6H2O (0.09 mmol) took the place of pure Er(NO3)3·6H2O.The erbium content in the final crystalline product was 8.55% determined by ICP-AES.Elemental Anal.Calcd.for C50H52Er0.1Y0.9N5O11P2S2(%):C, 53.54; H, 4.67; N, 6.24.Found(%): C, 53.18; H,4.84;N,6.39.

1.6 Synthesis of Nd0.1Y0.9 (SCN)2(NO3)(depma)2(4-hpy)2(2Nd@Y)

The synthetic procedure of 1Nd@Y was similar to 1Er@Y except that Er(NO3)3·6H2O was replaced by Nd(NO3)3·6H2O.The neodymium content in the final crystalline product was 8.47% determined by ICP -AES.Elemental Anal.Calcd.for C50H52Nd0.1Y0.9N5O11P2S2(%): C, 53.65; H, 4.68; N, 6.26.Found(%): C,53.37;H,4.34;N,6.44.

1.7 Crystallographicdatacollectionand refinement

Suitable single crystals were selected and mounted on a glass rod for X -ray measurements.Single crystals were used for data collections on a Bruker D8 Venture diffractometer using graphite-monochromated MoKαradiation(λ=0.071 073 nm).The data were integrated using the Siemens SAINT program[26], with the intensities corrected for Lorentz factor, polarization, air absorption, and absorption due to variation in the path length through the detector faceplate.Empirical absorption corrections were applied using the SADABS program[27].The structures were solved by direct method and refined onF2by full-matrix least squares using Olex2[28].All the non - hydrogen atoms were refined anisotropically.All hydrogen atoms were either put in calculated positions or found from the difference Fourier maps and refined isotropically.Details of the crystal data and refinements are summarized in Table S1(Supporting information).Selected bond lengths and angles are given in Table S2.

CCDC: 2303850, 1Er; 2303856, 2Nd; 2303854,3Y.

2 Results and discussion

2.1 Crystal structure

Single crystal X-ray diffraction analyses revealed that 1Er-3Y are isomorphic (Table S1).The purity of these phases was confirmed by PXRD and IR measurements (Fig.S1 and S2).So, 1Er is selected as a representative for a detailed structure description.Compound 1Er crystallizes in the monoclinic system with theP21/mspace group.As shown in Fig.1, the asymmetric unit contains one ErⅢion, two SCN-ion, one NO3-ion, two depma ligands and two 4-hpy molecules.Each ErⅢis eight-coordinated by two oxygen atoms(O1, O4) from two depma ligands, two oxygen atoms(O9 and O10) from NO3-, two axial oxygen atoms (O7 and O8) from another two 4-hpy ligands and two nitrogen atoms (N1 and N2) from two SCN-ions.The 4-hpy is a neutral zwitterionic ligand with a protonated pyridine N atom and a formally negative O atom.Two 4-hpy ligands occupy the axial positions providing the axial Er—O bond lengths of 0.221 4(3) and 0.223 2(3)nm which are shorter than the other Er—O(N) bonds(0.228 7(3)-0.256 6(3)nm).The axial O—Er—O angle is 163.54(13)° (Table S2).According to the continuous shape analysis[29], the [ErO6N2] polyhedron is mostly close to a geometry of snub diphenoid J84 inD2dsymmetry(CShM=1.807,Table S3).

Fig.1 (a-c)Molecule structures of compounds 1Er(a),2Nd(b)and 3Y(c);(d-f)2D supramolecular framework layer viewed along the a-axis for 1Er(d),2Nd(e),and 3Y(f)

As shown in Fig.1d, the adjacent monomers are connected along theb-axis with hydrogen bond interactions between the coordinated NO3-or SCN-anions and the 4-hpy through N4—H4…S1,N3—H3…O11,N4—H4…O11 and C42—H42…O10 contacts (Table S4),forming a one-dimensional (1D) supramolecular chain.The chains are stacked and further stabilized by faceto-faceπ-πinteraction between anthracene moieties in 1Er (center-to-center(dcc): 0.378 0 nm and nearest C2—C9′ (dC—C) distances: 0.379 0 nm).The dcc is much smaller than the limiting distance of Schmidt′s rule (0.42 nm) for solid-phase photodimerization to occur, so it is possible that anthracene-based photodimerization may occur in 1Er.

Compounds 2Nd and 3Y are isostructural to 1Er except for the slight difference in their unit cell volumes.The cell volume of 2Nd (2.556 7 nm3) is much larger than that of 1Er (2.508 5 nm3) and 3Y (2.520 5 nm3),attributed to the increased ionic radii of NdⅢcompared to ErⅢand YⅢ.As a result, the axial Nd—O(0.231 6(6), 0.233 7(6) nm) and equatorial Nd—O(N)(0.240 7(6)-0.263 6(6) nm) bond lengths are slightly longer than the corresponding Er—O(N) and Y—O(N)distances (Y—Oaxial: 0.221 1(7), 0.223 1(8) nm; Y—O(N)equatorial: 0.230 2(9)-0.255 7(8) nm) (Table S2).It is worth mentioning that theR1andwR2values are much higher for 3Y due to the poor single-crystal quality.We have tried several times to select new crystals for structural analysis and the result was not improved.

2.2 Optical properties

The UV-Vis absorption spectra, translated from the diffuse reflectance spectra using the Kubelka -Munk equationF(R)=(1-R)2/(2R), were recorded in the range of 200-800 nm in the solid state at room temperature (Fig.S4).All showed strong and broad absorption bands between 200-450 nm, corresponding to theπ→π* transition of the anthracene ligand.Thef-ftransitions of LnⅢwere also observed for 1Er, 2Nd and their diluted samples.For 1Er, the peaks at 487, 520, 544 and 651 nm are assigned to thef-ftransitions of ErⅢfrom4I15/2to4F7/2,2H11/2(2G11/2),4S3/2and4F9/2.For 2Nd,the peaks at 511, 525, 582 and 736 nm are ascribed to thef-ftransitions of NdⅢfrom4I9/2to4G9/2,4G7/2,2G7/2(2G5/2) and4F7/2[30].The diluted samples 1Er@Y and 2Nd@Y displayed similar spectra to 1Er and 2Nd,respectively, except that the peak intensities of thef-ftransitions were much weaker.

The photoluminescence (PL) properties of these compounds were further investigated in the solid state at room temperature.As shown in Fig.2a, the PL spectrum of depma exhibited a broad band peaking at 500 nm, attributed to the excimer formation[31].The broad band was also observed in 3Y but the peak was redshifted to 548 nm with a lifetime of 57.8 ns, attributed to the excimer formation of face-to-faceπ-πstacking of the anthracene groups.The quantum yield was 7.48%.Interestingly, the luminescence of 1Er and 2Nd was extremely weak in the range of 400-700 nm and was invisible to the naked eye,unlike their Dy and Y analogues.This is also supported by the undetectable emission life-time and low quantum yields in the visible region(0.43% for 1Er and 0.76% for 2Nd).The result indicates that an efficient energy transfer from the ligand to the lanthanide ion may occur in the cases of 1Er and 2Nd.

Fig.2 (a)Emission spectra for compounds depma,1Er,2Nd and 3Y excited at 365 nm;(b)NIR emission spectra for 1Er,2Nd,1Er@Y and 2Nd@Y excited at 365 nm

Noting that the maximum emission peak of the excimer in 3Y was 548 nm (18 248 cm-1), which was close to the excited states of NdⅢions (2H11/2and4F9/2)and ErⅢions(4F9/2)[15],we expected a favorable intramolecular energy transfer from the excimer of anthracene pairs to ErⅢor NdⅢ, thus sensitizing the luminescence of the lanthanide ions in 1Er and 2Nd.Indeed, upon excitation at 365 nm, 1Er emitted NIR emission at 1 530 nm, assigned to the4I13/2→4I15/2transition of ErⅢ,while 2Nd emitted NIR luminescence at 860 and 1 054 nm, assigned to the4F3/2→4I9/2and4F3/2→4I11/2transitions of NdⅢ,respectively[32-35](Fig.2b).

2.3 Photocycloaddition reaction

To investigate whether the photocycloaddition reaction occurs in complexes 1Er, 2Nd, and 3Y, we irradiated the samples under 395 nm UV light for 30 min.For 1Er and 2Nd, their PL and IR spectra and PXRD patterns showed no significant difference before and after light irradiation (Fig.S6-S8), indicating that the photochemical reaction did not occur.This is in accordance with the fact that an efficient energy transfer from the anthracene ligand to ErⅢor NdⅢion is present in these complexes,leading to the quenching of the cycloaddition reaction.By contrast, 3Y showed a continuous decline of the PL peak intensity at 547 nm upon 365 nm light irradiation, concomitant with the emergence and increase of new peaks at 422, 445 and 485 nm (Fig.S7).These new peaks are assigned to theπ*→πtransitions of the di-anthracene units.After 30 min of irradiation,the PL profile of 3Y did not change.Obviously, the face-to-faceπ-πstacked anthracene groups in 3Y underwent photocycloaddition to form a new phase [Y(SCN)2(NO3)(depma2)(4-hpy)2]n(3Y-UV),like the Dy analogue[18].The formation of dianthracene in 3Y-UV was supported by IR spectra, which showed a new peak at 687 cm-1characteristic of C—H vibration of dianthracene.Meanwhile the peak at 1 240 cm-1characteristic of C—H bending vibration of anthracene was almost disappeared (Fig.S8).The PXRD measurements suggest that 3Y-UV is isostructural to 2Dy but the structural transformation is incomplete(Fig.S9).

To examine the completeness of the photocycloaddition reaction, we measured the1H NMR spectra of 3Y and 3Y-UV in DMSO-d6solution.Compared with 3Y, compound 3Y-UV showed three additional signals with lower chemical shifts at 6.7-7.3 in addition to the four signals at 7.5-8.6(Fig.S10).The latter are assigned to the hydrogen atoms of the anthracene unit, while the former belong to the hydrogen atoms of the dianthracene unit.The photochemical reaction yield can be calculated according to the integral areas of the two sets of signals.For 3Y-UV, the photocycloaddition reaction has a yield of 86.1%.

Interestingly, the diluted samples 1Er@Y and 2Nd@Y can also undergo photocycloaddition reaction under UV light irradiation.Fig.3 shows the PL spectra of 1Er@Y and 2Nd@Y upon 395 nm light irradiation as a function of time.In both cases, the peak intensity at 547 nm decreased with increasing irradiation time and new peaks corresponding to dianthracene emerged at 422, 445 and 485 nm.The IR spectra showed a monotonous decreasing of the peak intensity at 1 240 cm-1and an increase of the peak intensity at 687 cm-1(Fig.S11), confirming the formation of photodimerized products.Samples obtained after 395 nm UV light irradiation of 1Er@Y and 2Nd@Y (20 mg) for 24 h were used for further measurements, named as 1Er@Y-UV and 2Nd@Y-UV.By comparing their PXRD patterns with that of 2Dy, it is clear that 1Er@Y-UV and 2Nd@Y-UV contain unreacted 1Er@Y and 2Nd@Y,respectively(Fig.S12),implying that the photodimerization reaction of the diluted samples was incomplete.The yields were 62.4% for 1Er@Y and 64.0% for 2Nd@Y based on the measurements of1H NMR spectra(Fig.S10).

Fig.3 Time-dependent photoluminescence spectra upon 395 nm UV-light irradiation of 1Er@Y(top)and 2Nd@Y(bottom)

2.4 Magnetic properties

The temperature dependence of magnetic susceptibility was carried out under 1 000 Oe dc field in the temperature range of 2-300 K.TheχMTvalue at 300 K for 1Er and 2Nd were 11.43, 1.47 cm3·K·mol-1,respectively, consistent with the theoretical values of 11.48 cm3·K·mol-1for ErⅢ(3H4,S=1,L=5,g=4/5), and 1.64 cm3·K·mol-1for NdⅢ(4I9/2,S=3/2,L=6,g=8/11).On cooling, theχMTvalue for 1Er decreased gradually in the range of 300-100 K, and then rapidly reached the minimum value of 5.01 cm3·K·mol-1at 2 K (Fig.4).This tendency can be attributed the progressive depopulation of theMJlevels of ErⅢand/or low-lying excited states.For 2Nd, theχMTvalue upon cooling remained almost constant down to 1.8 K.The field dependence of magnetization at 2 K reveals that the magnetizations at 70 kOe were 5.27Nβfor 1Er, which were obviously lower than the theoretical saturation values of 9Nβ(Fig.S13a), suggesting the presence of magnetic anisotropy.For 2Nd,isothermal magnetization curves at 2 K exhibited values of 1.23Nβat 70 kOe,which was close to the expected value of 1.27Nβfor one NdⅢion(Fig.S13b).

Fig.4 (a)Temperature dependence of χMT on cooling in a dc field of 1 kOe for 1Er and 2Nd;(b,c)Frequency dependence of the out-of-phase(χ″)signals for compounds 1Er(b)and 2Nd(c)at the depicted temperatures under a 500 Oe(for 1Er)or 1 kOe(for 2Nd)dc field;(d)Plot of τ vs T on the log-log scale for 2Nd

To probe the slow magnetization relaxation dynamics, ac magnetic susceptibilities were measured for 1Er and 2Nd.There were no out-of-phase (χ″) signals at 2 K for both compounds under zero dc field (Fig.S14), which can be attributed to the fast zero-field quantum tunneling of magnetization (QTM)[36-38].Since the QTM effect can be suppressed by applying an external dc field, ac magnetic data for both compounds were performed under an optimum dc field.For 1Er, theχ″ac susceptibility data exhibited the frequency dependence below 3.2 K under an optimum field of 500 Oe,but no maximum appeared even the frequency reached 1 000 Hz, implying a faster relaxation process in 1Er(Fig.4b and S15).

For 2Nd, both theχ′ andχ″ signals showed obvious frequency and temperature dependence under an optimal 1 kOe dc field (Fig.4c and S15), suggesting a field-induced slow relaxation of magnetization.The Cole-Cole plots can be fitted by the generalized Debye model[39]to extract the relaxation time (τ) (Fig.S16a).The distribution coefficient (α) values fall in the range of 0.23~0.30 (Table S5).We first attempted to fit the relaxation times by Direct and Orbach processes using Eq.(1), obtaining the parametersUeff=18.6(6) K,τ0=2.3(4)×10-6s,A=142(6) K-1·s-1for 2Nd (Fig.S16b).However, the smallUeffand largeτ0values suggest that other relaxation pathways were dominant.

An almost linear dependency was observed when theτvsTwas plotted in a log-log scale (Fig.4d), indicating that dominant relaxation mechanism may be a direct or Raman process.The fitting was fairly well by using the combination of Raman and direct process using Eq.2, giving parameters ofC=4.5(11) K-4.9·s-1,n=4.9(2).

To study the photo-controllable magnetic behaviour, the powdered samples 1Er@Y and 2Nd@Y and corresponding photoreaction products 1Er@Y - UV and 2Nd@Y-UV were used to measure the dc and ac susceptibilities in detail.The dc magnetic behaviours of 1Er@Y and 1Er@Y-UV were similar to those for 1Er (Fig.S17 and S18).Slow magnetization relaxation was also found for 1Er@Y and 1Er@Y-UV under 500 Oe dc field and had a slight difference between them (Fig.5 and S19).In addition, 2Nd@Y and 2Nd@Y-UV also showed similar behaviour to 2Nd(Figs.5 and S20,Tables S6 and S7).The ac data can also be fitted by Eq.2, giving the parametersn=3.8(9), andC=10(2)K-3.8·s-1for 2Nd@Y andn=5.2(8),andC=27(6)K-5.15·s-1for 2Nd@Y-UV.Compared to 2Nd@Y, thenvalue of the Raman process for 2Nd@Y-UV is increased, attributed to the significant change in the photon structure[40], which is reasonable owing to the formation of partial coordination polymer in 2Nd@YUV after photocycloaddition[41-42].

Table 1 The magnetic parameters for 1Er,1Er@Y,1Er@Y-UV,2Nd,2Nd@Y and 2Nd@Y-UV

Fig.5 Frequency dependence of the out-of-phase(χ″)signals at the depicted temperatures under a 500 Oe dc field for 1Er@Y(a),1Er@Y-UV(b)and 1 kOe dc field for 2Nd@Y(c)and 2Nd@Y-UV(d),respectively;Plots of τ vs T for 2Nd@Y(e)and 2Nd@Y-UV(f)

3 Conclusions

In summary, we have synthesized three isostructural lanthanide complexes with formulae Ln(SCN)2(NO3)(depma)2(4-hpy)2(Ln=Er(1Er),Nd(2Nd),Y(3Y)).1Er and 2Nd show a field-induced slow relaxation of the magnetization and near-infrared (NIR) luminescence but fail to undergo photocycloaddition reaction under 395 nm UV light irradiation.In contrast, the diamagnetic complex 3Y shows excimer emission and can undergo photocycloaddition.By doping ErⅢor NdⅢinto the 3Y lattice, we obtained Ln0.1Y0.9(SCN)2(NO3)(depma)2(4-hpy)2(Ln=Er (1Er@Y), Nd (2Nd@Y)),which not only can undergo partial photocycloaddition but also are accompanied by changes in magnetic and optical properties.Further work is in progress to explore new lanthanide systems with tunable magnetooptical properties.

Supporting information is available at http://www.wjhxxb.cn

猜你喜欢
化工学院南京大学鹏飞
我校党委书记柴林一行赴南京大学交流学习
国家开放大学石油和化工学院学习中心列表
《南京大学学报数学半年刊》征稿简则
【链接】国家开放大学石油和化工学院学习中心(第四批)名单
Quality Control for Traditional Medicines - Chinese Crude Drugs
为了避嫌
惩“前”毖“后”
举贤
《化工学报》赞助单位
Comprendre et s'entendre