亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜    

查看: 4784  |  回復(fù): 9

pariah

銅蟲 (小有名氣)

[求助] 如何讓Pwscf找出全部對稱性

BaMgF4晶體在用Xcrysden建模導(dǎo)出坐標(biāo)后,進(jìn)行scf計算,輸入文件如下:
CODE:
&control
    calculation='scf'
    restart_mode='from_scratch',
    tstress = .true.
    tprnfor = .true.
    prefix='bamgf',
    pseudo_dir = '/home/Xiaole/Desktop/Pwscf/espresso-5.0/pseudo/',
    outdir='/home/Xiaole/tmp/'
/
&system
    nosym = .false.,
    ibrav=  0,
    celldm(1) = 1.889725988,
    nat = 12 ,
    ntyp = 3,
    ecutwfc = 50.0
/
&electrons
    conv_thr =  1.0d-8
    mixing_beta = 0.7
/
ATOMIC_SPECIES
Ba 137.33  Ba.pz-n-nc.UPF
Mg  24.31  Mg.pz-n-nc.UPF
  F  19.00   F.pz-n-nc.UPF
ATOMIC_POSITIONS
Ba      0.0000000000    -2.1481299780    -2.6947971870
Mg     -2.0594950000    -1.2510149000     0.0000000000
  F     -2.0594950000    -2.7999593600    -1.1448674700
  F     -2.0594950000    -2.3950065600     1.5144795060
  F     -2.0594950000    -0.3948289800    -1.7911073820
  F      0.0000000000    -1.1295290600    -0.0586962510
Ba      0.0000000000     2.1481299780     0.2109578130
Mg      0.0000000000    -5.9802851000    -2.9057550000
  F      0.0000000000    -4.4313406400     1.7608875300
  F      0.0000000000    -4.8362934400    -1.3912754940
  F      0.0000000000    -6.8364710200     1.1146476180
  F      0.0000000000     1.1295290600     2.8470587490
K_POINTS {automatic}
6 6 6  1 1 1
CELL_PARAMETERS
    2.0594950000    7.2313000000    0.0000000000
   -2.0594950000    7.2313000000    0.0000000000
    0.0000000000    0.0000000000    5.8115100000

可是輸出中只找到了4個對稱操作(應(yīng)該是8個)。
附上BaMgF4的cif文件。
回復(fù)此樓

» 本帖附件資源列表

  • 歡迎監(jiān)督和反饋:小木蟲僅提供交流平臺,不對該內(nèi)容負(fù)責(zé)。
    本內(nèi)容由用戶自主發(fā)布,如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題,其責(zé)任在于用戶本人,如對版權(quán)有異議,請聯(lián)系郵箱:xiaomuchong@tal.com
  • 附件 1 : BaMgF4.cif
  • 2012-05-22 10:42:30, 10.55 K

» 收錄本帖的淘帖專輯推薦

pwscf 建模 第一性原理

» 猜你喜歡

» 本主題相關(guān)價值貼推薦,對您同樣有幫助:

已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

pariah

銅蟲 (小有名氣)

ATOMIC_POSITIONS {angstrom} 這一行是這樣的
2樓2012-05-22 10:44:43
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

just_play

至尊木蟲 (正式寫手)

【答案】應(yīng)助回帖

感謝參與,應(yīng)助指數(shù) +1
輸出文件中有Warning嗎?是否是分?jǐn)?shù)平移和FFT格點不匹配導(dǎo)致的?
So Trivial !
3樓2012-05-22 11:43:29
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

pariah

銅蟲 (小有名氣)

引用回帖:
3樓: Originally posted by just_play at 2012-05-22 11:43:29:
輸出文件中有Warning嗎?是否是分?jǐn)?shù)平移和FFT格點不匹配導(dǎo)致的?

輸出中并沒有warning。
4樓2012-05-22 12:29:49
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

pariah

銅蟲 (小有名氣)

這是輸出文件
CODE:
     Program PWSCF v.5.0        starts on 22May2012 at 12:28:34

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
          URL http://www.quantum-espresso.org",
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/quote.php

     Parallel version (MPI), running on     1 processors

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
     Reading input from standard input

     Subspace diagonalization in iterative solution of the eigenvalue problem:
     a serial algorithm will be used


     G-vector sticks info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Sum        1689    1689    517                55793    55793    9395



     bravais-lattice index     =            0
     lattice parameter (alat)  =       1.8897  a.u.
     unit-cell volume          =    1168.1339 (a.u.)^3
     number of atoms/cell      =           12
     number of atomic types    =            3
     number of electrons       =        64.00
     number of Kohn-Sham states=           32
     kinetic-energy cutoff     =      50.0000  Ry
     charge density cutoff     =     200.0000  Ry
     convergence threshold     =      1.0E-08
     mixing beta               =       0.7000
     number of iterations used =            8  plain     mixing
     Exchange-correlation      = LDA ( 1 1 0 0 0)
     EXX-fraction              =        0.00

     celldm(1)=   1.889726  celldm(2)=   0.000000  celldm(3)=   0.000000
     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000

     crystal axes: (cart. coord. in units of alat)
               a(1) = (   2.059495   7.231300   0.000000 )  
               a(2) = (  -2.059495   7.231300   0.000000 )  
               a(3) = (   0.000000   0.000000   5.811510 )  

     reciprocal axes: (cart. coord. in units 2 pi/alat)
               b(1) = (  0.242778  0.069144  0.000000 )  
               b(2) = ( -0.242778  0.069144  0.000000 )  
               b(3) = (  0.000000  0.000000  0.172072 )  


     PseudoPot. # 1 for Ba read from file:
     /home/Xiaole/Desktop/Pwscf/espresso-5.0/pseudo/Ba.pz-n-nc.UPF
     MD5 check sum: 0f03eb6889fec9f655a22505e1463c59
     Pseudo is Norm-conserving + core correction, Zval =  2.0
     Generated using "atomic" code by A. Dal Corso  (Quantum ESPRESSO distribution)
     Using radial grid of 1251 points,  1 beta functions with:
                l(1) =   1

     PseudoPot. # 2 for Mg read from file:
     /home/Xiaole/Desktop/Pwscf/espresso-5.0/pseudo/Mg.pz-n-nc.UPF
     MD5 check sum: 163d6ef2207de0fd4d03343463d2e2db
     Pseudo is Norm-conserving + core correction, Zval =  2.0
     Generated using "atomic" code by A. Dal Corso  (Quantum ESPRESSO distribution)
     Using radial grid of 1129 points,  2 beta functions with:
                l(1) =   0
                l(2) =   1

     PseudoPot. # 3 for  F read from file:
     /home/Xiaole/Desktop/Pwscf/espresso-5.0/pseudo/F.pz-n-nc.UPF
     MD5 check sum: 54fb74d952e105b2c7cdd3db7cc48c5b
     Pseudo is Norm-conserving, Zval =  7.0
     Generated using "atomic" code by A. Dal Corso  (Quantum ESPRESSO distribution)
     Using radial grid of 1105 points,  1 beta functions with:
                l(1) =   1

     atomic species   valence    mass     pseudopotential
        Ba             2.00   137.33000     Ba( 1.00)
        Mg             2.00    24.31000     Mg( 1.00)
        F              7.00    19.00000      F( 1.00)

      4 Sym. Ops. (no inversion) found ( 2 have fractional translation)



   Cartesian axes

     site n.     atom                  positions (alat units)
         1           Ba  tau(   1) = (   0.0000000  -2.1481301  -2.6947974  )
         2           Mg  tau(   2) = (  -2.0594952  -1.2510150   0.0000000  )
         3           F   tau(   3) = (  -2.0594952  -2.7999596  -1.1448676  )
         4           F   tau(   4) = (  -2.0594952  -2.3950067   1.5144796  )
         5           F   tau(   5) = (  -2.0594952  -0.3948290  -1.7911075  )
         6           F   tau(   6) = (   0.0000000  -1.1295291  -0.0586963  )
         7           Ba  tau(   7) = (   0.0000000   2.1481301   0.2109578  )
         8           Mg  tau(   8) = (   0.0000000  -5.9802856  -2.9057552  )
         9           F   tau(   9) = (   0.0000000  -4.4313410   1.7608877  )
        10           F   tau(  10) = (   0.0000000  -4.8362938  -1.3912756  )
        11           F   tau(  11) = (   0.0000000  -6.8364715   1.1146477  )
        12           F   tau(  12) = (   0.0000000   1.1295291   2.8470590  )

     number of k points=    52
                       cart. coord. in units 2pi/alat
        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0092593
        k(    2) = (   0.0000000   0.0000000   0.0286787), wk =   0.0185185
        k(    3) = (   0.0000000   0.0000000   0.0573574), wk =   0.0185185
        k(    4) = (   0.0000000   0.0000000  -0.0860362), wk =   0.0092593
        k(    5) = (  -0.0404630   0.0115240   0.0000000), wk =   0.0370370
        k(    6) = (  -0.0404630   0.0115240   0.0286787), wk =   0.0740741
        k(    7) = (  -0.0404630   0.0115240   0.0573574), wk =   0.0740741
        k(    8) = (  -0.0404630   0.0115240  -0.0860362), wk =   0.0370370
        k(    9) = (  -0.0809260   0.0230480   0.0000000), wk =   0.0370370
        k(   10) = (  -0.0809260   0.0230480   0.0286787), wk =   0.0740741
        k(   11) = (  -0.0809260   0.0230480   0.0573574), wk =   0.0740741
        k(   12) = (  -0.0809260   0.0230480  -0.0860362), wk =   0.0370370
        k(   13) = (   0.1213890  -0.0345719   0.0000000), wk =   0.0185185
        k(   14) = (   0.1213890  -0.0345719   0.0286787), wk =   0.0370370
        k(   15) = (   0.1213890  -0.0345719   0.0573574), wk =   0.0370370
        k(   16) = (   0.1213890  -0.0345719  -0.0860362), wk =   0.0185185
        k(   17) = (   0.0000000   0.0230480   0.0000000), wk =   0.0185185
        k(   18) = (   0.0000000   0.0230480   0.0286787), wk =   0.0370370
        k(   19) = (   0.0000000   0.0230480   0.0573574), wk =   0.0370370
        k(   20) = (   0.0000000   0.0230480  -0.0860362), wk =   0.0185185
        k(   21) = (  -0.0404630   0.0345719   0.0000000), wk =   0.0370370
        k(   22) = (  -0.0404630   0.0345719   0.0286787), wk =   0.0740741
        k(   23) = (  -0.0404630   0.0345719   0.0573574), wk =   0.0740741
        k(   24) = (  -0.0404630   0.0345719  -0.0860362), wk =   0.0370370
        k(   25) = (   0.1618520  -0.0230480   0.0000000), wk =   0.0370370
        k(   26) = (   0.1618520  -0.0230480   0.0286787), wk =   0.0740741
        k(   27) = (   0.1618520  -0.0230480   0.0573574), wk =   0.0740741
        k(   28) = (   0.1618520  -0.0230480  -0.0860362), wk =   0.0370370
        k(   29) = (   0.1213890  -0.0115240   0.0000000), wk =   0.0370370
        k(   30) = (   0.1213890  -0.0115240   0.0286787), wk =   0.0740741
        k(   31) = (   0.1213890  -0.0115240   0.0573574), wk =   0.0740741
        k(   32) = (   0.1213890  -0.0115240  -0.0860362), wk =   0.0370370
        k(   33) = (   0.0809260   0.0000000   0.0000000), wk =   0.0185185
        k(   34) = (   0.0809260   0.0000000   0.0286787), wk =   0.0370370
        k(   35) = (   0.0809260   0.0000000   0.0573574), wk =   0.0370370
        k(   36) = (   0.0809260   0.0000000  -0.0860362), wk =   0.0185185
        k(   37) = (   0.0000000   0.0460959   0.0000000), wk =   0.0185185
        k(   38) = (   0.0000000   0.0460959   0.0286787), wk =   0.0370370
        k(   39) = (   0.0000000   0.0460959   0.0573574), wk =   0.0370370
        k(   40) = (   0.0000000   0.0460959  -0.0860362), wk =   0.0185185
        k(   41) = (   0.2023150  -0.0115240   0.0000000), wk =   0.0370370
        k(   42) = (   0.2023150  -0.0115240   0.0286787), wk =   0.0740741
        k(   43) = (   0.2023150  -0.0115240   0.0573574), wk =   0.0740741
        k(   44) = (   0.2023150  -0.0115240  -0.0860362), wk =   0.0370370
        k(   45) = (   0.1618520   0.0000000   0.0000000), wk =   0.0185185
        k(   46) = (   0.1618520   0.0000000   0.0286787), wk =   0.0370370
        k(   47) = (   0.1618520   0.0000000   0.0573574), wk =   0.0370370
        k(   48) = (   0.1618520   0.0000000  -0.0860362), wk =   0.0185185
        k(   49) = (   0.0000000  -0.0691439   0.0000000), wk =   0.0092593
        k(   50) = (   0.0000000  -0.0691439   0.0286787), wk =   0.0185185
        k(   51) = (   0.0000000  -0.0691439   0.0573574), wk =   0.0185185
        k(   52) = (   0.0000000  -0.0691439  -0.0860362), wk =   0.0092593

     Dense  grid:    55793 G-vectors     FFT dimensions: (  64,  64,  50)

     Largest allocated arrays     est. size (Mb)     dimensions
        Kohn-Sham Wavefunctions         3.42 Mb     (   7004,   32)
        NL pseudopotentials             4.06 Mb     (   7004,   38)
        Each V/rho on FFT grid          3.12 Mb     ( 204800)
        Each G-vector array             0.43 Mb     (  55793)
        G-vector shells                 0.06 Mb     (   7586)
     Largest temporary arrays     est. size (Mb)     dimensions
        Auxiliary wavefunctions        13.68 Mb     (   7004,  128)
        Each subspace H/S matrix        0.25 Mb     ( 128, 128)
        Each matrix      0.02 Mb     (     38,   32)
        Arrays for rho mixing          25.00 Mb     ( 204800,   8)

     Initial potential from superposition of free atoms

     starting charge   63.96943, renormalised to   64.00000
     Starting wfc are   48 randomized atomic wfcs

     total cpu time spent up to now is       25.5 secs

     per-process dynamical memory:    44.1 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.00E-02,  avg # of iterations =  2.0

     total cpu time spent up to now is       79.5 secs

     total energy              =    -393.94513839 Ry
     Harris-Foulkes estimate   =    -394.46673437 Ry
     estimated scf accuracy    <       0.79873667 Ry

     iteration #  2     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.25E-03,  avg # of iterations =  2.0

     total cpu time spent up to now is      140.8 secs

     total energy              =    -394.07814873 Ry
     Harris-Foulkes estimate   =    -394.30450242 Ry
     estimated scf accuracy    <       0.41838968 Ry

     iteration #  3     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  6.54E-04,  avg # of iterations =  2.0

     total cpu time spent up to now is      201.8 secs

     total energy              =    -394.17392996 Ry
     Harris-Foulkes estimate   =    -394.18054387 Ry
     estimated scf accuracy    <       0.01433538 Ry

     iteration #  4     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.24E-05,  avg # of iterations =  2.0

     total cpu time spent up to now is      256.7 secs

     total energy              =    -394.17599680 Ry
     Harris-Foulkes estimate   =    -394.17622291 Ry
     estimated scf accuracy    <       0.00046669 Ry

     iteration #  5     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  7.29E-07,  avg # of iterations =  2.0

     total cpu time spent up to now is      314.5 secs

     total energy              =    -394.17610831 Ry
     Harris-Foulkes estimate   =    -394.17610648 Ry
     estimated scf accuracy    <       0.00000775 Ry

     iteration #  6     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  1.21E-08,  avg # of iterations =  2.0

     total cpu time spent up to now is      374.8 secs

     total energy              =    -394.17610903 Ry
     Harris-Foulkes estimate   =    -394.17610974 Ry
     estimated scf accuracy    <       0.00000159 Ry

     iteration #  7     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  2.48E-09,  avg # of iterations =  2.0

     total cpu time spent up to now is      435.0 secs

     total energy              =    -394.17610974 Ry
     Harris-Foulkes estimate   =    -394.17610953 Ry
     estimated scf accuracy    <       0.00000033 Ry

     iteration #  8     ecut=    50.00 Ry     beta=0.70
     Davidson diagonalization with overlap
     ethr =  5.14E-10,  avg # of iterations =  2.0

     total cpu time spent up to now is      492.7 secs

     End of self-consistent calculation

          k = 0.0000 0.0000 0.0000 (  6975 PWs)   bands (ev):

   -20.5506 -20.2556 -19.8993 -19.7958 -19.5932 -19.5287 -19.3678 -19.3159
    -2.6982  -2.2455  -2.2078  -2.0403  -1.7435  -1.6916  -1.6607  -1.1478
    -0.9031  -0.8913  -0.8572  -0.5373  -0.4265  -0.4234  -0.4152  -0.3082
    -0.2797  -0.1921   0.0895   0.0982   0.3130   0.3201   0.4171   0.4982

          k = 0.0000 0.0000 0.0287 (  6969 PWs)   bands (ev):

   -20.5300 -20.2765 -19.8811 -19.7897 -19.5911 -19.5254 -19.3814 -19.3300
    -2.6702  -2.2524  -2.2201  -1.9974  -1.7601  -1.6789  -1.5680  -1.1160
    -1.0173  -0.8726  -0.7924  -0.5966  -0.4513  -0.4040  -0.3421  -0.3013
    -0.2925  -0.1513  -0.0134   0.1216   0.2168   0.2652   0.4337   0.4446

          k = 0.0000 0.0000 0.0574 (  6971 PWs)   bands (ev):

   -20.4753 -20.3307 -19.8389 -19.7832 -19.5781 -19.5248 -19.4061 -19.3655
    -2.5766  -2.3515  -2.1656  -1.9205  -1.8838  -1.4848  -1.4746  -1.2487
    -1.0834  -0.7754  -0.7699  -0.6713  -0.4926  -0.4753  -0.2866  -0.2633
    -0.2329  -0.1953  -0.1213   0.0009   0.1505   0.2229   0.4053   0.4251

          k = 0.0000 0.0000-0.0860 (  6956 PWs)   bands (ev):

   -20.4026 -20.4026 -19.7989 -19.7989 -19.5482 -19.5482 -19.4007 -19.4007
    -2.4563  -2.4563  -2.0457  -2.0457  -1.6534  -1.6534  -1.3228  -1.3228
    -0.9495  -0.9495  -0.6838  -0.6838  -0.5005  -0.5005  -0.2459  -0.2459
    -0.2377  -0.2377  -0.1331  -0.1331   0.2037   0.2037   0.4123   0.4123

          k =-0.0405 0.0115 0.0000 (  6949 PWs)   bands (ev):

   -20.5115 -20.2533 -19.8763 -19.7985 -19.5797 -19.5244 -19.3801 -19.3379
    -2.9215  -2.3550  -2.1706  -1.9920  -1.7266  -1.4658  -1.3735  -1.2931
    -1.0453  -0.8642  -0.8127  -0.6964  -0.5654  -0.4648  -0.2685  -0.1901
    -0.1677  -0.0173   0.0379   0.1648   0.1723   0.3262   0.4101   0.4711

          k =-0.0405 0.0115 0.0287 (  6960 PWs)   bands (ev):

   -20.4942 -20.2718 -19.8616 -19.7941 -19.5777 -19.5256 -19.3872 -19.3476
    -2.8882  -2.4259  -2.1005  -1.9679  -1.6933  -1.4700  -1.3797  -1.2831
    -1.1377  -0.8704  -0.8199  -0.6932  -0.5325  -0.4537  -0.3247  -0.2558
    -0.1878  -0.0335   0.0261   0.1567   0.1848   0.2703   0.4069   0.4520

          k =-0.0405 0.0115 0.0574 (  6966 PWs)   bands (ev):

   -20.4472 -20.3202 -19.8281 -19.7887 -19.5690 -19.5327 -19.3982 -19.3711
    -2.7884  -2.5453  -1.9659  -1.8744  -1.6901  -1.4415  -1.3957  -1.2926
    -1.2203  -0.9135  -0.8454  -0.7275  -0.5182  -0.4149  -0.3520  -0.2725
    -0.2360  -0.0928  -0.0294   0.0307   0.2215   0.2540   0.4091   0.4275

          k =-0.0405 0.0115-0.0860 (  6978 PWs)   bands (ev):

   -20.3839 -20.3839 -19.7987 -19.7986 -19.5507 -19.5507 -19.3935 -19.3935
    -2.6690  -2.6690  -1.8335  -1.8335  -1.6578  -1.6578  -1.2999  -1.2999
    -1.1096  -1.1096  -0.8258  -0.8258  -0.4482  -0.4482  -0.3050  -0.3050
    -0.1847  -0.1847  -0.0728  -0.0728   0.2485   0.2485   0.4082   0.4082

          k =-0.0809 0.0230 0.0000 (  6966 PWs)   bands (ev):

   -20.4357 -20.2506 -19.8265 -19.8079 -19.5519 -19.5220 -19.4006 -19.3804
    -3.2542  -2.3850  -2.1164  -1.9371  -1.6788  -1.4837  -1.4178  -1.2307
    -1.0208  -0.8338  -0.8241  -0.7690  -0.6926  -0.6129  -0.4271  -0.2571
    -0.0399   0.0508   0.0891   0.2339   0.2626   0.3611   0.3693   0.4449

          k =-0.0809 0.0230 0.0287 (  6968 PWs)   bands (ev):

   -20.4241 -20.2643 -19.8202 -19.8047 -19.5606 -19.5389 -19.3859 -19.3725
    -3.2052  -2.5017  -2.0060  -1.9159  -1.6717  -1.4975  -1.4012  -1.2114
    -1.0387  -0.9302  -0.8677  -0.7786  -0.6627  -0.5780  -0.4007  -0.2531
    -0.0879   0.0532   0.0855   0.1663   0.2932   0.3419   0.3707   0.4422

          k =-0.0809 0.0230 0.0574 (  6974 PWs)   bands (ev):

   -20.3922 -20.3004 -19.8070 -19.7988 -19.5759 -19.5664 -19.3636 -19.3587
    -3.0839  -2.7013  -1.8821  -1.8191  -1.6057  -1.5929  -1.3955  -1.1979
    -1.1416  -0.9868  -0.8870  -0.7750  -0.6693  -0.5677  -0.3998  -0.2493
    -0.0952  -0.0172   0.0434   0.1257   0.2990   0.3182   0.3943   0.4351

          k =-0.0809 0.0230-0.0860 (  6954 PWs)   bands (ev):

   -20.3471 -20.3471 -19.7979 -19.7979 -19.5805 -19.5805 -19.3533 -19.3533
    -2.8969  -2.8969  -1.7849  -1.7849  -1.5895  -1.5895  -1.3385  -1.3385
    -1.0608  -1.0608  -0.7774  -0.7774  -0.6552  -0.6552  -0.3147  -0.3147
    -0.0580  -0.0580   0.0836   0.0836   0.3174   0.3174   0.4192   0.4192

          k = 0.1214-0.0346 0.0000 (  6968 PWs)   bands (ev):

   -20.3989 -20.2501 -19.8162 -19.7969 -19.5366 -19.5277 -19.4064 -19.4004
    -3.3873  -2.3070  -2.0696  -1.8576  -1.7398  -1.6883  -1.3488  -1.3065
    -1.1676  -0.8769  -0.7936  -0.7711  -0.6071  -0.5298  -0.5108  -0.2170
     0.0039   0.0610   0.1633   0.2346   0.2632   0.3352   0.4270   0.4381

          k = 0.1214-0.0346 0.0287 (  6980 PWs)   bands (ev):

   -20.3899 -20.2616 -19.8127 -19.7957 -19.5587 -19.5534 -19.3798 -19.3765
    -3.3405  -2.4776  -1.9502  -1.8165  -1.7924  -1.7125  -1.2841  -1.2403
    -1.1719  -0.9489  -0.8163  -0.7514  -0.5998  -0.5346  -0.4982  -0.3000
    -0.0409  -0.0071   0.1738   0.2435   0.2743   0.2792   0.4339   0.4479

          k = 0.1214-0.0346 0.0574 (  6960 PWs)   bands (ev):

   -20.3649 -20.2910 -19.8042 -19.7943 -19.5906 -19.5888 -19.3418 -19.3408
    -3.1953  -2.7286  -1.9172  -1.8146  -1.6771  -1.6445  -1.2886  -1.1503
    -1.1089  -1.0283  -0.8196  -0.6874  -0.6240  -0.6030  -0.3757  -0.3547
    -0.1383  -0.0941   0.1780   0.2096   0.2968   0.3170   0.4294   0.4734

          k = 0.1214-0.0346-0.0860 (  6980 PWs)   bands (ev):

   -20.3293 -20.3293 -19.7969 -19.7969 -19.6034 -19.6034 -19.3264 -19.3264
    -2.9931  -2.9931  -1.9109  -1.9109  -1.5421  -1.5421  -1.3294  -1.3294
    -0.9902  -0.9902  -0.7942  -0.7942  -0.6360  -0.6360  -0.3429  -0.3429
    -0.1756  -0.1756   0.1607   0.1607   0.2961   0.2961   0.4493   0.4493

          k = 0.0000 0.0230 0.0000 (  6961 PWs)   bands (ev):

   -20.5430 -20.2596 -19.9099 -19.7890 -19.5940 -19.5286 -19.3708 -19.3116
    -2.7043  -2.2635  -2.2334  -2.0249  -1.7820  -1.6515  -1.5987  -1.1017
    -0.9161  -0.8945  -0.8215  -0.5334  -0.4723  -0.4570  -0.3261  -0.2862
    -0.1892  -0.1847   0.0815   0.0866   0.2995   0.3161   0.3474   0.4898

          k = 0.0000 0.0230 0.0287 (  6973 PWs)   bands (ev):

   -20.5234 -20.2795 -19.8904 -19.7846 -19.5918 -19.5255 -19.3848 -19.3259
    -2.6799  -2.2908  -2.2265  -1.9892  -1.7387  -1.7167  -1.5346  -1.0971
    -0.9882  -0.8726  -0.7954  -0.6050  -0.5063  -0.4220  -0.3530  -0.2748
    -0.1791  -0.1562  -0.0305   0.1101   0.2119   0.2768   0.3573   0.4242

          k = 0.0000 0.0230 0.0574 (  6975 PWs)   bands (ev):

   -20.4709 -20.3316 -19.8444 -19.7814 -19.5793 -19.5255 -19.4084 -19.3618
    -2.5862  -2.3661  -2.1686  -1.9066  -1.8699  -1.5462  -1.4051  -1.2578
    -1.0377  -0.7801  -0.7702  -0.6687  -0.5398  -0.4272  -0.3669  -0.2761
    -0.2022  -0.1926  -0.1006   0.0136   0.1489   0.2373   0.3550   0.3850

          k = 0.0000 0.0230-0.0860 (  6970 PWs)   bands (ev):

   -20.4009 -20.4009 -19.8003 -19.8003 -19.5504 -19.5504 -19.3993 -19.3993
    -2.4677  -2.4677  -2.0516  -2.0516  -1.6540  -1.6540  -1.3117  -1.3117
    -0.9297  -0.9297  -0.6957  -0.6957  -0.5392  -0.5392  -0.2678  -0.2677
    -0.2413  -0.2413  -0.0922  -0.0922   0.1992   0.1992   0.3609   0.3609

          k =-0.0405 0.0346 0.0000 (  6989 PWs)   bands (ev):

   -20.5012 -20.2586 -19.8906 -19.7901 -19.5835 -19.5239 -19.3849 -19.3311
    -2.9494  -2.3723  -2.2489  -1.9734  -1.7046  -1.5293  -1.3114  -1.2922
    -1.0630  -0.8917  -0.8498  -0.6688  -0.5703  -0.3881  -0.2901  -0.2499
    -0.2369  -0.0920   0.0172   0.1368   0.1832   0.2965   0.3552   0.4333

          k =-0.0405 0.0346 0.0287 (  6989 PWs)   bands (ev):

   -20.4850 -20.2758 -19.8736 -19.7873 -19.5820 -19.5252 -19.3918 -19.3408
    -2.9089  -2.4537  -2.1547  -1.9380  -1.6748  -1.5292  -1.3734  -1.2453
    -1.1429  -0.8831  -0.8691  -0.6514  -0.5298  -0.4671  -0.3260  -0.2649
    -0.1956  -0.0853   0.0184   0.1144   0.1734   0.2783   0.3510   0.4118

          k =-0.0405 0.0346 0.0574 (  6972 PWs)   bands (ev):

   -20.4410 -20.3213 -19.8346 -19.7853 -19.5743 -19.5344 -19.4004 -19.3649
    -2.8000  -2.5522  -1.9951  -1.8363  -1.6767  -1.4810  -1.3635  -1.3253
    -1.1755  -0.9718  -0.8137  -0.6658  -0.5909  -0.4584  -0.3090  -0.2288
    -0.1939  -0.0842  -0.0320   0.0351   0.1795   0.2558   0.3480   0.3805

          k =-0.0405 0.0346-0.0860 (  6962 PWs)   bands (ev):

   -20.3813 -20.3813 -19.7993 -19.7993 -19.5557 -19.5557 -19.3904 -19.3904
    -2.6691  -2.6691  -1.8346  -1.8345  -1.6305  -1.6305  -1.3119  -1.3119
    -1.1098  -1.1098  -0.7449  -0.7449  -0.5275  -0.5275  -0.2414  -0.2414
    -0.1518  -0.1518  -0.0249  -0.0249   0.2237   0.2237   0.3517   0.3517

          k = 0.1619-0.0230 0.0000 (  6973 PWs)   bands (ev):

   -20.4300 -20.2530 -19.8325 -19.8039 -19.5593 -19.5203 -19.4030 -19.3742
    -3.2637  -2.3744  -2.1709  -1.8925  -1.6347  -1.4600  -1.3789  -1.3141
    -1.0750  -0.9273  -0.8301  -0.7732  -0.7012  -0.5359  -0.4134  -0.2130
    -0.0032   0.0501   0.1132   0.1891   0.2317   0.3097   0.3552   0.4594

          k = 0.1619-0.0230 0.0287 (  6979 PWs)   bands (ev):

   -20.4190 -20.2663 -19.8250 -19.8010 -19.5677 -19.5388 -19.3872 -19.3673
    -3.2189  -2.4968  -2.0488  -1.8611  -1.6731  -1.4791  -1.3788  -1.2928
    -1.0657  -0.9665  -0.8932  -0.7848  -0.6594  -0.4925  -0.3884  -0.2608
    -0.0341   0.0174   0.0775   0.1659   0.2683   0.2961   0.3488   0.4438

          k = 0.1619-0.0230 0.0574 (  6988 PWs)   bands (ev):

   -20.3885 -20.3010 -19.8088 -19.7958 -19.5824 -19.5689 -19.3632 -19.3553
    -3.0924  -2.6969  -1.8872  -1.8243  -1.6374  -1.5777  -1.3912  -1.2513
    -1.1309  -1.0168  -0.8830  -0.7677  -0.6251  -0.5303  -0.3674  -0.2484
    -0.1166  -0.0620   0.0395   0.1072   0.2926   0.3030   0.3602   0.4110

          k = 0.1619-0.0230-0.0860 (  6980 PWs)   bands (ev):

   -20.3457 -20.3457 -19.7968 -19.7968 -19.5855 -19.5855 -19.3516 -19.3516
    -2.9000  -2.9000  -1.8044  -1.8044  -1.5917  -1.5917  -1.3608  -1.3608
    -1.0944  -1.0944  -0.7379  -0.7379  -0.6460  -0.6460  -0.2812  -0.2812
    -0.1226  -0.1226   0.0469   0.0469   0.3057   0.3057   0.3937   0.3937

          k = 0.1214-0.0115 0.0000 (  6989 PWs)   bands (ev):

   -20.3989 -20.2502 -19.8163 -19.7969 -19.5370 -19.5277 -19.4062 -19.4007
    -3.3974  -2.3098  -2.0791  -1.8894  -1.7192  -1.6910  -1.3384  -1.2965
    -1.2210  -0.8966  -0.7858  -0.7597  -0.6620  -0.5906  -0.4449  -0.1723
    -0.0515  -0.0339   0.2058   0.2327   0.2679   0.3342   0.3835   0.4460

          k = 0.1214-0.0115 0.0287 (  6988 PWs)   bands (ev):

   -20.3900 -20.2615 -19.8127 -19.7958 -19.5592 -19.5531 -19.3797 -19.3767
    -3.3425  -2.4766  -1.9451  -1.8892  -1.7265  -1.7030  -1.3040  -1.2691
    -1.1681  -0.9566  -0.7920  -0.7501  -0.6450  -0.5929  -0.4317  -0.2522
    -0.1295   0.0314   0.1773   0.2052   0.2847   0.3233   0.3783   0.4544

          k = 0.1214-0.0115 0.0574 (  6965 PWs)   bands (ev):

   -20.3650 -20.2911 -19.8042 -19.7944 -19.5911 -19.5884 -19.3418 -19.3408
    -3.1922  -2.7257  -1.9217  -1.8207  -1.6936  -1.6151  -1.2941  -1.1782
    -1.1270  -1.0288  -0.7730  -0.6766  -0.6347  -0.6068  -0.4316  -0.3324
    -0.1847  -0.0251   0.1778   0.2281   0.2528   0.3066   0.4298   0.4659

          k = 0.1214-0.0115-0.0860 (  6960 PWs)   bands (ev):

   -20.3292 -20.3292 -19.7968 -19.7968 -19.6033 -19.6033 -19.3263 -19.3263
    -2.9831  -2.9831  -1.9004  -1.9004  -1.5322  -1.5322  -1.3122  -1.3122
    -1.0368  -1.0368  -0.7134  -0.7134  -0.6309  -0.6309  -0.4128  -0.4128
    -0.1186  -0.1186   0.1990   0.1990   0.2798   0.2798   0.4578   0.4578

          k = 0.0809 0.0000 0.0000 (  6979 PWs)   bands (ev):

   -20.4386 -20.2496 -19.8218 -19.8121 -19.5475 -19.5235 -19.3983 -19.3847
    -3.2448  -2.4001  -2.0818  -1.9690  -1.6883  -1.4849  -1.4381  -1.1920
    -1.0011  -0.8579  -0.8090  -0.8023  -0.6626  -0.6596  -0.4857  -0.1666
    -0.1414  -0.0162   0.1036   0.2283   0.3190   0.3337   0.3897   0.4230

          k = 0.0809 0.0000 0.0287 (  6954 PWs)   bands (ev):

   -20.4266 -20.2633 -19.8161 -19.8083 -19.5564 -19.5395 -19.3841 -19.3761
    -3.1909  -2.5003  -1.9686  -1.9500  -1.6498  -1.4995  -1.4042  -1.1762
    -1.0312  -0.9450  -0.8197  -0.7773  -0.6738  -0.6045  -0.4111  -0.2130
    -0.0995   0.0304   0.0805   0.1835   0.3177   0.3594   0.3931   0.4392

          k = 0.0809 0.0000 0.0574 (  6975 PWs)   bands (ev):

   -20.3940 -20.3002 -19.8052 -19.8012 -19.5725 -19.5655 -19.3631 -19.3611
    -3.0812  -2.7012  -1.9067  -1.7991  -1.6026  -1.5743  -1.3979  -1.1624
    -1.1565  -1.0037  -0.8588  -0.8149  -0.6797  -0.5937  -0.3551  -0.3449
    -0.0590   0.0005   0.0541   0.1284   0.2970   0.3160   0.4051   0.4220

          k = 0.0809 0.0000-0.0860 (  6984 PWs)   bands (ev):

   -20.3481 -20.3481 -19.7987 -19.7987 -19.5782 -19.5782 -19.3545 -19.3544
    -2.9045  -2.9045  -1.7842  -1.7842  -1.5989  -1.5989  -1.3396  -1.3396
    -1.0680  -1.0680  -0.8535  -0.8535  -0.6352  -0.6352  -0.3751  -0.3751
    -0.0681  -0.0681   0.1124   0.1124   0.2881   0.2881   0.4108   0.4108

          k = 0.0000 0.0461 0.0000 (  6973 PWs)   bands (ev):

   -20.5270 -20.2688 -19.9281 -19.7784 -19.5959 -19.5292 -19.3758 -19.3044
    -2.7201  -2.3522  -2.2515  -1.9855  -1.8385  -1.6226  -1.5498  -0.9912
    -0.9528  -0.8743  -0.7442  -0.6223  -0.4825  -0.4496  -0.3000  -0.2176
    -0.1781  -0.0769   0.0030   0.0600   0.0617   0.2812   0.3411   0.4488

          k = 0.0000 0.0461 0.0287 (  6980 PWs)   bands (ev):

   -20.5094 -20.2865 -19.9068 -19.7763 -19.5932 -19.5257 -19.3902 -19.3187
    -2.6917  -2.3631  -2.2318  -1.9567  -1.7821  -1.7119  -1.4755  -1.0588
    -0.9183  -0.8728  -0.7472  -0.6381  -0.5233  -0.4458  -0.3240  -0.2218
    -0.2002  -0.0862  -0.0743   0.0983   0.1056   0.2559   0.2886   0.3756

          k = 0.0000 0.0461 0.0574 (  6987 PWs)   bands (ev):

   -20.4617 -20.3337 -19.8547 -19.7784 -19.5816 -19.5271 -19.4118 -19.3553
    -2.5954  -2.4106  -2.1716  -1.8879  -1.8393  -1.6090  -1.3323  -1.2534
    -0.9330  -0.7822  -0.7465  -0.6632  -0.5895  -0.4771  -0.4274  -0.3180
    -0.2023  -0.1672  -0.0512   0.0527   0.1571   0.2216   0.2374   0.2968

          k = 0.0000 0.0461-0.0860 (  6984 PWs)   bands (ev):

   -20.3975 -20.3975 -19.8030 -19.8030 -19.5546 -19.5546 -19.3964 -19.3964
    -2.4895  -2.4895  -2.0584  -2.0584  -1.6488  -1.6488  -1.3030  -1.3030
    -0.8814  -0.8814  -0.6928  -0.6928  -0.6046  -0.6046  -0.3084  -0.3084
    -0.2613  -0.2613  -0.0036  -0.0036   0.2051   0.2051   0.2723   0.2723

          k = 0.2023-0.0115 0.0000 (  6995 PWs)   bands (ev):

   -20.4901 -20.2644 -19.9029 -19.7835 -19.5867 -19.5234 -19.3884 -19.3253
    -2.9653  -2.3823  -2.2953  -1.9461  -1.6420  -1.5707  -1.3036  -1.2231
    -1.0914  -0.9055  -0.7714  -0.6110  -0.5748  -0.4583  -0.3338  -0.3068
    -0.1801  -0.0915  -0.0247   0.0992   0.1948   0.2785   0.3387   0.4042

          k = 0.2023-0.0115 0.0287 (  6979 PWs)   bands (ev):

   -20.4750 -20.2804 -19.8841 -19.7818 -19.5855 -19.5249 -19.3951 -19.3349
    -2.9192  -2.4664  -2.1813  -1.9119  -1.6150  -1.5616  -1.3675  -1.1738
    -1.1310  -0.9278  -0.8217  -0.6705  -0.5153  -0.4432  -0.2999  -0.2234
    -0.1807  -0.0864  -0.0243   0.0781   0.2042   0.2764   0.3252   0.3668

          k = 0.2023-0.0115 0.0574 (  6971 PWs)   bands (ev):

   -20.4343 -20.3226 -19.8404 -19.7827 -19.5790 -19.5363 -19.4015 -19.3597
    -2.8099  -2.5655  -2.0087  -1.8102  -1.6582  -1.5059  -1.3415  -1.3132
    -1.1403  -0.9905  -0.8284  -0.7434  -0.5307  -0.4511  -0.2923  -0.1654
    -0.1502  -0.0854  -0.0154   0.0082   0.1911   0.2395   0.2907   0.3122

          k = 0.2023-0.0115-0.0860 (  6974 PWs)   bands (ev):

   -20.3788 -20.3788 -19.8002 -19.8002 -19.5606 -19.5606 -19.3877 -19.3877
    -2.6816  -2.6816  -1.8439  -1.8439  -1.6095  -1.6095  -1.3324  -1.3324
    -1.1165  -1.1165  -0.7232  -0.7232  -0.5816  -0.5816  -0.2105  -0.2105
    -0.1089  -0.1089  -0.0127  -0.0127   0.1936   0.1936   0.2899   0.2899

          k = 0.1619 0.0000 0.0000 (  6967 PWs)   bands (ev):

   -20.4270 -20.2542 -19.8348 -19.8024 -19.5628 -19.5197 -19.4036 -19.3717
    -3.2692  -2.3587  -2.1844  -1.8869  -1.5811  -1.4140  -1.3760  -1.3683
    -1.0991  -1.0120  -0.8355  -0.7737  -0.6202  -0.5881  -0.3569  -0.1081
    -0.0217   0.0348   0.0840   0.2155   0.2590   0.2604   0.3369   0.4669

          k = 0.1619 0.0000 0.0287 (  6984 PWs)   bands (ev):

   -20.4165 -20.2673 -19.8268 -19.7996 -19.5711 -19.5389 -19.3874 -19.3651
    -3.2275  -2.4917  -2.0691  -1.8541  -1.6374  -1.4757  -1.3641  -1.3460
    -1.0662  -1.0349  -0.8972  -0.7855  -0.6040  -0.5347  -0.3636  -0.1970
    -0.0243  -0.0159   0.0683   0.1859   0.2507   0.2675   0.3368   0.4488

          k = 0.1619 0.0000 0.0574 (  6984 PWs)   bands (ev):

   -20.3867 -20.3012 -19.8092 -19.7947 -19.5854 -19.5700 -19.3628 -19.3538
    -3.0896  -2.6952  -1.8944  -1.8170  -1.6227  -1.5789  -1.3905  -1.2766
    -1.1128  -1.0518  -0.8719  -0.7751  -0.5585  -0.5100  -0.3300  -0.3088
    -0.1237  -0.0487   0.0362   0.1257   0.2821   0.2829   0.3625   0.4138

          k = 0.1619 0.0000-0.0860 (  6976 PWs)   bands (ev):

   -20.3449 -20.3449 -19.7961 -19.7961 -19.5879 -19.5879 -19.3507 -19.3507
    -2.8972  -2.8972  -1.8212  -1.8212  -1.5721  -1.5721  -1.3700  -1.3700
    -1.1106  -1.1106  -0.7553  -0.7553  -0.5862  -0.5862  -0.2661  -0.2661
    -0.1410  -0.1410   0.0696   0.0696   0.2867   0.2867   0.3961   0.3961

          k = 0.0000-0.0691 0.0000 (  6970 PWs)   bands (ev):

   -20.5184 -20.2741 -19.9366 -19.7738 -19.5970 -19.5292 -19.3777 -19.3012
    -2.7273  -2.3871  -2.2500  -1.9495  -1.8842  -1.6080  -1.5398  -0.9689
    -0.9421  -0.8641  -0.6914  -0.6570  -0.4404  -0.3529  -0.3023  -0.2665
    -0.1331  -0.1059  -0.0243   0.0055   0.0662   0.2260   0.3573   0.4368

          k = 0.0000-0.0691 0.0287 (  7004 PWs)   bands (ev):

   -20.5020 -20.2906 -19.9146 -19.7730 -19.5940 -19.5259 -19.3926 -19.3155
    -2.7091  -2.3992  -2.2424  -1.9388  -1.8389  -1.7037  -1.4583  -1.0574
    -0.8917  -0.8834  -0.6831  -0.6757  -0.4742  -0.4561  -0.4217  -0.2671
    -0.1720  -0.0585  -0.0485  -0.0257   0.0997   0.2220   0.2806   0.3329

          k = 0.0000-0.0691 0.0574 (  6984 PWs)   bands (ev):

   -20.4568 -20.3350 -19.8595 -19.7771 -19.5826 -19.5279 -19.4131 -19.3524
    -2.6006  -2.4303  -2.1699  -1.8863  -1.8252  -1.6253  -1.2990  -1.2498
    -0.8765  -0.7295  -0.7177  -0.6974  -0.5874  -0.5713  -0.4550  -0.2817
    -0.2088  -0.1720   0.0248   0.0398   0.1637   0.2121   0.2309   0.2322

          k = 0.0000-0.0691-0.0860 (  6996 PWs)   bands (ev):

   -20.3958 -20.3958 -19.8045 -19.8045 -19.5567 -19.5567 -19.3950 -19.3950
    -2.5004  -2.5004  -2.0587  -2.0587  -1.6438  -1.6438  -1.2934  -1.2934
    -0.8462  -0.8462  -0.6930  -0.6930  -0.6460  -0.6460  -0.3273  -0.3273
    -0.2832  -0.2832   0.0401   0.0401   0.1961   0.1961   0.2037   0.2037

!    total energy              =    -394.17610980 Ry
     Harris-Foulkes estimate   =    -394.17610981 Ry
     estimated scf accuracy    <          3.9E-09 Ry

     The total energy is the sum of the following terms:

     one-electron contribution =    -247.08953131 Ry
     hartree contribution      =     160.96618263 Ry
     xc contribution           =     -81.44896836 Ry
     ewald contribution        =    -226.60379277 Ry

     convergence has been achieved in   8 iterations

     Forces acting on atoms (Ry/au):

     atom    1 type  1   force =     0.00000000    0.05430351    0.00545503
     atom    2 type  2   force =     0.00000000    0.00746900    0.00737313
     atom    3 type  3   force =     0.00000000    0.00798189    0.03300703
     atom    4 type  3   force =     0.00000000   -0.00245482   -0.06208810
     atom    5 type  3   force =     0.00000000    0.03895402   -0.00090912
     atom    6 type  3   force =     0.00000000    0.01524582    0.01716203
     atom    7 type  1   force =     0.00000000   -0.05430351    0.00545503
     atom    8 type  2   force =     0.00000000   -0.00746900    0.00737313
     atom    9 type  3   force =     0.00000000   -0.00798189    0.03300703
     atom   10 type  3   force =     0.00000000    0.00245482   -0.06208810
     atom   11 type  3   force =     0.00000000   -0.03895402   -0.00090912
     atom   12 type  3   force =     0.00000000   -0.01524582    0.01716203

     Total force =     0.142465     Total SCF correction =     0.000121


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=-1222.83
  -0.00797297   0.00000000   0.00000000      -1172.86      0.00      0.00
   0.00000000  -0.00872263   0.00000000          0.00  -1283.14      0.00
   0.00000000   0.00000000  -0.00824223          0.00      0.00  -1212.47


     Writing output data file bamgf.save

     init_run     :     25.02s CPU     25.37s WALL (       1 calls)
     electrons    :    461.64s CPU    467.20s WALL (       1 calls)
     forces       :      5.35s CPU      5.40s WALL (       1 calls)
     stress       :      6.76s CPU      6.82s WALL (       1 calls)

     Called by init_run:
     wfcinit      :     23.93s CPU     24.26s WALL (       1 calls)
     potinit      :      0.37s CPU      0.38s WALL (       1 calls)

     Called by electrons:
     c_bands      :    416.41s CPU    421.37s WALL (       8 calls)
     sum_band     :     44.95s CPU     45.50s WALL (       8 calls)
     v_of_rho     :      0.13s CPU      0.14s WALL (       9 calls)
     mix_rho      :      0.09s CPU      0.10s WALL (       8 calls)

     Called by c_bands:
     init_us_2    :      2.12s CPU      2.17s WALL (     988 calls)
     cegterg      :    413.78s CPU    417.70s WALL (     416 calls)

     Called by *egterg:
     h_psi        :    299.46s CPU    300.84s WALL (    1301 calls)
     g_psi        :      2.66s CPU      2.68s WALL (     833 calls)
     cdiaghg      :      3.72s CPU      3.74s WALL (    1249 calls)

     Called by h_psi:
     add_vuspsi   :     16.07s CPU     16.10s WALL (    1301 calls)

     General routines
     calbec       :     28.13s CPU     28.26s WALL (    1405 calls)
     fft          :      0.18s CPU      0.18s WALL (      43 calls)
     fftw         :    275.07s CPU    276.40s WALL (   93446 calls)
     davcio       :      0.05s CPU      1.57s WALL (    1404 calls)

     Parallel routines
     fft_scatter  :     33.84s CPU     34.01s WALL (   93489 calls)
     EXX routines

     PWSCF        :  8m19.02s CPU     8m25.07s WALL


   This run was terminated on:  12:37: 0  22May2012            

=------------------------------------------------------------------------------=
   JOB DONE.
=------------------------------------------------------------------------------=

5樓2012-05-22 13:38:50
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

just_play

至尊木蟲 (正式寫手)

【答案】應(yīng)助回帖

引用回帖:
5樓: Originally posted by pariah at 2012-05-22 13:38:50:
這是輸出文件

     Program PWSCF v.5.0        starts on 22May2012 at 12:28:34

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; plea

剛發(fā)現(xiàn)你建模用的是Primitive cell,元胞的對稱性本來就只有4個啊。。。
So Trivial !
6樓2012-05-22 13:59:03
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

pariah

銅蟲 (小有名氣)

引用回帖:
6樓: Originally posted by just_play at 2012-05-22 13:59:03:
剛發(fā)現(xiàn)你建模用的是Primitive cell,元胞的對稱性本來就只有4個啊。。。

對稱性操作的個數(shù)不只跟空間群有關(guān)嗎?跟我去的原胞和單胞也有關(guān)系么?
我用單胞建模卻只找到了2個對稱操作。。
7樓2012-05-22 15:37:12
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

just_play

至尊木蟲 (正式寫手)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
pariah: 金幣+20 2012-05-23 12:26:37
引用回帖:
7樓: Originally posted by pariah at 2012-05-22 15:37:12:
對稱性操作的個數(shù)不只跟空間群有關(guān)嗎?跟我去的原胞和單胞也有關(guān)系么?
我用單胞建模卻只找到了2個對稱操作。。

空間群中包含三種操作:
1.沿基矢整數(shù)倍的平移
2.純點操作
3.螺旋轉(zhuǎn)動和滑移(點操作+分?jǐn)?shù)平移)
第一性軟件輸出的對稱操作個數(shù)是后兩種的和。元胞和晶胞屬于同一個空間群,這三種類型操作的組合是一樣的,但它們對應(yīng)的平移操作不同,從而后兩種操作也未必相同吧?

你的輸出文件里顯示你建模用的是元胞(12個原子)而且也找到所有4個對稱操作了啊,可以用MS驗證下
So Trivial !
8樓2012-05-22 16:09:32
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

pariah

銅蟲 (小有名氣)

內(nèi)容已刪除
9樓2012-05-23 12:26:28
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

gswylq

木蟲 (正式寫手)

引用回帖:
9樓: Originally posted by pariah at 2012-05-23 12:26:28
多謝耐心的回復(fù),在聲子計算中的確找到了對的點群。但對于scf這步實在有些猶疑。我再好好琢磨琢磨。謝謝。...

看來這個5.0的版本改進(jìn)了,正好解決我的帖子的問題。http://www.gaoyang168.com/bbs/viewthread.php?tid=6173480&fpage=1
10樓2013-07-27 17:57:38
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 pariah 的主題更新
最具人氣熱帖推薦 [查看全部] 作者 回/看 最后發(fā)表
[考研] 352求調(diào)劑 +3 大米飯! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 324求調(diào)劑 +6 lucky呀呀呀鴨 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 269專碩求調(diào)劑 +6 金恩貝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 求調(diào)劑 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素顏傾城1988
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 085601調(diào)劑 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 303求調(diào)劑 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 晨昏線與星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 295求調(diào)劑 +4 一志愿京區(qū)211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機(jī)專業(yè)調(diào)劑) +3 123阿標(biāo) 2026-03-17 3/150 2026-03-20 15:58 by 飛行琦
[考研] 0856調(diào)劑,是學(xué)校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 無懈可擊111
[考研] 344求調(diào)劑 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 308求調(diào)劑 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 11408 一志愿西電,277分求調(diào)劑 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
請?zhí)钐幚硪庖?/div>
亚洲一区二区在线激情| 中文字幕一区二区人妻视频| 亚洲自拍偷拍av在线| 在线免费观看视频18| 92麻豆一区二区三区| 天天看片天天摸天天操| jiee日本美女视频网站| www国产亚洲精品久久久| 丰满少妇_区二区三区| 欧美激情视频第一页| 伊人网在线免费观看| 亚洲女人自熨在线视频| 日本国产亚洲欧美色综合| a级黄片免费观看| 欧美日本在线免费视频| 99国产精品国产精品毛片19| 色999日韩偷自拍拍免费| 自拍丝袜国产欧美日韩| 18福利视频在线观看| 国产激情视频在线观看的| 日本四十路人妻熟女| 91九色91在线视频| 日本免费人爱做视频在线观看不卡 | 最新国产午夜激情视频| 日本高清激情乱一区二区三区| 国产精品蝌蚪自拍视频| 午夜8050免费小说| 国产成人91色精品免费看片| 成人午夜麻豆大胆视频| 可以免费观看日韩av| 91精品国产人妻麻豆| 美女把逼扒开让男人桶| 蜜臀久久精品久久久久久av| 欧美操大黑鸡巴视频在线观看| 成人大片男人的天堂| 狠狠操狠狠操狠狠插| 在线观看中文字幕少妇av| 猫咪亚洲中文在线中文字幕| 国产黄色主播网址大全在线播放| 欧美插插插插插插| 2020国产成人精品视频| 丰满少妇人妻一区二区三区蜜桃| 开心激情五月天作爱片| 2020年亚洲男人天堂网| 四虎国产精品国产精品国产精品| 美女露阴道让男人捅| 日本高清在线观看不卡视频| 亚洲欧美成人午夜一区二区| 欧美黄色一区二区三区视频| 不卡一二三区别视频| 国产探花自拍亚洲av| 在线能看视频你懂的| 黄色av 在线观看| 偷拍欧美日韩另类图片| 91香蕉国产亚洲一二三区| 91精品国产91久久久久久密臀| 亚洲欧美日韩电影一区| 69国产精品成人aaaaa片| 天天插天天干天天狠| av天堂新资源在线| 不卡视频在线 欧美日韩| 国产激情一区二区视频| 天天摸天天舔天天操天天日| 欧美亚洲精品色图网站| 午夜福利在线不卡视频| 蜜臀一区二区日韩美女少妇视频| 熟女俱乐部jukujoclub| 熟女一区二区三区综合| 超碰在线观看97资源| 免费24小时人妻视频| 免费中文字幕a级激情| 欧美精品999不卡| 亚洲免费在线不卡视频| 国产精品亚洲精品亚洲| 顶级欧美色妇xxxx| 99久久国产精品免费热| 亚洲欧洲无码一区2区无码| 久久中文字幕av一区二区| 秋霞成人午夜鲁丝一区二区三区| 在线免费观看欧美小视频| 瑟瑟干视频在线观看| 国产美女视频带a∨黄色片| 天天碰天天摸天天搞| 在线免费观看a视频免费| 天天干夜夜撸天天操| 午夜美女福利视频在线| 日本小视频一区二区| 麻豆午夜激情在线观看 | 国产精品久久久久久成人久| 欧美日本国产一区二区| 欧美最新一区二区三区| 人人妻人人澡人人爽97| 色欲天天媓色媓香视频综合网| 日本老熟老熟妇七十路| 91精品视频在线观看视频| 中文字幕观看中文字幕免费 | 久久av色噜噜ai换脸| 亚洲在线观看中文字幕av| 真人一进一出抽搐大尺度视频| 69xx精品久久久久| 老熟女xxxⅹhd老熟女性| 久久99久久99久久97的人| 高清国产美女a一级毛片| 亚洲成人激情在线综合| 狠狠操狠狠操狠狠插| 夜色福利视频免费观看| 午夜8050免费小说| 丰满放荡熟妇在线播放| 欧美日韩成人高清中文网| 亚洲欧美小说中文字幕| 97精品久久久久久无码人妻| 免费看日韩黄视频在线观看| 亚洲国产美女主播在线观看| 视频免费在线观看网站| tobu8日本高清| 亚洲欧美精品日韩偷拍| 91福利高清在线播放| 男人的天堂aⅴ在线| 成人18禁高潮片免费日本| 欧美精品999不卡| 亚洲午夜高清在线观看| 久久sm人妻中出精品一区二区| 夜色17s精品人妻熟女av| 国产精品视频网站污污污| 嗯~嗯~啊啊啊~高潮了软件| 亚洲成年人精品国产| 亚洲三级综合在线观看| 欧美丝袜亚洲国产日韩| 77亚洲视频在线观看| 2020年亚洲男人天堂网| 欧美大胆a级视频秒播| 一区二区三区国产精华液区别大吗| 人妻系列中文字幕大乳丰满人妻| 天天操天天射天天操天天日 | 91系列视频在线播放| 欧美在线视频不卡一区| 欧美日韩亚洲tv不卡久久| 伊人网在线欧美日韩在线| 亚洲日本欧美韩国另类综合| 久久午夜免费鲁丝片| 桃色成人开心激情网| 人妻少妇的va视频| 青娱乐不卡视频在线| 日本高清激情乱一区二区三区| 91超碰国产在线观看| 天天日天天玩天天摸| 亚洲男人的天堂最新网址| 亚洲成人三级黄色片| 欧洲亚洲一区二区三区四区| 十八禁黄色免费污污污亚洲| 亚洲国产美女主播在线观看| 九九热在线精品播放| 国产在线小视频一区二区| 国产资源网站在线播放| 欧美第一激情综合网欧美激情| 亚洲av在线免费播放| 五月婷婷激情视频网| 在线中文字幕人妻av| 国产极品气质外围av| 久久免费视频ww一区| 九色91操最新在线观看网址| 午夜精品老牛av一区二区三区| 欧洲亚洲一区二区三区四区| 性感人妻 中文字幕| 成人精品动漫一区二区| 五月的婷婷综合视频| 人妻女侠被擒受辱记| 91人妻人人做人人爽高清| 久久久久久久久久久久久国产| 精品av天堂毛片久久久| 自拍偷拍 国产激情| 在线观看网站伊人网| 香港日本台湾经典三级| 精久久久久久久久久久久| 国产亚洲综合5388| 网站在线观看蜜臀91| 欧美亚洲另类精品第一页| 亚洲乱码av一区二区蜜桃av | 欧美久久一区二区伊人| 日韩欧美一区二区三区免费看 | 人妻免费视频黄片在线视频| 日韩一级欧美一级片| 亚洲国产日韩a在线欧美| 福利小视频免费在线| 中文字幕一区二区人妻视频| 正在播放麻豆精品一区二区| 日本丰满熟妇浓密多毛| 天天综合久久无人区| 天天干夜夜爽狠狠操| 男人的天堂av中文字幕| 丰满人妻熟女aⅴ一区| 91精品国产欧美在线| 成人av中文字幕在线看| 人妻中文字幕亚洲在线| 国产精品视频网站污污污| 欧美vr专区日韩vr专区| 日日夜夜免费视频精品| 亚洲午夜国产末满十八岁勿进网站| 福利视频免费在线播放| 看女人大BB群伦交| 天天色天天射天天日天天干| 亚洲经典av中文字幕| 亚洲人妻系列在线视频| 自拍丝袜国产欧美日韩| 男人av一区二区三区| 日本少妇人妻中文在线| 天天碰天天摸天天搞| 天天干夜夜操夜夜骑| 夜夜骚av一二三区| 午夜精品久久秘?18免费观看| av 一区二区三区 熟女| 日本东京热最新中文字幕| 亚洲精品激情视频在线观看| 午夜久久久久欠久久久久| 操死你美女在线视频| 99精品久久99久久久久一| 视频在线+欧美十亚洲曰本| 五月的婷婷综合视频| 欧美精品一区二区三区观看| 精品日本少妇久久久| 亚洲美女a级黄色在线播放| 亚洲成人激情在线综合| 99re这里是国产精品首页 | 久久国产精品久精国产爱| 三区美女视频在线观看| 欧美大胆a级视频秒播| 360偷拍蜜桃臀69式| 中国精品人妻一区二区| 日韩人妻精品久久久久| 天天操天天射天天操天天日 | 中文字幕久久久国产| 五十岁熟妇高潮喷水| 日韩无码国产一区二区| 开心五月综合激情婷婷| av无限看熟女人妻另类av| 国产黑色丝袜 在线日韩欧美| 国产精品美女免费视频观看| 极品风骚人妻3p视频| 美女欧美视频在线观看免费| 色噜噜噜噜色噜噜色合久一| 美女精品久久久久久久久| 日韩黄色在线观看网站上| 国产成人情侣激情视频| 伊人情人成综合视频| 日本四十路人妻熟女| 加勒比东京热绿帽人妻多人操| 亚洲综合在线视频在线播放| 91精产国品一二三产区区别网站| 日本男女免费福利视频| 亚洲色视频在线播放网站| 欧美大胆a级视频秒播| 全球高清中文字幕av| 欧美一级特黄大片在线| 国产黄色主播网址大全在线播放 | 久久无码高清免费视频| 欧美一区二区三区视频看| 亚州av嫩草av极品在线观看| 国产极品气质外围av| 亚洲av三级电影在线观看| 亚洲无人区乱码中文字幕一区| 欧美日韩黄片免费在线观看| 超碰在线观看97资源| 亚洲成人激情在线综合| 97人妻av人人澡人人爽| 搞乱在线在线观看视频| 强乱人妻中文字幕日本| 欧美视频亚洲视频在线| 小妹妹爱大棒棒免费观看视频| jizzjizz国产精品传媒| 国产精品网站的黄色| 91精品91久久久久| 国际精品熟女一区二区| 亚洲天堂av最新在线| 国产午夜在线播放视频| 日韩成人在线电影首页| 日韩黄色在线观看网站上| 人妻视频网站快射视频网站| 亚洲宅男噜噜噜66在线观看| 91精品夜夜夜一区二区| av一区二区三区四区五区在线| 91精品国产成人久久久久久| 自拍偷自拍亚洲精品10p| 自拍偷拍 国产激情| 日韩男女视频网站在线观看| 国产av嗯嗯啊啊av| 日韩美精品成人一区二区三区四区 | 成人午夜av电影网| 日本福利网站一区二区| 人妻少妇的va视频| 日本成人福利电影网| 黑人大吊大战亚洲女人。| 国产精品蝌蚪自拍视频| 一区二区在线观看视频网站| 5566熟女人妻人妻| 亚洲全国精品女人久久久| 黑鸡巴肏少妇逼视频| 亚洲色大WWW永久网站| 性感美女极品18禁网站在线| yellow在线亚洲精品一区| 97香蕉久久国产超碰| 宅男噜噜噜66国产在线观看| 中文字幕观看中文字幕免费| 99久9在线视频播放| 18禁男女啪啪啪无遮挡| 国产农村乱子伦精精品视频| 亚洲成人自拍图片网站 | 欧美精品激情在线不卡| 一级毛片特级毛片免费的| 男女69视频在线观看免费| 大屁股熟女一区二区视频| 嗯~嗯~啊啊啊~高潮了软件 | 日韩欧美黄色免费网站| 美女妩媚午夜诱惑网站| 加勒比东京热绿帽人妻多人操| 大成色亚洲一二三区| 天天曰天天摸天天爽| 青娱乐免费视频一二三| 日本一区二区高清av中文| 亚洲精品久久久人妻| 黄在线看片免费人成视频| 国产农村乱子伦精精品视频| 午夜偷拍的视频久久久免费大全| 欧美日韩亚洲tv不卡久久| 日本不卡 中文字幕| 久久视频 在线播放| 青青青国产精品视频| 国产在线小视频一区二区| 美女妩媚午夜诱惑网站| 天天弄天天草天天日天天| 亚洲乱熟女一区二区三区山| 亚洲激情噜噜噜久久久| 亚洲男人天堂最新网址大全| 人妻系列在线免费视频| 一区二区欧美 国产日韩| 97精品人妻免费视频| 夜夜骚av一二三区| 欧美黑人性猛交小矮人| 午夜精品老牛av一区二区三区| 欧美久久一区二区伊人| 69久久夜色精品国产69乱电影| 欧美一区二区三区爽爽| 久久久久国产精品二区| 亚洲精品色图1234| 国产精品福利久久久久| 国语精品视频自产自拍| 日本老女人日比视频| 成人18禁高潮片免费日本| 1级黄色片在线观看| 日本欧美亚洲国产啊啊啊| 欧美日本国产一区二区| 亚洲欧美另类丝袜另类自拍| 美女黄色啊啊啊啊视频| 99精品久久99久久久久一| 欧美情色av在线观看| 精品国产污污污免费入口| 欧美日韩高清片在线观看| 九九六视频,这里只有精品 | 美女妩媚午夜诱惑网站| 亚洲成a人77777| 亚洲综合一区二区三区四区| 黄很色很在线免费视频网站| 久久亚洲国产成人精品麻豆| 大香焦一道本一区二区三区| ysl蜜桃色7425| 国产青青青青草免费在线视频| 欧美日本亚欧在线观看| 天天日天天干天天日天天干天天| 亚洲欧美激情久久久| 99精品久久一区二区| 68福利精品在线视频| 人妻熟女 亚洲 一页二页| 精产国品一二三产品区别97| 亚洲午夜精品一级毛片app| 久久av色噜噜ai换脸| 伊人情人成综合视频| 另类欧美激情校园春色| 亚洲女人自熨在线视频| 国产男女无套?免费网站下载| 亚洲AV无码久久精品国产一区老| 91精品91久久久久| 国产精品无码无卡免费观| 50熟妇一区二区三区| 亚洲国产精品青青草| 亚洲欧洲一区二区三区在线| 99久久久久久亚洲精品免费| 中文字字幕在线精品乱码| 日日夜夜免费视频精品| 日本韩国福利在线播放| 国产一区两区三区福利小视频| 日本香港韩国三级黄色| 日韩免费黄色片在线观看| 亚洲春色av中文字幕| 黑人侵犯人妻森泽佳奈| 亚洲欧美一级特黄大片| 91色哟哟视频在线观看| 精产国品一二三产品区别97| 制服丝袜中文字幕熟女人妻| 丝袜美腿日韩av一区| 午夜亚洲国产精品中字| 久久久精品人妻无码专区不卡| 老司机伊人99久久精品| 中文字幕麻绳捆绑的人妻| 蜜桃tv一区二区三区| 亚欧洲乱码视频一二三区| 91精品麻豆91夜夜骚| 在线成人教育平台排名| 91精品视频在线观看视频| jiee日本美女视频网站| 欧美国产精品久久久免费| 丰满少妇_区二区三区| 91精品视频在线观看视频| 在线视频国产精品欧美| 亚洲宅男噜噜噜66在线观看| 中文字幕在线免费观看人妻| 亚洲天堂色综合久久| av一区二区三区四区五区在线| 青青青国产精品视频| 先锋人妻啪啪中文字幕| 国产 亚洲 欧美 自拍| 顶级欧美色妇4khd| 91激情四射婷婷综合| 美女网站视频久久精品| 亚洲制服丝袜美腿在线| 岛国av成人午夜高清| 福利视频导航在线观看| 韩国在线播放一区二区三区| 加勒比东京热绿帽人妻多人操| 婷婷一区二区三区五月丁| 日韩激情亚洲国产欧美另类激情| 国际日韩日韩日韩日韩日韩 | 日产国产欧美精品另类| 老鸭窝在线毛片观看免费播放| 综合激情网,激情五月| 91色乱一区二区三区| jandara在线观看| 亚洲韩精品一区二区三区| 大香焦一道本一区二区三区| 中文字幕在线免费观看人妻| 美女av色播在线播放| 婷婷综合缴情亚洲五月伊人 | 国产精品午夜无码AV体验区| 国产视频成人自拍蝌蚪视频| 自拍丝袜国产欧美日韩| 久久久人妻免费视频| 38av一区二区三区| 亚洲综合第一区二区| 亚洲制服丝袜资源网| 视频免费在线观看网站| 女人的天堂av在线网| 日韩国产欧美一区二区三区粉嫩| 天天日天天干天天日天天干天天| 在线视频国产精品欧美| 东京热日韩av影片| 中文字幕在线字幕乱码怎么设置| 亚洲人精品午夜射精日韩| 美女黄色啊啊啊啊视频| 情趣视频在线观看91| 日本黄页在线观看视频| 日韩免费黄色片在线观看| 亚洲AV无码一二三四区在线播放| xxxx69在线观看视频| 91日本精产品一区二区三区| 久久精品久久久久观看99水蜜桃| 丝袜美腿日韩av一区| 日本a级2020在线观看| 岳母的诱惑电影在线观看| 亚洲av毛片一区二区三区网| 午夜偷拍的视频久久久免费大全| 欧美熟女xx00视频| 人妻少妇精品二三区| 骚穴被阴茎插免费视频| 天天爱天天日天天爽| 一区二区九日韩美女| 国产极品气质外围av| 538欧美在线观看一区二区三区 | 激情九月天在线视频| 大片a免费观看在线视频观看| 亚洲精品综合欧美精品综合| 西野翔人妻中文字幕中字在| 熟妇高潮久久久久久久| 99久久国产精品免费热| 久久久久久久岛国免费观看| 黄色av 在线观看| 亚洲国产美女主播在线观看| 五十岁熟女高潮喷水| 色网站在线观看免费| 人人妻人人爽人人爽欧美一区| 2021国产在线视频| 蜜桃臀av在线一区二区| 亚洲综合一区二区三区四区| 亚洲av毛片一区二区三区网| 杜达雄啪啪毛片视频| xxnxx国产美女| 亚洲一区在线视频观看地址| 97cao在线视频| 在线免费视频999| 亚洲综合熟女乱中文| 亚洲av网站一区二区三区| 婷婷色综合五月天视频| 夜夜操夜夜爱夜夜摸| 黑人侵犯人妻森泽佳奈| 国产视频成人一区二区| 黑吊操欧美极品美女| 狠狠干狠狠操免费视频| 国产激情免费在线视频| 少妇被中出一区二区| 精久久久久久久久久久久| 91超精品碰国产在线观看| 日本亚洲精品视频在线观看| 亚洲免费午夜污福利| 超碰在线免费观看视频97| 2020国产成人精品视频| 麻豆国产91制片厂| 日本高清 中文字幕| 国产激情免费在线视频| 豆豆专区操逼性视频在线| 九九九九九久久久国产| avtt中文字幕手机版| 亚洲国产日韩精品在线| 蜜臀一区二区日韩美女少妇视频| 国产人妻777人伦精品hd超碰| 欧美性感美女热舞视频| 欧美日本国产一区二区| 鸡巴在里面福利视频在线观看| 中文字幕人妻精品精品| 国产成人情侣激情视频| 4438x亚洲最大的成人| 亚洲一区视频中文字幕在线播放| 中文字幕麻绳捆绑的人妻| aaaa级少妇高潮在线观看| 欧美黄色一区二区三区视频| 超peng视频在线免费播放97| 国产精美视频精品视频精品| 亚洲精品一区二区gif| 国产成人深夜福利短视频99| 一区二区三区免费版在线| 伊人久久综合国产精品| 精品高潮呻吟久久av| 日本久久久久久黄色| 欧美大胆a级视频秒播| 得得爱在线视频观看| 亚洲熟女一区二区三区250p| 亚洲欧美激情国产综合久久久| 青青青在线视频免费播放| 欧美精品激情在线不卡| 妈妈的朋友2中文字幕在线| av成人三级高清日韩| 在线视频自拍第三页| 开心激情五月天作爱片| 国产精品久久人人添| 福利小视频免费在线| 老色鬼精品视频在线观看播放| 不卡在线一区二区三区| 国产福利一区二区三区在线观看| 亚欧洲乱码视频一二三区| 欧美日韩福利视频网| 亚洲成人中文无码在线| 欧美大胆a级视频秒播| 亚洲午夜高清在线观看| 免费看一级高潮喷水片| 大乳丰满人妻中文字幕韩国hd| 中文字幕在线免费观看成人| 欧美亚洲愉拍一区二区三区| 国内自拍第一区二区三区| 在宿舍强奷两个清纯校花| 婷婷六月天在线视频| 18福利视频在线观看| 天天干夜夜撸天天操| 中文字幕亚洲乱码精品无限| 人妻少妇视频系列视频在线| 啊不行啊操逼好爽大鸡吧视频| 国色天香一二三期区别大象| 91久久精品美女高潮喷水白浆| 91中文字幕视频网站| 欧美情色av在线观看| 91色老久久精品偷偷蜜臀| 操烂你的骚逼天天欧美| 国产青青青青草免费在线视频| 在线观看免费啪啪啪| 99国产精品久久99久久久| 视频免费在线观看网站| 免费在线小视频你懂的| 久久无码高清免费视频| 国语精品视频自产自拍| 国内精品一区二区2021在线| 人妻少妇视频系列视频在线| 视频在线 一区二区| ass亚洲熟女ass| 亚洲一级熟妇丰满的女人| 天天操天天射天天操天天日| 日本成人福利电影网| 日韩久久不卡免费视频 | av男人站在线观看| 欧美亚洲精品色图网站| 91麻豆精品国产在线| 亚洲理论在线a中文字幕97| 黑川堇人妻88av| 免费在线小视频你懂的| av激情四射五月婷婷| 国产肥胖熟女又色又爽免费视频| 国际日韩日韩日韩日韩日韩| 情趣视频在线观看91| 美女福利网站在线播放| 人妻在线中文视频视频| 亚洲人妻系列在线视频| 网友自拍第一页99热| 美女av色播在线播放| 4438全国成人免费视频| 天海翼亚洲一区在线观看| 亚洲乱码国产乱码精品精视频| 欧美巨大另类极品video| 交换的一天中文字幕在线视频| 亚洲av毛片一区二区三区网| 公侵犯人妻中文字幕巨| 国产91免费在线观看| 欧美成人性生活视频播放| 亚洲一区二区偷拍女厕所| 亚洲一区视频中文字幕在线播放| 亚洲少妇色小说综合| 欧美亚洲愉拍一区二区三区| 亚洲欧洲一区二区三区在线| 成人午夜麻豆大胆视频| 国产精品蝌蚪自拍视频| 爱搞视频在线观看视频91| 天天操天天干天天舔天天| 中文字幕在线观看av观看| 交换的一天中文字幕在线视频| 丰满人妻被猛烈进入中文字幕| 午夜在线观看一级毛| 五月天男人的天堂中文字幕| 69精品人妻久久久久久久久久久 | 日韩最近中文在线观看| 在线免费视频999| 伊人网国产在线播放| 男女插鸡巴视频软件| 免费观看在线中文字幕视频| 人人妻人人爽人人摸| 啊~插得好快别揉我胸了视频| 人妻激情偷乱一区二区三区av| 狂操鸡巴小骚逼视频免费观看| 天天操天天日天天插天天舔| 岳的大肥屁熟妇五十路| 91佛爷视频在线观看| av天堂新资源在线| 自拍偷自拍亚洲精品10p| 狠狠操深爱婷婷综合一区| 熟妇人妻丰满久久久久久久| 日本亚洲精品视频在线观看| 午夜8050免费小说| 港台美女明星av天堂| 亚洲熟女在线免费观看| 久久99嫩草99久久精品| 久久久西西gogo日本美女人体| 亚洲人妻系列在线视频| 亚洲欧美日韩电影一区| 欧美一区二区播放视频| 日韩激情文学在线视频| 中文字幕 人妻 熟女| 91激情四射婷婷综合| 亚洲成人av在线一区二区| 亚洲人成小说网站色| 熟妇人妻丰满久久久久久久| 欧美一区二区播放视频| 96在线观看免费播放| 啪啪啪网站免费看视频| 国产肥胖熟女又色又爽免费视频| 91佛爷视频在线观看| 亚洲成人,国产精品| 亚洲欧美不卡专业视频| 中文字幕在线免费观看人妻 | 亚洲少妇视频在线观看| 中文字幕中文字幕在线中…一区| 亚洲自拍偷拍av在线| 天天透天天舔天天操| 在线观看网站伊人网| 免费成人av麻豆| 开心激情五月天作爱片| 亚洲AV无码一二三四区在线播放 | 久久久久性感美女偷拍视频| 大尺度久久久久久久| 日本久久久久久黄色| 夜夜爽夜夜操夜夜爱| 亚洲中文字幕最新地址| www国产亚洲精品久久久| 成年人黄色日本视频| 亚洲妹妹我爱你在线观看| 福利视频导航在线观看| 丝袜美腿日韩av一区| 99色在线观看免费观看| 人人妻人人爽人人摸| 青青青青午夜手机国产视频| 夜夜操天天干夜夜操| 欧美老熟妇xxoo老妇| 日本一区二区三区区别| 黑人大巨屌操美女逼| 天天干夜夜操91视频网站| avtt中文字幕手机版| 亚洲成人偷拍自拍在线| 亚洲免费在线不卡视频| 妈妈的朋友中字在线免费观看| 亚洲精品一区二区gif| 美女一区二区四区六区八区| 亚洲女人自熨在线视频| 福利视频导航在线观看| 国产精品午夜无码AV体验区| 日韩人妻一区二区三区在线观看 | 成人免费电影二区三区| av 一区二区三区 熟女| 中文字幕在线观看av观看| 99久久99九九九99九| 亚洲精品中文字幕手机在线免费看 | 亚洲乱码av一区二区蜜桃av| 污视频在线观看地址| 第一福利视频在线观看| 日韩最近中文在线观看| 大屁股熟女一区二区视频| 欧美男女一区二区三区| 在线免费观看a视频免费| 国产美女主播av在线| 中文字幕在线免费观看成人| 欧美一区二区三区视频看| avtt中文字幕手机版| 啪啪啪网站免费在线看| 天天透天天舔天天操| 亚洲黄色成人一级片| 亚洲av手机免费在线| 日韩黄色在线观看网站上 | 精久久久久久久久久久久 | 夜夜爽夜夜操夜夜爱| 久草视频在线看免费| 美女把逼扒开让男人桶| av在线播放观看h| 亚洲激情视频在线观看免费| 国产 亚洲 欧美 自拍| 少妇被粗大的猛进69视频| 大鸡扒操大逼大片免费关看| av男人站在线观看| 亚洲国产精品久久久久久无码| 亚洲avav天堂av在线网毛片| 午夜精品久久久久久久精品乱码| av丝袜免费在线观看| 一区二区三区资源视频| 天天操天天干天天舔天天| 人妻被强av系列一区二区| 9久re热视频在线精品| 天天色 天天操 天天好逼| 亚洲精品一区二区gif| 又爽又粗又猛又色又黄视频| 自拍偷拍 国产激情| 久久av色噜噜ai换脸| 福利美女视频在线观看| 不卡高清一区二区三区| av天堂新资源在线| 午夜久久久久欠久久久久| 黄色av 在线观看| 人妻超清中文字幕在线乱码| 久久99嫩草99久久精品| 亚洲成年人精品国产| caopeng97在线观看视频| 91麻豆精品国产在线| 国产福利一区二区三区在线观看 | 亚洲欧洲一区二区三区在线| 中文字幕欧美一区二区视频| 欧美性受黑人猛交裸体视频 | 中文字幕人妻精品精品| 欧美猛少妇色ⅹⅹⅹⅹⅹ猛叫| 欧美一级日韩一级亚洲一级va| 中文字幕 中文字幕 亚洲| 看女人大BB群伦交| 99久久99九九九99九| 亚洲美女露隐私av一区二区精品 | 四虎精品久久免费最新| 亚洲熟女在线免费观看| 天天看片天天摸天天操| 人妻少妇的va视频| 亚洲激情噜噜噜久久久| —区二区三区女厕偷拍| 美利坚合众国av天堂| 视频免费在线观看网站| 综合久久伊人久久88| 午夜偷拍的视频久久久免费大全| 亚洲欧美另类丝袜另类自拍| 一区二区三区国产在线成人av| 91九色尤物无套内射| 美女张开腿给男人桶爽的软件 | 美女福利视频一区二区三区四区| 91大神福利视频网| 奇米网首页神马久久| 蜜桃tv一区二区三区| 91超碰九色porny| 亚洲无人区乱码中文字幕一区| 91亚洲精品久久蜜桃| 国产激情免费在线视频| 亚洲精品乱码久久久久app | 99久久国产精品免费热| 黄片操操操操操操c| 青青操久久综合激情| 亚洲制服丝袜美腿在线| 天天干天天操天天要| 美女网站福利在线观看| 手机看电影一区二区三区| 91超碰九色porny| 日本欧美高清在线观看视频| 亚洲一区二区精品三区视频| 嗯~嗯~啊啊啊~高潮了软件| 精产国品一二三产品区别97 | 人人妻人人澡人人爽97| 色老头一区二区三区四区五区| 亚洲一区视频中文字幕在线播放 | 9420高清视频在线观看国语版| 后入日韩翘臀蜜桃臀美女| 蜜乳视频一区二区三区| 亚洲va999天堂va| 国产一级一国产一级毛片| 少妇被粗大的猛进69视频| 欧美日韩黄片免费在线观看| 麻豆白洁少妇在线播放| av成人三级高清日韩| 上床啪啪啪免费视频| 成人av中文字幕在线看| 视频在线+欧美十亚洲曰本| 91麻豆精品国产在线| 人妻少妇精品二三区| 午夜亚洲国产精品中字| 天天操天天舔天天射天天日天天干| 日韩一区二区在线播放观看| 国产做A爱免费视频在线观看| 午夜国产一区二区三区| 天天躁狠狠躁狠狠躁性色| 欧美大胆a级视频秒播| 国产av高清二区三区| 高清欧美色欧美综合网站| av 一区二区三区 熟女| 97精品国产91久久久| 欧美vs亚洲vs日韩| 欧美男男在线观看视频网站| 伊人精品久久一区二区| 女女抠逼白虎白丝袜| 欧美肥妇久久久久久| 欧美三区四区在线视频| 日韩一区二区在线播放观看| 国产不卡免费在线观看| 人妻少妇精品二三区| 七色福利视频在线观看| 亚洲av手机免费在线| 内地精品毛片在线观看| 亚洲高清一区二区三区久久| 国产大桥未久一区二区| 久久久久久a女人处女| 亚洲av网站一区二区三区| 99福利一区二区视频| 欧美久久蜜臀蜜桃资源吧| 久久久视频在线播放| 亚洲综合第一区二区| 99久久国产精品免费热| 91精品91久久久久| 精产国品一二三产品区别97| 女生抠逼自慰啊啊啊啊啊啊啊下载| 日本美女爱爱视频网站| 搞乱在线在线观看视频| 快色视频在线观看免费| 69国产精品成人aaaaa片| 妈妈的朋友中字在线免费观看| 色999日韩偷自拍拍免费 | 老熟妇一区二区三区v∧88| 国产,亚洲,欧美综合| 超peng视频在线免费播放97| 日本欧美国产在线一区| 91污污在线观看视频| 在线观看2022av| yellow在线亚洲精品一区| 男人的天堂av中文字幕| 一区二区三区内射美女| 五月天色婷婷狠狠爱| 国产成人情侣激情视频| 91青青青国产免费高清| 亚洲美女a级黄色在线播放| 91色哟哟视频在线观看| 中文乱码字幕人妻熟女人妻| 日本亚洲精品视频在线观看| 五月的婷婷综合视频| av里面的动作是真进去吗| 亚洲 偷拍 自拍 欧美| 成人黄色录像在线观看| 第一福利视频在线观看| 51vv精品视频在线观看| 色视频免费观看网址| 98热视频精品在线观看| 国产伦理二区三区在干嘛呢| 天天干天天操天天日天天日| 天天躁狠狠躁狠狠躁性色| 三区美女视频在线观看| 久久久视频在线播放| 午夜野花视频在线观看| 亚洲av三级电影在线观看| 一级毛片特级毛片免费的| 蜜桃tv一区二区三区| 男人和女人的逼视频| 自拍偷拍 国产激情| 韩日一级人添人人澡人人妻精品| 97人妻人人揉人人躁人人夜夜爽 | 中文字幕福利视频第四页| 精品美女洗澡一区二区| 啊不行啊操逼好爽大鸡吧视频 | av里面的动作是真进去吗| 91激情四射婷婷综合| 杜达雄啪啪毛片视频| 伊人综合在线视频免费观看| 97精品人妻免费视频| 亚洲av激情综合网| 91精品视频在线观看视频| 久久国产精品久精国产爱 | 天天操天天舔天天射天天日天天干| 伊人情人成综合视频| 日本特级黄片免费观看| 上床啪啪啪免费视频| 福利小视频免费在线| 97香蕉久久国产超碰| 在线观看网站伊人网| 亚洲最强的25个城市| 亚欧洲乱码视频一二三区| 丰满放荡熟妇在线播放| 色网站在线观看免费| 女同大尺度视频网站在线观看| 精产国品一二三产品区别97| 成人午夜麻豆大胆视频| 美女扒开逼逼给你看| 1区3区4区产品乱入视频| 美利坚合众国av天堂| 男插女视频大全免费| 一区二区在线观看视频观看| 在线免费视频999| 国产乱码有码一区二区三区| 91九色尤物无套内射| 亚洲少妇视频在线观看| 欧美久久蜜臀蜜桃资源吧| 亚洲综合天堂av网站在线观看| 中文字幕国产一区在线视频| av丝袜免费在线观看| 国产91九色视频在线观看| 五月的婷婷综合视频| 国产乱码有码一区二区三区| 亚洲全国精品女人久久久| 91麻豆精品国产在线| 国产成人在线观看视频播放| 不用付费特黄特色亚洲特级黄色片| 亚洲最大的自拍偷拍网| 亚洲一区二区三区四区入口| 亚洲av日韩久久网站| 青青草一个释放的网站| 91porny九色视频偷拍| 免费观看在线中文字幕视频| 午夜久久久久久av五月| 少妇被粗大的猛进69视频| 亚洲女人自熨在线视频| 国产一区二区手机在线观看| 青青操久久综合激情| 少妇熟女天堂网av| 精品久久久久久久久久久久久| 国产igao激情在线视频入口| 操操操操操操操操操网| 99久久99九九九99九| 在线成人教育平台排名| 欧美大鸡吧男操女啊啊啊视频| 青娱乐这里只有精品| 天天想要天天操天天干| 久久精品国产亚洲av清纯| 区一区二区三免费观看视频| 神马不卡视频在线视频| 午夜国产免费视频亚洲| 婷婷色九月综合激情丁香| 午夜野花视频在线观看| 麻豆午夜激情在线观看| 农村大炕有肉大屁股熟妇| 亚洲午夜高清在线观看| 99999久久久精品| 羞羞漫画无限免费观看秋蝉| 久久免费视频ww一区| 亚洲精品一区二区gif| 亚洲成人欧洲成人在线| av在线播放观看h| 女人扒开逼让男人操| 国产激情免费在线视频| 乱子伦国产一区二区三区| 亚洲综合第一区二区| 一区二区三区四区久久久久韩日| 中文字幕国产一区在线视频| 高潮喷水一区二区三区| 亚洲成人三级黄色片| 亚洲人成大片在线观看| 午夜国产免费视频亚洲| 琪琪日本福利伦理视频| 亚洲欧美另类校园春色| 亚洲国产精品一区二区第二页| 国产漂亮白嫩美女在线图片 | 中文字幕在线字幕乱码怎么设置 | 美国十次了亚洲天堂网国产| 5566熟女人妻人妻| 亚洲一区在线视频观看地址| 日本韩国福利在线播放| 国产欧美福利在线观看| 亚洲国产精品自产拍在线观看| 亚洲综合天堂av网站在线观看| 午夜福利国产精品久久久久| 久久99精品热在线观看| 黄色片黄色片黄色片黄色片黄色| 天天天天天天天天日日日| 黄色网络中文字幕日本| 免费的啪啪视频软件| 亚洲欧美一级特黄大片| 人人妻人人澡人人爽97| 免费的啪啪视频软件| 久久久久国产精品二区| 91porny九色视频偷拍| 午夜在线成人免费电影| 国产男人的天堂一区| 国产女人18毛片水真多精选| 在线视频自拍第三页| 午夜国产成人精品视频观看| 啊不行啊操逼好爽大鸡吧视频| 人妻女侠被擒受辱记| 色老头一区二区三区四区五区 | 中文字幕亚洲无线乱码| 久久99嫩草99久久精品| 人妻系列中文字幕大乳丰满人妻 | 亚洲自拍偷拍av在线| av里面的动作是真进去吗| 中文字幕 首页 人妻| 中文字幕综合网91| 亚洲妹妹我爱你在线观看| 欧美日本在线免费视频| 国产一级一国产一级毛片| 亚洲少妇色小说综合| 日本福利片在线播放| 男人电影天堂在线观看| 日韩A级毛片免费视频| 天天插天天透天天爽| 欧美日韩福利视频网| 亚洲自拍偷拍av在线| 天天操天天舔天天做| 69国产在线视频网站| 99热99这里免费的精品| 瑟瑟干视频在线观看| 黑人黄色免费一级av| 日本男女免费福利视频| 18岁禁一二三区免费体验| 96在线观看免费播放| 桃色成人开心激情网| 在线中文字幕人妻av| 亚洲美女色www色| 182tv精品免费在线观看| 中文字幕在线观看亚洲情色| 99久久久久久久久久久久久| 日本老熟老熟妇七十路| 人妻在线中文视频视频| 中文字幕熟女乱一区二区| 熟妇人妻av无码中文字幕| 天天干夜夜操夜夜骑| 97视频538在线观看| yy4080黄色片| 裸露视频免费在线观看| 日韩无码国产一区二区| 天天干天天操天天日天天日| 丰满放荡熟妇在线播放| 国产主播诱惑毛片av| 国产女主播在线观看一区| 一区二区三区国产在线成人av| 人妻激情综合久久久久蜜桃| 美女把逼扒开让男人桶| 9662av在线视频| 一区二区三区五区六区| 污视频在线观看地址| 精品人妻人人做人人爽| 欧美猛少妇色ⅹⅹⅹⅹⅹ猛叫| 人人妻人人爽人人爽欧美一区| 99久久久久久久久久久久久| 亚洲熟女乱一区二区精品成人| 日本东京热最新中文字幕| 日本久久久久久黄色| 亚洲一区二区精品三区视频| 亚洲蜜桃久久久久久| 久久一级片三上悠亚| 亚洲国产日韩欧美一区二区三区, 精久久久久久久久久久久 | 西野翔人妻中文字幕中字在| 天天干夜夜操夜夜骑| 91超精品碰国产在线观看| jizzjizz国产精品传媒| 亚洲综合在线视频在线播放| julia人妻av一区二区三区| 午夜偷拍的视频久久久免费大全| 亚洲美女午夜激情视频在线观看| 男插女视频大全免费| 一区二区三区av免费天天看| 69国产精品成人aaaaa片| 2018中文字字幕人妻| 天天在线播放日韩av| 性色蜜桃臀x88av天美传媒| 亚洲成人激情在线综合| 精品国产污污污免费入口| 欧美丝袜亚洲国产日韩| 中国精品人妻一区二区| 亚洲人人爽人人澡起碰av| 国产自拍偷拍在线精品| 天天日天天干天天日天天干天天| 亚洲免费在线不卡视频| 97精品久久久久久无码人妻| 亚洲av毛片在在线播放| 538欧美在线观看一区二区三区| 欧美aaaa性bbbbaaaa| 国产一级一国产一级毛片| 亚洲精品1卡2卡3卡| 天天操,天天射,天天爽| 亚洲综合首页综合在线观看| 欧美视频免费观看777| 亚洲国产精品青青草| 99精品久久精品一区二区| 国产资源网站在线播放| 亚洲唯美激情综合四射| 秋霞成人午夜鲁丝一区二区三区| 国产视频1区2区3区| 18禁男女啪啪啪无遮挡| 天天爱天天日天天爽| 日韩黄色在线观看网站上| 69视频在线精品国自产拍| 老司机伊人99久久精品| 男插女视频大全免费| a级片特黄免费看| 精产国品一二三77777| 人妻女侠被擒受辱记| 精品美女洗澡一区二区| 天天日天天干天天日天天干天天| 亚洲欧美激情久久久| 成人资源中文在线观看| 日韩久久九九精品视频| 国产夫妻视频在线观看免费| 天天色天天射天天日天天干| 亚洲欧美另类校园春色| 不卡视频在线 欧美日韩| 91激情四射婷婷综合| 日韩欧美国产一区二区在线观看| 天天干夜夜操夜夜骑| 天天干天天弄天天日| 日本一区二区三区调教性奴视频| 中出小骚货在线观看| 高潮喷水一区二区三区| 亚洲欧美激情国产综合久久久| 乱子伦国产一区二区三区| 日本小视频一区二区| 人妻系列在线免费视频| 天堂一区二区三区在线等| 亚洲国产日韩a在线欧美| 天天日天天玩天天摸| 两个奶被揉得又硬又翘怎么回事| 美女一区二区四区六区八区| 亚洲国产日韩精品在线| 可以免费观看日韩av| 神马不卡视频在线视频| 高潮喷水一区二区三区| av在线中文字幕在线| 日韩久久九九精品视频| 污视频在线观看地址| 国产午夜羞羞一区二区三区| 天天操天天日天天碰| 可在线免费观看av| 黄色av日韩在线观看| 久久国产精品久精国产爱| 天天干天天色综合久久| 亚洲欧洲一区二区三区在线| 18福利视频在线观看| 亚洲欧美成人午夜一区二区| 青娱乐免费最新视频| 午夜国产精品免费视频| 中文字幕国产一区在线视频| 亚洲美女a级黄色在线播放| 人人妻人人狠人人爽| 久久国产半精品99精品国产| 国产精品igao为爱寻找激情| 97成人老师在线视频| 亚洲日本欧美韩国另类综合| 欧美黄色性视频网站| 4438全国成人免费视频| 亚洲 自拍 激情 另类| 2018中文字字幕人妻| 天天干天天色综合久久| 欧美亚洲愉拍一区二区三区| 欧美日韩高清片在线观看| 国产精品久久久久久成人久| 亚洲 偷拍 自拍 欧美| 亚洲成人激情在线综合| 69xx精品久久久久| 好看的日本中文字幕在线观看二区| 在线 激情 亚洲 视频| 超peng视频在线免费播放97| 十八禁黄色免费污污污亚洲| 婷婷综合缴情亚洲五月伊人| 可以直接看av网站| www,日韩av,com| 亚洲欧美另类校园春色| 日本成年视频在线免费观看| 亚洲av毛片在在线播放| 国产视频成人自拍蝌蚪视频 | 色欲AV蜜桃一区二区三| 68福利精品在线视频| 最新中文字幕久久久久| 大秀成年人国产精品视频| 欧美日韩高清片在线观看| 视频在线 一区二区| 又粗又长又硬又黄又爽| 最新国产精品综合网高清| 国产午夜在线播放视频| 妈妈的朋友2中文字幕在线| 快进来插我的逼嗯啊视频| 五月天色婷婷狠狠爱| av日韩视频在线观看| 夜夜操夜夜爱夜夜摸| 一区二区三区观看在线| 激情九月天在线视频| 欧美视频免费观看777| 欧美成人久久久桃色aa| 欧美大鸡吧男操女啊啊啊视频| 美女露阴道让男人捅| 激情久久在线免费观看视频| 天天综合久久无人区| 高潮喷水一区二区三区| 国产视频成人自拍蝌蚪视频| 亚洲国产日韩欧美一区二区三区,| 中文字幕 一区二区在线观看| 老熟女xxxⅹhd老熟女性| 狠狠干狠狠操免费视频| 国产女人18毛片水真多精选| 日本少妇三级交换做爰做| 亚洲精品色图1234| 老熟妇一区二区三区v∧88| 欧美成人短视频在线播放| 日本一道中文字幕99| aa福利影视在线观看| aaaa级少妇高潮在线观看| 荣立三等功退休有什么待遇| 欧美成人短视频在线播放| 国产一区两区三区福利小视频| av无限看熟女人妻另类av| 猫咪亚洲中文在线中文字幕| 亚洲a级视频在线播放| 伊人久久综合国产精品| 50熟妇一区二区三区| 天天插天天透天天爽| 国产av高清二区三区| 精品一区二区三区免费毛片W| 蜜乳视频一区二区三区| 视频自拍偷拍视频自拍| 日韩人妻一区二区三区在线观看| 国产成人av在线你懂得| 一区二区三区免费版在线| 福利美女视频在线观看| 美女网站视频久久精品| 夜色福利视频免费观看| 91色乱一区二区三区| 激情九月天在线视频| 精品国产污污污免费入口| 日本小视频一区二区| 色就色综合偷拍区欧美在线| 4438全国成人免费视频| 国内销魂老女人老泬| 天天想要天天操天天干| 欧美色视频网址大全| 5d蜜桃臀女无痕裸感| 182tv精品免费在线观看| 亚洲国产美女主播在线观看| 蜜乳视频一区二区三区| 久久精品四虎夜夜拍拍拍| 日本亚洲精品视频在线观看| 免费的啪啪视频软件| 五月婷婷伊人久久中文字幕| 伊人情人成综合视频| 熟妇高潮久久久久久久| 日本在线免费观看国产精品| 中文字幕人妻一区色偷偷久久| 91精品视频在线观看视频| xxxx69在线观看视频| 亚洲欧美国产一本综合首页| 黄版视频在线免费观看| 久久人妻人人草人人爽| 女人高潮潮呻吟喷水网站| 欧美精品999不卡| 18岁禁一二三区免费体验| 蜜桃臀av在线一区二区| 少妇熟女天堂网av| 国产美女主播av在线| av在线中文字幕在线| 日本午夜福利免费在线播放| 免费看日韩黄视频在线观看| 中出小骚货在线观看| 亚洲色大WWW永久网站| 国产精品亚洲精品亚洲| 一区二区三区免费版在线| 国模伊人久久精品一区二区三区| 99热在线只有的精品| 亚洲综合色一区二区三区| 婷婷色综合五月天视频| 欧美成人少妇人妻精品| 特级aaaaa黄色片| 成人资源中文在线观看| 黄色网络中文字幕日本| 亚洲少妇色小说综合| 中文字幕日本一二三区| 日本久久久久久黄色| 97视频538在线观看| 99久久国语露脸国产精品| 青青操久久综合激情| 精品国产无乱码一区二区三区 | 91久久久精品成人国产| 欧美国产精品久久久免费| 高潮喷水在线视频观看| 欧美视频免费观看777| 天天色天天射天天日天天干| 午夜一区二区三区视频在线观看| 美女网站视频久久精品| julia人妻av一区二区三区| 国产精品免费看一区二区三区| 一区二区在线观看视频网站| 亚洲资源在线免费观看| 2020精品视频在线| 日本熟妇乱妇熟色视频| 91偷拍被偷拍在线播放| xxoo福利视频导航| 精品久久久久久久久久久久久| 911美女片黄在线观看| 自拍偷拍亚洲综合第一页| 两个奶被揉得又硬又翘怎么回事| 天天插天天干天天狠| 欧美在线观看视频欧美| 国产美女高潮精品视频| 日本不卡视频一二三区| 婷婷色综合五月天视频| 丰满少妇_区二区三区| 天天色天天射天天日天天干| 欧美情色av在线观看| 国产大桥未久一区二区| 久久99嫩草99久久精品| 亚洲欧洲无码一区2区无码| 亚洲成人av在线一区二区| 狠狠操深爱婷婷综合一区| 手机视频在线观看一区| 又爽又粗又猛又色又黄视频| 麻豆国产精品777777在| 欧美成人性生活视频播放| 国色天香一二三期区别大象| 欧美一级特黄大片在线| 久久av色噜噜ai换脸| 熟女人妻少妇一区二区| 大乳丰满人妻中文字幕韩国hd| 精品国产污污污免费入口| 精品国模一区二区三区欧美| 91久久久久久最新网站| 夜夜爽夜夜操夜夜爱| 中文字幕福利视频第四页| 黄片视频免费观看视频| 精品人妻人人做人人爽| 一区二区三区不卡免费视频网站| 最新国产精品综合网高清| julia人妻av一区二区三区| 亚洲女人自熨在线视频| 久操资源在线免费播放| 日本不卡视频一二三区| 人妻中文字幕亚洲在线 | 七色福利视频在线观看| 不卡视频在线 欧美日韩| jizzjizz国产精品传媒| 亚洲最大的自拍偷拍网| 精品不卡一区二区三区| 丰满少妇_区二区三区| 欧美老熟妇xxoo老妇| 国产精品中文字幕丝袜| 亚洲熟妇在线视频观看| 最新激情中文字幕视频| 欧美日本亚欧在线观看| 夜夜躁婷婷av蜜桃妖| 开心五月综合激情婷婷| 成人av中文字幕在线看| 日韩男女视频网站在线观看| 午夜精品视频免费观看| 大鸡扒操大逼大片免费关看| 欧美黑人性猛交小矮人| 亚洲综合在线视频在线播放| 久操资源在线免费播放| 68福利精品在线视频| 视频自拍偷拍视频自拍| 日韩激情文学在线视频| 日本少妇三级交换做爰做| 91精品国产91久久久久久密臀| 99久久碰碰人妻国产| 国产精品网站亚洲发布| 真人一进一出抽搐大尺度视频| 亚洲一区二区偷拍女厕所| 汤姆提醒30秒中转进站口| 国产av剧变态维修工虐杀美女| 丰满少妇人妻一区二区三区蜜桃| 美女扒开逼逼给你看| 午夜精品久久秘?18免费观看| 日本人妻少妇xxxxxxx| 亚洲乱熟女一区二区三区影片| jiee日本美女视频网站| 美利坚合众国av天堂| 亚洲日本欧美韩国另类综合 | 美国伦理片午夜理论片| 美女妩媚午夜诱惑网站| 午夜五十路久久福利| 亚洲成人三级黄色片| 免费在线观看黄色小网站| 黑人爆操女人免费视频| 午夜福利在线不卡视频| 九九热视频1这里只有精品| 在线看的免费网站黄| 绿巨人浩克在线视频观看| 天天天天天天天天干夜夜| 在线观看网站伊人网| 婷婷色九月综合激情丁香| 一区二区三区 国产日韩欧美| 国产自拍偷拍视频在线免费观看| 久久午夜免费鲁丝片| caopeng97在线观看视频| 东京热日韩av在线| 天天躁狠狠躁狠狠躁性色| 亚洲全国精品女人久久久| 91 精品视频在线看| 成人免费电影二区三区| 亚洲熟女一区二区三区250p| 大成色亚洲一二三区| 国产福利小视频在线观看网站 | 精品不卡一区二区三区| 99热这里只有精品免费播放| 一区二区三区资源视频| 成人午夜av电影网| 日本五六十路熟女视频| 4438x亚洲最大的成人| 亚洲一区二区精品在线播放| 一区二区在线观看视频观看| 午夜在线观看一级毛| 久久久久久久岛国免费观看 | 亚洲熟妇在线视频观看| 人妻激情综合久久久久蜜桃 | 免费成人av麻豆| 天天操,天天射,天天爽| 中文字幕观看中文字幕免费 | 精品国产无乱码一区二区三区| 黄色片黄色片黄色片黄色片黄色| 97香蕉久久国产超碰| 中文字幕欧美一区二区视频| 一区二区三区高清视频3| 日本不卡 中文字幕| 九九六视频,这里只有精品| 久久久久国产精品二区| 最新日韩中文字幕啪啪啪| 青青国产95免看视频| 欧美日韩精品aaa| 亚洲一区二区在线激情| 偷拍欧美日韩另类图片| 国产激情免费在线视频| 国内自拍第一区二区三区| 正在播放麻豆精品一区二区| 天天操天天日天天插天天舔| 日韩少妇免费在线播放| 午夜国产免费视频亚洲| 中文字幕丰满子伦无码专区| 日本电影一级人妻在线播放四区| caopeng97在线观看视频| 2020精品视频在线| 妈妈的朋友中字在线免费观看| 国产一区两区三区福利小视频| 亚洲熟女人妻自拍在线视频 | 99精品久久精品一区二区| 大尺度av毛片在线网址| 亚欧洲乱码视频一二三区| 99久久国产精品免费热| 伦理在线观看未删减中文字幕| 亚洲在线观看中文字幕av| 熟女人妻aⅴ一区二区三| 伊人网在线观看 视频一区| 少妇精品视频一区二区免费看| 中文字幕人妻一区色偷偷久久| 欧美黑人性猛交小矮人|