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

    Theoretical and experimental study of phase optimization of tapping mode atomic force microscope

    2022-08-01 06:02:14ZhengWei魏征AnJiePeng彭安杰FengJiaoBin賓鳳姣YaXinChen陳亞鑫andRuiGuan關(guān)睿
    Chinese Physics B 2022年7期
    關(guān)鍵詞:魏征

    Zheng Wei(魏征), An-Jie Peng(彭安杰), Feng-Jiao Bin(賓鳳姣), Ya-Xin Chen(陳亞鑫), and Rui Guan(關(guān)睿)

    College of Mechanical and Electrical Engineering,Beijing University of Chemical Technology,Beijing 100029,China

    Keywords: TM-AFM,phase image,excitation frequency,energy dissipation

    1. Introduction

    Since its advent,the atomic force microscopes(AFM)has become one of the nanoscale surface metrology tools used most widely.[1,2]An atomic force microscope is similar to an atomic-scale phonograph. Its core component is a microcantilever(phonographic arm)with a sharp tip(phonographic stylus).[3]The surface information can be obtained by transforming the interaction force between the tip and a sample into the displacement of the micro-cantilever tip.[4]As one of the most commonly used operational modes in atomic force microscope,the tapping mode(TM)can provide high-resolution images of a variety of samples,as well as greatly reducing the wear of the tip and damages to the samples in comparison with contact modes.[5]Therefore, tapping-mode atomic force microscopes(TM-AFM)are widely used in biomolecules,polymers,and nanoscale structures.[6–16]

    As an element of TM-AFM, phase-imaging has greatly promoted the establishment of key nanoscale characterization methods of force microscopy in various material science fields.[4,17]By recording the phase shifts of microcantilever oscillations relative to the driving signals, the interaction characteristics of a heterogeneous sample can be plotted.[18]Phase-imaging has been used to study various material properties, such as capillary forces,[19]wettability,[20]and viscoelasticity.[21]In addition,owing to its wide and powerful applications, phase-imaging has consistently been a research hotspot in the relevant atomic force microscope.[22,23]

    Phase images are used as supplementary information for topography in TM-AFM.In operational processes,various parameters in the control panels are adjusted to ensure the integrity and contrasts of the phase images. However,there are currently few theories to guide the adjustments of those parameters. Therefore,it is significant to study the factors affecting phase-imaging theoretically and experimentally.[21,24–28]The dependence of phase contrast on elastic interaction and inelastic interaction in TM-AFM have been investigated.[21]Some results have revealed the relationships between the phase images and the amplitude values of the free oscillationA0and the setpoint amplitudeASP.[24–26]Fasolket al. found that the selecting of the appropriate temperature levels for the phase scanning according to the characteristics of the samples can improve the phase contrast results.[27]Furthermore, a large number of studies have shown that energy dissipation is the key factor affecting the phase imaging of TM-AFM.[29–32]

    Generally, there are various energy dissipations in TMAFM systems, including contact dissipation between the tip and sample, air viscous dissipation, squeeze film dissipation,and liquid bridge dissipation.[29,34–36]All of these dissipations will affect the phases of a system.In the present study,the contact dissipation between the tip and sample is of direct interest since it is considered to accurately reflect the properties of the samples. Meanwhile, the background dissipation masked the key dissipation, so that the real phase information of the sample is masked. Therefore, it is determined in this study that phase imaging can be successfully optimized by enhancing key dissipation or eliminating other irrelevant dissipations.

    In this work, starting from the energy dissipation of the TM-AFM system, a theoretical model of the key factors affecting the phases is established. Then, through the analysis of the theoretical model, two ways to optimize the phase imaging are identified. The effects of the excitation frequency on a two-component blend are first studied both theoretically and experimentally. In this study it is proven both theoretically and experimentally that there exists an optimal excitation frequency that maximizes the phase contrast. In addition,by eliminating the background dissipation via a self-excitation method, the system is able to only reflect the key dissipation between the tip and the sample. As a result,only the real sample information is obtained and the phase contrast is improved to a certain extent.Furthermore,it is considered to be most advantageous to optimize the phase imaging of a two-component blend by considering a combination method of identifying an optimal excitation frequency and eliminating the background dissipation.

    2. Factors affecting phases of TM-AFM

    2.1. Phases and phase contrasts of TM-AFM

    The key component of TM-AFM is a micro-cantilever with a sharp tip at the end which is very sensitive to small interaction forces. During its application,the micro-cantilever is oscillated at a frequency close to a resonance frequency while scanning a sample.[1]The tip contacts the sample intermittently, and information about the sample surface is obtained by monitoring the changes in the tip’s dynamic responses. In Fig. 1(a),ldenotes the length of the micro-cantilever,bthe width, anddthe thickness. The differential equation of the micro-cantilever under the interactions between the tip and sample is expressed as follows:

    wherewrepresents the absolute displacement of the microcantilever,uthe relative displacement of the micro-cantilever relative to the substrate,EandIare the Young’s modulus and section moment of inertia of the microcantilever,respectively;ρindicates the density of the micro-cantilever,cthe equivalent damping coefficient per unit length of the micro-cantilever,m′the mass of the tip,F(xiàn)tsthe interaction forces between the tip and the sample,andδ(xl)the Dirac function.

    The term with equivalent damping coefficient (i.e., the last term on the left-hand side) in Eq. (1) is the dissipation term, which reflects the energy dissipation in the microcantilever system. As shown in Fig.1(a),the micro-cantilever has various energy dissipation paths when working in an atmospheric environment. From the perspective of phase imaging, there are two types of energy dissipation in the system.The first type is the contact dissipation between the tip and the sample, which can directly reflect the sample surface information,and this type of energy dissipation is referred to as the key dissipation. The other type of dissipation which does not reflect the surface information of the sample is termed the background dissipation.

    Fig.1. (a)Schematic diagram of energy dissipation process in TM-AFM and(b)simplified model of micro-cantilever system in TM-AFM.

    In the tapping mode, the movement of the microcantilever in the air will drag some of the medium with it.Part of the energy of the micro-cantilever is converted into the kinetic energy of the medium, and this part of energy loss is referred to as the air viscous dissipation.[33,37]Generally speaking, the quality factors of the micro-cantilever under high vacuum condition can reach hundreds of thousands.Meanwhile, in atmospheric environments, the quality factors of the micro-cantilever tend to be only several hundred.[38,39]Consequently, air viscous dissipation is one of the most important energy dissipations in the TM-AFM. For example, if the gaps between the micro-cantilever and the sample surfaces become smaller,the micro-cantilever will need to squeeze the medium out of the gaps.[40,41]The energy loss in these types of processes is named the squeeze film dissipation. In addition, when the gap between the tip and sample surface is smaller than a critical distance (approximately the length of the Kelvin radius), liquid bridge will be formed between the tip and the sample surface. When the tip leaves the sample surface, the liquid bridge breaks, during which time the liquid bridge fracture energy will be a dominant part of the energy dissipation.[2,34]In addition, a micro-cantilever system also has support dissipation due to the strain concentrating at the attachment points of the resonators, thermoelastic damping dissipation caused by phonon–phonon interactions, surface dissipation,etc.[42–45]However,owing to the fact that the above-mentioned dissipations exist in the system itself and do not reflect the sample information, they are put into the category of background dissipation in this study.The effects of the background dissipation on the phase and phase contrast results will be analyzed later.

    In the TM-AFM, the piezoelectric actuator at the fixed end of the micro-cantilever drives the micro-cantilever through displacement excitation, which is expressed aszf=Dsin?t.The absolute displacement of the micro-cantilever isw(x,t)=u(x,t)+Dsin?t, and the differential equation of the microcantilever under displacement excitation is as follows:

    The boundary conditions of the micro-cantilever are given by

    Therefore, according to the method of separating variables,there exists the following equation:

    Substituting Eq. (4) into Eq. (2), for a linearized system ignoring damping and excitation, the general solution of the modal function can be obtained from the boundary conditions as follows:

    where

    whereωiis thei-th natural frequency of the micro-cantilever.

    Subsequently, the following equation is obtained in this study by substituting Eq.(4)into Eq.(2):

    where

    and the relative damping coefficients of each orderξi=ci/2ρbdωi. The values ofφi(l),Mi, andmiare obtained from the modal functions and the parameters of the microcantilever. In this work, the time solution of the relative displacement can be solved by using Eq. (8), and the absolute displacement of the micro-cantilever is given as follows:

    Therefore,by considering the first-order mode,whenx=l,the absolute displacement of the tip is obtained below

    In the TM-AFM,the micro-cantilever vibrates at about its first resonant frequency.The sample surface is then imaged by monitoring the displacements at the tip of the cantilever. The steady-state responses of the tip are shown in Eq.(10). Therefore, the position of the tip and the first-order mode can be used to simplify the micro-cantilever. The interaction forces between the tip and the sample include nonlinear elastic forces and dissipative forces. As shown in Fig.1(b),for the tip equation of motion, a micro-cantilever which has been excited by displacement at its first natural frequency or close to it can be obtained from the following equation:[2,46]

    wheremerepresents the effective mass, which is approximately one quarter of the total mass of the micro-cantilever;keis the force constant of the micro-cantilever;c0denotes the equivalent damping coefficient describing the background dissipation;zrepresents the instantaneous deflection of the tip;zf=Dsin?tis the displacement excitation;Dindicates the excitation amplitude;and?is the excitation frequency;hrepresents the instantaneous distance between the tip and the sample;Ftsrepresents the tip-sample interaction forces. The force between the tip and the sample is further simplified into a nonlinear elastic force and a viscous damping force,and the equation is expressed as follows:

    wherectsrepresents the equivalent damping coefficient representing the contact dissipation between the tip and the sample;ktsis the elastic coefficient representing the conservative force between the tip and the sample. Previous studies have shown that the force between the tip and the sample has little effect on the natural frequency of the system.[41,47]Therefore,ktsis ignored in the present investigation. The equation which describes the TM-AFM is given as follows:w

    hereQ=meω/(c0+cts)is the total quality factor of the system,andω=(ke/me)-1/2represents the angular frequency of the first resonance. In this study,the quality factor represents the ratio of the total energy to the energy lost per oscillation of the micro-cantilever and is expressed as follows:[33,48,49]

    whereWindicates the total energy stored in the microcantilever system, and ΔWrepresents the energy dissipated in each cycle of the micro-cantilever. Therefore, the energy dissipation of the system can be expressed by quality factor.As shown in Fig. 1(a), the total quality factor of the microcantilever system is expressed as follows:

    whereQkeyis the quality factor representing the key dissipation between the tip and the sample;Qbgis the quality factor representing the background dissipation of the microcantilever system;Qviscousis the quality factor representing the viscous dissipation;Qsqueezeis the quality factor representing the squeeze film dissipation;Qbridgeis the quality factor representing the liquid bridge dissipation. Subsequently,the phase lag angle of the steady-state response obtained from Eq.(13)is

    wheres=?/ω0denotes the frequency ratio. Therefore, according to Eq.(15),the phase in Eq.(16)reflects not only the dissipation between the tip and the sample,but also the background dissipation of the micro-cantilever system. However,it is determined that the phase does not reflect the real phase information of the sample,but the distorted phase masked by the background dissipation.In this study,in order to obtain the real phase information of the sample,it is necessary to eliminate the background dissipation. The specific method adopted for the elimination of the background dissipation is presented in Subsection 2.3.

    2.2. Effects of excitation frequencies on phase contrast results

    In the TM-AFM, the phase images are obtained by detecting the phase lag angles of the tip at each position on the sample surface. The difference in phase lag angle between two different components on the sample surface is referred to as the phase contrast, which reflects the difference in sample surface property.The greater the phase contrast,the easier it is to distinguish the properties of a sample’s surface area. In the present investigation,according to Eq.(16),the phase contrast between two areas with different properties can be expressed as

    whereQ1andQ2are the quality factors of the system when the tip is located at two areas with different properties,respectively.In this study,Q1is less thanQ2by default.Therefore,it can be seen from the above formula that whenQ1andQ2are determined,the phase contrast will be related to the excitation frequency. As a result, the phase contrast can be controlled and operated by using the excitation frequency.

    In the TM-AFM operation process, the excitation frequency of the piezoelectric actuator is one of the most basic and easiest parameters to change. For example, it can be expected that a maximum phase contrast can be identified by selecting the optimal excitation frequency. As can be seen from Eq. (17) when the quality factors of the microcantilever system working in two areas with different properties are determined, the phase contrast will only be related to the excitation frequency. The excitation frequency corresponding to the maximum phase contrast can be obtained by calculating the extreme value of the above-mentioned expression of the phase contrast. Since an excitation frequency slightly smaller than the natural frequency is generally used in the operations, the corresponding frequency when the phase contrast is taken as a maximum value is

    It can be seen from Eq.(18)that an excitation frequency to maximize the phase contrast exists in theory,which also explains why the excitation frequency is generally adjusted in order to obtain the best phase contrast in experimental processes.

    In the present experiments,a Bruker dimension XR AFM is used to scan the surfaces of a blend of polystyrene and lowdensity polyethylene(PS-LDPE).The experiments are carried out by using a Tap AN-NSC10 probe at room temperature and a relative humidity of approximately 30%. According to the results of frequency sweep experiments,the natural frequency of the micro-cantilever system is approximately 188 kHz. The frequency ratio increases from 0.95 to 1 with equal intervals,and 11 groups of experiments are performed. The phase images corresponding to the different excitation frequencies are detailed in Fig.2.

    Fig.2. Phase images of PS-LDPE with different excitation frequencies,where white scale bar is 1 μm;total scan size is 4.5 μm×4.5 μm;sizes of the yellow box and the blue box are both 1 μm×1 μm;sample corresponding to the yellow box area is LDPE,and sample corresponding to blue box area is PS.Inset in lower right corner of each phase image shows distribution density of the data in the corresponding specific phase image.

    Fig. 3. Normalized theoretical and experimental phase contrast results of PS-LDPE.

    As shown in Fig.2,the characteristics of the sample can obviously be seen in the phase images, where the LDPE is represented by a black circular area, and the PS is located in the yellowish-brown background area. Owing to the fact that the properties of the PS and LDPE are very different, excellent contrast results are obtained in phase images at all frequency ratios.It is observed that when the frequency ratios are 0.955,0.995,and 1,there occur obvious scratched areas in the phase images due to the instability or interference effects of the micro-cantilever system. In Fig.2,the insets in the lower right corners of phase images show the corresponding distribution densities of the data in those specific phase images. In each inset,the horizontal axis represents the phase,and vertical axis the density of the phase data. The phases of the start position and end position of the horizontal axis correspond to the minimum phase and maximum phase on the scale on the right side of the figure, respectively. Generally speaking, the sharper the two peaks in the binormal distribution,and the farther apart from each other the two peaks,the greater the phase contrast is. However, since the phase contrast of each image cannot be perceived intuitively,it is necessary to quantify the phase contrast in each phase image.

    It is determined from Eq. (18) that there is an optimal excitation frequency for maximizing the phase contrast. However,in order to prove this experimentally,the phase images at different frequency ratios are quantified. In this study the data in the blue box and yellow box,which are shown in Fig.2,are used respectively,and the method proposed in Refs.[50–52]is adopted to quantify the phase contrast. The statistical formula is given as follows:

    where Ashman’sDrefers to the phase image contrast;μ1andμ2represent the normal distribution means of the phases of the two components; andσ1andσ2are the standard deviations. The quantization results of the phase images with different frequency ratios are detailed in Fig.3.

    It can be seen in Fig.3 that the theoretical curve fitted by Eq.(17)coincides exactly with the experimental results,which indicates that the theoretical description of the phase contrast in this study is in line with reality and reflects the essence of the phases in the TM-AFM. The experimental results are fitted using Eq. (17). The values ofQ1andQ2are 18 and 54,respectively. These represent the quality factors of the system when the tip contacts the PS and LDPE,respectively.Then,by substitutingQ1andQ2into Eq. (18), the theoretical optimal frequency ratio is determined to be 0.984. The experimental results show that there is an optimal frequency ratio(approximately 0.985)available to maximize the phase contrasts of the two areas with different properties, which is almost exactly consistent with the theoretical results.

    Therefore, this study is able to select a certain number of frequency ratios near the natural frequency for the phase scanning processes,and then fit the experimental results using Eq. (17) in order to obtain the correspondingQ1(Qkey1) andQ2(Qkey2) of the two components of the blend, respectively.Subsequently, the optimal excitation frequency of the microcantilever system is theoretically determined. Then, the optimal excitation frequency obtained by the theoretical fitting is selected for the phase imaging process,and the phase images with the best contrast results are selected.

    2.3. Acquisition of real phase information of samples

    The phase image reflects all the energy dissipations of the TM-AFM system, including the background dissipation and the key dissipation. As shown in Fig. 1(a) and Eq. (15),in the atmospheric environment of this study, there are many background dissipations of micro-cantilever system, such as air viscous damping,squeeze film damping,liquid bridge dissipation. However, only the contact dissipation between the tip and the sample can accurately reflect the real information about the sample surface.It is found that the background dissipation can potentially reduce the phase contrast or even mask or distort the effective information in the phase images.Therefore,it is necessary to determine a method which only reflects the key dissipation and eliminates the effect of the background dissipation on phase contrast result.

    Fig.4. Schematic diagram of self-excitation method for TM-AFM.

    Figure 4 shows the self-excitation method (also known as the Q-control)[53]in the TM-AFM.The micro-cantilever is excited jointly by the transformed tip response signals and the original excitation. The goal of this study is to offset the effect of the background dissipation on the phase contrast by using an effective compensation method. The general equation of the tip motion with self-excitation gain in the TM-AFM is expressed as[53,54]

    whereGis the gain factor of the self-excitation,andφdenotes the phase shift between the self-excitation signal and the instantaneous deflection. Therefore,whenφis taken asπ/2,the phase lag angle can be expressed as

    It can be seen from Eq. (21) that the background dissipation in the system can be reduced or eliminated by adjusting the gain factorG. When the gain factorG=s/Qbg,the background dissipation can be effectively eliminated. Figure 5 shows the phases at different gain factors,obtained using Eq.(21),whereQkey1andQkey2represent the key dissipation of the systems on two samples of different properties,respectively. In the figure, whenG=0, the relationship between the phase and the key dissipation without self-excitation is indicated. In addition,G=0.00197 shows that the background dissipation is just offset by the gain factor. WhenG=0.00788 andG=0.01970,it is determined that the gain factor has offset the background dissipation and reduced the key dissipation to a certain extent. In addition, Δφ1, Δφ2, Δφ3, and Δφ4represent the phase contrasts between two samples with different properties under different gain factors.

    Fig.5.Curves of phase versus background dissipation for different gain factors,where frequency ratio is s=0.985,and quality factor Qbg representing background dissipation is 500, with different colors denoting phases under different gain factors.

    The curves of phaseversusbackground dissipation for different gain factors are compared and the phase change rates increase with gain factorGincreasing, which are shown in Fig. 5. Therefore, whenQkey1andQkey2reflecting the key dissipations between the tip and the sample are determined,it can be observed that the phase contrasts of two components with different properties increase with gain factor increasing.That is to say, Δφ4>Δφ3>Δφ2>Δφ1. WhenG=0, the phases undergo the effect of background dissipation. WhenG=0.00197,it is indicated that the effects of the background dissipation on the phases are eliminated,and the phases reflect only the real core dissipation between the tip and the samples.Furthermore, whenG=0.00788 andG=0.01970, the gain factor offsets part of the key dissipation, and the phases and phase contrasts are both distorted. Incidentally, in the liquid environment,the quality factor of micro-cantilever is very low,usually only tens or even single digit,and the conventional tapping mode is very difficult to image in this environment.[30]Therefore,if the self-excitation method is used to eliminate the influence of background dissipation on phase, it is of significance in implementing the phase imaging of some living biological cells and biological macromolecules.Generally speaking, the eliminating of the background dissipation can eliminate its masking effects on the key dissipation and improve the phase contrast results to a certain extent. However, it should be noted that if the gain factor has an unreasonable value and thus offsets only part of the key dissipation, then the phase information will also be distorted. Therefore, based on this study’s results, it is necessary to select an appropriate value for the gain factor.

    The quality factor of a micro-cantilever system oscillating in air can be obtained generally using frequency sweep in experiment. Figure 6 shows the curves of amplitudeversussweep frequency for different distances between the tip and the sample. It can be seen in the figure that the maximum amplitude of the curve decreases with the relative distance between the needle tip and the sample decreasing. In addition,the quality factor of the system also decreases with the relative distance decreasing. It is indicated that the decreases in the quality factor are caused by such factors as the squeeze film dissipation of the micro-cantilever and the liquid bridge dissipation of the tip when the micro-cantilever is close to the sample surface. It is observed that when the distance between the tip and the sample is further reduced, the frequency response curve is truncated.At that time,the quality factor of the system cannot be obtained by using frequency sweep in experiment. Since the truncation of the sweep frequency curve represents the contact between the tip and the sample,the sweep frequency curve before truncation reflects the background dissipation.

    Fig.6.Curves of amplitude versus frequency for different distances between tip and sample.

    A new method is proposed in this study in order to obtain the quality factor representing the background dissipation of the system. Multiple frequency sweeping experiments are performed near the sample surface for obtaining the approximate range of the quality factor representing the background dissipation. Figure 6 shows the experimental results about frequency sweep at different initial distances between the tip and the sample surface. As can be seen in the figure, when the tip is in“virtual”contact with the sample surface at the maximum displacement(for example, the tip is infinitely close to the sample surface but does not contact the sample surface or interaction between the tip and sample surface is very inconsiderable), the quality factor obtained from the frequency sweep experiment can be utilized as the reference value of theQbg.

    In summary, the background dissipation is found to reduce the phase contrasts of the samples or even distort the phase information. Therefore, through the above-mentioned self-excitation methods,this study can obtain the phase information which includes only the key dissipation between the tip and the sample (in other words, the real phase information),thereby improving the phase contrast results of TM-AFM to a certain extent.

    3. Conclusions

    In this study, based on the dynamic model of TM-AFM,two factors affecting phase contrast are determined, and two optimization methods of the phases are proposed. An optimal excitation frequency proves both theoretically and experimentally to be available for maximizing the phase contrasts between the two components of the blend.Therefore,the optimal excitation frequency can be selected to optimize the phase contrast between the two components. Furthermore,the masking effect of the background dissipation on the phase information of sample surface is analyzed theoretically. A self-excitation method is utilized in this study to reduce or eliminate the effect of background dissipation on the phase, thereby reaching the goal of accurately reflecting the real phase information of each component. The results show that the elimination of the background dissipation can reflect the real phase information of the sample and improve the phase contrast between each component of the sample. The proposed method is applicable to optimizing not only the phases of the two-component blend,but also the phases of each component within multi-component sample.In addition,in this study a frequency sweep method of obtaining the background dissipation is proposed,thereby providing a reference for eliminating the background dissipation.The aforementioned method is considered to have both theoretical and experimental significance in optimizing the phase of TM-AFM.

    Acknowledgement

    Project supported by the National Natural Science Foundation of China(Grant No.11572031).

    猜你喜歡
    魏征
    唐太宗與魏征
    中國史上最能杠的“找碴兒王”
    Design of a novel correlative reflection electron microscope for in-situ real-time chemical analysis?
    魏征身后的落寞
    中外文摘(2021年6期)2021-11-11 17:43:34
    唐太宗與魏征
    讀書文摘(2020年10期)2020-11-23 06:30:39
    魏征身后的落寞
    讀書文摘(2020年10期)2020-11-23 06:30:39
    唐太宗因何會砸了魏徵的墓碑
    北廣人物(2018年17期)2018-05-09 06:42:48
    唐太宗的“砸”與“立”
    文存閱刊(2016年1期)2016-03-17 04:06:52
    唐太宗為何不殺魏征
    魏征進(jìn)諫
    前線(2015年7期)2015-05-14 11:00:08
    国产老妇伦熟女老妇高清| 国产爱豆传媒在线观看| 一级爰片在线观看| 久久久亚洲精品成人影院| 国内精品一区二区在线观看| 欧美潮喷喷水| av黄色大香蕉| 男人狂女人下面高潮的视频| 少妇猛男粗大的猛烈进出视频 | 白带黄色成豆腐渣| 日本一本二区三区精品| 国产精品蜜桃在线观看| 久久久色成人| 亚洲国产最新在线播放| 69av精品久久久久久| 亚洲最大成人av| 嫩草影院精品99| 高清av免费在线| 高清日韩中文字幕在线| 日韩伦理黄色片| 亚洲自拍偷在线| 精品国产一区二区三区久久久樱花 | 欧美潮喷喷水| 别揉我奶头 嗯啊视频| 美女高潮的动态| 亚洲欧美清纯卡通| 2021少妇久久久久久久久久久| 国产av国产精品国产| 欧美zozozo另类| 嫩草影院精品99| 中国国产av一级| 国产精品一区二区三区四区久久| 日韩制服骚丝袜av| 亚洲综合色惰| 欧美日韩在线观看h| 日韩欧美 国产精品| 日本黄大片高清| 99视频精品全部免费 在线| 国产激情偷乱视频一区二区| 免费看不卡的av| 精品久久久噜噜| 日日摸夜夜添夜夜添av毛片| 亚洲最大成人中文| 国产一级毛片在线| 一级黄片播放器| 麻豆成人午夜福利视频| videossex国产| 欧美激情在线99| 久久久精品欧美日韩精品| 日本猛色少妇xxxxx猛交久久| 美女国产视频在线观看| 国产伦精品一区二区三区视频9| 日韩伦理黄色片| av播播在线观看一区| 亚洲精品视频女| 亚洲精品视频女| 亚洲精品视频女| 久久国产乱子免费精品| 日韩不卡一区二区三区视频在线| 国产 亚洲一区二区三区 | 久久国产乱子免费精品| freevideosex欧美| 国产极品天堂在线| 亚洲精品,欧美精品| 三级男女做爰猛烈吃奶摸视频| 熟女电影av网| 久久精品国产亚洲av涩爱| 国产色婷婷99| 免费黄频网站在线观看国产| 在线观看人妻少妇| 欧美成人精品欧美一级黄| 啦啦啦啦在线视频资源| 亚洲一区高清亚洲精品| 能在线免费观看的黄片| 亚洲精品第二区| 最近的中文字幕免费完整| 秋霞在线观看毛片| 嘟嘟电影网在线观看| 中文字幕免费在线视频6| 国产色爽女视频免费观看| 欧美高清成人免费视频www| 国产淫语在线视频| 内地一区二区视频在线| 老师上课跳d突然被开到最大视频| 一级爰片在线观看| 免费看不卡的av| 国产成人freesex在线| 久久久久久久久久久免费av| 成人亚洲精品一区在线观看 | 国产亚洲5aaaaa淫片| 日韩欧美 国产精品| 久久久亚洲精品成人影院| av免费观看日本| 中文字幕av在线有码专区| 男人和女人高潮做爰伦理| 久久97久久精品| 亚洲va在线va天堂va国产| 国产精品久久久久久精品电影小说 | 国内精品宾馆在线| 国产精品一区www在线观看| 少妇人妻精品综合一区二区| 婷婷色av中文字幕| 欧美极品一区二区三区四区| 久久久久性生活片| 国产v大片淫在线免费观看| 日本三级黄在线观看| 亚洲三级黄色毛片| 亚洲在久久综合| 日韩,欧美,国产一区二区三区| 看十八女毛片水多多多| 免费观看av网站的网址| 成年女人看的毛片在线观看| 精品熟女少妇av免费看| 身体一侧抽搐| 成人毛片a级毛片在线播放| 久久热精品热| 久久99蜜桃精品久久| 欧美潮喷喷水| 又爽又黄a免费视频| 99热全是精品| 特大巨黑吊av在线直播| 国产大屁股一区二区在线视频| 亚洲天堂国产精品一区在线| 成人美女网站在线观看视频| 国产伦精品一区二区三区四那| 国产激情偷乱视频一区二区| 国产伦在线观看视频一区| 久久久久久久久久人人人人人人| 两个人的视频大全免费| 日本色播在线视频| 黑人高潮一二区| 亚洲精品国产成人久久av| 日本av手机在线免费观看| 春色校园在线视频观看| 午夜激情福利司机影院| 亚洲av不卡在线观看| 麻豆av噜噜一区二区三区| 国产一级毛片在线| 又黄又爽又刺激的免费视频.| 日韩欧美一区视频在线观看 | 国产大屁股一区二区在线视频| 日韩国内少妇激情av| 一个人看视频在线观看www免费| 亚洲av免费高清在线观看| 在线免费十八禁| 非洲黑人性xxxx精品又粗又长| 精品久久久噜噜| 成人毛片60女人毛片免费| 亚洲国产欧美人成| 六月丁香七月| 午夜激情欧美在线| 亚洲av不卡在线观看| 欧美日韩国产mv在线观看视频 | 精品人妻一区二区三区麻豆| 高清欧美精品videossex| 大香蕉97超碰在线| 超碰97精品在线观看| 亚洲在线自拍视频| 搞女人的毛片| 纵有疾风起免费观看全集完整版 | 国国产精品蜜臀av免费| 一夜夜www| 国产精品美女特级片免费视频播放器| 中文字幕制服av| a级一级毛片免费在线观看| 99久久九九国产精品国产免费| 成人毛片60女人毛片免费| 又大又黄又爽视频免费| 国产精品三级大全| 男女视频在线观看网站免费| 国产精品一区www在线观看| 亚洲av电影不卡..在线观看| 亚洲精品一二三| 三级经典国产精品| 2021天堂中文幕一二区在线观| 久久国产乱子免费精品| 久久人人爽人人爽人人片va| 亚洲国产最新在线播放| 蜜桃亚洲精品一区二区三区| 国产激情偷乱视频一区二区| 波野结衣二区三区在线| freevideosex欧美| 欧美精品国产亚洲| 日本免费a在线| 亚洲精品第二区| 男女那种视频在线观看| 欧美变态另类bdsm刘玥| 日本黄大片高清| 精品久久久久久电影网| 舔av片在线| 日韩亚洲欧美综合| 麻豆精品久久久久久蜜桃| 日韩视频在线欧美| 白带黄色成豆腐渣| ponron亚洲| 色吧在线观看| 日韩av在线大香蕉| 免费看av在线观看网站| 久久这里有精品视频免费| 亚洲电影在线观看av| 亚洲乱码一区二区免费版| 欧美精品国产亚洲| 日本av手机在线免费观看| 色视频www国产| a级一级毛片免费在线观看| 国产av码专区亚洲av| 亚洲欧洲国产日韩| 亚洲av电影在线观看一区二区三区 | 日本-黄色视频高清免费观看| 美女大奶头视频| 深爱激情五月婷婷| 午夜爱爱视频在线播放| 色尼玛亚洲综合影院| 伊人久久精品亚洲午夜| 午夜福利高清视频| 一级毛片 在线播放| 国产精品久久久久久av不卡| 免费看日本二区| 亚洲性久久影院| 久久久久久九九精品二区国产| 日韩大片免费观看网站| 在线观看免费高清a一片| 狂野欧美白嫩少妇大欣赏| 国产日韩欧美在线精品| 一个人看视频在线观看www免费| 久久久久久九九精品二区国产| 国产精品av视频在线免费观看| 嫩草影院新地址| 欧美丝袜亚洲另类| 国产单亲对白刺激| 婷婷六月久久综合丁香| 日日干狠狠操夜夜爽| 男的添女的下面高潮视频| 你懂的网址亚洲精品在线观看| 99久国产av精品| 国产毛片a区久久久久| 春色校园在线视频观看| 在线免费十八禁| 真实男女啪啪啪动态图| 精品国产露脸久久av麻豆 | 久久人人爽人人爽人人片va| 亚洲四区av| 老师上课跳d突然被开到最大视频| 夫妻性生交免费视频一级片| 黄片wwwwww| 欧美激情在线99| 亚洲国产成人一精品久久久| 十八禁网站网址无遮挡 | 成人国产麻豆网| 两个人视频免费观看高清| 中文字幕免费在线视频6| av在线亚洲专区| 2021天堂中文幕一二区在线观| 久久久久久国产a免费观看| 精品一区二区三区视频在线| 精品一区二区三卡| 日韩av不卡免费在线播放| a级毛片免费高清观看在线播放| 日日摸夜夜添夜夜爱| 国产精品美女特级片免费视频播放器| 午夜精品一区二区三区免费看| 欧美日韩一区二区视频在线观看视频在线 | 国产男女超爽视频在线观看| 天天躁夜夜躁狠狠久久av| 黑人高潮一二区| 亚洲最大成人手机在线| 国产综合懂色| 欧美不卡视频在线免费观看| 26uuu在线亚洲综合色| 久久精品久久精品一区二区三区| 麻豆精品久久久久久蜜桃| 丝瓜视频免费看黄片| .国产精品久久| 一级毛片 在线播放| 一区二区三区乱码不卡18| 国产成人a区在线观看| 亚洲欧美精品自产自拍| 色吧在线观看| 日韩亚洲欧美综合| 联通29元200g的流量卡| av.在线天堂| 亚洲精品色激情综合| 亚洲欧美清纯卡通| 别揉我奶头 嗯啊视频| av福利片在线观看| 久久99热这里只有精品18| 少妇猛男粗大的猛烈进出视频 | 国产精品一二三区在线看| 91狼人影院| 黄片无遮挡物在线观看| 欧美日韩视频高清一区二区三区二| 人妻系列 视频| 亚洲国产日韩欧美精品在线观看| 男的添女的下面高潮视频| 777米奇影视久久| 搞女人的毛片| 男的添女的下面高潮视频| 亚洲国产欧美人成| 久久99热6这里只有精品| 国产成人freesex在线| 男人舔女人下体高潮全视频| 亚洲精品影视一区二区三区av| 免费黄网站久久成人精品| 少妇高潮的动态图| 菩萨蛮人人尽说江南好唐韦庄| 最后的刺客免费高清国语| 精品久久久噜噜| 男人舔奶头视频| 国产精品女同一区二区软件| 欧美3d第一页| 国产精品一区二区三区四区免费观看| 婷婷六月久久综合丁香| 大香蕉97超碰在线| 欧美97在线视频| 成人毛片60女人毛片免费| 亚洲av成人av| 日韩 亚洲 欧美在线| 亚洲精华国产精华液的使用体验| 久久99精品国语久久久| 天美传媒精品一区二区| 日本三级黄在线观看| 精品国产一区二区三区久久久樱花 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av成人精品| 成人综合一区亚洲| 亚洲精品乱码久久久久久按摩| 国产亚洲最大av| 女人被狂操c到高潮| 岛国毛片在线播放| 免费播放大片免费观看视频在线观看| 欧美成人a在线观看| 高清毛片免费看| 国产成人精品久久久久久| 高清av免费在线| 麻豆精品久久久久久蜜桃| 岛国毛片在线播放| 国产精品av视频在线免费观看| 性色avwww在线观看| 欧美潮喷喷水| 久久99精品国语久久久| 成人午夜高清在线视频| 丝瓜视频免费看黄片| 免费观看av网站的网址| 男女边吃奶边做爰视频| a级毛色黄片| 伦理电影大哥的女人| 搡女人真爽免费视频火全软件| 性色avwww在线观看| 在线观看人妻少妇| 人妻夜夜爽99麻豆av| 成人漫画全彩无遮挡| 国产在线一区二区三区精| 午夜福利网站1000一区二区三区| 亚洲国产欧美人成| 亚洲欧美中文字幕日韩二区| 尾随美女入室| 婷婷色综合大香蕉| 最近手机中文字幕大全| 欧美高清成人免费视频www| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| 色视频www国产| 不卡视频在线观看欧美| 免费av毛片视频| 在线免费观看不下载黄p国产| 99re6热这里在线精品视频| 日韩精品有码人妻一区| 国产精品久久视频播放| 又黄又爽又刺激的免费视频.| 蜜臀久久99精品久久宅男| 久久6这里有精品| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 久久99热这里只频精品6学生| 午夜久久久久精精品| 成人毛片60女人毛片免费| 舔av片在线| 午夜精品国产一区二区电影 | 大又大粗又爽又黄少妇毛片口| 亚洲自拍偷在线| 好男人视频免费观看在线| 欧美日本视频| av播播在线观看一区| 亚洲高清免费不卡视频| 国产色爽女视频免费观看| freevideosex欧美| 毛片一级片免费看久久久久| 夜夜爽夜夜爽视频| 日产精品乱码卡一卡2卡三| 午夜福利视频精品| 成年女人在线观看亚洲视频 | 在线观看av片永久免费下载| 乱人视频在线观看| 亚洲国产精品国产精品| av网站免费在线观看视频 | 成年人午夜在线观看视频 | 日产精品乱码卡一卡2卡三| 日日干狠狠操夜夜爽| 精品久久久噜噜| 婷婷色麻豆天堂久久| 亚洲欧美日韩卡通动漫| 国产日韩欧美在线精品| 日日啪夜夜撸| 成人无遮挡网站| 国产午夜精品论理片| 亚洲美女视频黄频| 人妻少妇偷人精品九色| h日本视频在线播放| 国国产精品蜜臀av免费| 蜜桃亚洲精品一区二区三区| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 亚洲精品国产av蜜桃| 国内少妇人妻偷人精品xxx网站| 久久97久久精品| 国产精品麻豆人妻色哟哟久久 | 好男人在线观看高清免费视频| 高清视频免费观看一区二区 | 国产色爽女视频免费观看| 国产亚洲5aaaaa淫片| 免费观看的影片在线观看| 伊人久久精品亚洲午夜| 床上黄色一级片| 别揉我奶头 嗯啊视频| 婷婷色麻豆天堂久久| 午夜老司机福利剧场| 久久久久免费精品人妻一区二区| 最近最新中文字幕免费大全7| 亚洲一区高清亚洲精品| 日本-黄色视频高清免费观看| 一二三四中文在线观看免费高清| 国产激情偷乱视频一区二区| 日韩国内少妇激情av| 97在线视频观看| 日本一本二区三区精品| 精品国产三级普通话版| 久热久热在线精品观看| 国产成人精品久久久久久| 丝袜美腿在线中文| 亚洲欧美成人精品一区二区| 在线免费观看的www视频| 晚上一个人看的免费电影| 国产爱豆传媒在线观看| 欧美成人a在线观看| 国产亚洲精品av在线| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕制服av| 在线观看av片永久免费下载| 大片免费播放器 马上看| 在线播放无遮挡| 国产亚洲午夜精品一区二区久久 | 国产成人精品一,二区| 少妇丰满av| 久久99热6这里只有精品| 纵有疾风起免费观看全集完整版 | 熟妇人妻不卡中文字幕| 汤姆久久久久久久影院中文字幕 | 久久久久久久久久久丰满| 久久久午夜欧美精品| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 亚洲精品国产成人久久av| 午夜精品一区二区三区免费看| 热99在线观看视频| 免费看av在线观看网站| 欧美日韩一区二区视频在线观看视频在线 | 一级毛片黄色毛片免费观看视频| 国产黄频视频在线观看| 一本一本综合久久| 欧美不卡视频在线免费观看| 噜噜噜噜噜久久久久久91| av在线老鸭窝| 成人国产麻豆网| 欧美3d第一页| 日韩av在线免费看完整版不卡| 日韩亚洲欧美综合| av天堂中文字幕网| 青青草视频在线视频观看| av专区在线播放| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品一二三| 日韩三级伦理在线观看| 乱码一卡2卡4卡精品| 男人舔女人下体高潮全视频| 人妻系列 视频| 纵有疾风起免费观看全集完整版 | av在线老鸭窝| 日韩制服骚丝袜av| 久热久热在线精品观看| 国产黄频视频在线观看| 水蜜桃什么品种好| 亚洲第一区二区三区不卡| 在线观看免费高清a一片| 熟妇人妻久久中文字幕3abv| 国产视频首页在线观看| 久久精品熟女亚洲av麻豆精品 | 成人无遮挡网站| 国产黄片美女视频| 日韩成人av中文字幕在线观看| 尾随美女入室| 亚洲最大成人av| 日韩一区二区视频免费看| 丰满人妻一区二区三区视频av| 精品99又大又爽又粗少妇毛片| 国产真实伦视频高清在线观看| 韩国av在线不卡| 蜜臀久久99精品久久宅男| 麻豆精品久久久久久蜜桃| 国产乱人视频| 国产成人aa在线观看| 大话2 男鬼变身卡| 午夜精品在线福利| 亚洲电影在线观看av| 性色avwww在线观看| 在线a可以看的网站| 亚洲国产精品成人综合色| 99re6热这里在线精品视频| 美女xxoo啪啪120秒动态图| 国产色爽女视频免费观看| 七月丁香在线播放| 成年av动漫网址| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 最近的中文字幕免费完整| 日韩av不卡免费在线播放| 久久久久网色| 身体一侧抽搐| 免费观看的影片在线观看| 你懂的网址亚洲精品在线观看| 欧美成人a在线观看| 岛国毛片在线播放| 国产高清有码在线观看视频| 亚洲在久久综合| 久久综合国产亚洲精品| 欧美高清成人免费视频www| 91久久精品国产一区二区三区| 午夜免费观看性视频| 高清av免费在线| 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| 久久久久久久久中文| 国产一区二区亚洲精品在线观看| 日韩一区二区视频免费看| 国产亚洲午夜精品一区二区久久 | 午夜视频国产福利| 丰满乱子伦码专区| 国产又色又爽无遮挡免| 搡老乐熟女国产| 久久亚洲国产成人精品v| 少妇的逼好多水| 亚洲丝袜综合中文字幕| 精品久久国产蜜桃| 国产午夜精品一二区理论片| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 亚洲精品国产成人久久av| 久久久欧美国产精品| 亚洲最大成人av| 国产成人午夜福利电影在线观看| 91久久精品国产一区二区三区| 国产又色又爽无遮挡免| 亚洲av电影在线观看一区二区三区 | 老司机影院成人| 人人妻人人澡人人爽人人夜夜 | 熟妇人妻不卡中文字幕| 看非洲黑人一级黄片| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 色尼玛亚洲综合影院| 精品一区二区免费观看| 最近最新中文字幕大全电影3| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 日韩三级伦理在线观看| 亚洲在线观看片| 国产国拍精品亚洲av在线观看| 青春草视频在线免费观看| 国产黄a三级三级三级人| 亚洲人成网站在线观看播放| 中国国产av一级| 亚洲欧美日韩东京热| 免费观看无遮挡的男女| .国产精品久久| 一区二区三区乱码不卡18| 在线a可以看的网站| 三级毛片av免费| 22中文网久久字幕| 非洲黑人性xxxx精品又粗又长| 男人舔奶头视频| 成人二区视频| 婷婷六月久久综合丁香| 中文字幕人妻熟人妻熟丝袜美| 亚洲av免费在线观看| 成年女人在线观看亚洲视频 | 一级毛片电影观看| 色5月婷婷丁香| 丝袜美腿在线中文| 成人亚洲精品一区在线观看 | 日韩一区二区视频免费看| 夫妻午夜视频| 亚洲性久久影院| 日韩av在线大香蕉| 九九在线视频观看精品| 中文乱码字字幕精品一区二区三区 | xxx大片免费视频| 99久久精品一区二区三区| 中国国产av一级| www.av在线官网国产| 插阴视频在线观看视频| 边亲边吃奶的免费视频| 青春草国产在线视频| 黄片wwwwww| 乱码一卡2卡4卡精品| 九色成人免费人妻av|