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

    A comparison of transient reflection and transmission techniques for electron and lattice dynamics in metal nanofilms

    2021-12-28 11:24:22WUWenZhiRENShiWeiKONGDeGuiCHAIZhiJun

    WU Wen-Zhi, REN Shi-Wei, KONG De-Gui, CHAI Zhi-Jun

    (College of Electronic Engineering, Heilongjiang University, Harbin 150080, China)

    Abstract:The electron and lattice dynamics of gold and copper nanofilms are measured with femtosecond transient reflection and transmission techniques. Two-temperature model (TTM) and Crude-model approximation are applied to study the mechanism of heat transfer for metal nanofilms theoretically. Electron and lattice temperatures are calculated by TTM, which assume electron-lattice coupling coefficient as a constant for gold and copper. Compared the results of reflection with transmission, the difference of electron relaxation is negligible at the first two picoseconds, but electron-lattice coupling effect in the transmission method is stronger and more sensitive than that of the reflection method under the identical experimental conditions. Change of temperature gradient along with the direction of film thickness and electron-lattice coupling resistance due to the boundary scattering are responsible for the difference between transient reflection and transmission data. Both transient reflection and transmission of semitransparent films should be considered together for the investigation on the mechanism of ultrafast dynamics.

    Key words:two-temperature model (TTM); electron-lattice coupling; ultrafast pulsed laser heating; substrate effect; Crude-model approximation

    Some ultrathin semiconductor and metal nanofilms are semitransparent at the range of UV-Visible and near infrared wavelength, the steady and excited-state absorption spectra of them can be obtained after the combination with reflection and transmission spectra[1]. Optical transmission is the characteristics of a semitransparent sample to permit the passage of light, with the incident light being absorbed partly. In the past twenty years, semitransparent metal[2-4]or semiconductor nanofilms[5]is applied extensively in photoelectric devices for their excellent photoelectric and thermal properties, and these properties of semiconductor nanofilms take an important role in determining the heat dissipation and response capability of micro-structure devices. With the development of molecular beam epitaxy and magnetron sputtering techniques, semiconductor, and metal nanofilms can be prepared with high smoothness. Relaxation mechanism on free carriers between the electron-electron[6]and electron-lattice[7-8]remains a hot topic of solid-state physics. In semiconductor or metal thin films, the time-resolved measurements of the ultrafast dynamics show important temperature dependence on the electron-lattice coupling, which in the nonequilibrium transient heat transfer process of the highly excited electron-lattice system[9]. Due to the high density of electrons in metals, the initial stages of electron relaxation are connected with electron-electron interactions and their electronic effects, which can be illuminated by relatively simple models[10].

    Extensive investigations on the ultrafast dynamics of gold[11]and copper films[12]have been conducted previously. Hopkins and Norris[13]studied in the nonequilibrium heat transfer of Au film, the substrate of the gold film is effected by the electron-lattice coupling, it must exist the energy loss and transferred from the electron system to the metal lattice. The 5 nm thick gold film on BK7 substrate as a sample, and the ultrafast electron dynamics of it have been investigated as a function of film morphology, through the transient reflectivity measurements results, a sub-100 fs rise time and a complex relaxation profile on picosecond timescale can be observed[14]. Guo, et al.[15]investigated the interfaces heat transfer between metal and dielectric after the ultrafast-laser acting by transient reflection technique. The temperature nonequilibrium between electrons and lattices exist in gold film obviously after utilizing ultrafast-laser heating, transient thermo-reflectance measurement using femtosecond laser pulses on Au-Si samples and the simulation based on a two-temperature model are compared with the measured data. The results illustrate the direct coupling between electrons in metal and lattice in dielectric is strong when the contact resistance is sufficiently small. A molecular dynamics method based on TTM for predicting and numerical modeling the thermal resistance between metal and nonmetal interface[16]. The effective nonlinear optical absorption coefficient is measured for 20-nm-thick Au films as a function of pulse width by use of Z-scan measurement. Differential transmissivity and reflectivity measurements are performed at 775 nm pump and 630 nm probe pulses. By simulating the experiments which based on TTM and the analysis results of electrons and lattice can be obtained[17]. However, the detailed interaction mechanism between film and substrate is still in debate. Comparison between transmission and reflection experiments of metal nanofilms based on the change of dielectric function can provide abundant information to understand substrate effect of metal nanofilms.

    Two typical metal nanofilms and measure the ultrafast dynamics using transient reflection and transmission technique is choosed under the same experimental conditions, respectively. TTM is used to simulate the temperature of electron and lattice at the surface of nanofilms. Time dependence on electron and lattice temperatures shows an obvious deviation between transient reflection and transmission results. Crude-model approximation is used to explain the noticeable difference of dynamic profiles between transient reflectivity and transmissivity measurements. Fitting of experimental data shows combining the results of predicted lifetimes with the observation are reasonably close, and the time dependence of the relaxation rate cannot be neglected when analyzing a systematic deviation between transient reflectivity and transmissivity measurements. Considering interface thermal resistance and gradient change of temperature along with the direction of film thickness, transient reflectivity and transmissivity can be derived, and the equilibrium state of both electrons and lattice could be affected by the substrate.

    1 Experimental section

    1.1 Gold and copper nanofilms preparation and characterization

    BK7 glass substrates with the thickness of 170 μm are purchased from Agar Scientific Ltd. The substrates are first cleaned with acetone, ethanol then soaked in a detergent solution, and finally rinsed with de-ionized water and dried with nitrogen. Gold and copper nanofilms are fabricated by electron-beam evaporation deposition technique[18], which is the most consistent method. The gold film (60.0 nm) and Cu film (57.5 nm) are prepared by evaporating 99.999% pure gold and copper onto the substrates respectively. The thickness of metal nanofilms is measured by use of atomic force microscope (Nanosurf-Easyscan). The substrates are not intentionally heated by pump laser pulse. The evaporation rates and the thicknesses are monitored by an IC/5 thin film deposition controller (INFICON) with a 6 MHz Au-coated crystal[19]. For gold nanofilm, the obvious absorption peak locates at 510 nm, and the thickness of 60 nm can be evaluated as shown in Fig.1(a). All the measurements are performed quickly after the evaporation of thin film in order to decrease the oxidation of thin film.

    1.2 Experimental section

    The home-built transient reflection/transmission setup is demonstrated as shown in Fig.1(b).The source of ultrashort laser pulses consists of a mode-locked Ti: sapphire regenerative amplifier oscillator generating 120 fs pulses at the wavelength of 800 nm (Coherent, Legend), with pulse energy close to 1 mJ at a repetition rate of 1 kHz. The ultrafast pulses are used in collinear pump-probe geometry, where the beams are frequency doubled at 400 nm using a nonlinear beta-barium borate (BBO) crystal. The collinear pump and probe pulses are focused onto the sample surface to a spot size of 100 μm and 80 μm respectively. The radii of the pump and probe beams are measured with a sweeping knife edge method. The reflection and transmission data collected with a photo detector is locked into the pump modulation frequency to get the induced reflectivity ΔR/Rand transmissivity ΔT/Tat the same wavelength using a time-delayed probe pulse. Though 400 nm excitation light penetrates depth of gold and copper are 16 nm and 13 nm respectively, nanoscale gold film possess nonlinear absorption after 400 nm femtosecond laser excitation through the focal lens[17]. So, the substrate and boundary effect should be considered in the process of energy transfer between electrons and surrounding material systems.

    Fig.1 (a) Absorption spectrum of translucent gold film, insets are the morphology of cross section and surface of golf nanofilms with the thickness of 60.0 nm, respectively; (b) Schematic diagram of transient reflection and transmission experimental setup

    2 Results and discussions

    2.1 Two-temperature model and numerical analysis

    TTM can be solved numerically to predict the time dependence of electron temperature and lattice temperature[20-24], which commonly describes the evolution of the temperature changes due to the absorption of a laser pulse within the solid in spatial and temporal domain. TTM has been applied to simulate the process of heat transfer in metal thin films by femtosecond laser pulses[25-26]. According to TTM, the evolution of electron temperature Te and lattice temperature TL can be expressed as:

    g[TL(x,t)-Te(x,t)]+Q(x,t)

    (1)

    (2)

    where the specific heat capacity of the electronCeis linearity with the electron temperature,CLrepresents the specific heat capacity of the lattice and electronic heat conductivityκe, which is expressed byκeTe/TLandκeis the equilibrium electron thermal conductivity measured at room temperature. Temperature dependence on metal films is shown to be inaccurate for electron temperature beyond 1×104K, but at the low pump fluence the electron-lattice coupling factor is a constant. Detailed parameters are listed in Table 1.

    Table 1 TTM parameters for gold and copper nanofilms

    The pump laser pulse is assumed to be Gaussian in shape with a pulse width of 0.2 ps. Eq.(1) and (2) can be numerically solved as shown in Fig.2, which shows depth dependence on the electron and lattice temperatures at a fixed value of the electron-lattice coupling constant. Because of the electron heat capacity is less than lattice heat capacity, electrons are heated to a high temperature in a short time. We can see that the temperature of electron rises rapidly within the range of 0~500 fs at the film surface (x=0), and the maximum temperature reaches 2 240 K. When an ultrafast laser pulse heats the electrons in the metal, and for a short time (<1 ps), the electrons are driven out of equilibrium with the host lattice because the lattice temperature cannot rise too high in a short time to maintain equilibrium. Consequently, the electron temperature rises far above that of the lattice. There occurs non-equilibrium thermal phenomenon of electron and lattice due to the presence of the temperature difference between the electron and the lattice. Subsequent to the excitation of the electrons which they transfer energy to the lattice through the lattice scattering, and the electron and lattice temperatures equilibrium is the rate limiting step by the strength of the electron-lattice interaction.

    Fig.2 (a)Electron temperature gradient change of 60.0 nm thickness gold nanofilm relates to the time and the depth along with film thickness direction; (b)Lattice temperature of 60.0 nm thickness gold nanofilm connected with the depth along with the film thickness direction and the time.

    In Eq.(1),Q(x,t) describes the absorbed energy, which has a Gaussian temporal and spatial shape alongtandxwith the film depth coordinate:

    (3)

    where the fluenceFand pulse widthtpare parameters of incident laser pulse, the reflectivityRand transmissivityTrelate to thickness and property of substrate material, the magnitude ofRis roughly eight times than that ofT, and d denotes film thickness. Moreover,xis the distance from the inside of the film to the surface, andtdenotes the time delay after the laser reaches the surface of the material. The reflectivityRand transmissivityTat normal incidence on plane-parallel surfaces absorbing film bounded by air and glass substrate which both are transparent and non-absorbing media are given[27]by,

    (4)

    T=

    (5)

    whereε=ε1+iε2is the dielectric function of the metal film,ε1andε2represent the real and imaginary parts of dielectric function respectively. In addition,εs=ε1s+iε2sis the dielectric function of BK7 glass substrate,εwhich is transparent for the 800 nm spectral region and thenε2s=0.α1is the absorption coefficient of the metal film. From the equation (5), we can see that the transmissivityTis related to the dielectric function of BK7 glass substrateεs. Three kinds of coupling exist at the interface which between the gold nanofilm and substrate as follows: The electron and lattice coupling within gold nanofilm, the electron in gold nanofilm coupling with the lattice in substrate, and the lattice in gold nanofilm coupling with the lattice in substrate, and the lattice-lattice coupling give rise to the resistance on the interface. The main reasons which causing the decrease of the temperature are the electron-lattice coupling and heat diffusion, which is due to the existence of temperature gradient that energy spreads by diffusion in two materials here. Moreover, the reason why there exists heat diffusion is the temperature difference between gold nanofilm and substrate, for BK7 is dielectric, there are no free electrons in it, when the electron-lattice coupling reaches equilibrium, the temperature in gold nanofilm is different with it on the substrate surface all the time, so a part of electrons which have reached the equilibrium coupling with lattice in substrate and the temperature decreased. The dielectric constants at specific wavelength are obtained from the database. The reflection and transmission change are related to the alteration of the dielectric functionε[28]as,

    ΔR/R=(Δε1?R/?ε1+Δε2?R/?ε2)/R

    (6)

    ΔT/T=(Δε1?T/?ε1+Δε2?T/?ε2)/T

    (7)

    For a thin film, where ?R/?ε1, ?R/?ε2, ?T/?ε1, and ?T/?ε2can be derived by Eq.(4)-(5) respectively. Δε1and Δε2are the change of the real and imaginary part of dielectric function, which can be calculated by integrating all the available transitions from d-band to conduction band, located atLpoint of the band structure by Rosei’s dielectric function model[29].The difference between reflection and transmission can be discussed theoretically. For opaque film, the effect of substrate needed to be considered due to the change of refraction index.

    2.2 Gold and copper nanofilms using TTM analysis

    The lowest energy d-band to available Fermi level transition is small for gold (2.40 eV) and copper (2.15 eV), it’s necessary to consider the effects of interband and intraband transitions on the transient reflection and transmission due to the excitation photon energy of 3.1 eV[30]. When the excitation photon energy is larger than interband transition threshold, electrons in the d band can also absorb the energy and transit over Fermi energy. Here, the output at 800 nm (1.55 eV) from the laser which is below interband transition threshold, the excitation fluence is low enough to make the number of s/p band is much more than that of d-band electrons, under this condition the effect of the d band perturbation to the transient reflectance is ignorable[30].

    Electron relaxation in metal exists at the first two picoseconds, the transmission and reflection results are similar at the first two picoseconds as reported, so inter- or intraband transition of electrons is not directly due to the pump pulses after the initial fast transient.

    (8)

    where e-(t/τ)2is the instrument response function, ΔTeand ΔTLare the variation of electron and lattice temperature, respectively. According to the literature[30], at early times after 150 fs laser excitation, the normalized reflectivity change ΔR/Ror the normalized transmission change ΔT/Tis proportional to the normalized temperature change of the electron gas and lattice, hence, the surface reflectivity kinetics is obtained and assigned to the surface temperature variations. But for ΔT/T, probe and pump pulses transmitted through the glass substrate, thermal diffusion and interface resistance should be considered. The classical TTM is hard to estimate the dynamics measured by transient transmission method. With decreasing film thickness to the order of the thermal penetration depth, the boundary effect is essential to be solved during the electron-lattice coupling equilibrium process. The effect of substrate is negligible for TTM. From TTM, the thickness dependence has been considered but substrate dependence has not been considered. Through modify electron-lattice constant to optimize the measured dynamics with least square errors from TTM fitting. The magnitude of ΔTis roughly three times higher than that of ΔR. Consequently, the total change in the absorption for the Au film is about three times than previously obtained from the analysis of the transmissivity.

    At the stage of rise time and longer delay time, the experimental results are not satisfied with the simulation of TTM because of the limitation of time resolution of the setup, the coevolution is applied to fitting. The initial stage of eight picoseconds is discussed with the classical TTM, which is shown above in Fig.3. The agreement between the fitted results and the experimental data is generally favorable from zero to eight picoseconds. The small discrepancy between the experimental and the simulation results should stem from substrate effect and thermal diffusion. The change rate of gold film reflectivity is slower than the change rate of the copper film reflectivity. Firstly, because the electron-lattice coupling coefficient of gold is 2.3×1016Wm-3K-1, under the same laser energy density and film thickness condition, the electron-lattice coupling coefficient of copper is 1.0×1017Wm-3K-1, and the electron-lattice coupling coefficient is used to represent the energy transfer rate from electron to lattice. The larger the electron lattice coupling coefficient, the shorter time of the electron and the lattice to reach the same temperature, and the temperature change faster. Secondly, the thermal conductivity of copper is larger than that of gold at room temperature, the heat conductivity capability of copper is greater than that of gold, so the temperature of copper changes faster than that of gold.

    Fig.3 (a) Time-dependent differential reflection curves for Au nanofilm pumped with pulses of 200 fs, 400 nm and energy from 50 to 110 nJ; (b) Time-dependent differential reflection curves for Cu nanofilms pumped with pulses of 200 fs, 400 nm pulse with pulse and energy from 50 to 110 nJ

    2.3 Electron and lattice dynamics using Crude-model approximation

    In order to investigate the influence of the data acquisition mode on ultrafast dynamics, we chose gold and copper semi-transparent metal nanofilms to measure the transient signal. For most semiconductor at high laser excitation, the local dielectric function is modified and creates a discontinuity at their surface. So, we can observe the change of the sign of transient reflectivity. However, for gold and copper thin film we measure the reflectivity and transmissivity signals using transient reflection and transmission technique at the same experimental conditions, the sign of transient transmission results varies with the delay time. Transient transmissivity is pushed down compared with transient reflectivity.

    Since the change of the transmission, the reflection, and the absorption must sum to unity, the change in the total absorbance ΔAequals-(ΔT+ΔR).

    (9)

    whereI0is the intensity of incident laser pulse, ΔTand ΔRare the intensity of transmission and reflection under steady-state condition. For Au and Cu nanofilms the ratio of intensity between reflectivity and transmissivity is about 20 at steady-state measurement. Comparing the transient reflectivity changes with transient transmissivity changes shown in Fig.4, the time-resolved experiments reveal that both ΔTand ΔRare negative, leading to an increase of the absorption, but they are different for the profile of ultrafast dynamics. The transient reflectivity and transmission data of gold and copper nanofilms are presented for the case of incident laser pulse energy from 50 to 110 nJ. The Crude-model approximation which separates the electron-gas fluctuation into thermalized and non-thermalized components, and a further analysis of the transient reflectivity and transmissivity changes which can be observed[9]. In all cases the data show a sub-picosecond rise time and a picosecond timescale relaxation. These data are initially fit with the function[19, 20]as,

    e-(t-t′)/τep+B(1-e-(t-t′)/τep)]e-(t′/τ)2dt′

    (10)

    whereτrrepresents the sub-picosecond rise time and it is related to the relaxation dynamics of the non-thermal initial distributions, which giving a general description of the build-up of the Fermi distribution temperature through energy transfer. Andτepreflects the electron-lattice relaxation time related to thermal equilibrium.τis the instrumental response function, the value is about 800 fs. The final rising component is introduced to reflect the long-lived offset[9, 31]. This function gives an adequate fit to the data with rise times of 0.8 ps and about 1.1 ps for gold and copper nanofilms, respectively. With the energy of pump laser pulse increasing, rise time and transient transmission data are comparable with that of transient reflection data. Although the energy dependence on the electron relaxation can be ignored in this phenomenological approach, it provides a simple and qualitative description of the electron dynamics. This model results are almost agreement with the experimental results without considering small discrepancies. A better fitting is achieved by adding an additional term, introducing a second component in the electron-lattice relaxation as shown in Fig.4. In all cases the data show a sub-picosecond rise time and a picosecond timescale relaxation. The influence of the thermal resistance of the glass substrate also reduces the temperature of the backside of the gold film. In the ultra-fast process testing of metallic semipermeable membrane materials, the measurement of transient transmission or reflection may bring about experimental deviation, and transient transmission and reflection signals should be measured at the same time. Transmission data of them is semblable by comparing the reflection data of Au and Cu nanofilm at the first ~2 ps time range, but after ~2 ps, the change of sign for transmission dynamics exists for Cu and Au nanofilms at various pump laser excitations. The electron-lattice interaction of Au (or Cu) and lattice-dielectric interaction between Au (or Cu) and glass should be attributed to the difference between ΔR/Rand ΔT/T. At first, the electron-lattice and lattice-dielectric interaction happens at ~2 ps after the pump laser excitation. With the pump laser excitation, the increment exterior and interior temperature of metal film is different, the influence of gradient temperature is dominant to ΔT/Tbecause of the penetration of the whole film. We can also get pump power dependence of rise time and electron-lattice relaxation time.

    Fig.4 Transient reflectivity and transmission data and fitting for gold nanofilm (a) and copper nanofilm, (b) recorded with different energies of pump laser pulse (50, 70, 90, 110 nJ)

    Under the nonequilibrium thermal condition, the electron-lattice coupling of metal-nonmetal interface cannot be ignored for transient transmission. The higher the temperature is, the larger the electron-lattice coupling resistance[32]. From Fig.2(a), we can see that the highest electron temperature is 2 240 K, therefore, the electron-lattice coupling resistance is also large. Considering the electron-lattice coupling resistance, the total thermal conductivityκdecreases and the relaxation time of electron-latticeτepincreases, so with the laser energy density becomes higher, the recovery time to equilibrium is longer. Compared with the unexcited initial state, the total thermal conductivity of the electron-lattice system decreases and overcomes the larger electron-lattice coupling resistance, the temperature is lower than the initial state.

    The combination between transient reflectivity and transmissivity for Au and Cu nanofilms are fitted using Crude model as shown in Fig.5.As the size of the devices drops to microscale and nanoscale, the role of Au-glass interface becomes increasingly important in determining the thermal and electrical performance. Electron-lattice nonequilibrium near the interface leads to lattice generation, which can disturb the original lattice scattering, and hence affects lattice transport, we would also consider the penetration of electrons into the nonmetal side. Reflectivity is a function of the surface, and transmission is a function affected by the interface, so we can see that transmissivity is more susceptible to the substrate than reflectivity. A linear increasing from time zero to rough 2.5 picoseconds, then keeping a constant of 0.3 for gold nanofilm and 0.4 for copper nanofilm. At first two picoseconds we can see no obvious change of relaxation rate between the transient and transmission signals because the effect of electron-lattice coupling exists. However, after two picoseconds, lattice thermal diffusion and interface thermal resistance should be considered, in addition, an account of the interaction process might be associated with interface thermal resistance and has a relaxation time with tens of picoseconds[33].

    Fig.5 Combination between transient reflectivity and transmissivity for Au (a) and Cu (b) nanofilms at the same experimental conditions

    The rising time and electron-lattice relaxation time of gold and copper nanofilms are shown in Fig.6. Because of the thicknesses is on the order of the thermal penetration depth, the result of transient reflection and transmission measurements are influenced by nonequilibrium between the electrons in the metal films and the lattices at the boundary of metal and substrate, and this phenomenon caused by electron energy loss to both the lattices in the metal film and those in the underlying substrate. Using Crude-model the electron-lattice coupling factor can be obtained and the excess energy loss from the electron system is observed for ~60 nm gold and copper nanofilms on glass substrates[9, 14]. With the pump pulse energy increasing, electron-lattice relaxation time of gold and copper nanofilms have a weak increment. Non-thermal rise time of transient transmission result is shorter than that of transient reflection result under the same experimental conditions. Lattice thermal diffusion and interface thermal resistance should be considered because the interaction process might be associated with interface thermal resistance and has a relaxation time with tens of picoseconds[34].

    Fig.6 Rise time and electron-lattice relaxation time under various pump laser excitation based on therecombination between transient reflectivity and transmissivity for Au and Cu nanofilms

    In the high temperature limit[13]:

    (11)

    where me is the effective mass, s is the speed of sound, ne is the electron density, andτis the relaxation time[35]. At high ultrafast laser excitation, electron-lattice relaxation time is not a constant; Using the Debye speed of sound,s=2.65×103ms-1for Cu ands=1.40×103ms-1for Au are calculated. From Eq.(11) we can plot the relationship between pump pulse energy and electron-lattice relaxation time. Energy transfer mechanism of opaque nanofilms during thermal nonequilibrium states between electrons and the surrounding material systems has a significant application to semiconductor devices[16].

    3 Conclusions and perspectives

    In summary, we investigate the ultrafast electron and lattice dynamics in nanofilms of gold and copper deposited on a glass substrate using transient reflectivity and transmissivity method. The experimental data are further analyzed in terms of TTM and Crude-model approximation which accounts for the relaxation of the initial optically excited and electron-lattice interaction, and we show the reason why the sign for transmission changes by analyzing the substrate effect. From these experimental and fitting results, the sub-picosecond relaxation of the optically excited electrons shows a relatively weak dependence on pump pulse energy, while the electron-lattice relaxation time which governs the lifetime of the thermalized electrons depends on the pump pulse energy. Electron-lattice coupling effect in transmission method is stronger and more sensitive than that of reflection method under the same experimental conditions. Therefore, through the transient transmission results, the effect of substrate should be considered more in the view of electron-lattice coupling and interface thermal resistance. These conclusions are being useful for investigating relaxation process of opaque nanofilms further.

    久久久久久久精品吃奶| 久久久久久久久中文| 成人高潮视频无遮挡免费网站| 免费在线观看影片大全网站| 丰满的人妻完整版| 99riav亚洲国产免费| 日本 欧美在线| 一边摸一边抽搐一进一小说| 首页视频小说图片口味搜索| 欧美绝顶高潮抽搐喷水| 女人高潮潮喷娇喘18禁视频| 观看美女的网站| 成人性生交大片免费视频hd| 国产伦精品一区二区三区四那| 精品国产超薄肉色丝袜足j| 美女免费视频网站| 免费看美女性在线毛片视频| 免费在线观看成人毛片| 久久久久免费精品人妻一区二区| 大型黄色视频在线免费观看| 床上黄色一级片| 亚洲av免费高清在线观看| 成人亚洲精品av一区二区| 国产伦一二天堂av在线观看| 亚洲欧美日韩东京热| 人妻丰满熟妇av一区二区三区| 国产三级黄色录像| 国产精品三级大全| 桃红色精品国产亚洲av| 天堂网av新在线| 精品一区二区三区视频在线观看免费| 欧美黑人欧美精品刺激| 在线观看免费午夜福利视频| netflix在线观看网站| 亚洲 国产 在线| 久久香蕉国产精品| 亚洲欧美激情综合另类| 中文字幕人妻熟人妻熟丝袜美 | 国产成人系列免费观看| 欧美在线黄色| 色精品久久人妻99蜜桃| 日本在线视频免费播放| 日本成人三级电影网站| 99精品欧美一区二区三区四区| 俺也久久电影网| 国产精品 欧美亚洲| 老熟妇乱子伦视频在线观看| 亚洲激情在线av| 亚洲精华国产精华精| 国产久久久一区二区三区| 欧美不卡视频在线免费观看| 国产不卡一卡二| 亚洲欧美日韩卡通动漫| 中文在线观看免费www的网站| 三级男女做爰猛烈吃奶摸视频| 国产精品av视频在线免费观看| 亚洲av熟女| 色av中文字幕| 成人国产综合亚洲| 高清日韩中文字幕在线| 中亚洲国语对白在线视频| 久久伊人香网站| 国语自产精品视频在线第100页| 久久精品夜夜夜夜夜久久蜜豆| 91久久精品国产一区二区成人 | 亚洲精品亚洲一区二区| 国产伦人伦偷精品视频| svipshipincom国产片| 国产黄片美女视频| 五月玫瑰六月丁香| tocl精华| 亚洲午夜理论影院| 国产一区二区在线观看日韩 | 国内久久婷婷六月综合欲色啪| 亚洲一区高清亚洲精品| 99久久精品热视频| 国产私拍福利视频在线观看| 真人一进一出gif抽搐免费| 免费观看精品视频网站| 日本五十路高清| 深爱激情五月婷婷| 99热这里只有精品一区| 日本成人三级电影网站| 国产一区二区在线av高清观看| 在线观看日韩欧美| 亚洲成av人片在线播放无| 亚洲成人精品中文字幕电影| 欧美日韩乱码在线| 亚洲欧美精品综合久久99| 国产精品一区二区三区四区免费观看 | 99国产精品一区二区三区| 此物有八面人人有两片| а√天堂www在线а√下载| 精品午夜福利视频在线观看一区| 国产视频一区二区在线看| 黄片小视频在线播放| 美女cb高潮喷水在线观看| 日本一本二区三区精品| 97碰自拍视频| 亚洲黑人精品在线| 精品人妻一区二区三区麻豆 | 男女视频在线观看网站免费| 日本在线视频免费播放| 成人一区二区视频在线观看| 啦啦啦观看免费观看视频高清| 中文字幕人妻丝袜一区二区| 国产欧美日韩一区二区三| 久久久久久久久久黄片| a级毛片a级免费在线| www.www免费av| 欧美日韩亚洲国产一区二区在线观看| 91九色精品人成在线观看| 老熟妇仑乱视频hdxx| 丰满人妻一区二区三区视频av | 亚洲人与动物交配视频| 精品一区二区三区av网在线观看| 日韩国内少妇激情av| 国产亚洲精品久久久久久毛片| 动漫黄色视频在线观看| 亚洲国产高清在线一区二区三| 国产野战对白在线观看| 亚洲国产欧美网| 亚洲av五月六月丁香网| 国产高清视频在线播放一区| 男女下面进入的视频免费午夜| 99国产极品粉嫩在线观看| 一进一出好大好爽视频| 18美女黄网站色大片免费观看| 婷婷丁香在线五月| 国产精品国产高清国产av| 日韩成人在线观看一区二区三区| 亚洲熟妇熟女久久| 免费人成视频x8x8入口观看| 少妇人妻精品综合一区二区 | 国产私拍福利视频在线观看| 观看美女的网站| 精品熟女少妇八av免费久了| 国产久久久一区二区三区| 国产精品野战在线观看| 亚洲中文字幕日韩| 一进一出好大好爽视频| av视频在线观看入口| 国产日本99.免费观看| 亚洲国产中文字幕在线视频| netflix在线观看网站| 在线a可以看的网站| 欧美在线黄色| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 一本综合久久免费| 三级国产精品欧美在线观看| 亚洲,欧美精品.| 国产精品嫩草影院av在线观看 | 成年女人永久免费观看视频| 天天躁日日操中文字幕| 很黄的视频免费| 看免费av毛片| 欧美最黄视频在线播放免费| 欧美日韩瑟瑟在线播放| 一级黄色大片毛片| 人人妻,人人澡人人爽秒播| 国产精品一区二区三区四区久久| 国产精品女同一区二区软件 | 国产亚洲精品久久久久久毛片| 成年女人看的毛片在线观看| 亚洲国产精品999在线| 亚洲熟妇中文字幕五十中出| 国产97色在线日韩免费| 亚洲第一电影网av| 好看av亚洲va欧美ⅴa在| 国内精品美女久久久久久| 久久精品人妻少妇| 国产不卡一卡二| 成人特级av手机在线观看| 母亲3免费完整高清在线观看| 最近在线观看免费完整版| 色av中文字幕| 国产亚洲精品一区二区www| 欧美一区二区精品小视频在线| 美女大奶头视频| 久久久久久国产a免费观看| 久久久色成人| 丰满人妻熟妇乱又伦精品不卡| 老司机在亚洲福利影院| 久久久久久九九精品二区国产| 精华霜和精华液先用哪个| 国产精品三级大全| 日韩高清综合在线| 最新在线观看一区二区三区| 韩国av一区二区三区四区| 国产午夜精品久久久久久一区二区三区 | 夜夜躁狠狠躁天天躁| 免费看光身美女| 91av网一区二区| 久久久久久久精品吃奶| 99热6这里只有精品| 国产又黄又爽又无遮挡在线| 亚洲国产日韩欧美精品在线观看 | 女人十人毛片免费观看3o分钟| 免费观看的影片在线观看| 国产日本99.免费观看| 免费观看精品视频网站| 亚洲成人久久性| 国产精品一区二区三区四区久久| 久久6这里有精品| 精品无人区乱码1区二区| 国产在视频线在精品| 一本综合久久免费| 日韩欧美 国产精品| 国产欧美日韩精品亚洲av| 亚洲av免费高清在线观看| 黄色丝袜av网址大全| 国产黄片美女视频| 国产私拍福利视频在线观看| 99久久无色码亚洲精品果冻| 国产精品久久久久久久久免 | 美女 人体艺术 gogo| 中文在线观看免费www的网站| 国产极品精品免费视频能看的| 国产一区二区激情短视频| 免费观看的影片在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲人成网站在线播| 欧美性猛交╳xxx乱大交人| 亚洲七黄色美女视频| 狂野欧美白嫩少妇大欣赏| 国产精品电影一区二区三区| 成人特级黄色片久久久久久久| 精品人妻1区二区| 高清日韩中文字幕在线| 国模一区二区三区四区视频| 黄色丝袜av网址大全| 夜夜看夜夜爽夜夜摸| 18禁国产床啪视频网站| 一个人免费在线观看电影| 亚洲人成网站在线播| 岛国在线观看网站| 国产高潮美女av| 中文字幕人妻丝袜一区二区| 久久九九热精品免费| 蜜桃久久精品国产亚洲av| 在线国产一区二区在线| 国产高潮美女av| 亚洲精品国产精品久久久不卡| 欧美成人一区二区免费高清观看| 日韩欧美在线二视频| 三级男女做爰猛烈吃奶摸视频| 国产私拍福利视频在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品美女特级片免费视频播放器| 丁香六月欧美| 97碰自拍视频| 免费看美女性在线毛片视频| 亚洲精品456在线播放app | 搡老妇女老女人老熟妇| 久久精品91无色码中文字幕| 国产亚洲欧美98| 看片在线看免费视频| 久久久久亚洲av毛片大全| 国产精品自产拍在线观看55亚洲| 亚洲五月婷婷丁香| 午夜精品一区二区三区免费看| 波多野结衣高清无吗| 国产精品久久久久久人妻精品电影| 九色国产91popny在线| 午夜福利在线在线| 欧美一区二区亚洲| 狂野欧美激情性xxxx| 国产aⅴ精品一区二区三区波| 色av中文字幕| a级一级毛片免费在线观看| 日韩中文字幕欧美一区二区| avwww免费| 欧美黄色片欧美黄色片| 午夜福利在线观看吧| 一本久久中文字幕| 中文字幕久久专区| av中文乱码字幕在线| 成人午夜高清在线视频| 日本三级黄在线观看| 日本黄色视频三级网站网址| 99精品欧美一区二区三区四区| 99精品久久久久人妻精品| 国产精品永久免费网站| 国产伦在线观看视频一区| 无遮挡黄片免费观看| 啪啪无遮挡十八禁网站| 最近最新中文字幕大全电影3| 精品一区二区三区视频在线观看免费| 99精品欧美一区二区三区四区| 欧美丝袜亚洲另类 | 久久天躁狠狠躁夜夜2o2o| 国产在视频线在精品| 亚洲国产精品合色在线| 深夜精品福利| 国产成人aa在线观看| a级一级毛片免费在线观看| 久久精品国产自在天天线| 国产精品 国内视频| 国内久久婷婷六月综合欲色啪| 在线观看av片永久免费下载| 露出奶头的视频| 在线免费观看不下载黄p国产 | 嫁个100分男人电影在线观看| 一本精品99久久精品77| 不卡一级毛片| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合久久99| 欧美日韩综合久久久久久 | av片东京热男人的天堂| 无限看片的www在线观看| 国产蜜桃级精品一区二区三区| 天天躁日日操中文字幕| 怎么达到女性高潮| 精品人妻偷拍中文字幕| 天天添夜夜摸| 操出白浆在线播放| 亚洲欧美日韩东京热| 亚洲av免费在线观看| 99在线视频只有这里精品首页| 一区二区三区国产精品乱码| 女人高潮潮喷娇喘18禁视频| 国产黄色小视频在线观看| 很黄的视频免费| 免费av观看视频| 最好的美女福利视频网| 久久精品91无色码中文字幕| 国产在线精品亚洲第一网站| 日韩精品中文字幕看吧| 热99在线观看视频| 可以在线观看的亚洲视频| www日本在线高清视频| 美女黄网站色视频| 成年免费大片在线观看| 99久久九九国产精品国产免费| 美女高潮喷水抽搐中文字幕| 狂野欧美激情性xxxx| 亚洲最大成人中文| 香蕉av资源在线| 国产单亲对白刺激| 变态另类成人亚洲欧美熟女| www.999成人在线观看| 蜜桃亚洲精品一区二区三区| 一级黄色大片毛片| 亚洲av美国av| 18禁黄网站禁片午夜丰满| 婷婷亚洲欧美| 久久久精品欧美日韩精品| 亚洲天堂国产精品一区在线| 亚洲 国产 在线| 日韩高清综合在线| 麻豆一二三区av精品| 日韩人妻高清精品专区| 最近最新中文字幕大全免费视频| 免费一级毛片在线播放高清视频| eeuss影院久久| 午夜老司机福利剧场| 午夜视频国产福利| 18禁黄网站禁片免费观看直播| 在线观看免费午夜福利视频| 又黄又爽又免费观看的视频| 深夜精品福利| 亚洲av电影不卡..在线观看| 欧美三级亚洲精品| 免费电影在线观看免费观看| 国产成人aa在线观看| 免费av不卡在线播放| a级毛片a级免费在线| svipshipincom国产片| 日韩欧美 国产精品| 久久久久国内视频| 脱女人内裤的视频| 2021天堂中文幕一二区在线观| 美女cb高潮喷水在线观看| 又黄又粗又硬又大视频| 国产伦精品一区二区三区视频9 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 美女被艹到高潮喷水动态| 日韩欧美精品v在线| 免费观看的影片在线观看| 99久久九九国产精品国产免费| 色综合站精品国产| 美女高潮的动态| 手机成人av网站| 天堂√8在线中文| 中文字幕人妻熟人妻熟丝袜美 | 国产麻豆成人av免费视频| 白带黄色成豆腐渣| 久久久久久国产a免费观看| 亚洲成人久久性| 亚洲国产欧美人成| 久久草成人影院| 成人国产一区最新在线观看| 母亲3免费完整高清在线观看| 国产熟女xx| 久久久国产精品麻豆| 久久人妻av系列| 亚洲专区中文字幕在线| 久久久久国内视频| 无限看片的www在线观看| 岛国视频午夜一区免费看| 亚洲欧美精品综合久久99| 国产精品98久久久久久宅男小说| 久久久精品欧美日韩精品| 天堂网av新在线| 欧美日韩黄片免| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久久黄片| 亚洲一区高清亚洲精品| 久久香蕉精品热| 18禁在线播放成人免费| 欧洲精品卡2卡3卡4卡5卡区| 日本免费一区二区三区高清不卡| 男女视频在线观看网站免费| 99精品在免费线老司机午夜| 精品国产美女av久久久久小说| 欧美色欧美亚洲另类二区| 999久久久精品免费观看国产| 亚洲内射少妇av| 欧美日韩亚洲国产一区二区在线观看| 欧美一区二区精品小视频在线| 久久久久久久精品吃奶| 色综合婷婷激情| 国产精品免费一区二区三区在线| 在线免费观看不下载黄p国产 | 国产av不卡久久| 欧美国产日韩亚洲一区| 亚洲乱码一区二区免费版| 嫩草影院入口| 一个人看视频在线观看www免费 | 成人亚洲精品av一区二区| 免费电影在线观看免费观看| 国产精品久久久久久精品电影| 国产欧美日韩一区二区精品| 国产精品98久久久久久宅男小说| 日本a在线网址| 在线播放无遮挡| 国产私拍福利视频在线观看| 色综合欧美亚洲国产小说| 人人妻人人澡欧美一区二区| 欧美色视频一区免费| av福利片在线观看| 日本五十路高清| 俺也久久电影网| 内射极品少妇av片p| 国产一区二区亚洲精品在线观看| 内地一区二区视频在线| 免费看a级黄色片| 亚洲五月婷婷丁香| 99久久九九国产精品国产免费| 俄罗斯特黄特色一大片| 天堂网av新在线| www.色视频.com| 欧美成人免费av一区二区三区| 午夜影院日韩av| 国产亚洲欧美98| 亚洲专区中文字幕在线| 色精品久久人妻99蜜桃| 成人永久免费在线观看视频| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 亚洲国产精品999在线| 在线观看日韩欧美| 九九在线视频观看精品| 观看美女的网站| 久久精品国产综合久久久| 99精品欧美一区二区三区四区| 欧美极品一区二区三区四区| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 日本熟妇午夜| 深爱激情五月婷婷| 国产午夜精品论理片| 免费高清视频大片| 午夜激情欧美在线| 欧美最新免费一区二区三区 | 精品乱码久久久久久99久播| 国产精品 欧美亚洲| 亚洲18禁久久av| 看黄色毛片网站| 两个人视频免费观看高清| 国产精品亚洲一级av第二区| 狂野欧美白嫩少妇大欣赏| 18禁黄网站禁片免费观看直播| 级片在线观看| 久久精品综合一区二区三区| 亚洲一区二区三区不卡视频| 在线观看午夜福利视频| 成人国产一区最新在线观看| 亚洲avbb在线观看| 长腿黑丝高跟| 日韩国内少妇激情av| 俺也久久电影网| 伊人久久大香线蕉亚洲五| 精品午夜福利视频在线观看一区| 亚洲精品色激情综合| 在线播放无遮挡| 夜夜爽天天搞| 久久久久久九九精品二区国产| 国产单亲对白刺激| 岛国视频午夜一区免费看| 国产亚洲精品综合一区在线观看| 国产精品久久久久久人妻精品电影| 国产午夜福利久久久久久| 免费在线观看成人毛片| 国产精品久久久久久精品电影| 丁香六月欧美| 欧美在线一区亚洲| 久久久久久人人人人人| 精品日产1卡2卡| 天堂动漫精品| 欧美乱色亚洲激情| 高清毛片免费观看视频网站| 人人妻人人澡欧美一区二区| 国产在视频线在精品| 露出奶头的视频| 欧美日韩精品网址| 亚洲人成电影免费在线| 国产一区二区在线观看日韩 | av欧美777| 99热精品在线国产| 国产高清有码在线观看视频| 不卡一级毛片| 中出人妻视频一区二区| 亚洲国产欧洲综合997久久,| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲av嫩草精品影院| 午夜福利18| 九色国产91popny在线| 叶爱在线成人免费视频播放| 精品不卡国产一区二区三区| 久99久视频精品免费| av欧美777| 国产伦精品一区二区三区四那| 亚洲国产精品成人综合色| 不卡一级毛片| 亚洲欧美精品综合久久99| 国产欧美日韩精品亚洲av| 99国产精品一区二区三区| 国产精品嫩草影院av在线观看 | 国产毛片a区久久久久| 亚洲欧美日韩高清专用| 欧美在线黄色| 99久久精品一区二区三区| 久久久久久久精品吃奶| 日本免费一区二区三区高清不卡| 日韩欧美在线二视频| 色综合站精品国产| 在线a可以看的网站| 熟女人妻精品中文字幕| 搡女人真爽免费视频火全软件 | 五月玫瑰六月丁香| 免费看光身美女| 麻豆国产av国片精品| 一a级毛片在线观看| 99久久精品一区二区三区| 欧美黑人巨大hd| 嫩草影视91久久| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 午夜两性在线视频| 狂野欧美激情性xxxx| 十八禁网站免费在线| 亚洲欧美日韩东京热| 亚洲国产色片| 热99在线观看视频| 久久精品国产综合久久久| 嫁个100分男人电影在线观看| 在线视频色国产色| 淫妇啪啪啪对白视频| 一夜夜www| 国产爱豆传媒在线观看| svipshipincom国产片| 国产精品精品国产色婷婷| 亚洲成人免费电影在线观看| 少妇的丰满在线观看| 免费高清视频大片| 亚洲av五月六月丁香网| 亚洲欧美激情综合另类| 亚洲av美国av| 成人国产综合亚洲| 51国产日韩欧美| 国产一区在线观看成人免费| 午夜久久久久精精品| 国产黄a三级三级三级人| 国产黄片美女视频| 一进一出抽搐动态| 免费av毛片视频| 久久精品亚洲精品国产色婷小说| 久久久色成人| 搡女人真爽免费视频火全软件 | 麻豆成人av在线观看| 麻豆国产av国片精品| 久久人人精品亚洲av| 国产精品一区二区免费欧美| 日本a在线网址| 日韩人妻高清精品专区| 日韩欧美国产在线观看| 亚洲av第一区精品v没综合| 久久久久久九九精品二区国产| 成人18禁在线播放| 欧美日韩乱码在线| 一进一出抽搐gif免费好疼| 久久久久精品国产欧美久久久| 国产一区二区在线av高清观看| 搞女人的毛片| 18禁裸乳无遮挡免费网站照片| 少妇的逼水好多| 观看免费一级毛片| 日韩 欧美 亚洲 中文字幕| 淫妇啪啪啪对白视频| 一个人免费在线观看电影| 国产精品亚洲一级av第二区| 好男人在线观看高清免费视频| 亚洲欧美日韩无卡精品|