• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Excited-state conformation capture by supramolecular chains towards triplet-involved organic emitters

    2021-11-06 03:18:12HoLiuNokiAnoShigehiroYmguchiPnceNumovHongyuZhng
    Chinese Chemical Letters 2021年5期

    Ho Liu,Noki Ano,Shigehiro Ymguchi,c,Pnˇce Numov,Hongyu Zhng,*

    a State Key Laboratory of Supramolecular Structure and Materials, College of Chemistry, Jilin University, Changchun 130012, China

    b Department of Chemistry, Graduate School of Science, and Integrated Research Consortium on Chemical Sciences (IRCCS), Nagoya University, Furo Chikusa,Nagoya 464-8602, Japan

    c Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo Chikusa, Nagoya 464-8602, Japan

    d New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates

    ABSTRACT Nowadays, the development of triplet-involved materials becomes a hot research topic in solid-state luminescence fields.However,the mechanism of triplet-involved emission still remains some mysteries to conquer.Here,we proposed a new concept of excited-state conformation capture for the constructions of different types of triplet-involved materials.Firstly, excited-state conformation could be trapped by supramolecular chains in crystal and form a new optimum excited-state structure which is different from that in solution or simple rigid environment,leading to bright thermally activated delayed fluorescence(TADF) emission.Based on excited-state conformation capture methodology, next, we obtained roomtemperature phosphorescence (RTP) by introducing Br atoms for the enhancement of intersystem crossing.It could be concluded from experimental results that TADF may originate from aggregate effect while RTP may derive from monomers.Finally, heavy-atom free RTP and ultra RTP were achieved by eliminating aggregate effect.This work could not only extend the design methodology of triplet-involved materials but also set clear evidences for the mechanism of triplet-involved emissions.

    Keywords:Room-temperature phosphorescence Supramolecular chains Conformation capture Organic emitters Triplet-involved emission

    Solid-state luminescent materials have been always a desirable cutting-edge research field on account of the potential values in optoelectronic technologies [1-8].The emergence of tripletinvolved emitters thrives many advanced material fields such as electroluminescence, molecular sensors, and time-resolved bioimaging [9-11].Emission mechanism of these materials involves triplet state so that quantum efficiency could surpass the limitation of the efficiency of traditional fluorescent materials and/or materials with afterglow emission which is used to exclude the influence of fluorescence in some fields[12-14].How to design triplet-involved materials is always an important subject in materials research.The design method for constructing Thermally Activated Delayed Fluorescence(TADF)had already been proposed by reducing the energy gap between singlet and triplet states[15].However, the suitable design approaches for Room-Temperature Phosphorescent (RTP) materials are still under developing.

    Initially, the construction of RTP materials has been largely focused on metal-containing inorganic materials on the basis of transition metal complexes [16,17], because metal atoms could promote the spin-orbit coupling(SOC)and boost the intersystem crossing (ISC) [18].In 2011, Kim et al.designed a doped organic crystal with persistent RTP characteristic, demonstrating the importance of heavy atom effect in the generation of RTP [19].The proposal of crystalline-induced phosphorescence (CIP) made an enormous splash on single molecular RTP materials,indicating the imporance of crystalline state [20,21].From then on, several types of crystal structures, including H-aggregation [22-24],intermolecular electronic coupling [25], π???π interactions[26,27],halogen bonds[28],and so on[29,30],have been proposed for illustrating the mechanism of RTP.Even though, the known mechanisms for triplet-involved emitters are still underdeveloped.Further investigations on the mechanism of triple-invovled materials are crucial for molecular design and the fabrication of materials.

    Most of RTP organic crystalline materials could be assigned to crystallization-induced emission, namely emission by the restriction of intramolecular motions (RIM) [20,21].It means emission behavior originates from the limitation of molecular rotation and vibration.Molecules could rotate and vibrate freely in solution,where non-radiative transitions are usually dominant.However,in a rigid environment (for example, frozen solution), molecules cannot rotate and vibrate freely, leading to RTP emission.In this work, we propose a new concept of excited-state conformation capture to construct crystallization-induced emission RTP materials which are different from the aforementioned RTP materials.The optimum excited states of monomers in solution or rigid enviornment have totally distinct conformations when compared with those in crystal structure.A simple restriction of molecular movements is not enough to catch the optimum excited-state conformations with bright RTP emissions.However,because of the strong interactions in crystalline state, excited-state structures would be changed, forming new optimum excited-state conformations with RTP emission behaviors.In other words, supramolecular interactions could capture new optimum excited-state conformations for bright RTP emissions in crystal.The design concept, emission behavior, and intermolecular interactions will be discussed in this work.

    Herein, we reported a series of molecules lack of emissions in solution and rigid environment but have different triplet-involved emission behaviors in crystalline forms.This article will be divided into four steps to demonstrate our diverse design concepts on the basis of excited-state conformation capture.Two crystallineinduced TADF molecules are designed and synthesized (step 1).Those molecules could form supramoelecular chains by strong hydrogen bonds which could capture excited-state conformations.Starting with those chains,efficient RTP(step 2),heavy-atom-free RTP (step 3) and ultra RTP materials (step 4) are attained by seperating aggregate effectsand ehnhancing ISC process.Our work harvests the triplet energies of pure organic materials efficiently by capturing excited-state structures,and investigates the essence of partial crystalline-induced triplet-involved emitters in-depth.

    Step 1: Excited-state conformation capture for TADF.

    Compounds 1a and 1b were obtained by one step synthesis(details see Supporting information)[31].As shown in Fig.1a,the crystals 1a and 1b demonstrate bright emissions with peak maxima at 427 and 463 nm,respectively.The decay spectra further revealed that the crystals have long-lived luminescent behaviors besides prompt luminescence(Figs.S1a-d in Supporting information).In order to confirm the origin of such emission component,temperature-dependent decay spectroscopy (Figs.1b and c) was performed.The ratio of delayed component increased with the temperature rising from 250 K to 330 K, which determined the behavior of thermal activation energy for TADF.For further investigation of TADF emission behavior, monomer’s photophysical properties of 1 were employed.It should be noted that 1a and 1b are all non-emissive in deoxygenated solution, on the plates of thin-layer chromatography (TLC plate, aggregate state),and on the spin-coated films (Figs.S2—S4 in Supporting information), which is totally different from tranditional TADF emitters reported in the literatures, indicating the existance of triplet-involved crystallization-induced emission(CIE)behavior.In general,the solutions of triplet-involved molecules emit phosphorescence under frozen conditions on account of RIM [20,21].In contrast, at 77 K, the solutions of 1a and 1b could not emit luminescence (Fig.S5 in Supporting information).Such crucial difference reveals that molecular structures and the rigidity of environment are not the most important factors for the TADF properties and there must be some other impacts in crystal structures which may be momentous for the observed TADF behavior.

    Fig.1.Emission spectra of crystals 1,insert:photograph of 1a(upper)and 1b(lower)under 365 UV light(a);temperature-dependent decay spectra of crystals 1a(b)and 1b(c)measured at 420 and 460 nm,respectively;adiabatic ΔEst between optimum singlet excited state and optimum triplet excited state in crystal structure(d);excited-state conformations of 1a and 2b calculated in THF(e);excited-state conformation(red)were calculated by freezing the upper and lower molecules(green)based on S1 state(f).

    For the sake of investigation on such significant emission between monomer and crystal, single-crystal X-ray diffraction analysis was performed on 1a and 1b.Both crystals adopt the same crystal systems (monoclinic) and space groups (P21/c).Molecules 1a and 1b take molecular conformations with dihedral angles of 76.66°and 64.45°between the central phenyl ring and the fluorophenyl ring, respectively (Figs.S7a and b in Supporting information).

    Each molecule in both 1a and 1b connects others through four robust N--H???O--C interactions (2.10 ? for 1a and 2.27 ? for 1b),forming infinite supramolecular chains (Figs.2a and b).These chains connect with other chains through weak C--H???π (2.60 ? for 1a and 2.65 ? for 1b)interactions(Figs.S7c and d in Supporting information).In the structures of these two crystals, the main packing modes are supramolecular chains formed by hydrogen bonds.This unique chain-like structure may be essential for triplet-involved emission.

    To verify that supramolecular chain is the key factor of tripletinvolved emission, compounds 2a and 2b were synthesized in which two NH moieties were replaced with O atoms to eliminate the supramolecular chains.The monomers of 2 could not emit luminescence in solution, on the plates of thin-layer chromatography (TLC plate, aggregate state), and on the spin-coated films,even at a extreme low temperature (77 K) (Figs.S2—S5).Interestingly, in crystalline forms, the crystals 2 still could not emit luminescence (Fig.S8 in Supporting information), ensuring that the rigidity of crystals is not the main reason for TADF emission behavior of 1.For a further comparison,single-crystal Xray diffractions of 2 were investigated in comparison with those of 1.For 2,they inherit similar chain-like packing modes,but without intermolecular interactions in chains (Figs.2c and d).The lack of intermolecular interactions in packing modes prevents the formation of supramolecular chains.Considering such important differences in the crystal structures of 1 and 2, it could be concluded that the unjoined structures of 2a and 2b handicaps the formation of triplet-involved emissions,and that the formation of supramolecular chains must be crucial for triplet-involved emissions.

    The above experimental results proved the key role of supramolecular chains in emission behaviour.However, two problems still exist: (1) How do supramolecular chains influence the emission behavior between monomers(in solution or in rigid environment) and molecules in crystal? (2) What is the origin of TADF emission?

    In order to solve these two problems, theoretical calculations using time-dependent density functional theory (TD-DFT) in solution and crystal state were performed by Gaussian 16 program[32].Taking 1a and 2b as examples (Fig.1e), in solution, excitedstate conformations are significantly twisted because of the lack of sufficient constraint in molecular conformation, leading to a decreation of molecular energy and ensuring a spin-forbidden transition process(oscillator strength of single state=0).Thus such twist excited-state structure leads to a non-radiative transition process which prohibits the emission of monomers.In order to investigate emission behaviors in crystal state, a segment of a supramoleculat chain with three molecules were employed to provide similar crystal environment in which upper and lower molecules are frozen to provide interactions, while the central molecule could act freely(Fig.1f).The calculation results were all without imaginary frequencies which means that real local minima were reached.In crystal structure, molecules 1a locked by strong hydrogrn bonds in supramolecular chain could not rotate or vibrate frequently, forming a new type of exicted state conformation.In crystals 2, molecules could tranform to the twisting conformation as in solution because of the absence of interactions in packing mode, though crystal has a rigid enviroment.Such evidences proved that the supramolecular chains influence emission behaviors by capturing a new optimum excited state conformation, which leads to the emission in crystalline form.

    As for the origin of TADF, the adiabatic energy gap in crystal structure between the optimum singlet excited state and the optimum triplet exicited state was calculated to be 0.4 eV,which is suitable for TADF emission (Fig.1d) [15].Simultaneously, the energy gap could be reduced with the extension of chain lengths based on the crystal structure (Table S2 in Supporting information).

    Fig.2.Supramolecular chains of 1a(a),1b(b),forming by hydrogen bonds.Packing modes of 2a (c), and 2b (d).

    In this part, a new optimum excited-state conformation was obtained by using supramolecular chains as a driving force of excited-state conformation capture technique.In order to utilize such method sufficiently, we tend to design materials with RTP emission and investigate the mechanism of RTP.

    Step 2: Spin-orbit coupling enhancement for RTP.

    TADF emitters could be required by excited-state conformation capture.How to foster RTP materials is the main topic in this step.Heavy atoms (for example, Br atom) could enhance SOC and accelerate ISC process between singlet and triplet states,thus they were constantly selected to attain RTP materials[33,34].Similarly,compound 3 was designed by introducing two Br atoms on the central phenyl ring to enhance SOC process.Compound 3 does not emit in frozen solution (77 K), on TLC plate, and on spin-coating film (Figs.S3—S5), implying that the molecular structure and the rigidity of environment are not the most crucial factors for emission.However, the crystals obtained by sublimation demenstrate green emission.The emission maximum was located at 501 nm as shown in Fig.3a.Temperature-dependent decay spectra indicated that the emission of 3 could be assigned to RTP(Fig.3b).

    The space group of 3 is P-1, confirmed by Single crystal X-ray diffraction analysis.Each molecule connecting with others through four C=O???H--N hydrogen bonds(2.08 ?)forms the supramolecular chains (Fig.3e).The Br atoms in one molecule could provide efficient external heavy atom effects due to the appropriate distances (3.66 ?) and tremendous overlap between Br atoms in one molecule and phenyl rings in other molecules (Fig.S10 in Supporting information).

    Additionally, a pseudo-polymorph (3′) of compound 3 was obtained.Crystal 3′was grown in DMF at a chilling temperature.In contrast to the bright emission of crystal 3, 3′is non-emissive under UV irradiation though it possesses suitable crystalline morphology(Fig.S11 in Supporting information).The unit cell of 3′contains one molecule 3 and one DMF molecule.The isolated molecule forms N--H???O--C interactions(2.09 ?)with two other DMF molecules, prohibiting the construction of supramolecular chains(Fig.3d).The hydrogen bonds,which are formed with DMF molecules,could not lock excited-state conformation,resulting in free rotation and vibration of molecule 3.Such results further proved the importance of supramolecular chain and its capability in excited-state conformations capturing.

    The theoretical calculations of 3 were carried out in the same way with 1a and 2b.Single molecule was calculated by TD-DFT and trimer state was calculated by UDFT method.In solution, the excited-state conformation is also a twisting structure with the oscillator strength of single state=0 (Fig.S9b in Supporting information).However, in crystalline state, excited-state conformation could be captured by supramolecular chains as an optimum structure with radiative transition (Fig.3f).Excited-state conformation capture method constitutes the cornerstone of RTP emission behavior of 3.For a further investigation of the difference of emission behavior between 1 and 3, natural transition oribital(NTO) analysis was performed.As shown in Fig.3g, Figs.S9c and S9d (Supporting information), singlet excited-state transition forms of 1 and 3 both possess aggregate effect while triplet excited-state transition forms are both attributed to the influence of monomers whose conformations are captured by supramolecular chains.Similar transitions mean transition forms are not the momentous factors for RTP emissions.Hence, distinct emission behavior could be assigned to the existence of Br atoms.Br atoms could enhance ISC process, which is proved by the large SOC constants shown in Fig.3c and Fig.S13a(Supporting information).SOC constants were calculated by ORCA [35,36], and corresponding diagrams were drawn by VMD and Multiwfn [37,38].Most of the SOC constants of monomers are larger than those of trimers,indicating that RTP should be mainly attributed to the influence of monomers.Br atoms could facilitate the RTP emissions by reducing reverse intersystem crossing(RISC)process as well as accelerating ISC process.

    Fig.3.Emission spectrum of crystal 3, insert: photograph of 3 under 365 UV light (a); temperature-dependent decay spectra of 3 (b); energy level diagram and the corresponding SOC constants of the monomer of 3 (c); crystal structure of 3′ (d); supramolecular chain of 3 (e); excited-state conformation (red color) was calculated by freezing the upper and lower molecules (green color) based on T1 state (f); NTO isosurfaces of the trimer of 3 with different transition types (g).

    Based on the NTOs results, it could be proposed that TADF emissions may originate from aggregate state while RTP emissions may derive from the monomer's excited-state conformation captured by supramolecular chains.Therefore, in the next step,we manage to weaken the hydrogen bonds in supramolecular chains to impair the status of aggregate state and the construction of heavy-atom-free RTP materials.

    Step 3:Elimination of aggregate effect for heavy-atom-free RTP.

    To obtain efficient and heavy-atom-free RTP materials, the influence of aggregate effect is intended to undermine and thus the effect of monomer's excited-state conformation capture would be enhanced.The hydrogen atoms bonding to the nitrogen atoms are replaced with phenyl rings in order to reduce the aggregate state.This design concept would enlarge the distance of two neighbouring molecules,leads to the weakening of intermolecular-aggregation effect.Simultaneouly, the extension of intermolecular distances may impair the hydrogen bonds strength, which may also handicap the transition process between molecules.Based on these considerations,compound 4 was designed and synthesized.

    The emission behavior of 4 in solution is identical to compouds 1-3.Compound 4 does not emit in solution (77 K), on TLC plate,and on spin-coating film (Figs.S2—S5).In crystalline state, 4 exhibits intensive emission in air with the peak maxima at 512 nm(Fig.4a).Temperature-dependent decay spectra from 77 K to 328 K confirmed that the emission of 4 is undoubtedly ascribed to RTP(Fig.4b).Space group of compound 4 is P-1.As shown in Fig.4d,although H atoms in N--H moieties are replaced with phenyl rings,packing mode is quite similar to those of 1-3(molecules are locked by four hydrogen bonds).Each molecule connects with adjacent molecules through four C=O???H hydrogen bonds, by which supramolecular chains are formed.The phenyl rings forming hydrogen bonds with carbonyl groups are almost parallel to the hydrogen bonds,which indicates that the replacement with phenyl groups does not disturb the formation of supramolecular chains.Compared with those of 1-3, the length of hydrogen bonds of 4 increases to 2.71 ?, indicating an extension of intermolecular distance.Both of the experimental results and single crystal structure analysis prove the reliability of our hyothesis that TADF emissions originate from aggregate states while RTP emissions derive from the monomer's excited-state conformation captured by supramolecular chains.

    Fig.4.Emission spectrum of crystal 4.Insert:photograph of 4 under 365 UV light(a); temperature-dependent decay spectra of 4 (b); energy level diagram and the corresponding SOC constants of the monomer of 4 (c); supramolecular chain of 4(d);excited-state conformation(red color)was calculated by freezing the upper and lower molecules(green color)based on T1 state(e);NTO isosurfaces of the trimer of 3 with different transition types (f).

    Theoretical calculations were carried out in the same way of 3.In crystalline state (Fig.4e), the excited-state conformation could still be captured bystronghydro genbonds in supramolecularcha ins,forming the cornerstone of triplet-involved emissions.To investigate the origin of RTP,ISC process was investigated at first.As can be seen in SOC constants(Fig.4c),compared with 3(Fig.3c),SOC process is not enhanced dramatically,indicating that RTP could be caused by otherreasons.NTOs indicatethat the transitiontypes bothofsinglet excited state and triplet excited state convert to monomer’s forms,demonstrating the weakening of aggregate state which is caused by the extension of intermolecular distances and the reduction of hydrogen bonds(Fig.4f).Besides,the energy gaps between singlet and triplet state are 0.62 eV for monomer and 0.68 eV for trimer,respectively, which prohibit RISC process.All above mentioned facts contribute in the promoting of the RTP emission.

    It is concluded from this step that RTP emission originates from monomer’s excited-state confomration capture.In the next step,we intend to design ultralong RTP by summarising the concepts of Steps 1-3.

    Step 4: Stabilization of triplet state for URTP.

    In Step 1,we figured out that excited-state conformation could be captured by supramolecular chains, and constructs tripletinvolved materials.In Step 2, the introduction of Br atoms enhances spin-orbit coupling and facilitates the ISC process.In Step 3,heavy-atom-free RTP material was obtained by eliminating aggregate effect and achieving monomer’s excited-state confomration capture.Based on these results,we tend to design ultralong RTP(URTP)materials by using other interactions as implements for supramolecular chains.

    Four conditions should be satisfied in order to obtain URTP materials by supramolecular chains: (1) excited-state conformation capture, (2) monomer’s transition formation, (3) enhancement of ISC process, and (4) stability of triplet excitions.Hence,π???π stacking which is widely used to stablise triplet excitions[22,23] and Br atom which is usually introduced for an enhancement of ISC process [33-35] were employed in the construction of URTP material.In this step, π???π stacking could also serve as an implement of excited-state conformation capture for the formation of supramolecular chains.At the same time, it was found that the range of π???π stacking is around 3.50 ?(face to face), which could undermine molecular aggregation.Therefore,compound 5 was designed and synthesized.

    Crystals of compound 5 were obtained with good crystalline morphology from methanol/DCM solution.In crystalline state,compound 5 displays fluorescence peaking at 430 nm and excellent RTP with the emission maximum at 500 nm (Fig.5a).The decay spectrum of RTP component is shown in Fig.5b,and the lifetime of 5 was 99 ms,which is indicative of a good performance of RTP characteristic.In order to find the origin of RTP emission,the emission behavior of monomers was investigated.The solution of 5 at room temperature is non-emissive (Fig.S2).However, under frozen solution, phophorescent emission could be achieved,though intensity is feebly detected by naked eyes (Fig.S6b in Supporting information).Such weak emission with peak maximum around 500 nm proved that the origin of RTP should be attributed to the influence of monomers, but such restriction of molecular movements is not enough for RTP emission as in crystals.The interactions in crystals must play key role in the RTP behavior.

    In the crystal structure, the unit cell possesses two molecules.These two molecules form supramolecular chains through π???π packing.The distances between adjacent phenyl rings are 3.57 ?(Fig.5e) and 3.62 ? (Fig.S12 in Supporting information), and the centroid-centroid distances of adjacent phenyl rings are both 3.98 ?.It could be concluded that, compared with the hydrogen bonds in 1-4, π???π stacking interaction increases the distances among molecules,which reduce the influence of aggregate effects.Besides,the existence of π???π stacking may not only capture a new optimum excited-state conformation but also stabilize triplet excitions while Br atom could benefit the ISC process, both of which are responsible for the URTP.Such crystal structure and emission under frozen conditions preliminarily proved the efficacy of our design methods mentioned previously.

    Theoretical calculation were performed on one of the conformations in crystal form as a further evidence.The excited-state conformation is captured by π???π stacking(Fig.5f),while Br atom benefits the ISC process as shown by SOC constants in Fig.5c.NTO analysis (Fig.5g) showed that, because of the extension of intermolecular distances by π???π stacking, transition forms of both singlet and triplet excited states are located on a monomer,eliminating aggregate effect.

    Fig.5.Emission spectra of crystals 5,insert:photograph of 5 under daylight and with afterglow(a);decay spectrum of 5(b);energy level diagram and the corresponding SOC constants of the monomer of 5 (c); time-resolved indentification pattern of 5 (d); supramolecular chain of 5 (e); excited-state conformation (red color) was calculated by freezing the upper and lower molecules (green color) based on T1 state (f); NTO isosurfaces of the trimer of 5 with different transition types (g).

    In the summary of this step, (1) the introducation of π???π stacking, which is an implement for excited-state conformation capture in this step,could more effectively restrict the movements of molecules than a simple rigid environment for RTP emission;(2)π???π stacking eliminates aggregate effect by extending the distance among molecules, which may hinder the RISC process,and stabilizes triplet excitons; (3) Br atom benefits spin-orbit coupling and enhance intersystem crossing process.It could be concluded that RTP emission behavior should originate from monomer’s excited-state conformation capture.

    The feature of URTP is persistent emission after the removal of excitation sources.Consequently, the URTP materials generally function as color-encoding or secure patterns.As shown in Fig.5d,we demonstrated the potential application of 5 in time-resolved identification as an illustration.The pattern was consisted of two components.The moon was constructed by 5 with URTP and the stars were comprised of 1a.The moon and stars displayed cyan emission under 365 UV light.While only the moon exhibited pure green emission after the removal of UV light and the emission could be detected by naked eyes, indicating a marvelous application in secure protections.

    In this work, a new concept of excited-state conformation capture by supramolecular chains was proposed for the explanation of mechanism of triplet-involved emissions and for the construction of different types of triplet-involved emitters(TADF,RTP,heavy-atom-free RTP,and URTP).First of all,we designed two molecules with crystalline-induced TADF.It could be summarised that, different from the optimum excited-state structure in solution, excited-state conformation in crystal state is captured by supramolecular chain, forming a new optimum sturcture with bright emission.Next, bright RTP emitter was designed by introducing Br atoms based on the similar supramolecular chains.The introduction of Br atoms enhances spin-orbit coulping and improves ISC process.From a comparison of TADF and RTP emitters, it could be proposed that TADF may originate from aggregate effect.Then aggregate effect was eradicated by extending the distances between molecules in supramolecular chains for heavy-atom-free RTP.Finally, π???π stacking was employed as an interaction in supramolecular chains for excited-state conformation capture and for the elimination of aggregate effects, constructing URTP materials.There methods not only enrich the mechanism of triplet-involved materials but also provide potient ways, which intramolecular interactions are used to capture excited-state conformations, for constructing amorphous phosphorescent organic materials.Further investigations are undergoing.

    Declaration of competing interest

    The authors report no declarations of interest.

    Acknowledgment

    This work was supported by the National Natural Science Foundation of China (No.51773077).

    Appendix A.Supplementary data

    Supplementary material related to this article can be found,inthe online version,at doi:https://doi.org/10.1016/j.cclet.2020.12.013.

    老汉色av国产亚洲站长工具| 久久国产精品影院| 亚洲性夜色夜夜综合| 亚洲中文日韩欧美视频| 成人手机av| 人人妻人人爽人人添夜夜欢视频| 老熟妇乱子伦视频在线观看 | 最近中文字幕2019免费版| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 精品一品国产午夜福利视频| 91精品国产国语对白视频| 久久久久久久久久久久大奶| 青春草视频在线免费观看| 色老头精品视频在线观看| 999久久久精品免费观看国产| 欧美黄色片欧美黄色片| 亚洲国产精品成人久久小说| 一本色道久久久久久精品综合| 男女下面插进去视频免费观看| 丁香六月天网| 黄色视频在线播放观看不卡| 青春草亚洲视频在线观看| 亚洲色图 男人天堂 中文字幕| 国产伦理片在线播放av一区| 又紧又爽又黄一区二区| 国产成人a∨麻豆精品| 久久人人爽av亚洲精品天堂| 色老头精品视频在线观看| kizo精华| 欧美久久黑人一区二区| 成年人黄色毛片网站| 国产亚洲精品久久久久5区| 欧美变态另类bdsm刘玥| 亚洲精品在线美女| 91九色精品人成在线观看| 欧美另类一区| 日韩制服骚丝袜av| 国产精品自产拍在线观看55亚洲 | 亚洲精品美女久久久久99蜜臀| 国产精品一二三区在线看| 啦啦啦在线免费观看视频4| 在线观看免费午夜福利视频| 欧美激情高清一区二区三区| 午夜福利,免费看| 亚洲精品中文字幕在线视频| 美国免费a级毛片| 亚洲精品国产区一区二| 99热国产这里只有精品6| 最近最新免费中文字幕在线| 日韩三级视频一区二区三区| 久久亚洲精品不卡| 国产精品一二三区在线看| 婷婷色av中文字幕| 欧美精品啪啪一区二区三区 | av国产精品久久久久影院| 久久久久国产一级毛片高清牌| 1024视频免费在线观看| 19禁男女啪啪无遮挡网站| 亚洲成人免费av在线播放| 久久久精品区二区三区| 热99国产精品久久久久久7| 国产精品欧美亚洲77777| 在线观看免费午夜福利视频| 免费看十八禁软件| 久久女婷五月综合色啪小说| 精品一品国产午夜福利视频| 国产日韩欧美亚洲二区| 国产一级毛片在线| 脱女人内裤的视频| 欧美中文综合在线视频| 久久精品久久久久久噜噜老黄| 极品人妻少妇av视频| 中文字幕另类日韩欧美亚洲嫩草| 极品少妇高潮喷水抽搐| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区字幕在线| 国产精品一区二区免费欧美 | 久久久久久免费高清国产稀缺| 欧美人与性动交α欧美软件| 制服人妻中文乱码| 国产一区二区在线观看av| 亚洲五月婷婷丁香| 交换朋友夫妻互换小说| 一二三四社区在线视频社区8| 免费黄频网站在线观看国产| 婷婷成人精品国产| 一级黄色大片毛片| 欧美精品啪啪一区二区三区 | 国产在视频线精品| 搡老熟女国产l中国老女人| 老司机福利观看| 91国产中文字幕| 亚洲九九香蕉| 久久天堂一区二区三区四区| 国产高清视频在线播放一区 | 亚洲精品美女久久av网站| 亚洲精品一区蜜桃| www日本在线高清视频| av天堂久久9| 大香蕉久久成人网| 中文字幕人妻丝袜制服| 99re6热这里在线精品视频| 亚洲国产成人一精品久久久| 一本久久精品| 色综合欧美亚洲国产小说| 国产黄频视频在线观看| av有码第一页| 久久精品熟女亚洲av麻豆精品| 国产欧美日韩一区二区三区在线| 咕卡用的链子| 久久青草综合色| 最新在线观看一区二区三区| av网站免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 国产一区二区三区综合在线观看| 精品视频人人做人人爽| 男人舔女人的私密视频| 99久久综合免费| 国产国语露脸激情在线看| 成人手机av| 成人黄色视频免费在线看| 69av精品久久久久久 | 成年动漫av网址| 午夜福利视频在线观看免费| 丝袜在线中文字幕| 久久国产精品影院| 国产精品1区2区在线观看. | 9191精品国产免费久久| 免费久久久久久久精品成人欧美视频| 99热国产这里只有精品6| 亚洲精品久久久久久婷婷小说| 精品人妻熟女毛片av久久网站| 久久精品国产a三级三级三级| 欧美黄色淫秽网站| 久久免费观看电影| 欧美日本中文国产一区发布| 国产不卡av网站在线观看| 又紧又爽又黄一区二区| 精品国产国语对白av| 色视频在线一区二区三区| 法律面前人人平等表现在哪些方面 | 我的亚洲天堂| 99热网站在线观看| 成年人黄色毛片网站| 母亲3免费完整高清在线观看| 日韩三级视频一区二区三区| 美女午夜性视频免费| 精品人妻1区二区| 飞空精品影院首页| 一级黄色大片毛片| 男人爽女人下面视频在线观看| 国产成人免费无遮挡视频| 精品欧美一区二区三区在线| 久久精品久久久久久噜噜老黄| 中国美女看黄片| 国产成人精品无人区| 久久中文字幕一级| 久久中文看片网| 精品一品国产午夜福利视频| 黄色视频在线播放观看不卡| √禁漫天堂资源中文www| 亚洲av男天堂| 美女脱内裤让男人舔精品视频| 最近中文字幕2019免费版| 人人妻,人人澡人人爽秒播| 又紧又爽又黄一区二区| 久久久久网色| 在线十欧美十亚洲十日本专区| 中文字幕高清在线视频| 中文字幕最新亚洲高清| 黑人猛操日本美女一级片| 成人三级做爰电影| 老司机午夜十八禁免费视频| 黄片大片在线免费观看| a级毛片在线看网站| 91av网站免费观看| 99re6热这里在线精品视频| 深夜精品福利| 国产免费av片在线观看野外av| 搡老岳熟女国产| 成人国产一区最新在线观看| 麻豆乱淫一区二区| 丝袜脚勾引网站| 欧美国产精品va在线观看不卡| 搡老熟女国产l中国老女人| 少妇猛男粗大的猛烈进出视频| 91精品国产国语对白视频| 亚洲情色 制服丝袜| 精品少妇久久久久久888优播| 男女国产视频网站| 国产一区二区 视频在线| 伦理电影免费视频| 伦理电影免费视频| 免费av中文字幕在线| 国产激情久久老熟女| 国产亚洲午夜精品一区二区久久| 国产男女内射视频| 欧美日韩亚洲国产一区二区在线观看 | 免费高清在线观看日韩| 视频区欧美日本亚洲| 国产精品二区激情视频| 亚洲视频免费观看视频| 一本一本久久a久久精品综合妖精| 久久中文看片网| 人妻一区二区av| 亚洲精品粉嫩美女一区| 五月开心婷婷网| 中文字幕人妻丝袜一区二区| 国产伦理片在线播放av一区| 日韩欧美一区二区三区在线观看 | 国产真人三级小视频在线观看| 久久久国产欧美日韩av| 精品人妻一区二区三区麻豆| 嫩草影视91久久| 免费少妇av软件| 美女中出高潮动态图| 国产欧美亚洲国产| 国产精品自产拍在线观看55亚洲 | 黄色毛片三级朝国网站| xxxhd国产人妻xxx| 啦啦啦在线免费观看视频4| 国产成人a∨麻豆精品| 免费在线观看视频国产中文字幕亚洲 | 丝袜脚勾引网站| 久久国产精品影院| 久久久久国产精品人妻一区二区| 一边摸一边抽搐一进一出视频| 亚洲欧美激情在线| 不卡一级毛片| 国产高清视频在线播放一区 | 在线av久久热| 我的亚洲天堂| 亚洲色图综合在线观看| 国产成人精品久久二区二区免费| 人妻人人澡人人爽人人| 在线观看人妻少妇| 激情视频va一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 老司机靠b影院| 999久久久国产精品视频| 女性被躁到高潮视频| 搡老熟女国产l中国老女人| av免费在线观看网站| 天天躁夜夜躁狠狠躁躁| 精品亚洲乱码少妇综合久久| av在线老鸭窝| av欧美777| www.999成人在线观看| 亚洲成av片中文字幕在线观看| 亚洲国产日韩一区二区| 免费观看人在逋| 纵有疾风起免费观看全集完整版| 大香蕉久久成人网| 黄色片一级片一级黄色片| 国产区一区二久久| av一本久久久久| 窝窝影院91人妻| 久久精品国产综合久久久| 男人爽女人下面视频在线观看| 9热在线视频观看99| 亚洲精品中文字幕在线视频| 亚洲国产毛片av蜜桃av| 99国产精品免费福利视频| 国产伦人伦偷精品视频| 久久天躁狠狠躁夜夜2o2o| 久9热在线精品视频| 18禁裸乳无遮挡动漫免费视频| 欧美久久黑人一区二区| 人妻人人澡人人爽人人| 在线观看www视频免费| 丁香六月天网| 欧美激情 高清一区二区三区| 99国产综合亚洲精品| 亚洲人成电影观看| 婷婷丁香在线五月| 黄片播放在线免费| 五月天丁香电影| 午夜视频精品福利| 午夜福利一区二区在线看| 午夜福利在线免费观看网站| 十八禁人妻一区二区| 桃花免费在线播放| 人人澡人人妻人| www日本在线高清视频| 12—13女人毛片做爰片一| 欧美午夜高清在线| 久久精品人人爽人人爽视色| 亚洲av欧美aⅴ国产| 18禁黄网站禁片午夜丰满| 亚洲国产精品一区三区| 中文字幕精品免费在线观看视频| 国产日韩欧美视频二区| 欧美精品高潮呻吟av久久| 国产精品二区激情视频| 天天操日日干夜夜撸| 日本撒尿小便嘘嘘汇集6| 久久人人爽av亚洲精品天堂| 欧美国产精品va在线观看不卡| 自线自在国产av| 捣出白浆h1v1| 麻豆av在线久日| 两人在一起打扑克的视频| 精品国产乱子伦一区二区三区 | 黄色视频不卡| a级毛片在线看网站| 可以免费在线观看a视频的电影网站| 精品国内亚洲2022精品成人 | 日韩电影二区| 亚洲精品粉嫩美女一区| 久久天躁狠狠躁夜夜2o2o| 久久久国产一区二区| 午夜福利影视在线免费观看| 丝袜在线中文字幕| 美女中出高潮动态图| av有码第一页| 欧美成狂野欧美在线观看| 色婷婷av一区二区三区视频| 久久久久久久久免费视频了| a级毛片黄视频| 国产精品影院久久| 亚洲成av片中文字幕在线观看| 天天操日日干夜夜撸| 亚洲国产av新网站| 日本av手机在线免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 久久亚洲国产成人精品v| 母亲3免费完整高清在线观看| 99热网站在线观看| 淫妇啪啪啪对白视频 | 精品人妻一区二区三区麻豆| 五月开心婷婷网| 亚洲美女黄色视频免费看| www.熟女人妻精品国产| 亚洲精品国产av成人精品| 天堂俺去俺来也www色官网| h视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 极品人妻少妇av视频| 两个人看的免费小视频| 老熟妇乱子伦视频在线观看 | 岛国毛片在线播放| 日韩大码丰满熟妇| 成人国语在线视频| av一本久久久久| 亚洲 欧美一区二区三区| 777久久人妻少妇嫩草av网站| 成年人免费黄色播放视频| 纵有疾风起免费观看全集完整版| 中文字幕人妻熟女乱码| 18禁裸乳无遮挡动漫免费视频| 日本av手机在线免费观看| 一边摸一边做爽爽视频免费| 成年动漫av网址| 国产一区二区三区在线臀色熟女 | 美女大奶头黄色视频| 桃花免费在线播放| 热re99久久国产66热| 亚洲 国产 在线| 亚洲第一青青草原| 亚洲一区二区三区欧美精品| 久久久国产一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 19禁男女啪啪无遮挡网站| 男女下面插进去视频免费观看| 少妇的丰满在线观看| 久久精品国产综合久久久| 另类精品久久| 久久精品国产亚洲av香蕉五月 | 久久久精品免费免费高清| 超碰成人久久| 99久久国产精品久久久| 一二三四社区在线视频社区8| 国产欧美日韩精品亚洲av| 麻豆av在线久日| 多毛熟女@视频| 女性被躁到高潮视频| 男人添女人高潮全过程视频| 免费观看av网站的网址| 久久久精品国产亚洲av高清涩受| 丁香六月欧美| 亚洲精品中文字幕在线视频| 精品人妻熟女毛片av久久网站| 免费人妻精品一区二区三区视频| 正在播放国产对白刺激| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 亚洲精品自拍成人| 一级毛片电影观看| 亚洲五月色婷婷综合| 国产淫语在线视频| 精品国产乱码久久久久久男人| 操美女的视频在线观看| 精品一区二区三区av网在线观看 | 亚洲精品国产一区二区精华液| 成年人午夜在线观看视频| 男女午夜视频在线观看| 成在线人永久免费视频| 成人国产av品久久久| 国产欧美日韩综合在线一区二区| 日日爽夜夜爽网站| 极品少妇高潮喷水抽搐| 欧美xxⅹ黑人| 亚洲欧洲日产国产| 又紧又爽又黄一区二区| 午夜福利,免费看| 美女主播在线视频| 91九色精品人成在线观看| 新久久久久国产一级毛片| 亚洲欧美一区二区三区黑人| 精品久久蜜臀av无| 女人爽到高潮嗷嗷叫在线视频| 女人高潮潮喷娇喘18禁视频| 老司机影院毛片| 色婷婷av一区二区三区视频| 窝窝影院91人妻| 婷婷成人精品国产| www.av在线官网国产| 18禁黄网站禁片午夜丰满| 丝袜美足系列| 亚洲精品一区蜜桃| 日韩欧美免费精品| 一本—道久久a久久精品蜜桃钙片| 9191精品国产免费久久| 国产日韩欧美亚洲二区| 天堂中文最新版在线下载| 国产亚洲欧美精品永久| 国产精品久久久久成人av| 亚洲午夜精品一区,二区,三区| 黄片小视频在线播放| 操出白浆在线播放| 欧美人与性动交α欧美软件| 午夜免费观看性视频| a在线观看视频网站| av国产精品久久久久影院| 国产亚洲av片在线观看秒播厂| 香蕉丝袜av| 色综合欧美亚洲国产小说| 下体分泌物呈黄色| 18禁国产床啪视频网站| 久久久国产成人免费| 一区二区三区乱码不卡18| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 搡老熟女国产l中国老女人| 精品国产一区二区三区四区第35| 永久免费av网站大全| 久久热在线av| 成人国产一区最新在线观看| av有码第一页| 99热网站在线观看| 成年人午夜在线观看视频| kizo精华| 在线观看免费午夜福利视频| 亚洲成人手机| 男女午夜视频在线观看| 成人免费观看视频高清| 大香蕉久久成人网| 精品卡一卡二卡四卡免费| av欧美777| 免费在线观看黄色视频的| 五月天丁香电影| 欧美老熟妇乱子伦牲交| 国产一区二区在线观看av| 九色亚洲精品在线播放| 美女高潮到喷水免费观看| 老司机福利观看| 欧美日韩亚洲高清精品| 最近最新中文字幕大全免费视频| 制服人妻中文乱码| 97在线人人人人妻| 亚洲精品日韩在线中文字幕| 亚洲av成人一区二区三| av欧美777| 一区二区av电影网| 欧美日韩精品网址| 成年人免费黄色播放视频| 久久久久视频综合| 久久久久久免费高清国产稀缺| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| 狠狠精品人妻久久久久久综合| 欧美在线一区亚洲| 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 免费高清在线观看视频在线观看| 热re99久久国产66热| 两个人免费观看高清视频| 日韩有码中文字幕| 水蜜桃什么品种好| www.av在线官网国产| 亚洲欧洲日产国产| 999久久久国产精品视频| 精品福利永久在线观看| 亚洲色图综合在线观看| avwww免费| 久久这里只有精品19| 一级a爱视频在线免费观看| 国产免费现黄频在线看| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 啪啪无遮挡十八禁网站| 啦啦啦视频在线资源免费观看| 精品一品国产午夜福利视频| 别揉我奶头~嗯~啊~动态视频 | 中文字幕人妻丝袜一区二区| 亚洲精品乱久久久久久| 不卡一级毛片| 亚洲av日韩在线播放| 日日夜夜操网爽| 久热这里只有精品99| 日本一区二区免费在线视频| 国产精品久久久久久精品电影小说| 国产三级黄色录像| 多毛熟女@视频| 亚洲少妇的诱惑av| 一区二区三区精品91| 国产精品熟女久久久久浪| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 少妇精品久久久久久久| 成人黄色视频免费在线看| av免费在线观看网站| 最近中文字幕2019免费版| 欧美另类一区| 波多野结衣av一区二区av| 91国产中文字幕| 97人妻天天添夜夜摸| 国产精品自产拍在线观看55亚洲 | av欧美777| 精品久久久久久电影网| 超碰97精品在线观看| 亚洲国产欧美一区二区综合| 99久久综合免费| 日本精品一区二区三区蜜桃| 亚洲欧洲日产国产| 国产精品麻豆人妻色哟哟久久| 亚洲一卡2卡3卡4卡5卡精品中文| 水蜜桃什么品种好| 久久中文字幕一级| 精品国内亚洲2022精品成人 | 俄罗斯特黄特色一大片| 91老司机精品| 精品一品国产午夜福利视频| 一本—道久久a久久精品蜜桃钙片| 老熟妇仑乱视频hdxx| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 男女下面插进去视频免费观看| 美国免费a级毛片| 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 91麻豆精品激情在线观看国产 | www.熟女人妻精品国产| h视频一区二区三区| 色综合欧美亚洲国产小说| 国产野战对白在线观看| 亚洲欧美成人综合另类久久久| 老司机靠b影院| 亚洲五月婷婷丁香| 久久人人爽人人片av| cao死你这个sao货| 男人爽女人下面视频在线观看| 肉色欧美久久久久久久蜜桃| netflix在线观看网站| 日本黄色日本黄色录像| 久久综合国产亚洲精品| 他把我摸到了高潮在线观看 | 亚洲精品乱久久久久久| videos熟女内射| 日韩大片免费观看网站| 国产在线一区二区三区精| 女警被强在线播放| 亚洲全国av大片| 99久久国产精品久久久| 日韩制服丝袜自拍偷拍| bbb黄色大片| 自线自在国产av| 日本wwww免费看| 黑人巨大精品欧美一区二区蜜桃| 国产欧美日韩综合在线一区二区| 自线自在国产av| 国产精品秋霞免费鲁丝片| 欧美+亚洲+日韩+国产| 国产成人一区二区三区免费视频网站| 18禁观看日本| 最黄视频免费看| 18禁观看日本| 亚洲国产av新网站| 日韩欧美国产一区二区入口| 亚洲国产av新网站| 97在线人人人人妻| 两人在一起打扑克的视频| 精品人妻1区二区| 十分钟在线观看高清视频www| 婷婷色av中文字幕| 精品人妻在线不人妻| 久久久国产成人免费| av片东京热男人的天堂| 嫩草影视91久久| 精品久久久久久电影网| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区大全| 一级毛片精品| 亚洲第一av免费看| 大片电影免费在线观看免费| 久久久久久人人人人人| 国产成人一区二区三区免费视频网站| 亚洲中文字幕日韩| 最近中文字幕2019免费版| 又黄又粗又硬又大视频| 岛国毛片在线播放| 久久青草综合色| 一个人免费看片子| 亚洲人成电影观看|