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

    Electromagnetically induced transparency via localized surface plasmon mode-assisted hybrid cavity QED

    2023-12-02 09:22:36XiaomiaoLi李曉苗FaminLiu劉發(fā)民ZigengLi李子更HongyanZhu朱虹燕FanWang王帆andXiaolanZhong鐘曉嵐
    Chinese Physics B 2023年11期
    關鍵詞:李子

    Xiaomiao Li(李曉苗), Famin Liu(劉發(fā)民), Zigeng Li(李子更), Hongyan Zhu(朱虹燕),Fan Wang(王帆), and Xiaolan Zhong(鐘曉嵐)

    School of Physics,Beihang University,Beijing 100191,China

    Keywords: quantum optics,electromagnetically induced transparency,all-optical switching,localized surface plasmon

    1.Introduction

    Electromagnetically induced transparency (EIT) is a quantum interference phenomenon where an electromagnetic field controls the response of a three-level quantum emitter (QE) to a probe field, i.e., the ensemble system can be switched from an initially opaque state to a transparent state depending on the control field.[1]The destructive interference strongly depends on the ground-state coherence and will result in the characteristic dark states.The characteristic super slow light velocity of the dark state makes many effects possible,such as slow light,[2,3]quantum memory[4,5]and quantum cooling.[6,7]

    In addition, through tuning the dark and bright states,more research is focused on the EIT effect which provides a highly efficient optical switch[8–15]based on a cavity quantum electrodynamics (QED) system in a quantum information network,[16]especially in long-distance transport.[17]As a high-performance optical switch, switching contrast is one of the most important factors in optical networks.A completely transparent or opaque result means that the system can selectively allow or reject transport for a specific signal.An important class of schemes for quantum communication and computing is based on one-sided and two-sided systems.On the one hand, highly efficient quantum memory and perfect photon absorption have been demonstrated in onesided cavity QED systems which meet an impedance-matched condition.[18–20]Nevertheless,high switching contrast has not been obtained in one-sided cavity QED systems for high cooperativity of the systemC=g2/(2κγs) (gis the QE–cavity coupling strength,κis the decay rate of the cavity andγsis the dephasing rate of the QE).Therefore, this one-sided system hinders multitasking devices including transistor action and quantum storage.[21]However, the elimination of perfect absorption in the two-sided system owing to the EIT effect has been achieved theoretically excited by two coherent light fields from two mirrors.[22]Notably,good transparency may turn off the control field fast enough,but make quantum storage inefficient.Hence,improving the high and stable switching contrast and manipulating the EIT effect remain to be elucidated in the one-sided cavity QED system.

    Then, exploring the system and how to manipulate the coupling strength and improve switching contrast in the expanded cavity QED is primary.In our previous work, we found that when the light is strongly coupled with two different emitters, a cascade frequency-dependent dual Rabisplitting can be achieved by adjusting the cavity mode reflecting the cavity-molecule detuning.[23]Importantly, Doelemanet al.and Hoanget al.have reported a great deal of promising phenomena in plasmon–photon cavity systems, such as cooperative Purcell enhancements and ultrafast spontaneous emission.[24,25]Thus, it provides the possibility that through manipulating the coupling strength in a hybrid plasmon–photon cavity QED system, the dark and bright states can be flexibly switched resulting in a variable output frequency.The nanometer-sized metal nanoparticle(MNP),as one of the potential candidates,has many unique advantages.

    Firstly, the localized surface plasmon (LSP) mode does not require wave vector compensation and generates a subwavelength localized mode because of the much smaller mode volume with deep subwavelength order.[26–29]Secondly, the switching time is a potential challenge in the traditional cavity QED-based optical switch limited by the QE–cavity interaction strength via tuning their spectral overlaps.[14]By means of adding an MNP structure, a three-spectral overlap can be achieved because the coupling strength of the system is boosted greatly dependent on the more confined field of the LSP mode and then the radiative rate of emitters is increased to a great extent.[30–32]Particularly,in the LSP modeprovided cavity QED, the stronger coupling strength of the MNP–QE interaction than the microcavity–QE interaction is shown in Refs.[33,34] due to the much lower mode volume of the LSP mode,which can achieve lower-threshold lasing in the latter.[33]Hence, quantum plasmonics, as a new research field, exhibits very attractive prospects in quantum information,plasmonic quantum lasers and highly sensitive detection,etc.[35]

    Here, we propose an all-optical switch designed by the LSP mode-assisted cavity QED system,which is composed of a three-level state(3LS)QE and a metal nanosphere confined in the center of a cavity.In this one-sided system,flexible and variable frequency of transmitted light is accessible by changing the intensity of the control field and the detuning of frequencies.In addition,high and stable switching contrast could be implemented showing the switching on/off state by the control field, which is vital for chip-based quantum communications and fiber-based long-distance transport.In particular,we employ a precise quantum interference explanation to demonstrate the occurrence of a converted dark state.Ascribing to multiple interactions of three modes, three systematic states are formed.Afterward, two transition pathways are formed primarily passing through one same systematic state, resulting in the destructive interference and a dark state without the control field.Then the dressed state of the QE is produced by the resonant control field and the dark state is eliminated.The evolution is analyzed based on the Heisenberg–Langevin (HL)equation and input-output relation which is well consistent with the master equation.

    2.General model

    We study the LSP mode-assisted hybrid cavity QED system interacting with a 3LS QE as shown in Fig.1(a).The system consists of a single-mode cavity and evanescent field provided by a metal nanosphere (MNP) which are tuned in resonance with each other.

    The QE is assembled onto the MNP surface with a fixed distancedby manipulating the optical trap with a highly focused laser.[36,37]The hybrid structures are adjustable in an optical microcavity.Here the radius of the MNP isR=10 nm to ensure that plasmonic resonance approximates the Lorentzian line shape.The mode volume of the MNP is much smaller than the cavity by an amount of order(~5)to break the diffraction limit.[26]The selection principle of the MNP surface–QE distance is a compromise between avoiding quenching and producing a more confined field.For this purpose, it is prudentially selected asd=10 nm referring to quenching distance,and provides the potential for ultrafast photon emission.[38]

    Under the dipole approximations,the Hamiltonian of the composite system can be written asH=H0+HI+Hdrivewhich is given by

    The ground state|1〉 and the excited state|3〉 with the excitonic transition frequencyω31are coupled to the microcavity and the LSP mode with the coupling strengthg1andg2,respectively.Mode enhanced factorJis the coupling strength due to the superposition of the vacuum and confined evanescent field between the microcavity and the MNP.

    The microcavity–LSP and QE–LSP (microcavity) couplings are described within the rotating-wave approximation(RWA) with the input operator defined by the unitary operator exp[-iω0H0t], i.e.,O(t)(O=a,b)is replaced byO(t)→Oeiω0t.The RWA is valid provided that the following criteria are satisfied:ωc-ωm?ωc+ωm,ω12(23)-ωc?ω12(23)+ωc,J ?ωc,ωm, andg1(2)?ωc,ωm,ωi j(i,j=1,2,3).The hybrid system is described by the Hamiltonian

    where the detunings are defined as the transition of QE–probe detuning?0=ω0-ω31,microcavity–QE detuning?c=ωcω31,and control frequency detuning?23=ωp-ω23.The frequency of the LSP mode and microcavity fields is regarded as the same all the time (ωc=ωm).Compared to the traditional cavity-induced EIT,[39]the major difference lies in the mode enhanced factorJbetween the microcavity and the plasmon mode through introducing the MNP into the microcavity.The frequency of the microcavity and transition|3〉?|1〉 of the QE is resonant, i.e., the detuning?c=0.Moreover, our model can provide a distinct method to control EIT and optical switching simultaneously.

    2.1.Steady-state analytical solution

    We introduce the significant damping Hamiltonian term of the open systems as presented in Ref.[40],which includes the heating loss of the bath and the interaction with considered system components.Then the terms are given as

    wheremandnare defined as the radiative decay due to the mode radiative output and the nonradiative decay due to ohmic losses of the open environment,which obey the commutation relation [m,n]=0, respectively.The reservoirsT1–T3represent the coupling strength between the reservoir and the microcavity, the MNP and the QE.Substituting Eq.(3) into the Heisenberg–Langevin equations of motion,the dynamic equations of a coupled microcavity-MNP-QE system are given as

    where operators denote the expectation valueO=〈O〉on the left side of the equation.κ1=πT21,κ2=πT22,γs=πT23=(γ32+γ31)/2 represent the decay of the microcavity mode,the sum of the ohmic loss and the dephasing rate of the MNP and the dephasing rate of the QE.Moreover,fa,fb,f13,f23,f12,f11,f22andf33are the noise operators which are derived from the reservoir interaction and reservoir modes which satisfy commutation relation at all times.σz=σ33-σ11is the population inversion of the QE.The Langevin noise operators have zero mean values by assumption in this treatment,as shown already in earlier works.[40]When the system is at the initial point,the great majority of the population of the QE distribute in the|1〉 state, i.e., (σ11≈1,σ22=σ33≈0).The steady-state solutions of Eq.(4) can be readily obtained analytically.The analytical solution of the cavity modeais given as

    2.2.Numerical description

    In general,the properties of the entire system can be obtained numerically by replacing Eq.(3)with the master equation as

    Then,we elucidate how the LSP mode-assisted one-sided system manipulates and changes the characteristics of the photon emitting cavity and switching contrast with a control field.We also quantitatively discuss the quantum interference effect of the variable frequency of output light (VF-EIT) and traditional EIT by the presence and absence of the MNP structure,respectively.

    3.Results and discussion

    3.1.Characteristic output properties

    Fig.2.For the (a) variable frequency EIT (VF-EIT) (red line) and (b) traditional EIT (blue line), the output intensity of the MNP-assisted system and traditional cavity QED system numerically (solid line) and analytically (dot line)based on the H-L and master equation with ? =0,? =1 meV and ? =2 meV (from top to bottom panels).(c) Selective frequency in the VF-EIT with tuning strength of control field ? =0–2 meV.A pair of bright modes are plotted(dark dashed line).Other parameters are{κ1,κ2,γs,γ21,η,g1,g2,J}=0.07,0.7,0.086,0.002,0.1,0.0206,1,2.9 meV for (a).Other parameters are{κ1,κ2,γs,γ21,η,g1,g2,J}=0.07,0,0.086,0.002,0.1,0.15,0,0 meV for(b).

    Figure 2(a)shows the output transmission as the expression|aout|/|ain| of the proposed scheme.The fixed parameters in this section are{κ1,κ2,γs,γ21,η}=0.07, 0.7, 0.086,0.002, 0.1 meV.The decay of the cavity is much smaller than the MNP with an order of around 10, and the weak probeηis smaller than the coupling strength of the system as well to satisfy the nonlinear effect at low optical intensity.The coupling strengths{g1,g2,J}=0.0206,1,2.9 meV of the interaction of the QE–cavity mode, the QE–MNP and the cavity–MNP are calculated referring to Ref.[26]but here providing the smaller dipole momentμ13=1e·nm.Due to the much smaller mode volume of the MNP than the microcavity, the coupling strengthg1?g2is reasonable in spite of the stronger loss of the latter.As we can see, the output transmission of the analytical H-L equation and the numerical master equation are accordant with each other as seen in the experimental results of Refs.[14,41].At?0=0, the narrow spectral peak representing the dark mode shows without the control field (?=0), while it disappears at?/=0.The detailed and quantitative description of the dark state is deduced using an interaction Hamilton later.

    Figure 2(b) plots the normalized transmitted light power versus?0without the assisted LSP mode (J=0,κ2= 0,g2= 0), which shows that the EIT effect makes the obvious conversion with the VF-EIT in the MNP-assisted cavity QED system.The system is composed of one microcavity and a three-level QE, and the parameters are chosen as{κ1,κ2,γs,γ21,η,g1,g2,J}=0.07,0,0.086,0.002,0.1,0.15,0,0 meV.When?= 0, the Rabi splittingThe narrow dip at?0=0 shows that the system is under strong coupling in accordance withg>κ ?γ.Conversely,when the control field is present,i.e.,?/=0,the narrow spectral peak at?0=0 shows that the control field induces destruction interference between excited states of the QE.The probe light is essentially transmitted in the cavity and has little absorption by the medium.It is ascertained that the control field plays a critical role to achieve switching performance from the above two results.In view of Figs.2(a)and 2(b),there is no deviation between the peaks of the analytical and numerical results,but the strength of transmission shows differences due to the static approximation of the H-L equation.

    Figure 2(c) depicts the VF-EIT with the wide range of strength of control field?=0–2 meV.This demonstrates the existence of the EIT effect and exhibits the switching response in our system from?= 0 to a certain strength.The proposed model could provide an easier operation of light switching which provides continuous transmission by manipulating the detuning?0and the intensity of the control field?in the change from blue to red lines shown in Fig.2.Additionally,a pair of bright modes show in the off-resonant frequencies,and total coupling strength is enhanced.

    Fig.3.The dependence of the EIT effect on detunings ?0 and ?p.(a)–(c)are the traditional cavity QED system,while(d)–(f)are the VF-EIT effects provided by the plasmon-assisted system,both including a 3LS QE for ? =0 meV,1 meV,2 meV.Parameters are the same as above.

    We consider the variation of output intensity if the probe light and control field are tuned off-resonance simultaneously in order to verify the dependence of additional conditions depicted in Fig.3.Here we discuss traditional and MNP-assisted systems as above with the variation of control field?=0–2 meV depicted in the first and second rows.The first panel is simplified as the cavity QED system with a three-level QE shown in Figs.3(a)–3(c).There are just two bright modes without the control field,while splitting once with the increasing control field.In the second panel,note that the intracavity dark mode splits into two modes which are enhanced by the increasing control field strength plotted in Figs.3(d)–3(f).However, the bright modes remain in the fixed detuning position?0=0.16 meV in spite of the varied?.The second anticrossing is more obvious at the position of the bright mode.When the applied control field boosts,the frequency of the second Rabi splitting center is resonant?p=0 meV.Thus,we can not only control the splitting position and strong coupling of the system by the control field but also the switch‘on/off’.Notably,the dark state apparently exists in our system without the need for the control field so the system shows a strong transmitted intensity.The bright modes remain almost unchanged in both results.Contrary to the EIT effect, the extra advantage of the proposed system is the shifting frequency of output light, showing that when the probe light is off-resonant with the cavity,the same transmitted efficiency can be exhibited by manipulating the intensity or power of the control field.This is useful for selective undistinguished frequencies of photons to enter and entangle in the entanglement swapping by means of tuning the strength of the coupling field.Referring to the results in Ref.[43],we have already given a general explanation of VF-EIT by adding an extra field so as to enrich the coupling process.This system of Guoet al.studied an N-type fourlevel QE and discussed the contribution of a microwave field added to couple the two ground states.This provides a helpful method to regulate the entry of photons and exclude the noise photons in the network node.The range of tuning frequency of probe light is ideally wide within the gap of bright modes as shown in Figs.3(e)and 3(f)(?0=-0.5 meV to 0.5 meV).Notably,the time for the steady state in the plasmon-assisted cavity QED system is<50 ns which is three times smaller than that of the traditional cavity QED(150 ns)given in Ref.[44].

    3.2.Switching contrast

    The switching contrast is defnied by the difference of normalize transmission of the resonant photonTcst=At the resonance?0=0, the contrast value of the LSP mode-assisted cavity QED system and conventional cavity QED has a significant discrepancy with 60%and<10% given in Figs.2(a) and 2(b).Further, the dependence of the switching contrast on the decoherence of the QE,cavity dissipation, imperfect mode matching and the LSP mode-QE cooperationC=(g2)2/(2κ2γs)is given in Figs.4(a)–4(d).Obviously,the mode-enhanced factor J is quite informational due to the existence of a maximum point in the switching contrast curve in the middle strengthJ ≈0.6 meV, which shows the critical dependence of a well-designed system.The FWHM corresponds to a range of coupling strengthJ(0.17–2.6 meV),which can be modulated by the spatial position and size of the MNP in the cavity.Parametersκ2andγsplay subtle roles in exploring the maxima or minima.Weak dependence of the decay rate of the LSP mode provides an advantage for eliminating cavity loss.Simultaneously, switching contrast diminishes exponentially with the increased dephasing rate of the excited state of the QE.The above discussion reveals that the moderate coupling strength of two fields and a small dephasing rate of the excited state of the QE are preferred for high switching contrast.In Fig.4(d), stable switching contrast (about 90%) is ensured onJ=0.6 meV even at the prior high cooperationCof Ref.[21].Conversely, due to the existence of the LSP mode, light is transmitted for?= 0 and absorbed for?/=0.So less dependence on the cooperativity is shown.Moreover, we also focus on the effect of QE–microcavity couplingg1in the lack of LSP mode-assisted system.For the bad cavityg13).[43]Wanget al.demonstrated that only at a certain probe light–QE detuning and coupling strength,the resonant photons are completely reflected and the cavity QED system could be regarded as aperfect absorber with cumbersome experimental operation.[22]This indicates that the proposed LSP mode-assisted system can serve as a high-contrast and stable photonic switch using a relatively simple setup for one-sided input light.

    3.3.Optical property of the QE

    Next we discuss qualitatively the linear property of the medium to explain this VF-EIT effect.The real and imaginary parts of the susceptibility reflect the dispersive and absorbing nature.Without the extra control field, there is an interaction between the cavity field and QE level transition|1〉→|3〉.According to the Heisenberg–Langevin analytical method, the significant transmission peak appears at the resonant frequency due to rare absorption of the medium shown in Fig.5(a)(red line),which corresponds to the dark mode in Fig.5(a)(?=0),manifesting that the incident probe photon is nearly all trapped in the microcavity.At off-resonant frequencies shown in Fig.5(b),the absorption has two peaks and also corresponds to the gap between bright and dark modes in Fig.2 (?=1,2 meV).However, with the control field, the complete reversal happens at?0=0 and exhibits the vanishing of the dark state at the fully resonant frequency.The real part of permittivity plotted in Fig.5(c) shows the anomalous dispersion from one to two contours corresponding to the central linewidth of the absorption profile.To make a contrast,the absorption and dispersion of a one QE–one cavity system are also depicted in Figs.5(a)–5(d) (blue line).Without the control field,the output transmission exhibits the strong peaks of Rabi splitting which illustrates the non-existent dark state,as shown in Fig.5(a)(blue line).Due to the manipulation of the control field, the strong absorption at?0/=0 corresponds to the weak bright mode,while the stronger EIT peak is shown at a resonant frequency in Figs.5(b) (blue line) showing the occurrence of the dark state.The reflective index of the LSP mode-assisted system is four times as large as that of the traditional cavity QED system shown in Figs.5(c)and 5(d)showing stronger manipulation of the structure.

    Fig.5.The absorption and dispersion of the medium of ? =1 meV(a)and(c),and ? =0(b)and(d)for VF-EIT of the MNP-assisted cavity QED system (red line) and EIT system (blue line).Other parameters are equal to the MNP-assisted and traditional cavity QED system as in the above discussions,respectively.

    3.4.Quantum interference explanation

    Fig.6.Energy-level and transition pathway diagram of the coupled cavity–QE(a)without and(b)with the control field in the LSP modeassisted system.The quantum interference effect between different transition pathways for the transition |3〉?|2〉 of the QE due to the presence of a control field leads to the counterintuitive EIT effect.The Dirac operator is a direct product of the Fork state of the cavity photon and the population of the QE.Here one Fork state is considered. g1 and g2 is the coupling between the transition of |3〉?|1〉 and the microcavity,and the LSP mode,respectively. J refers to the mode-enhanced factor between the microcavity and LSP mode.(c) Population evolution of mean photon number of LSP mode for VF-EIT versus different strength of control field ? =0 meV and ? =2 meV.

    Considering the Hamilton term of all interactions and the pump of the control field, the density matrix of the proposed system can be given as

    The eigenvalues of the matrix reflect the states of the scheme.

    Next, in the proposed MNP-assisted system, the middle EIT peak represents the intracavity dark state when the control field is absent.The quantum explanation is destructive interference due to mode interaction between the extra LSP mode and cavity field which plays the main role.The ground state basis of the coupled system|0,0,1〉 experiences the probe lightηof the cavity and yields|1,0,1〉.Two coherent transition pathways are plotted in Fig.6(a).One is coherent coupling|1,0,1〉→|0,0,3〉withg1,while the other is two steps|1,0,1〉→|0,1,1〉→|0,0,3〉 with coupling strengthJandg2.The systematic eigenstates can be written as

    During this process,the population on middle|1,0,1〉is suppressed, steering the system to a dark state as seen with two fields in Ref.[42].Therefore, there is no medium absorption and no effects on the cavity photons.Notably,variation of?has little effect on the frequency range of bright modes shown in Figs.3(b)and 3(c).The broad frequency of bright modes is located on the off-resonant frequency at?p=±0.5 meV.

    In the case of applying the control field,the excited state of the QE degenerates to a dressed state, and the whole systematic states include|1,0,1〉,|0,1,1〉,|0,0,±〉as shown in Fig.6(b).The destructive interference between two coherent excitation paths of the QE(|2〉±|3〉)seen in the conditional EIT effect is induced due to the same strength and opposite phase[43]at the frequency of probe resonance.However, the third coherent pathway is formed due to mode-enhanced factor J between|1,0,1〉 and|0,1,1〉.Thus, the formation of destructive interference and the dark state vanishes.Additionally,variation of?has proportional manipulation on Rabi splitting of dark modes.As shown in Fig.2,a single EIT peak splits into two peaks near resonance frequency and resonant frequency can be given as four eigenvalues and written asλ2=From the above expression,the decisive factor to prepare the dark state is whether or notJ2?2=0.Hence,reliable evidence to generate a dark state in the LSP mode-assisted system is again given without the control field.

    In Fig.6(c),it can be seen that the transient evolution of the population of the LSP cavity reflects the distribution of systematic states.For the dark stateΨdat?=0,the|0,1,1〉state is destructive interference and the number of LSP cavity photons is substantially reduced.Contrarily,the population shows an increasing trend at?/= 0 because of the vanished dark state and all contributions of the four eigenstates of the system.After destructive interference occurs,the electronic state is trapped in the excited state of the QE.This corresponds to the inversion of the dark state to manipulate the switch’on/off’state.

    4.Conclusion

    In conclusion, we have achieved the high and stable switching contrast and variable frequency of photon output from the one-sided LSP mode-assisted cavity QED system,even with smallJand high cooperativityCby tuning the intensity of the control field.A full analytical solution is given to describe transmission spectra and eigenstates for the evolution of the system on the basis of the Heisenberg–Langevin equation and input–output relation the solution of which is universal for the interaction of three modes on cavity QED systems.The quantum interference is presented to describe the evolution of systematic states by the absence or presence of the control field.Due to the destructive interference,a certain systematic state vanishes that renders one of the microcavity or LSP photons trapped without the control field.At macro scales,the electric susceptibility of the medium is significantly different between VF-EIT and EIT effects.In long-distance entanglement swapping,selecting a certain frequency of photons for Bell-state measurement(BSM)is a critical step to conduct the entangled photons.As the network nodes via the photonic link,the LSP mode-assisted system is embedded in the quantum circuit to help select a certain frequency of input photon.

    Here we propose several experimental realizations to achieve the LSP mode-assisted cavity QED to manipulate optical switching.Jauffredet al.optically trapped an individual airborne gold nanoparticle using a single infrared laser beam.[45]By constructing three chambers to trap a single particle easily, the measurement time is about 0.04 s, which is much longer than the switching process(the order of ns).TheΛ-type artificial QE can also be trapped by applying another laser with visible wavelength.Then ideally, through infinite approaching MNP and QE by modulating two lasers,the transient sub-nanometer distance between the MNP surface and the QE can be achieved.An order of GHz microwave drives the magnetic-dipole transition between two hyperfine ground states of the QE to form aΛ-type QE.Finally, by turning on/off the controlling field applying to the optical medium state-excited state, transmission of the probe light can be detected in the same direction in the one-sided system.

    However,the accurate control of the position of 10–20 nm between the QE and MNP in the air and the total resolution of the two particles is challenging.Hence,providing the fixed QE and movable MNP through the optical tip makes the theoretical system possible.Firstly,the QE can be fixed on the surface of molecules or spheres which cannot respond to the characteristics of interaction with light by a self-assembly technique.Then, the position of the MNP is controlled by an AFM tip with the precision of a few dozen nanometers.[46,47]The cavity could be provided by the waveguide.The whole experiment can be achieved with the on-chip device.

    Acknowledgements

    The authors acknowledge support from the National Natural Science Foundation of China (Grant Nos.62075004 and 11804018)and the Beijing Natural Science Foundation(Grant No.4212051).

    猜你喜歡
    李子
    我戰(zhàn)勝了黑暗
    小主人報(2022年7期)2022-08-16 06:59:28
    智能鼠捉鼠記
    謝謝你
    小主人報(2022年2期)2022-01-19 11:46:20
    一次難忘的生日
    秋天
    李子有多少
    給“小松”洗澡
    奔跑吧!李子柒
    海峽姐妹(2020年1期)2020-03-03 13:35:52
    我的糊涂媽媽
    Formation of electron depletion layer and parallel electric field in the separatrix region of anti-parallel magnetic reconnection?
    色94色欧美一区二区| av网站在线播放免费| 夜夜躁狠狠躁天天躁| 九色亚洲精品在线播放| 午夜福利一区二区在线看| 国产伦人伦偷精品视频| 最新在线观看一区二区三区| 无限看片的www在线观看| 好男人电影高清在线观看| 80岁老熟妇乱子伦牲交| 国产精品 国内视频| 日本vs欧美在线观看视频| 一级毛片女人18水好多| 曰老女人黄片| 久久精品aⅴ一区二区三区四区| 99精国产麻豆久久婷婷| 亚洲人成电影免费在线| 成人18禁在线播放| 国产精品一区二区精品视频观看| 国产不卡一卡二| 国产精品av久久久久免费| 人人妻人人澡人人看| 国产99久久九九免费精品| 精品国产亚洲在线| 日本精品一区二区三区蜜桃| 1024视频免费在线观看| www.精华液| 不卡av一区二区三区| 日本五十路高清| 女人久久www免费人成看片| 男人操女人黄网站| 中文字幕最新亚洲高清| 亚洲精品国产区一区二| 精品人妻熟女毛片av久久网站| 午夜福利在线观看吧| 免费久久久久久久精品成人欧美视频| 国产精品国产高清国产av | 中国美女看黄片| 亚洲免费av在线视频| 国产1区2区3区精品| 久久久精品区二区三区| 午夜福利,免费看| 久久久久国内视频| 我的亚洲天堂| 99久久人妻综合| 亚洲av电影在线进入| 69精品国产乱码久久久| 精品高清国产在线一区| 在线视频色国产色| 18禁美女被吸乳视频| 最新在线观看一区二区三区| 欧美日韩精品网址| 国产av一区二区精品久久| av网站免费在线观看视频| av网站在线播放免费| 色老头精品视频在线观看| 欧美在线黄色| 国产黄色免费在线视频| 丰满人妻熟妇乱又伦精品不卡| 九色亚洲精品在线播放| 久久精品aⅴ一区二区三区四区| 男女之事视频高清在线观看| 国产单亲对白刺激| 水蜜桃什么品种好| 亚洲男人天堂网一区| av超薄肉色丝袜交足视频| 国产精品久久电影中文字幕 | 亚洲男人天堂网一区| 超色免费av| 亚洲精品一二三| 国产亚洲精品久久久久5区| 别揉我奶头~嗯~啊~动态视频| 一级毛片女人18水好多| 飞空精品影院首页| 国产av精品麻豆| 国产成人免费观看mmmm| 精品福利观看| 老汉色av国产亚洲站长工具| 亚洲成a人片在线一区二区| 精品国产国语对白av| 亚洲美女黄片视频| 国产精品自产拍在线观看55亚洲 | 久久精品国产综合久久久| 日韩欧美在线二视频 | 搡老岳熟女国产| 91精品三级在线观看| 成人亚洲精品一区在线观看| 久久久久精品人妻al黑| 五月开心婷婷网| 中文字幕人妻丝袜一区二区| 亚洲第一青青草原| 国产精品国产av在线观看| 国产精品免费一区二区三区在线 | 久久精品91无色码中文字幕| 男人操女人黄网站| 国产单亲对白刺激| 50天的宝宝边吃奶边哭怎么回事| 一区二区三区国产精品乱码| 国产在视频线精品| 黄色成人免费大全| 成人av一区二区三区在线看| 免费黄频网站在线观看国产| 在线视频色国产色| www.999成人在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图 男人天堂 中文字幕| 欧美大码av| 日本黄色视频三级网站网址 | 美国免费a级毛片| 亚洲av成人不卡在线观看播放网| 国产伦人伦偷精品视频| 丰满迷人的少妇在线观看| 婷婷成人精品国产| 国产av精品麻豆| 久久久久国产一级毛片高清牌| a级毛片黄视频| 人妻 亚洲 视频| 18禁裸乳无遮挡免费网站照片 | 精品一区二区三卡| 国产色视频综合| av国产精品久久久久影院| 国产高清视频在线播放一区| 一夜夜www| 免费久久久久久久精品成人欧美视频| 妹子高潮喷水视频| 中文字幕另类日韩欧美亚洲嫩草| 美女扒开内裤让男人捅视频| 国产成人欧美| 久久中文字幕一级| 亚洲欧洲精品一区二区精品久久久| 国产亚洲一区二区精品| 黄色视频不卡| 首页视频小说图片口味搜索| 侵犯人妻中文字幕一二三四区| 亚洲精品美女久久久久99蜜臀| 亚洲人成电影免费在线| a级毛片黄视频| 他把我摸到了高潮在线观看| 别揉我奶头~嗯~啊~动态视频| 一级毛片女人18水好多| 日韩制服丝袜自拍偷拍| 人成视频在线观看免费观看| 国产激情欧美一区二区| videos熟女内射| 在线看a的网站| 国产又色又爽无遮挡免费看| 久久午夜亚洲精品久久| 中文字幕精品免费在线观看视频| 午夜福利,免费看| 99久久国产精品久久久| 国产精品99久久99久久久不卡| 亚洲成人国产一区在线观看| 国产成人av教育| 国产精品久久久久成人av| 久久久精品免费免费高清| 久久精品亚洲av国产电影网| 成人永久免费在线观看视频| 深夜精品福利| 少妇的丰满在线观看| 亚洲国产精品合色在线| 国产在线观看jvid| 欧美激情极品国产一区二区三区| 午夜福利免费观看在线| 久久久久久久精品吃奶| 少妇的丰满在线观看| 91精品国产国语对白视频| 国产亚洲av高清不卡| 天堂动漫精品| 91麻豆av在线| 午夜视频精品福利| 亚洲国产精品一区二区三区在线| 久久人妻av系列| 亚洲av成人不卡在线观看播放网| 国产精品久久视频播放| 日日夜夜操网爽| videos熟女内射| 久久天躁狠狠躁夜夜2o2o| 夜夜躁狠狠躁天天躁| 国产日韩欧美亚洲二区| 中文欧美无线码| 黑人操中国人逼视频| 日韩欧美一区视频在线观看| 亚洲人成伊人成综合网2020| 久久国产精品影院| 国产精品av久久久久免费| 国内久久婷婷六月综合欲色啪| 亚洲国产欧美日韩在线播放| 啪啪无遮挡十八禁网站| 精品卡一卡二卡四卡免费| 午夜福利欧美成人| 精品国产一区二区久久| 一区二区三区激情视频| 国产精品av久久久久免费| 欧美精品av麻豆av| 久久久国产欧美日韩av| 久9热在线精品视频| 99国产精品一区二区三区| 1024香蕉在线观看| 国产成人欧美在线观看 | 人人妻人人澡人人看| 成年版毛片免费区| 精品国产一区二区三区久久久樱花| 亚洲国产中文字幕在线视频| 新久久久久国产一级毛片| 久久精品91无色码中文字幕| 好看av亚洲va欧美ⅴa在| 久久精品aⅴ一区二区三区四区| av福利片在线| 午夜两性在线视频| 欧美在线一区亚洲| 无限看片的www在线观看| 亚洲,欧美精品.| 精品一区二区三区四区五区乱码| 亚洲欧美激情在线| 国产成人精品久久二区二区91| 亚洲欧美日韩高清在线视频| 欧美av亚洲av综合av国产av| 国产欧美日韩综合在线一区二区| 精品国产亚洲在线| 午夜亚洲福利在线播放| 啦啦啦视频在线资源免费观看| 一a级毛片在线观看| 69精品国产乱码久久久| 欧美成人午夜精品| 男女下面插进去视频免费观看| 亚洲 欧美一区二区三区| 亚洲国产精品sss在线观看 | 免费久久久久久久精品成人欧美视频| 人人妻,人人澡人人爽秒播| 亚洲五月天丁香| 99国产精品99久久久久| 亚洲男人天堂网一区| 精品欧美一区二区三区在线| 成人av一区二区三区在线看| 嫁个100分男人电影在线观看| 久久精品国产综合久久久| a级片在线免费高清观看视频| 日本撒尿小便嘘嘘汇集6| 亚洲五月色婷婷综合| 90打野战视频偷拍视频| 一级毛片女人18水好多| 1024香蕉在线观看| 国产单亲对白刺激| 99国产极品粉嫩在线观看| 欧美日韩成人在线一区二区| 超碰成人久久| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 精品久久久久久,| 中文欧美无线码| 色婷婷久久久亚洲欧美| 免费高清在线观看日韩| av国产精品久久久久影院| 成人av一区二区三区在线看| 99精国产麻豆久久婷婷| 在线观看免费高清a一片| 国产欧美日韩一区二区三区在线| 叶爱在线成人免费视频播放| 9色porny在线观看| 一进一出好大好爽视频| 精品欧美一区二区三区在线| videos熟女内射| 波多野结衣一区麻豆| 色老头精品视频在线观看| 视频区欧美日本亚洲| av一本久久久久| av线在线观看网站| 国产欧美日韩综合在线一区二区| 丝袜在线中文字幕| 日韩大码丰满熟妇| 美女国产高潮福利片在线看| 男女之事视频高清在线观看| 好看av亚洲va欧美ⅴa在| 久久久久精品人妻al黑| av天堂久久9| 欧美在线黄色| 18禁裸乳无遮挡动漫免费视频| 悠悠久久av| 国产精品 国内视频| 美国免费a级毛片| 19禁男女啪啪无遮挡网站| 色尼玛亚洲综合影院| 国产精品影院久久| 69精品国产乱码久久久| 身体一侧抽搐| 后天国语完整版免费观看| 欧美激情久久久久久爽电影 | 久久天堂一区二区三区四区| 婷婷丁香在线五月| 久久精品国产a三级三级三级| av超薄肉色丝袜交足视频| 国产成+人综合+亚洲专区| 国产欧美日韩精品亚洲av| 中出人妻视频一区二区| 午夜免费成人在线视频| 日本一区二区免费在线视频| 侵犯人妻中文字幕一二三四区| 天堂中文最新版在线下载| 啦啦啦视频在线资源免费观看| 宅男免费午夜| 亚洲综合色网址| 国产男靠女视频免费网站| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 亚洲欧美一区二区三区黑人| 老熟妇仑乱视频hdxx| 国产精品久久视频播放| 97人妻天天添夜夜摸| 18禁美女被吸乳视频| 亚洲成人手机| 亚洲中文av在线| 久久国产精品人妻蜜桃| 80岁老熟妇乱子伦牲交| 国产精品免费大片| 极品少妇高潮喷水抽搐| 午夜老司机福利片| 大型黄色视频在线免费观看| 国产av又大| 女同久久另类99精品国产91| av网站免费在线观看视频| 在线观看舔阴道视频| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 国产欧美日韩一区二区三| 亚洲成国产人片在线观看| 亚洲av成人不卡在线观看播放网| 日韩有码中文字幕| 国产99白浆流出| 色精品久久人妻99蜜桃| 亚洲一区二区三区不卡视频| 久久久国产精品麻豆| 欧美久久黑人一区二区| 91麻豆精品激情在线观看国产 | 欧美日韩视频精品一区| 欧美久久黑人一区二区| 在线观看66精品国产| 亚洲久久久国产精品| 夜夜夜夜夜久久久久| 一边摸一边做爽爽视频免费| 欧美午夜高清在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产亚洲av麻豆专区| 国产精品影院久久| 极品教师在线免费播放| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 国产精品久久久av美女十八| 欧美日韩一级在线毛片| 身体一侧抽搐| 欧美成狂野欧美在线观看| 国产精品99久久99久久久不卡| 亚洲一码二码三码区别大吗| 亚洲专区字幕在线| 国产精品九九99| 多毛熟女@视频| 手机成人av网站| 一进一出好大好爽视频| aaaaa片日本免费| 亚洲精品中文字幕在线视频| 一夜夜www| 国产精品欧美亚洲77777| 妹子高潮喷水视频| 久久国产精品大桥未久av| 日本撒尿小便嘘嘘汇集6| 国产精品乱码一区二三区的特点 | 婷婷丁香在线五月| 亚洲熟女毛片儿| 女人久久www免费人成看片| 国产亚洲欧美98| 变态另类成人亚洲欧美熟女 | 91在线观看av| 热99久久久久精品小说推荐| 中文字幕人妻熟女乱码| 久久精品亚洲熟妇少妇任你| 一进一出好大好爽视频| 极品人妻少妇av视频| 亚洲欧美日韩另类电影网站| 日韩视频一区二区在线观看| 色尼玛亚洲综合影院| 自线自在国产av| 18禁黄网站禁片午夜丰满| 三上悠亚av全集在线观看| tocl精华| 老熟妇仑乱视频hdxx| 免费在线观看黄色视频的| 亚洲全国av大片| 日韩有码中文字幕| 日韩制服丝袜自拍偷拍| 欧美黄色淫秽网站| tube8黄色片| 可以免费在线观看a视频的电影网站| 免费看a级黄色片| 色婷婷av一区二区三区视频| 久久国产乱子伦精品免费另类| 在线av久久热| 国产主播在线观看一区二区| 精品无人区乱码1区二区| 国产亚洲精品久久久久久毛片 | 午夜福利视频在线观看免费| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲精品一区二区精品久久久| 悠悠久久av| 狠狠婷婷综合久久久久久88av| 一级毛片女人18水好多| 黄色女人牲交| 欧美黄色淫秽网站| www.自偷自拍.com| 成年人黄色毛片网站| www日本在线高清视频| 国产欧美亚洲国产| 老司机午夜福利在线观看视频| 12—13女人毛片做爰片一| 电影成人av| 国产精品乱码一区二三区的特点 | 岛国在线观看网站| 啦啦啦 在线观看视频| 大片电影免费在线观看免费| 成人av一区二区三区在线看| 亚洲精品自拍成人| 久久精品亚洲精品国产色婷小说| 欧美激情极品国产一区二区三区| 国内久久婷婷六月综合欲色啪| 久久精品国产a三级三级三级| 久久精品成人免费网站| 美女国产高潮福利片在线看| 岛国毛片在线播放| 窝窝影院91人妻| 日本黄色日本黄色录像| 大片电影免费在线观看免费| 国产精品综合久久久久久久免费 | av天堂久久9| 亚洲欧美色中文字幕在线| 久久九九热精品免费| 亚洲五月色婷婷综合| 成在线人永久免费视频| 脱女人内裤的视频| 精品国产美女av久久久久小说| 性少妇av在线| 亚洲欧美日韩另类电影网站| 啪啪无遮挡十八禁网站| 亚洲av成人一区二区三| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 大型黄色视频在线免费观看| 18禁裸乳无遮挡免费网站照片 | 校园春色视频在线观看| 自线自在国产av| 超色免费av| av国产精品久久久久影院| 久久国产乱子伦精品免费另类| 国产成人一区二区三区免费视频网站| 成人av一区二区三区在线看| 两性夫妻黄色片| 久久国产精品影院| 天天操日日干夜夜撸| 无限看片的www在线观看| 亚洲精品中文字幕在线视频| 在线观看午夜福利视频| 曰老女人黄片| 国产一区二区激情短视频| 国内久久婷婷六月综合欲色啪| 国产亚洲精品第一综合不卡| 一区在线观看完整版| 黄色成人免费大全| 51午夜福利影视在线观看| 久久人妻av系列| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 成年版毛片免费区| 欧美激情高清一区二区三区| 黄色毛片三级朝国网站| 欧美黑人精品巨大| 欧美丝袜亚洲另类 | 韩国精品一区二区三区| 青草久久国产| 亚洲欧美精品综合一区二区三区| 最近最新中文字幕大全免费视频| 99精品在免费线老司机午夜| 一a级毛片在线观看| 亚洲成国产人片在线观看| 免费看a级黄色片| 午夜老司机福利片| 国产一区二区三区视频了| 亚洲aⅴ乱码一区二区在线播放 | 很黄的视频免费| 乱人伦中国视频| 亚洲成a人片在线一区二区| 国产精品久久久av美女十八| 大片电影免费在线观看免费| av视频免费观看在线观看| 最新美女视频免费是黄的| 久久久久精品人妻al黑| 亚洲在线自拍视频| 成年动漫av网址| 一进一出抽搐动态| 国产无遮挡羞羞视频在线观看| 亚洲五月婷婷丁香| 激情视频va一区二区三区| 18在线观看网站| 欧美日韩瑟瑟在线播放| 51午夜福利影视在线观看| 日日夜夜操网爽| 亚洲av成人av| av天堂久久9| 一本一本久久a久久精品综合妖精| www.精华液| 久久天堂一区二区三区四区| 91大片在线观看| 高清视频免费观看一区二区| 久久午夜综合久久蜜桃| 深夜精品福利| 国产日韩一区二区三区精品不卡| 黄片小视频在线播放| 最新美女视频免费是黄的| 在线免费观看的www视频| 欧美大码av| 国产xxxxx性猛交| 色综合欧美亚洲国产小说| 飞空精品影院首页| 制服诱惑二区| 久久国产精品影院| 国产精品成人在线| 亚洲精品国产区一区二| 免费看十八禁软件| 成年版毛片免费区| 日韩人妻精品一区2区三区| 日韩免费高清中文字幕av| 国产欧美日韩一区二区三| 欧美中文综合在线视频| 精品人妻在线不人妻| 多毛熟女@视频| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 天天躁狠狠躁夜夜躁狠狠躁| 99久久国产精品久久久| 亚洲精品在线美女| 99香蕉大伊视频| 嫩草影视91久久| 一进一出抽搐动态| 丝瓜视频免费看黄片| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品久久久久久毛片 | 岛国在线观看网站| 天堂中文最新版在线下载| 激情在线观看视频在线高清 | 亚洲av成人一区二区三| 国产精品九九99| 制服人妻中文乱码| av欧美777| www.精华液| av视频免费观看在线观看| 国产人伦9x9x在线观看| 国产不卡av网站在线观看| 欧美成人免费av一区二区三区 | 99riav亚洲国产免费| 国产精品综合久久久久久久免费 | 国产视频一区二区在线看| 中国美女看黄片| 亚洲熟女毛片儿| 久久九九热精品免费| 三上悠亚av全集在线观看| 国产精品秋霞免费鲁丝片| 国产男靠女视频免费网站| 亚洲成国产人片在线观看| 亚洲成人免费av在线播放| 久久久久国产一级毛片高清牌| 免费黄频网站在线观看国产| 久久久国产一区二区| 精品一品国产午夜福利视频| 久久ye,这里只有精品| 欧美激情 高清一区二区三区| 国产亚洲欧美精品永久| 在线观看一区二区三区激情| 18禁黄网站禁片午夜丰满| 久久人人97超碰香蕉20202| 欧美色视频一区免费| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看 | 最近最新中文字幕大全电影3 | 精品一品国产午夜福利视频| 成熟少妇高潮喷水视频| 精品一品国产午夜福利视频| 黄色怎么调成土黄色| 亚洲精品国产精品久久久不卡| 啦啦啦在线免费观看视频4| 色94色欧美一区二区| 亚洲自偷自拍图片 自拍| 久久久久国内视频| 亚洲精品在线美女| 三级毛片av免费| 亚洲国产精品sss在线观看 | 交换朋友夫妻互换小说| 精品福利观看| 又大又爽又粗| 久久婷婷成人综合色麻豆| 精品无人区乱码1区二区| 亚洲国产看品久久| 午夜久久久在线观看| 人妻 亚洲 视频| 制服人妻中文乱码| 成人免费观看视频高清| 午夜日韩欧美国产| 欧美在线一区亚洲| 12—13女人毛片做爰片一| svipshipincom国产片| 国产男女超爽视频在线观看| 国产99白浆流出| 法律面前人人平等表现在哪些方面| 国产男女超爽视频在线观看| 超碰97精品在线观看| 亚洲精品乱久久久久久| 国产一区二区三区综合在线观看| 亚洲欧美精品综合一区二区三区| 十分钟在线观看高清视频www| 日韩一卡2卡3卡4卡2021年|