999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

摻雜對鉺、釹蒽基配合物的光二聚反應及光磁性能調控的影響

2024-01-20 03:56:06麻秀芳秘鵬飛鮑松松鄭麗敏
無機化學學報 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

猜你喜歡
實驗室化學
電競實驗室
電子競技(2020年4期)2020-07-13 09:18:06
電競實驗室
電子競技(2020年2期)2020-04-14 04:40:38
電競實驗室
電子競技(2019年22期)2019-03-07 05:17:26
電競實驗室
電子競技(2019年21期)2019-02-24 06:55:52
電競實驗室
電子競技(2019年20期)2019-02-24 06:55:35
電競實驗室
電子競技(2019年19期)2019-01-16 05:36:09
奇妙的化學
奇妙的化學
奇妙的化學
奇妙的化學
主站蜘蛛池模板: 日韩精品亚洲一区中文字幕| 欧美区日韩区| 伊人色天堂| 福利在线不卡一区| 又黄又爽视频好爽视频| 国产制服丝袜91在线| 8090成人午夜精品| 无码 在线 在线| 国产主播福利在线观看| a级毛片免费播放| 国产在线八区| 第一页亚洲| 国产三级精品三级在线观看| 欧美国产日韩一区二区三区精品影视| 国产精品永久久久久| 高清久久精品亚洲日韩Av| 久久久久亚洲AV成人网站软件| 国产精品亚洲αv天堂无码| 五月天久久婷婷| 一级毛片免费高清视频| 青青青视频蜜桃一区二区| 亚洲三级影院| 国产日韩欧美精品区性色| 国产在线日本| 国产喷水视频| 久久成人国产精品免费软件 | 69视频国产| av尤物免费在线观看| 欧美精品影院| 日韩成人午夜| 欧美成人第一页| 自拍中文字幕| 国产一级毛片网站| 国产呦视频免费视频在线观看| 青青草原偷拍视频| 色丁丁毛片在线观看| 精品中文字幕一区在线| 怡红院美国分院一区二区| 亚洲欧美国产视频| 日韩精品视频久久| 日本人又色又爽的视频| 欧美α片免费观看| 亚洲国产成人麻豆精品| 国产亚洲日韩av在线| 国产一区二区精品福利| 永久免费AⅤ无码网站在线观看| 91最新精品视频发布页| 亚洲欧美精品在线| 国产精品露脸视频| 亚洲精品第五页| 99福利视频导航| 国产系列在线| 久久免费视频播放| 国产99视频在线| 亚洲视频免费播放| 欧美日韩精品一区二区在线线| 色噜噜狠狠狠综合曰曰曰| 久久91精品牛牛| 国产区在线看| 狠狠色丁香婷婷综合| 无码'专区第一页| 亚洲欧美日韩动漫| 国产日韩欧美精品区性色| 久久综合婷婷| 国产在线专区| 影音先锋亚洲无码| 久久毛片免费基地| 国产亚洲成AⅤ人片在线观看| 一本二本三本不卡无码| 福利小视频在线播放| 色综合久久久久8天国| 欧美日韩在线成人| 国产欧美视频综合二区| 在线视频精品一区| 手机精品视频在线观看免费| 在线观看无码a∨| 亚洲侵犯无码网址在线观看| 亚洲男人在线天堂| 国产精品女主播| 暴力调教一区二区三区| 色久综合在线| 538精品在线观看|