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

    In situ study of calcite-III dimorphism using dynamic diamond anvil cell

    2022-09-24 08:01:00XiaZhao趙霞ShengHuaMei梅升華ZhiZheng鄭直YueGao高悅JiangZhiChen陳姜智YueGaoLiu劉月高JianGuoSun孫建國YanLi李艷andJianHuiSun孫建輝
    Chinese Physics B 2022年9期
    關(guān)鍵詞:孫建國李艷升華

    Xia Zhao(趙霞) Sheng-Hua Mei(梅升華) Zhi Zheng(鄭直) Yue Gao(高悅) Jiang-Zhi Chen(陳姜智)Yue-Gao Liu(劉月高) Jian-Guo Sun(孫建國) Yan Li(李艷) and Jian-Hui Sun(孫建輝)

    1CAS Key Laboratory for Experimental Study under Deep-sea Extreme Conditions,Institute of Deep-sea Science and Engineering,Chinese Academy of Sciences(CAS),Sanya 572000,China

    2University of Chinese Academy of Sciences,Beijing 100049,China

    3Department of Mechanical Engineering,National University of Singapore(NUS),Singapore 117575,Singapore

    Keywords: calcite-III,dimorphism,dynamic diamond anvil cell,in situ Raman spectroscopy,pressurizing rate

    1. Introduction

    As the main carbon carriers, carbonate minerals control the migration of carbon in different geochemical reservoirs.[1-4]Calcite (CaCO3), the most common carbonate mineral, has over 12 proven structural polymorphs, such as calcite-I-calcite-VII(cc-I-cc-VII),aragonite,and pyroxene CaCO3.[2,5-12]Although it is accepted that calcite reacts with enstatite in the upper mantle to form dolomite,[13]the diamond inclusions and ancient subduction zone rocks suggest that calcite could exist at the top of the lower mantle.[14-16]Under the mantle pressure and temperature condition, calcite must undergo a substantial phase transition to adapt the changes in pressure and temperature with the depth.[2,5,6,11,17,18]Owing to its pressure sensitivity,structural phase transitions in calcite may also be used as geobarometers.[19]Thus,phase transition in calcite is of great importance in understanding the structure and dynamics of the Earth.

    The pressure-induced phase transition of cc-III dimorphism at room temperature has been extensively studied.[2,20-24]Despite a long history of investigation, the transition process among the dimorphisms of calcite,cc-III in triclinic and cc-IIIb in pseudomonoclinic triclinic structure,remains unclear. On the one hand, some researchers claimed that the cc-III and cc-IIIb can coexist in the same pressure range. For example, Liu and Mernagh identified the coexistence of cc-III and cc-IIIb at a pressure of 5.2 GPa, based on the observation of different cell parameters and Raman vibrational modes.[20]The coexistence was later confirmed by Merliniet al. but at a lower pressure of 3 GPa in polycrystalline samples.[21]Yuanet al. further determined the coexisting pressure range to be 1.97 GPa-3.38 GPa.[25]On the other hand,several researchers observed that each phase exists separately in different pressure ranges. For example,Pippingeret al. found that the pressure ranges for the existence of cc-IIIb and cc-III are 2.14 GPa-3.14 GPa and 3.42 GPa-5.42 GPa,respectively,at a temperature of 298 K,and the pressure ranges change to 2.2 GPa-3.71 GPa and 3.84 GPa-4.45 GPa,respectively,at a temperature of 283 K.[22]In both cases,the cc-IIIb and cc-III do not coexist. Furthermore, there are also studies suggesting that cc-IIIb may appear or disappear repeatedly with the increase of pressure,which complicates the issue even more. Koch-Muelleret al. found experimentally that the cc-IIIb first appears at 2.2 GPa, then disappears at~3 GPa due to the transition to cc-III,and finally reappears at~5 GPa.[24]The first-principle simulation also indicated that the transition between cc-IIIb and cc-III can oscillate as pressure increases from cc-IIIb to cc-III at 4.3 GPa and then from cc-III to cc-IIIb at 14.9 GPa.[23]These findings are summarized in Fig.1.In order to clarify the coexistence of cc-III and cc-IIIb, it is necessary to conduct further research.

    In addition,the Earth’s interior is persistently in the process of dynamic evolution, including seafloor spreading and plate subduction.[26-29]Some geological zones,such as those around fault zones during earthquakes or adjacent areas during volcanic eruptions, may be squeezed at a high pressurizing rate. Accordingly, the physical and chemical properties of the constituent minerals including calcite in the relevant Earth interior should be affected by these dynamic processes. For example, different carbonates (dolomite, aragonite,calcite,magnesite)were found in different high-pressure(HP) metamorphic rocks (metapelite, marble, eclogite) from the Western Tianshan subduction zone.[30]Rapid magma ascents contacted calcite inclusions within quartz phenocrysts at Santa Rita, New Mexico (USA).[31]However, up to now,the studies aiming at exploring the physical and chemical processes in the Earth’s interior,whether mineral physics or mineral chemistry,were mainly carried out under static conditions,and the dynamic research is almost blank. Undoubtedly, it is of great scientific significance to carry out dynamic research of the high-pressure phase transition of calcite under the conditions of temperature and pressure simulating mineral occurrence environment in geological processes.

    Fig.1. Summary of previous results of coexistence of cc-III and cc-IIIb at room temperature. Merlini et al. (2012) found that both cc-III and cc-IIIb can be stable in a pressure range of 2.5 GPa-15 GPa.[2] Later(2014), they observed the coexistence of these two phases at a pressure of ~3 GPa.[21]

    The spectral analysis provides useful information about the phase transitions in solids. Specifically, Raman scattering occurs when monochromatic light interacts with molecular vibrations inside a crystal. With the spectrum being able to record any change in bonding structure, Raman spectroscopy is used extensively in the study of phase transition in calcite.[5,7,20,32-40]The isolated carbonate ion[CO3]2-has four distinct vibration bands: symmetric stretching mode(ν1),out-of-plane bending mode(ν2),asymmetric stretching mode(ν3),and in-plane bending mode(ν4).Carbonate minerals also have two lattice vibration bands reflecting the lattice vibration(translation and oscillation) between Ca2+and [CO3]2-.[41]The Raman spectra of rhombohedral cc-I are characterized by five Raman-active peaks located at 156, 282, 712 (=ν4),1085 (=ν1), and 1434 cm-1(=ν3), respectively.[37]Owing to the high sensitivity in pressure, the change ofν4mode often receives special attention in the spectral analysis of phase changes in calcite. In calcite,the cc-I and cc-II have different crystal structures: the former is rhombohedral (R3c) and the latter is monoclinic(P21/m). The transition from cc-I to cc-II is achieved by shifting Ca2+cations and rotating the[CO3]2-group~11°.[25,42]Correspondingly, theν4mode in the Raman spectrum will split from one peak for cc-I into two peaks for cc-II. The structure of cc-III is less symmetrical than that of cc-II.The change from cc-II to cc-III results from the breaking of the triclinic(P)symmetry within the thermodynamically stable field of cc-III.[2]This change reflected in the Raman spectra is that theν4mode changes from two to more than five peaks, which has been confirmed by Pippingeret al.[22]and Yuanet al.[25]

    Diamond anvil cell (DAC) has been used extensively in high-pressure studies due to its capability of being incorporated with spectrum analysis techniques, and the understanding of the calcite phase has been greatly improved with the DAC technique.[6-10,43]However, the observation of detailed structure change is a challenge to the conventional DAC technique.In the conventional DAC,the increase of pressure in the experiment can only be achieved step-by-step mechanically.Because of a large pressure interval produced by discontinuous pressurizing,some transient processes,such as the coexistence of phases,could be easily overlooked. Dynamic loading DAC(dDAC)was developed from the conventional DAC.At a designated rate,it can provide continuous pressurizing. Several dDAC studies[44-48]on mineral phases(including Raman studies on carbon tetrachloride and ice, and x-ray diffraction studies on KCl and zirconium)have demonstrated the prospect of this new technology.

    In this study, using the dDAC combined with Raman spectral analysis,we realizein situmonitoring of calcite phase transition upon pressure increasing. Our results confirm the coexistence of cc-III and cc-IIIb,and the pressure range of the coexistence depends on the pressurizing rate. More importantly, our discovery reveals that kinetics plays a key role in controlling the transition, which may provide an insight into the structural transition in calcite.

    2. Materials and methods

    2.1. Sample composition

    All samples for DAC experiments in this study were from the same natural calcite crystal(see Fig.S1). The crystal was colorless and transparent, about 3 cm in diameter. Scanning electron microscopy (SEM) image and x-ray energy dispersive spectrum (EDS) of the crystal were taken on a Phenom ProX microscope (see Fig. S2). The EDS results confirmed the high purity of CaCO3crystal containing only three major elements: calcium,carbon,and oxygen.

    2.2. Experimental procedure

    High-pressure experiments of the single-crystal calcite were performed at the CAS Key Laboratory for Experimental Study under Deep-sea Extreme Conditions,Institute of Deepsea Science and Engineering, Chinese Academy of Sciences(Sanya, Hainan, China). Experiments were carried out on our dDAC system,which combines piezoelectric actuators and conventional DAC. Conventional DAC was located between pillars and pressure screws(see Fig.S3). The sample assembly was the same as that in a conventional DAC,equipped with two type-IIA low-fluorescence diamonds with 300-μm culets.The sample chamber was a hole of~100μm in diameter,produced by a Drilex-2000 laser punch,located in the center of a stainless-steel gasket of 100μm in thickness. Water was used as a pressure-transmitting medium. In the process of sample assembly,a ruby particle was enclosed with the sample in the sample chamber as a pressure indicator,and the pressure was measured byin situruby fluorescence.[49]

    The dDAC device provides controllable loading functions. By setting the target pressure and loading path (i.e.,the times for loading-stay-unloading)in advance through the signal generator, selected loading modes including the designated loading rate could be obtained. The correctness of the input signal was ensured by an oscilloscope(see Fig.S4).

    Combined with spectral analysis technology, our highpressure test system has the capability ofin situobservation. The Raman and fluorescence spectra of calcite and ruby were simultaneously collected by a charge coupled device(CCD) detector and an intensified CCD (ICCD) detector, respectively. The Raman spectra from 100 cm-1to 930 cm-1were recorded by an Andor Shamrock SR500i spectrometer equipped with a Leica DM2700 microscope, and a Plan Apo SL Infinity Corrected 20×objective with a numerical aperture value of 0.28. The excitation laser was a 532-nm frequencydoubled Nd:YAG laser. A 1800-gr·mm-1holographic grating coupled with a CCD detector(2000×256 pixels per inch,1 inch=2.54 cm) with a spectral resolution of 1.5 cm-1was used for the Raman spectrum measurement. Raman spectra of calcite at room temperature and atmospheric pressure were recorded to ensure a good instrument condition, and three standard Raman-active peaks were located at 156, 282, and 712 cm-1(=ν4), respectively. A 600-gr·mm-1holographic grating coupled with an ICCD detector(1024×1024 pixels per inch)was used to measure the ruby fluorescence spectra. The exposure time for fluorescence and Raman spectra were about 1 s and 4 s, respectively. The time-resolved Raman spectrogram in the dynamic loading process is shown in Fig.2.

    Fig.2. Typical time-resolved Raman spectra for calcite investigated in dynamic loading process.

    2.3. Ab initio simulation using density functional theory

    To corroborate the experimental Raman spectra,we used Materials Studio to calculate the Raman peaks of cc-III and cc-IIIb. The basic setup was similar to that in Bayarjargalet al.[39]with an 880-eV plane-wave cut-off and Monkhorst-Pack grids[50]of 1×1×1 grid parameters. We also used the same criterion of self-consistency as that in Bayarjargalet al.,[39]and the maximum displacement is 0.005 ?A.

    3. Results and discussion

    3.1. Experimental conditions

    Eight experiments from five samples (C1, C2, C3, C4,and C5) are performed with dDAC. Experiments are carried out at high pressures and different pressurizing rates. The detailed information for all runs is summarized in Table 1.

    Table 1. Pressure and loading rate settings for experiments.

    3.2. Raman spectrum analyses

    Previous studies revealed that in calcite, the phase transition leads the characteristic modeν4to split(e.g.,from one peak for cc-I to two peaks for cc-II, to more than five peaks for cc-III); the feature is therefore used to identify the phase state in calcite. While there is consensus on the feature ofν4mode of cc-I,cc-II,and cc-III,the feature ofν4for cc-IIIb remains uncertain. In this case, we try to find clues throughabinitiosimulation to reveal the relationship between the bonding structure and the feature of their characteristic modeν4for both cc-III and cc-IIIb.

    We perform first-principle simulations based on DFT for the characteristicν4mode of CaCO3at a spectra resolution of 4 cm-1separately at a pressure of 2.2 GPa and 7 GPa. As presented in Fig. 3, cc-III shows seven peaks at 2.2 GPa and eight peaks at 7 GPa,while cc-IIIb exhibits three peaks at both pressures. Therefore, we use the features ofν4mode: three peaks for cc-IIIb and more than seven for cc-III, as the criterion to identify cc-IIIb and cc-III in our analysis. The simulation results are in good agreement with previous observations.For example, Yuanet al. observed that theν4mode in calcite shows 3 sharp peaks in cc-IIIb,and splits into 5 peaks in cc-III.[25]Similarly, Pippingeret al. found that theν4mode in calcite changes from 4 to 7 peaks in the process of phase transition from cc-IIIb into cc-III.[22]

    Fig.3. Curves of characteristic ν4 mode of calcite,by first-principle simulation,for cc-III and cc-IIIb at P=2.2 GPa(a)and 7 GPa(b).

    3.3. Transition from cc-I to cc-II and from cc-II to cc-III

    Although this study focuses on the dimorphism of cc-III,the transition from cc-I to cc-II and from cc-II to cc-III,which are essential steps for the introduction of cc-III,need to be analyzed systematically.

    Figure 4 shows the Raman spectra recorded from run DD0012, indicating the transition of calcite from cc-I to cc-II. When the pressure rises to about 1.8 GPa, the vibration mode ofν4splits into two peaks(713.5 cm-1and 718.8 cm-1)from one(715.5 cm-1),which indicates the occurrence of the phase transition from cc-I to cc-II. The same phenomenon is observed in other runs, but the transition pressure is slightly different.

    The transition of calcite from cc-II to cc-III can be observed in the Raman spectra collected from run DD0006 as shown in Fig.5. In the phase of cc-II,the doubletν4peaks are located respectively at 712.3 cm-1and 721.3 cm-1at 2.2 GPa.When the pressure reaches 2.7 GPa,the doubletν4mode splits into eight peaks,located at 679.3,693.1,698.3,712.9,723.7,733.0,741.8,and 748.2 cm-1,respectively. This indicates the transition from cc-II to cc-III,and there is a transient process before the completion of phase transition, in which cc-II coexists with cc-III(e.g.,atP ≈2.4 GPa). Experimental results regarding the transitions, from cc-I to cc-II and from cc-II to cc-III,obtained from all runs are summarized in Table 2.

    Fig. 4. Raman spectra recorded from run DD0012 under dynamic loading pressurizing. At P ≈1.8 GPa,the vibration mode of ν4 splits from one peak to two peaks,indicating the transition of calcite from cc-I into cc-II.

    Fig.5. Raman spectra collected from run DD0006 under dynamic loading pressurizing.The vibration mode of ν4 splits from two peaks at P ≈2.2 GPa to eight peaks at 2.7 GPa, indicating the transition of calcite from cc-II to cc-III.There is a transient process before the completion of phase transition at ~2.4 GPa,in which cc-II coexists with cc-III.

    Table 2. Pressure conditions (P, in unit of GPa) for the stability of each phase and coexistence of cc-III dimorphism in calcite.

    Our results show that the transition pressures from cc-I to cc-II and from cc-II to cc-III are around 1.5±0.3 GPa and 2.5±0.3 GPa,respectively. Although the transition pressure from cc-I to cc-II is consistent with most of previous studies,[5,24,42,51,52]the transition pressure from cc-II to cc-III agrees only with some of the results[20,25]at a pressure of about 2.5 GPa. We attribute the discrepancy between this study and others in the transition pressure from cc-II to cc-III to the ambiguity of the cc-III and cc-IIIb dimorphisms. In some researches, cc-II is preferentially transformed into cc-IIIb rather than cc-III as pressure increases.[22-24]As a result,the transition pressure from cc-II to cc-III cannot be compared with their results,because the transition did not occur in their studies.

    3.4. Coexistence of cc-III and cc-IIIb

    As mentioned above,the dimorphism of cc-III is the focus of this study, and more attention is paid to the coexistence of cc-III and cc-IIIb. These two phases,based on the analyses of their vibrational characteristics, can be recognized from theirν4mode,with seven or more peaks for cc-III and three for cc-IIIb. Figure 6 shows the Raman spectra in the pressure range for the coexistence of cc-III and cc-IIIb. AtP ≈3.1 GPa,the characteristicν4mode exhibits eight peaks, located respectively at 679.3, 693.1, 699.5, 711.2, 724.5, 733.8, 742.2, and 747.9 cm-1, indicating that the calcite is in cc-III phase. As the pressure increases,theν4mode changes continuously. AtP ≈3.3 GPa, the peak located at 710.8 cm-1disappears, indicating the beginning of the phase transition from cc-III into cc-IIIb. At 6.7 GPa, theν4mode shows mainly three peaks at 699.1, 745.0, and 750.7 cm-1, marking the completion of the phase transition from cc-III to cc-IIIb. From that point on,a few weak external modes appear or disappear with pressure increasing, but theν4mode remains almost unchanged,indicating that calcite is stable in cc-IIIb phase.In the pressure range between 3.3 GPa and 6.4 GPa,however,the peaks representing cc-III diminish and the feature for cc-IIIb emerges,which can be considered as the coexistent pressure range of cc-III and cc-IIIb. Similar results are obtained in all other runs,which are summarized in Table 2.

    Fig. 6. Raman spectra recorded from run DD0006 under dynamic-loading pressurizing. The ν4 mode changes from eight peaks at 3.1 GPa into three peaks at 6.7 GPa, indicating the phase transition of calcite from cc-III to cc-IIIb. Between these two pressures,two phases coexist.

    Through continuously pressurizing andin situspectroscopic observations,we confirm the coexistence of cc-III and cc-IIIb. Although the coexistent pressure range varies slightly run by run,our observations are consistent with previous studies in which coexistence was also found.[2,20,21,25]

    3.5. Influence of pressurizing rates on phase transition

    In order to evaluate the influence of pressurizing rate on the phase transition of calcite, experiments are carried out at different pressurizing rates. The effect of pressurizing rate on the phase transition of calcite is displayed in Fig.7.

    As shown in Fig. 7, at different pressurizing rates, the transition pressures from cc-I to cc-II and from cc-II to cc-III are~1.5 GPa and~2.5 GPa,respectively,indicating that the transition pressures are not influenced by the pressurizing rate. In comparison,the transition pressures from cc-III to cc-IIIb are~9.9,~6.9,~6.1,and~4.5 GPa for the pressurizing rate of 0.006,0.012,0.021,and 0.056 GPa/s,respectively.The transition pressure decreases with the increase of pressurizing rate,showing a remarkable pressurizing rate dependence.

    Furthermore, it is found that the pressure range for the coexistence of cc-III and cc-IIIb is also dependent on the pressurizing rate. As shown in Fig. 7, as the pressurizing rate increases from 0.006 GPa/s to 0.056 GPa/s, the pressure range for the coexistence of cc-III and cc-IIIb decreases from~7 GPa to~1.5 GPa. The higher the pressurizing rate, the narrower the coexistence pressure range of cc-III to cc-IIIb is.

    Our results show that part of the phase transition process of calcite is related to the pressurizing rate,while the other part is independent of the pressurizing rate. This phenomenon can be explained from the perspective of energy.

    In thermodynamics, a system tends to be in an energy“well”(that is, the Gibbs energy is the lowest)to remain stable. In order to adapt to the changes of thermodynamic conditions(such as pressure),the system must find new energy well accordingly. Although thermodynamics provides the driving force to bring the system from one well to another, effective activation energy(Qeff)is required to overcome an energy barrier (Q) for completing the transition. In other studies,[53]it is found that theQeffdecreases exponentially with the pressurizing rate. The points mentioned here are schematically illustrated in Fig.S5 for later discussion.

    Comparing with the transition of calcite from cc-I to cc-II, and cc-III, the coexistence of cc-III and cc-IIIb can be regarded as a metastable state with a lower energy barrier for the next transition. If the transition between the two states with a highQ, even if theQis reduced toQeffdue to the influence of pressurizing rate, the reduced part of the energy,ΔQ=Q-Qeff?Q, is still too small to affect considerably the transition process. So,the effect of the pressurization rate is not obvious. However,for the metastable state of the coexistence of cc-III and cc-IIIb with a lowQ, any change of ΔQresulting from the influence of pressurizing rate(say,ΔQ~Qcan be significant. The higher the pressurizing rate,the lower the energy barrier is,and the easier it is to pull the system out of the metastable state. This is very consistent with our observations shown in Fig.7.

    Our results suggest that the coexistence of cc-III and cc-IIIb depends on the pressurizing rate and may provide clues to the explaining of the lack of observed coexistence in some static studies. In the conventional DAC, each stage of stepby-step pressurization is commonly completed in a very short time, that is, the pressurizing rate is very high. As shown in Fig. 7, the higher the pressurizing rate, the narrower the coexistence pressure range is. Since the pressurizing rate in stepwise-loading DAC can be several orders of magnitude higher than the fastest pressurizing rate in this study, it is not surprising that the coexistence, if any, of cc-III and cc-IIIb is overlooked in those researches carried out by the conventional DAC.

    Fig.7. Variations of transition pressures on calcite phase with the pressurizing rate for cc-I,cc-II,cc-III,and cc-IIIb. The portion of the coexistence of cc-III and cc-IIIb are highlighted and the pressurizing rate is log-scaled.

    4. Conclusions

    Using the d-DAC, the phase transitions in calcite, from cc-I to cc-II, cc-III, and cc-IIIb, are investigated based onin situexperiments. A number of new findings are revealed,especially about the coexistence of cc-III dimorphism,which are summarized as follows.

    (i)At room temperature, the phase transitions of CaCO3from cc-I to cc-II and from cc-II to cc-III occur at~1.5 GPa and~2.5 GPa,respectively.The pressures of these transitions are independent of the pressurizing rate.

    (ii)The transition pressure from cc-III to cc-IIIb is pressurizing rate dependent. As the pressurizing rate increases from 0.006 GPa/s to 0.056 GPa/s, the transition pressure decreases from 9.8 GPa to 4.5 GPa.

    (iii)The transition from cc-III to cc-IIIb is a sluggish process. Before the complete transition from cc-III to cc-IIIb,those two phases can coexist in a certain pressure range,which is dependent on the pressurizing rate. The higher the pressurizing rate,the narrower the pressure range for the coexistence of the two phases is.

    (iv)The dependence of the coexistence of cc-III to cc-IIIb on pressurizing rate may result from the influence of pressurizing rate on the activation process by reducing energy barrier.The higher the pressurizing rate, the lower the energy barrier is,and the easier it is to pull the system out of the metastable state.

    This study marks some of the first attempts ofin situRaman spectroscopy of mineral phase transition in dynamic loading DAC. It reveals that the phase transition of solid can be pressure-rate-dependent, providing new insights into the understanding of the structural changes of minerals in the dynamic environment.

    Acknowledgements

    We are indebted to Junwei Li for his technical help in the experiment. We want to thank High Performance Computing(HPC),NUS for providing technical support in part of our simulation work.

    Project supported by the Fund from the Chinese Academy of Sciences (Grant No. QYZDY-SSW-DQC029)and the National Natural Science Foundation of China(Grant No.41674097).

    猜你喜歡
    孫建國李艷升華
    歡喜? ?
    驚蟄,萬物復(fù)蘇萌人懂 ……
    貪吃蛇升華記
    小麥籽粒大小相關(guān)基因TaGS2克隆及功能分析
    “三法”化解表面積與體積問題
    6 Things You Must Know about Major Cold
    春天的電話
    胡娟:奉獻(xiàn)也是自我升華
    海峽姐妹(2018年1期)2018-04-12 06:44:09
    孫建國:不懈挑戰(zhàn),讓生命自由呼吸
    賣刀的人下落不明
    青春(2012年2期)2012-04-29 05:25:45
    夜夜夜夜夜久久久久| 中文字幕高清在线视频| 可以免费在线观看a视频的电影网站| 精品久久蜜臀av无| 亚洲精品国产色婷婷电影| 91九色精品人成在线观看| 亚洲专区国产一区二区| 好男人电影高清在线观看| 国产精品国产av在线观看| 国产男女超爽视频在线观看| 欧美日韩视频精品一区| 老司机福利观看| 高清av免费在线| 色婷婷av一区二区三区视频| 精品人妻在线不人妻| 国产xxxxx性猛交| 男人舔女人的私密视频| 国产无遮挡羞羞视频在线观看| 一区二区三区乱码不卡18| 欧美乱码精品一区二区三区| 色综合婷婷激情| 亚洲色图综合在线观看| 怎么达到女性高潮| 国产激情久久老熟女| 欧美日韩成人在线一区二区| 另类亚洲欧美激情| avwww免费| 欧美av亚洲av综合av国产av| 老汉色∧v一级毛片| 国产精品99久久99久久久不卡| 丁香六月欧美| 久久青草综合色| 99精品在免费线老司机午夜| 亚洲午夜理论影院| 高清欧美精品videossex| 18禁观看日本| 国产91精品成人一区二区三区 | 亚洲久久久国产精品| 国产在线免费精品| 国产av精品麻豆| 在线永久观看黄色视频| 一级毛片精品| 免费在线观看日本一区| 成人18禁在线播放| 大片电影免费在线观看免费| 免费在线观看黄色视频的| 欧美人与性动交α欧美精品济南到| 欧美乱妇无乱码| 波多野结衣一区麻豆| 超色免费av| 少妇的丰满在线观看| 国产1区2区3区精品| 国产无遮挡羞羞视频在线观看| 久久精品aⅴ一区二区三区四区| 少妇裸体淫交视频免费看高清 | 欧美成人午夜精品| 午夜福利乱码中文字幕| 亚洲av第一区精品v没综合| 99久久精品国产亚洲精品| 国精品久久久久久国模美| 国产亚洲一区二区精品| kizo精华| 亚洲精华国产精华精| 一本—道久久a久久精品蜜桃钙片| 日日摸夜夜添夜夜添小说| 天堂俺去俺来也www色官网| 久热这里只有精品99| 久久精品国产亚洲av高清一级| 极品人妻少妇av视频| 日日爽夜夜爽网站| 成人av一区二区三区在线看| 岛国毛片在线播放| 国产精品美女特级片免费视频播放器 | 成人三级做爰电影| 成人永久免费在线观看视频 | 久久久久国产一级毛片高清牌| 黑人操中国人逼视频| 999久久久国产精品视频| 一二三四社区在线视频社区8| 亚洲一区二区三区欧美精品| av天堂久久9| 一级a爱视频在线免费观看| av在线播放免费不卡| 人妻一区二区av| 成人亚洲精品一区在线观看| 成年人黄色毛片网站| 久久香蕉激情| 黄色毛片三级朝国网站| 午夜视频精品福利| av有码第一页| 亚洲精品乱久久久久久| 高清欧美精品videossex| 国产亚洲精品第一综合不卡| 中文字幕另类日韩欧美亚洲嫩草| 久久精品国产99精品国产亚洲性色 | 亚洲精品成人av观看孕妇| 精品视频人人做人人爽| 久久久久久久久久久久大奶| 国产高清videossex| 九色亚洲精品在线播放| 久久久精品免费免费高清| 久久99一区二区三区| 18禁裸乳无遮挡动漫免费视频| 午夜91福利影院| 18禁观看日本| 性高湖久久久久久久久免费观看| 精品久久久久久电影网| 日日爽夜夜爽网站| 精品人妻1区二区| videosex国产| 亚洲精品乱久久久久久| 国产精品国产高清国产av | 亚洲欧美色中文字幕在线| 精品欧美一区二区三区在线| 狂野欧美激情性xxxx| 亚洲国产中文字幕在线视频| 欧美大码av| 久久青草综合色| 亚洲va日本ⅴa欧美va伊人久久| 在线 av 中文字幕| 美女视频免费永久观看网站| 一本一本久久a久久精品综合妖精| 满18在线观看网站| 午夜精品久久久久久毛片777| 成人手机av| 国产一区有黄有色的免费视频| 亚洲五月色婷婷综合| 麻豆乱淫一区二区| 老汉色∧v一级毛片| 亚洲全国av大片| 国产成人系列免费观看| 免费观看av网站的网址| 亚洲人成电影观看| 咕卡用的链子| 亚洲国产欧美一区二区综合| 日本黄色日本黄色录像| 欧美+亚洲+日韩+国产| 国产精品1区2区在线观看. | 国产在线观看jvid| 午夜视频精品福利| 在线观看舔阴道视频| 一本—道久久a久久精品蜜桃钙片| 欧美激情高清一区二区三区| 国产免费av片在线观看野外av| 999精品在线视频| 精品国产亚洲在线| 亚洲伊人色综图| 久久精品国产亚洲av香蕉五月 | 亚洲精品美女久久av网站| 久久中文字幕一级| 国产亚洲欧美精品永久| 亚洲欧美日韩另类电影网站| 美女高潮到喷水免费观看| 国产极品粉嫩免费观看在线| 亚洲欧美一区二区三区久久| 国产高清激情床上av| 久久精品成人免费网站| av天堂久久9| 美女视频免费永久观看网站| 香蕉丝袜av| 亚洲成国产人片在线观看| 99国产精品一区二区三区| 天天影视国产精品| 又大又爽又粗| 成人精品一区二区免费| 一个人免费看片子| 一级,二级,三级黄色视频| 最近最新中文字幕大全电影3 | 亚洲精品在线美女| 国产三级黄色录像| 极品人妻少妇av视频| 午夜福利影视在线免费观看| 午夜福利欧美成人| 亚洲av国产av综合av卡| 高潮久久久久久久久久久不卡| 国产精品久久久久久精品古装| 精品少妇黑人巨大在线播放| 首页视频小说图片口味搜索| 精品少妇一区二区三区视频日本电影| 成年人免费黄色播放视频| 中文字幕人妻熟女乱码| 一区二区三区乱码不卡18| 亚洲第一av免费看| 一边摸一边抽搐一进一小说 | 9热在线视频观看99| 少妇被粗大的猛进出69影院| 如日韩欧美国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 亚洲美女黄片视频| av超薄肉色丝袜交足视频| 国产高清videossex| 超碰97精品在线观看| 大香蕉久久网| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成a人片在线一区二区| 欧美久久黑人一区二区| 大片免费播放器 马上看| 少妇裸体淫交视频免费看高清 | 9191精品国产免费久久| 老司机福利观看| av网站免费在线观看视频| 久久99一区二区三区| 久久午夜亚洲精品久久| 菩萨蛮人人尽说江南好唐韦庄| 国产激情久久老熟女| 亚洲色图综合在线观看| 老司机亚洲免费影院| 大型黄色视频在线免费观看| 国产精品二区激情视频| 成年人黄色毛片网站| 少妇猛男粗大的猛烈进出视频| 激情视频va一区二区三区| 在线天堂中文资源库| 国产成人一区二区三区免费视频网站| 亚洲av成人一区二区三| av欧美777| 日韩大码丰满熟妇| 男男h啪啪无遮挡| 国产成人系列免费观看| 大香蕉久久成人网| 国产免费视频播放在线视频| 亚洲国产看品久久| 国产成人影院久久av| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 亚洲性夜色夜夜综合| 如日韩欧美国产精品一区二区三区| 欧美中文综合在线视频| 欧美人与性动交α欧美软件| 丰满饥渴人妻一区二区三| 亚洲五月色婷婷综合| 90打野战视频偷拍视频| 国产老妇伦熟女老妇高清| 黑人欧美特级aaaaaa片| av一本久久久久| 成人三级做爰电影| 国产精品亚洲av一区麻豆| 十分钟在线观看高清视频www| 亚洲伊人久久精品综合| 日本vs欧美在线观看视频| 黑人欧美特级aaaaaa片| 激情视频va一区二区三区| 国产午夜精品久久久久久| 婷婷丁香在线五月| 悠悠久久av| 久久久久久久精品吃奶| 黄色毛片三级朝国网站| 国产免费av片在线观看野外av| 99热国产这里只有精品6| 激情在线观看视频在线高清 | 叶爱在线成人免费视频播放| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 欧美成狂野欧美在线观看| 亚洲欧美一区二区三区久久| www.999成人在线观看| 久久香蕉激情| 中文字幕高清在线视频| 国产日韩欧美视频二区| 美国免费a级毛片| 日韩制服丝袜自拍偷拍| 亚洲av日韩在线播放| 这个男人来自地球电影免费观看| 国产一区二区 视频在线| 岛国在线观看网站| 国产主播在线观看一区二区| 最新的欧美精品一区二区| 精品国内亚洲2022精品成人 | 欧美大码av| 黄色成人免费大全| 亚洲国产毛片av蜜桃av| 高清黄色对白视频在线免费看| 亚洲avbb在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产成人系列免费观看| 亚洲国产精品一区二区三区在线| 日韩有码中文字幕| 精品卡一卡二卡四卡免费| www.自偷自拍.com| 国产欧美日韩一区二区三| 久久久国产精品麻豆| 麻豆国产av国片精品| 精品国内亚洲2022精品成人 | 亚洲av片天天在线观看| 国产不卡一卡二| 国产精品美女特级片免费视频播放器 | 精品人妻在线不人妻| 三级毛片av免费| a级毛片黄视频| 免费高清在线观看日韩| 色视频在线一区二区三区| 国产日韩欧美在线精品| 老司机午夜十八禁免费视频| 狠狠精品人妻久久久久久综合| 久久免费观看电影| 国产成+人综合+亚洲专区| 美女高潮喷水抽搐中文字幕| 成人影院久久| 纯流量卡能插随身wifi吗| 欧美日韩亚洲综合一区二区三区_| 亚洲九九香蕉| 天堂中文最新版在线下载| 日韩一卡2卡3卡4卡2021年| 搡老熟女国产l中国老女人| 悠悠久久av| 少妇裸体淫交视频免费看高清 | 人成视频在线观看免费观看| 国产成人av激情在线播放| 嫩草影视91久久| 久久精品人人爽人人爽视色| 人成视频在线观看免费观看| 首页视频小说图片口味搜索| 又黄又粗又硬又大视频| 亚洲伊人色综图| 大陆偷拍与自拍| 国产成人免费观看mmmm| 777久久人妻少妇嫩草av网站| 久久国产精品大桥未久av| 国产男女内射视频| 亚洲成人免费电影在线观看| 亚洲少妇的诱惑av| 国产成人欧美| 少妇精品久久久久久久| 少妇裸体淫交视频免费看高清 | 人人妻人人爽人人添夜夜欢视频| 亚洲国产欧美一区二区综合| 丰满迷人的少妇在线观看| 久久99热这里只频精品6学生| svipshipincom国产片| 欧美日韩精品网址| 国产又色又爽无遮挡免费看| 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区 | 少妇 在线观看| 欧美黄色片欧美黄色片| 国产一区二区激情短视频| 久9热在线精品视频| 国产成人免费无遮挡视频| 免费看a级黄色片| 精品国产乱码久久久久久男人| 夜夜夜夜夜久久久久| 9热在线视频观看99| 日韩欧美三级三区| 成年人黄色毛片网站| 黄色毛片三级朝国网站| 日韩一区二区三区影片| 变态另类成人亚洲欧美熟女 | 国产精品一区二区精品视频观看| 日韩 欧美 亚洲 中文字幕| 91成人精品电影| 夜夜爽天天搞| 一本综合久久免费| 久久久久久免费高清国产稀缺| 91成年电影在线观看| 欧美日韩福利视频一区二区| 一区二区av电影网| 老司机深夜福利视频在线观看| 一进一出好大好爽视频| 亚洲av国产av综合av卡| a级片在线免费高清观看视频| 日韩欧美免费精品| 日日夜夜操网爽| 人妻 亚洲 视频| 免费在线观看黄色视频的| 国产日韩一区二区三区精品不卡| 韩国精品一区二区三区| 免费观看人在逋| 日韩大片免费观看网站| 一本大道久久a久久精品| 女人被躁到高潮嗷嗷叫费观| 午夜福利在线免费观看网站| 高清av免费在线| 午夜福利影视在线免费观看| 亚洲精品美女久久av网站| 首页视频小说图片口味搜索| 国产成人欧美| 一级片免费观看大全| 美国免费a级毛片| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 成年人午夜在线观看视频| 男女床上黄色一级片免费看| 如日韩欧美国产精品一区二区三区| 久久热在线av| 亚洲国产毛片av蜜桃av| 香蕉国产在线看| netflix在线观看网站| 两个人免费观看高清视频| 亚洲人成电影免费在线| 久久精品亚洲精品国产色婷小说| 国产不卡av网站在线观看| 热99re8久久精品国产| 美女国产高潮福利片在线看| 91麻豆精品激情在线观看国产 | 免费在线观看影片大全网站| 国产男女超爽视频在线观看| 亚洲 欧美一区二区三区| 丁香欧美五月| 国产成人影院久久av| netflix在线观看网站| 人妻 亚洲 视频| 91麻豆精品激情在线观看国产 | 日韩免费高清中文字幕av| 欧美变态另类bdsm刘玥| 性少妇av在线| 热re99久久国产66热| 99re在线观看精品视频| 肉色欧美久久久久久久蜜桃| 亚洲欧洲精品一区二区精品久久久| 自线自在国产av| av超薄肉色丝袜交足视频| 亚洲综合色网址| 亚洲精品在线观看二区| 1024视频免费在线观看| 香蕉久久夜色| 一二三四社区在线视频社区8| 91麻豆av在线| 日韩一区二区三区影片| 老汉色av国产亚洲站长工具| av网站在线播放免费| 久久人妻av系列| 女人被躁到高潮嗷嗷叫费观| 天天操日日干夜夜撸| 久久国产精品人妻蜜桃| 丝袜喷水一区| 黑丝袜美女国产一区| 日本一区二区免费在线视频| 亚洲成人手机| 久久天堂一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 久热爱精品视频在线9| 男女之事视频高清在线观看| 青青草视频在线视频观看| 久久精品国产亚洲av香蕉五月 | 久久九九热精品免费| 人人澡人人妻人| 日韩 欧美 亚洲 中文字幕| 国产精品久久久久久人妻精品电影 | 国产成人精品久久二区二区免费| 亚洲成av片中文字幕在线观看| 日韩熟女老妇一区二区性免费视频| 午夜福利乱码中文字幕| 无人区码免费观看不卡 | 中亚洲国语对白在线视频| av一本久久久久| 十八禁高潮呻吟视频| 久久精品成人免费网站| 91国产中文字幕| 午夜91福利影院| 免费在线观看完整版高清| 日韩欧美三级三区| 热99re8久久精品国产| 成年版毛片免费区| 岛国毛片在线播放| √禁漫天堂资源中文www| 国产亚洲精品第一综合不卡| 国产成人av教育| 久久婷婷成人综合色麻豆| 一边摸一边抽搐一进一小说 | 国产精品久久久久成人av| 久久精品人人爽人人爽视色| 亚洲欧美日韩高清在线视频 | 搡老乐熟女国产| 亚洲国产看品久久| 国产精品亚洲av一区麻豆| 新久久久久国产一级毛片| 桃红色精品国产亚洲av| 9热在线视频观看99| 亚洲av电影在线进入| 精品第一国产精品| 新久久久久国产一级毛片| 欧美人与性动交α欧美精品济南到| 日韩中文字幕欧美一区二区| 99精品欧美一区二区三区四区| www.精华液| 乱人伦中国视频| 欧美人与性动交α欧美精品济南到| 午夜两性在线视频| 一级片免费观看大全| 91九色精品人成在线观看| 岛国毛片在线播放| 一级毛片电影观看| 午夜福利一区二区在线看| 男女高潮啪啪啪动态图| 国产亚洲精品第一综合不卡| 久久九九热精品免费| 日韩中文字幕欧美一区二区| 国产成人精品久久二区二区91| 亚洲avbb在线观看| 欧美黑人精品巨大| av一本久久久久| 日韩视频在线欧美| 最新在线观看一区二区三区| 久久久久久久大尺度免费视频| 欧美日韩成人在线一区二区| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 色94色欧美一区二区| 国产精品一区二区在线不卡| 巨乳人妻的诱惑在线观看| 久久久精品国产亚洲av高清涩受| 麻豆成人av在线观看| 满18在线观看网站| 美女国产高潮福利片在线看| 狠狠婷婷综合久久久久久88av| 欧美 亚洲 国产 日韩一| 精品国产一区二区三区久久久樱花| 人人妻人人添人人爽欧美一区卜| 久久毛片免费看一区二区三区| 精品一区二区三区四区五区乱码| 午夜视频精品福利| 高清毛片免费观看视频网站 | 国产精品.久久久| 国产片内射在线| 久久久欧美国产精品| 国产人伦9x9x在线观看| 久久久精品免费免费高清| 欧美另类亚洲清纯唯美| 午夜免费成人在线视频| 69精品国产乱码久久久| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 精品亚洲成国产av| 99久久精品国产亚洲精品| 成人国产一区最新在线观看| 男人操女人黄网站| 在线观看www视频免费| 又大又爽又粗| 一边摸一边做爽爽视频免费| 一级,二级,三级黄色视频| 少妇 在线观看| 亚洲五月婷婷丁香| 欧美日韩黄片免| 欧美日韩一级在线毛片| 韩国精品一区二区三区| av片东京热男人的天堂| 国产精品久久久av美女十八| 一本—道久久a久久精品蜜桃钙片| 男女午夜视频在线观看| 精品久久蜜臀av无| 成年动漫av网址| av片东京热男人的天堂| 97在线人人人人妻| 亚洲成a人片在线一区二区| 国产一区有黄有色的免费视频| 国产男女内射视频| 两个人免费观看高清视频| 在线av久久热| 久久久国产欧美日韩av| 在线av久久热| 老司机靠b影院| 免费观看av网站的网址| 亚洲欧美一区二区三区久久| h视频一区二区三区| bbb黄色大片| 大型黄色视频在线免费观看| 日本一区二区免费在线视频| 久久久久久免费高清国产稀缺| 亚洲人成电影观看| 在线观看免费视频日本深夜| 亚洲黑人精品在线| 精品国产亚洲在线| 一区在线观看完整版| 韩国精品一区二区三区| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 精品一区二区三区av网在线观看 | 国产成人av教育| 久久久久久久国产电影| 久久这里只有精品19| 日韩制服丝袜自拍偷拍| 国产成人免费观看mmmm| 精品亚洲成国产av| 免费不卡黄色视频| 涩涩av久久男人的天堂| 首页视频小说图片口味搜索| 国产精品久久久av美女十八| 国产男女内射视频| 欧美午夜高清在线| 国产一区二区 视频在线| 在线天堂中文资源库| 一级毛片精品| 露出奶头的视频| 免费日韩欧美在线观看| 一级片免费观看大全| 日韩免费av在线播放| 久久久精品94久久精品| videosex国产| 成年版毛片免费区| 成人黄色视频免费在线看| 国产在线观看jvid| 淫妇啪啪啪对白视频| 精品免费久久久久久久清纯 | 免费久久久久久久精品成人欧美视频| 性高湖久久久久久久久免费观看| 亚洲专区国产一区二区| 性少妇av在线| 国产淫语在线视频| 亚洲av片天天在线观看| 久久性视频一级片| 欧美+亚洲+日韩+国产| 正在播放国产对白刺激| 中文字幕人妻丝袜制服| 侵犯人妻中文字幕一二三四区| 十分钟在线观看高清视频www| 精品人妻熟女毛片av久久网站| 国产老妇伦熟女老妇高清| 三级毛片av免费| 国产伦理片在线播放av一区| av不卡在线播放| 久久精品亚洲av国产电影网| 国产成人系列免费观看| 黄色毛片三级朝国网站| 99riav亚洲国产免费|