91视频下载软件-91视频亚-91视频亚洲-91视频亚洲电影-91视频亚洲精品-91视频一区

Technical Articles

技術文章

當前位置:首頁  >  技術文章  >  A study of adsorption of dodecylamine on quartz surface using quartz crystal

A study of adsorption of dodecylamine on quartz surface using quartz crystal

更新時間:2012-10-18      點擊次數:325

This article appeared in a journal published by Elsevier. The attached
copy is furnished to the author for internal non-commercial research
and education use, including for instruction at the authors institution
and sharing with colleagues.
Other uses, including reproduction and distribution, or selling or
licensing copies, or posting to personal, institutional or third party
websites are prohibited.
In most cases authors are permitted to post their version of the
article (e.g. in Word or Tex form) to their personal website or
institutional repository. Authors requiring further information
regarding Elsevier’s archiving and manuscript policies are
encouraged to visit:
http://www.elsevier.com/copyright
Author's personal copy
Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83
Contents lists available at ScienceDirect
Colloids and Surfaces A: Physicochemical and
Engineering Aspects
journal homepage: www.elsevier.com/locate/colsurfa
A study of adsorption of dodecylamine on quartz surface using quartz crystal
microbalance with dissipation
J. Koua,b, D. Taoa,∗, G. Xua
a Department of Mining Engineering, University of Kentucky, Lexington, KY 40506, USA
b School of Civil and Environment Engineering, University of Science and Technology Beijing, 30 Xueyuan Road, Haidian District, Beijing 100083, PR China
a r t i c l e i n f o
Article history:
Received 23 April 2010
Received in revised form 13 July 2010
Accepted 15 July 2010
Available online 23 July 2010
Keywords:
Adsorption
Quartz
QCM-D
Dodecylamine hydrochloride
Zeta-potential
FTIR
a b s t r a c t
In this study the adsorption characteristics of dodecylamine hydrochloride (DACL) on quartz surface have
been investigated using a high sensitivity surface characterization technique referred to as quartz crystal
microbalance with dissipation (QCM-D) technique in conjunction with zeta-potential and FTIR analyses.
The experimental results have demonstrated the versatility and accuracy of the QCM-D for surface
adsorption characterization and for the first time revealed the changes in structure and orientation of
the amine adsorption film on quartz surface during the adsorption process by real-time measurements
of frequency and dissipation shifts with quartz coated sensor. Five distinct adsorption behaviors were
identified from the different slopes of D−f plots at different concentrations of DACL at pH 6 and 9.5.
The physisorption, coadsorption of dodecylammonium and dodecylamine, and surface precipitation of
neutral amine molecules with variable conformation and orientation were revealed on the quartz surface
by FTIR and QCM-D. Physisorption of ammonium ion and coadsorption of dodecylammonium and dodecylamine
dominated at concentrations <0.11mMatpH6 and 9.5 forming a rigid and thin adsorption layer.
A compaction stage was present at pH 9.5 at concentrations lower than 1.13mM. Surface precipitation
of neutral molecules dominated at higher concentrations at pH 6 and 9.5 to form a thick but dissipated
adsorption layer. The adsorption density was calculated with Sauerbrey equation and Voigt model and
the results indicated the existence of a critical concentration of 0.45mM at pH 6 and 1.13mM at pH 9.5
which led to a significant increase in adsorption density and a structural change in adsorption layer.
© 2010 Elsevier B.V. All rights reserved.
1. Introduction
Long-chain primary alkylammonium salts are employed in the
reverse flotation of silicates as well as in the flotation of quartz,
dolomite and calcite from phosphate minerals [1]. Since theamount
of surfactant adsorbed, the orientation of the adsorbed molecules,
and the nature of the adsorbed layer have amajor influence on the
hydrophobicity of mineral surface, a good understanding of adsorption
mechanisms of long-chain primary alkylammonium salt on
quartz surface and properties of the adsorbed layer is crucial for a
variety of industrial applications [2].
The mechanism of amine-silicate interaction has been studied
extensively by indirect methods and ex situ measurements such as
flotation recovery response [3], adsorption experiments [4,5], AFM
[6], contact angle [7], zeta-potential [8,9] and FTIR spectroscopy
[10–13] in the last several decades. Many concepts have been proposed
based on these studies.
∗ Corresponding author. .: +1 859 257 2953; +1 859 323 1962.
address: dtao@engr.uky.edu (D. Tao).
Bijsterbosch [14] investigated the mechanism of cationic surfactant
adsorption on silica surfaces by electrophoresis experiments.
The results supported the theory of monolayer formation at low
surfactant concentrations due to interaction of opposite charges on
the silica surface and the surfactant ions. The formation of bilayers
at high surfactant concentration on quartz surface was demonstrated
by Menezes et al. [15] using contact angle and zeta-potential
measurements, which indicated that the decrease in contact angles
at concentration above the CMC can be attributed to the formation
of bilayers.
The adsorption mechanism of cationic surfactant dodecylpyridinium
chloride on quartz surface was also delineated through
measurements of adsorption isotherms, zeta-potentials, suspension
stability, contact angles, induction times, and flotation
response by Fuerstenau and Jia [9]. The results showed that
dodecylpyridinium adsorption on quartz occurred in four distinct
regions as the concentration of surfactant was increased. The
adsorption was controlled primarily by electrostatic interactions at
low concentrations, but at higher concentrations, adsorbed surfactant
ions began to associate at the interface forming hemimicelles.
At the CMC, the adsorbed surfactant ions existed as a bilayer or
its equivalent, with the charged head groups oriented towards
0927-7757/$ – see front matter © 2010 Elsevier B.V. All rights reserved.
doi:10.1016/j.colsurfa.2010.07.017
Author's personal copy
76 J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83
the solution phase. Churaev et al. [16] developed two models of
Langmuir-type adsorption of cationic surfactant on quartz surfaces,
including the adsorption of cations on both charged and neutral
sites which formed the monolayer and the consecutive adsorption
that first occurred on charged sites only, followed by adsorption
on hydrophobic tails of pre-adsorbed ions, forming two layers or
hemimicelles.
Novich and Ring [3] investigated the flotation behavior and
adsorption mechanism for the alkylamine-quartz flotation system.
They demonstrated the existence of three adsorption regions in
the adsorption isotherm. In the lower concentration range single
ion adsorbed on the quartz surface by electrostatic effect and the
adsorption density increased linearly with equilibrium concentration.
At the equilibrium concentration of approximay 0.1mM
(neutral pH) the surface micelle formed with the flattened adsorption
curves, which then abruptly steepened at an equilibrium
concentration closed to the CMC and continued to rise linearly with
the increase of concentration. Multilayer adsorption occurred at
high surfactant concentrations. These results are consistent with
G–F model [7,17] which also postulated the formation of surface
micelles and hemimicelle prior to monolayer coverage at HMC
(hemimicelle concentration).
The adsorption mechanisms of long-chain alkylamines and their
acetate salts on quartz and mixed cationic and anionic collectors
on feldspar were investigated by Vidyadhar and Rao [18,19] using
Hallimond flotation, zeta-potential, FTIR, and X-ray photoelectron
spectroscopy (XPS) at neutral and acidic pH. In differentiation
to the electrostatic adsorption theory of G–F model, their FTIR
results exhibited the strong hydrogen bonds between the amine
cations and surface silanol groups. The acetate counterions were
found to influence the amine adsorption. The presence of neutral
amine molecules together with protonated ammonium ions was
also revealed by the XPS spectra above the critical hemimicelle
concentration (CHC). Vidyadhar and Rao [18] also proved spectroscopically
that long-chain alcohols were coadsorbed along with
amine cations, which led to formation of a closely packed surface
layer, as compared to the case of adsorption of pure amine alone at
the same concentration.
The correlation between solid surface free energy and flotation
activity was verified by Chibowski and Holysz [20] with quartz conditioned
with DACL dissolved in methanol. They concluded that
at 0.25 statistical monolayer of DACL the sample of quartz floats
in 90%, which can be interpreted by the drastic reduction in the
polar component of quartz surface free energy, and the maximum
flotation activity appears at one statistical monolayer (20 ?A2 was
assumed for the molecule) of the amine.
The hydrophobic attractive forces between mica surfaces in
dodecylammonium chloride solution were reported by Yoon and
Ravishankar [21]. At concentration of 0.003mM and pH 5.7 where
the surfactant is in ionized form (dodecylammonium), only shortrange
hydrophobic forces were observed, which can be attributed
to the difficulty in forming close-packed monolayers on the mica
surface. At pH 9.5 where neutral molecule of dodecylamine was
formed as a result of hydrolysis long-range hydrophobic forces
were observed, which may be attributed to the coadsorption of
both dodecylammonium and dodecylamine, increasing the packing
density of hydrocarbon chains on the surface. They also concluded
that the pH at which the long-range hydrophobic force appeared
corresponded to the maximum quartz flotation.
McNamee et al. [22] investigated the adsorption of the cationic
surfactant octadecyl trimethylammoniumchloride (C18TAC) at low
concentrations to negatively charged silica surfaces in water using
atomic force microscopy (AFM). From the AFM images of silica surface,
C18TAC was seen to form islands of bilayers and patches on
mica at a concentration of 0.03mM, which was identified as partial
surfactant bilayers or hemimicelles. The electrostatic attraction
between the bilayer islands of surfactant and bare areas of the
substrate was observed at low surfactant concentrations.
Most of these concepts about the adsorption of amines on
silicate minerals were based on ex situ studies and indirect methods
such as measurement of contact angle, zeta-potential, surface
forces, and recovery response, which unfortunay cannot monitor
the formation process of the adsorbed layer. The objective of
the present study was to perform in situ investigation of the behavior
of the adsorbed layer of dodecylamine on quartz surface by
use of the quartz crystal microbalance with dissipation technique
(QCM-D), which is capable of providing real-time information on
the adsorption layer density, strength, and structural change under
various process conditions. To our knowledge, this technique has
not been employed previously for studying the adsorption of amine
on quartz.
Quartz crystal microbalance (QCM) measures the mass change
per unit area by monitoring the change in frequency of a quartz
crystal resonator. It offers an opportunity to study molecular interactions
and conformation changes of adsorbed layer on many
different types of surfaces with acute sensitivity [23]. QCM-D is the
second generation of QCM, which cannot only determine the mass
of surface bound layers, but also simultaneously give information
about their structural (viscoelastic) properties based on the data of
dissipation factor shift. A brief discussion of this technique is provided
later in this paper.QCM-Dhas recently been used to study the
viscoelastic properties of protein adsorbed on biosensors surface
[24] and the adsorption behavior of surfactants and polymers from
aqueous solutions [25–28]. In this work, in situ adsorption behavior
of dodecylammonium onto the quartz coated sensor at various
concentrations and pH was studied for the first time using QCMD.
The obtained data was fitted with Voigt model to get physical
properties and mechanical properties of the adsorbed layer.
2. Materials and methods
2.1. Materials
Dodecylamine hydrochloride (CH3(CH2)11NH2·HCl) and anhydrous
ethanol (C2H5OH) with 99% and 99.5% purity, respectively
were acquired from Acros Organics. Sodium hydroxide, which was
purchased from Fisher Scientific, contained 99.8% NaOH. Deionized
water was used throughout the experiments.
The quartz coated sensor used in QCM-D analysis was an ATcut
quartz disc (14mm in diameter and 0.3mm in thickness with
an active sensor crystal area of 0.2cm2) with 5 nm Cr, 100nm Au,
50nm Ti and 50nm SiO2 sputter-coated onto the crystal surface
successively. The sensors and Q-sense E4 system were supplied by
Q-Sense Co.
Zeta-potential and FTIR tests were conducted with quartz powder
(50% <16.93m, 90% <85.90m) made of pulverized pure
quartz crystals purchased from Ward’s Natural Science. XRD analyses
with the sample did not show any impurity.
2.2. Quartz crystal microbalance with dissipation
All QCM-D measurements in this study were conducted at 25? C
(±0.02 ?C). The stock solution was prepared by dissolving dodecylamine
hydrochloride (DACL) in deionized water in the presence of
ethanol. The solutionpHwas adjusted to 6.00±0.05 and 9.50±0.05
using 1% NaOH and 1% HCl solution. To ensure dissolution and
degassing, the solutions were left in an ultrasonic bath for 5–10 min.
For each experiment, the data generated with the solvent only
(deionized water) was accepted as baseline when it became stable.
Following this, the DACL solutions were injected into the flow
module by a chemical feeding pump capable of precise control of
Author's personal copy
J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83 77
flow rate. The flow rate in the experiment was kept at 0.50 mL/min.
Prior to each test, the flow modules, crystals and tubings were
cleaned with 2% Hellmanex II solution (an alkaline cleaning liquid
from Hellma GmbH & Co. KG, Germany) for 0.5–1 h and deionized
water for 0.5 h. Then the crystal was removed and rinsed with more
water or ethanol. After drying with nitrogen gas, the crystal was
remounted into the module and deionized water was injected into
the system.
Software QTools 3.0 was used for data modeling and analysis.
For a rigid, thin, uniform film the change in dissipation factor
D<1×10−6 for 10 Hz frequency change [29] and the Sauerbrey
equation (Eq. (1)) was used for the mass calculation:
m = −qtqf
f0n
= −qvqf
2f 2
0 n
= −Cf
n
(1)
where the constant C has a value of 17.8 ngcm−2 Hz−1 and n is
the harmonic number (when n=1, f0 = 5MHz). If the adsorbed film
is “soft” (viscoelastic), it will not fully couple to the oscillation of
the crystal [30], causing energy dissipation of the system. The dissipation
factor D is proportional to the power dissipation in the
oscillatory system (Eq. (2)) and can give valuable information about
the rigidity of the adsorbed film [30]:
D =
Edissipated
2Estored
(2)
where Edissipated is the energy dissipated during one oscillation and
Estored is the energy stored in the oscillating system [30]. Therefore,
the changes in the property of adsorbed layer will result in the
measured change in D. During the tests, QCM-D simultaneously
measures the changes in resonance frequency f and dissipation
D as a result of adsorption on a crystal surface. The quantitative
information on the adsorbed layer on the mineral surface can be
obtained by data modeling.
When the adsorption caused greater shift in D value, i.e.,
D>1×10−6, as a result of the viscous and soft layer, Voigt modeling
(Eqs. (3) and (4)) [28] was used:
f ≈ − 1
20h0
??
?
3
?3
+

j=1,2

hjjω−2hj

3
?3
2 jω2
2j
+ 2j
ω2
 ??
? (3)
D ≈ 1
2f0h0
??
?
3
?3
+

j=1,2

2hj

3
?3
2 jω2
2j
+ 2j
ω2
 ??
? (4)
According to Voigt model for viscous adsorbed layer, f and D
depend on the density (), thickness (h), elastic shear modulus ()
and shear viscosity () of adsorption layer (j: number of adsorbed
layer). The Sauerbrey equation and Voigt model are the theoretical
basis for data modeling using software QTools 3.0 (Q-Sense Co.).
2.3. Zeta-potential measurements
The zeta-potential measurements were made with Zeta-plus
analyzer of Brook Haven Instruments Corporation. All experiments
were conducted with DACL solutions with 1.0mMKCl under ambient
conditions. 1.0 g quartz powder was conditioned in 20mL stock
solution with amagnetic stirrer for 1 h. The mineral suspension was
filtered using Whatman filter paper (pore size 25m) and then
poured into the rectangular cell for zeta-potential measurements.
2.4. FTIR analysis
The infrared transmission spectra were recorded on a Thermo
Nicolet Nexus 470 FTIR spectrometer. The quartz powders were
conditioned with 20mL DACL solution with ethanol at different
Fig. 1. Real-time experimental data of frequency shifts (a) and dissipation shifts (b)
for the third overtone (15 MHz) of QCM-D resonator for different concentrations of
dodecylamine hydrochloride adsorption on quartz surface atpH6. Arrow 1 indicates
the injection of dodecylamine hydrochloride solution.
pH’s and different concentrations while being agitated with a magnetic
stirrer for 0.5 h to make 1% suspension. The suspension was
then filtered with Whatman filter paper (pore size 25m) and the
solids were air-dried overnight at the room temperature. The samples
were prepared by dispersing 0.025 g air-dried powder in 5 g
KBr followed by pressing into a transparent tablet for scanning.
The untreated (initial) quartz powder was used as reference. Each
spectrum is the average of 250 scans.
3. Results
3.1. Adsorption behavior of DACL on quartz surface
Fig. 1 displays the real-time experimental data of the frequency
shiftf (Fig. 1(a)) and dissipation shiftD (Fig. 1(b)) from the third
overtone (15 MHz) associated with DACL adsorption onto quartz
surface at different concentrations at pH 6 (the same trend was
observed at 0.11 and 6.76mM as at 0.07 and 2.25mM, respectively,
and thus not shown in Fig. 1). Immediay after the injection
of DACL at arrow 1, a gradual decrease in f (adsorption density
increase) and a slight increase in D appeared with 0.07 and
0.11mM DACL. At the steady state D was <0.3×10−6 and f
was −2.5 Hz at 0.07 and 0.11mM, which indicated the formation of
thin and rigid adsorption film. After the injection of 0.45mMDACL,
f decreased gradually until it reached the steady state at about
−5HzwhileD remained <0.4×10−6, suggesting the formation of
a thicker adsorption layer than at 0.11 and 0.07mM DACL.
At concentrations of 1.13, 2.25 and 6.76mM, the injection of
DACL caused an initial rapid decrease in f and a significant
Author's personal copy
78 J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83
Fig. 2. Real-time experimental data of frequency shifts (a) and dissipation shifts (b)
for the third overtone (15 MHz) of QCM-D resonator for different concentrations
of dodecylamine hydrochloride adsorption on quartz surface at pH 9.5. Arrow 1
indicates the injection of dodecylamine hydrochloride solution.
increase inDfollowed by a gradual decrease inf and an increase
inD which became stable after 1.7 h. At the steady state,f of−9,
−12 and −17 Hz of was achieved with 1.13, 2.25 and 6.76mMDACL
solution, respectively and theD was higher than 2×10−6 at these
concentrations, revealing the formation of a thick but less rigid
adsorption layer with adsorbed molecules not densely compacted.
Fig. 2 shows the f and D as a function of time when the
quartz surface was exposed to different dosages of DACL at pH 9.5
(the same trend was observed at 0.07mMas at 0.11mMand therefore
is not shown in Fig. 2). The first observation was the gradual
decrease in f and increase in D after the injection of ≤0.45mM
solutions (Fig. 2, arrow 1). The steady states of f and D were
reached with≤0.45mMpH9.5 solutions 1.2 hafter the solutionwas
injected, or 0.5 h earlier than with the same concentration at pH 6.
The value of D smaller than 1×10−6 indicates that the adsorbed
layer was non-dissipative. The stable D values were almost the
same as those observed with the adsorption at the same concentration
of ≤0.45mM at pH 6 but the f was much lower at pH 9.5,
which suggests an increase in adsorption density with increasing
solution pH value.
After the injection of 1.13, 2.25 and 6.76mM DACL solutions
(Fig. 2a, arrow 1), a sharp decrease in f was observed to a steady
state value of−13,−16 and−24 Hz, respectively, which were much
lower and appeared much earlier than at pH 6. It should be noted
that the steady state value of 0.8×10−6 for D at 1.13mM was
much lower at pH 9.5 than 2.6×10−6 produced at pH 6, suggesting
that a thicker and more rigid adsorption layer was formed in
1.13mM DACL solutions at pH 9.5 than at pH 6. The D’s at 2.25
and 6.76mMwere also lower at pH 9.5 than at pH 6, but still higher
than 1×10−6.
Fig. 3. Real-time experimental data of frequency and dissipation shifts for the third
overtone (15 MHz) of QCM-D resonator for DACL adsorption and rinsing with deionized
water from quartz surface. The straight lines are the f and D of 2.25mM
DACL at pH 6; the dash lines are the f and D of 1.13mM DACL at pH 6; the
dash-dotted curves are the f and D of 6.76mM DACL at pH 9.5; the dotted
curves are the f and D of 0.45mM DACL at pH 9.5. Arrows 1 and 2 indicate
the injection of dodecylamine hydrochloride solution and reagent-free water,
respectively.
Fig. 3 shows the f and D as a function of time when the
quartz surface was exposed to 1.13 and 2.25mM DACL solutions
at pH 6 and 0.45 and 6.76mM DACL solutions at pH 9.5 at
arrow 1 followed by rinsing with water of the same pH value but
without DACL at arrow 2. It can be observed that after the injection
of water at arrow 2, f and D at all these concentrations
changed significantly. With a 0.45mM solution at pH 9.5, the D
and f were immediay back to zero after rinsing with water,
implying that the desorption occurred readily at lower concentrations
at pH 9.5. At concentrations of 1.13 and 2.25mM at pH
6 and of 6.76mM at pH 9.5, D changed to about 0–0.5×10−6
with f increased to about −2 Hz, which corresponded to
0.5×10−10 mol/cm2 adsorption density. This adsorption density
was much lower than the statistical monolayer, which was
6.64×10−10 mol/cm2 on the crystal sensor (the active sensor crystal
area was 0.2cm2 and 25 ?A2 was assumed for the molecule) [3].
It can be concluded that even at higher pH or higher concentration,
almost all adsorbed DACL can be washed off from the quartz
surface.
3.2. Adsorption density of DACL
The measured data of f and D from different overtone
(from 3rd to 13th) was fitted using the Sauerbrey equation for
D<10−6 or Voigt model for D>10−6. The calculated adsorption
density at concentrations from 0.07 to 6.76mM is shown
in Fig. 4. It can be seen that the adsorption density increased
with increasing concentration at both pH values. In pH 6 solutions
the adsorption density was lower than 5×10−10 mol/cm2 at
concentrations ranging from 0.07 to 0.45mM. The adsorption density
jumped to 19.4×10−10 mol/cm2 at 1.13mM and continued to
increase to about 56.7×10−10 mol/cm2 with increasing the concentration
to 6.76mM. A thicker adsorption layer appeared at all
concentrations at pH 9.5 than at pH 6. For example, the adsorption
density was 19.4×10−10 mol/cm2 and 47.3×10−10 mol/cm2
at pH 6 and 9.5, respectively in 1.13mM solutions. The dissipation
shift D was only 0.8×10−6 with 1.13mM DACL at pH 9.5,
which was <2.5×10−6 observed at pH 6. It can be concluded that
the adsorption layer was thicker and less dissipated at higher
pH.
Author's personal copy
J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83 79
Fig. 4. Adsorption density for different concentrations of dodecylamine hydrochloride
on quartz surface at pH 6 and 9.5.
3.3. Kinetics of DACL adsorption
To characterize the adsorption behavior of DACL on quartz, the
measured data of f and D from third overtone (15 MHz) in
Figs. 1 and 2 were correlated in D−f plots in Fig. 5. D−f
plots serve as the fingerprint of an adsorption process and provide
the information on the energy dissipation per unit mass added
to the crystal sensor [29]. The slope of the plots is defined as K
(K =D/f, absolute value) which is indicative of kinetic and structural
alternation during adsorption process [29]. Small value of K
indicates the formation of a compact and rigid layer, while a high
value indicates a soft and dissipated layer [23,29].
Fig. 5 shows the D−f plots for the adsorption of DACL on
the quartz surface at (a) 0.11mM at pH 6, (b) 2.25mM at pH 6, (c)
0.45mM at pH 9.5, and (d) 2.25mM at pH 9.5. The data in Fig. 5
(a) shows only one phase (one slope) and Fig. 5(b) and (d) shows
two distinct phases (two slopes) of kinetics whereas Fig. 5(c) shows
three phases (three slopes). According to the previous studies conducted
by Rodahl et al. [23], Hook et al. [31], and Paul et al. [29], a
morerigid and compact adsorption mass is expected to yield a small
K value and a soft and dissipated layer is associated with a higher
K value. One slope only indicates a direct adhesion on the surface.
The presence of multiple slopes suggests direct adhesion and orientational
changes associated with hydrodynamically coupled water
[29,23]. Data shown in Fig. 5 reveals that DACL adsorption took
Fig. 5. D−f plots for DACL adsorbed on the quartz surface at (a) 0.11mM and pH 6; (b) 2.25mM and pH 6; (c) 0.45mM and pH 9.5; (d) 2.25mM and pH 9.5.
Author's personal copy
80 J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83
Table 1
Slopes in D−f plots for adsorption of different concentrations of dodecylamine hydrochloride at pH 6 and 9.5 on quartz.
Concentration of DACL (mM) 0.07 0.11 0.45 1.13 2.25 6.76
pH 6 slope K (×10−6)Hz−1
0.0313 0.102
K1 = 0.2495 K1 = 0.3609 K1 = 0.4904 K1 = 0.2516
K2 = 0.0288 K2 = 0.093 K2 = 0.1409 K2 = 0.0026
pH 9.5 slope K (×10−6)Hz−1 K1 = 0.1038 K1 = 0.0438 K1 = 0.1417
K1 = 0.1549 K1 = 0.0467 K1 = 0.0387
K2 = 0.0051 K2 = 0.0027 K2 = 0.0072
K3 = 0.1035 K3 = 0.0314 K3 = 0.0672 K2 = 0.006 K2 = 0.1544 K2 = 0.4146
place with constant kinetics at 0.11mM in pH 6 solution. However,
the adsorption in pH 9.5 solutions or in pH 6 solutions with
DACL concentration higher than 0.11mM resulted in changes in
adsorption characteristics.
Table 1 shows calculated K values for adsorption processes at
various DACL concentrations and pH’s represented in Fig. 4. Four
distinct adsorption processes can be identified: (1) the adsorption
layer was thin and rigid at lower DACL concentrations (0.07
and 0.11mM) with only one K value of <0.15 at pH 6, which indicates
the direct adhesion with no kinetic changes; (2) two slopes
(K1 > 0.15 > K2) were observed with DACL concentrations of greater
than 0.45mM at pH 6 and 1.13mM at pH 9.5, which indicates the
formation of an initial flexible and water rich adsorption layer, followed
by a structural and orientational changes as evidenced by a
decrease in energy dissipation per unit adsorption mass K1 values
increased with increasing concentration of DACL at pH 6, as shown
in Table 1, which indicates that the adsorption layer became more
and more dissipated. (3) Three slopes (0.15 > K1 ≈K3 > 0.01 > K2)
were observed with lower DACL concentrations (0.07, 0.11 and
0.45mM) at pH 9.5, which shows the same slopes at the beginning
and the end of the adsorption process separated by a stage with
a much lower slope of <0.01. Small values of three slopes (<0.15)
indicate the formation of compact and rigid layers as evidenced by
D<10−6 and orientational changes during the adsorption. (4) At
higher DACL concentrations of 2.25 and 6.67mM and pH 9.5, two
slopes (K2 > 0.15 > K1) were observed, indicating the rapid adhesion
followed by conformational changes or water molecules trapped in
the adsorption layer with low structural stability. Comparing the K
values in Table 1 also shows that a smaller K1 value was observed
at pH 9.5 than at pH 6 with the same DACL concentrations, which
indicates that although the initial adsorption of DACL on quartz was
rapid at pH 9.5, as shown in Figs. 1 and 2, the adsorbed layer was
more rigid and less dissipated than at pH 6. When the adsorption
became stable, the adsorbed layer of DACL had higher adsorption
density and lower D at pH 9.5 than at pH 6, which indicated better
organized structure of DACL on the quartz surface at higher pH.
These adsorption results can be used to explain why quartz flotation
at pH 9.5 can achieve better performance than at pH 6 when
amine is used as collector. According to Mielczarski et al. [32] the
hydrophobicity of mineral is closely related to the structure of the
adsorbed surfactant on the surface. The highly packed adsorption
layer of surfactant imparts strong hydrophobicity to solid surface
whereas a poorly organized structure does not.
3.4. Zeta-potential measurement
To better understand the adsorption mechanism of DACL in
weakly acidic solutions, the relationship between the initial DACL
concentration and the zeta-potential of quartz particles at pH 6 is
shown in Fig. 6. It is evident that increasing the concentration to
about 0.36mMsharply increased the surface charge of quartz from
−110 to 0mV. Novich and Ring [3] suggested that the quartz surface
will be negatively charged if the suspension pH is above the point of
zero charge (pH 2.0) and cationic amine surfactant will absorbwhen
it is placed in solution with the negatively charged quartz surface.
According to G–F model, alkylammonium ions specifically adsorb
at negative sites on the quartz surface in the low-concentration
region, oriented with the charged head toward the surface and the
hydrocarbon tail into the solution, which accounted for the increase
in zeta-potential. The QCM-D results also showed a low adsorption
density of 1.8×10−10 mol/cm2 and only one slope of K≤0.15 on
the D−f plot at 0.07 and 0.11mM at pH 6. It can be concluded
that a rigid and well-ordered monolayer was present on the quartz
surface at a concentration lower than 0.36mM at pH 6.
A further increase in DACL concentration from 0.36 to 6.76mM
led to positively charged phase although the rate of increase in surface
charge slowed with increasing DACL concentration. It can be
postulated based on the G–F model that after the monolayer coverage,
a second layer forms with the charged amine heads toward
the solution and the hydrocarbon tails oriented toward the surface.
This is responsible for the positive zeta-potential at the beginning
of the higher concentration region and agrees well with the QCMD
results that showed two slopes (K1 > 0.15 > K2) on the D−f
plot and significantly increased adsorption density (decreased f)
at higher DACL concentrations at pH 6.
3.5. FTIR analysis
Infrared transmission spectra of quartz after adsorption in DACL
solution at different pH’s and concentrations were investigated to
confirm the adsorption mechanism of amine on the quartz surface.
Fig. 7 shows the infrared transmission spectra of dodecylamine
hydrochloride (DACL) as reference. The bands characteristic of
alkyl chains for dodecylamine are identified at 2950, 2920, and
2850cm−1, which were assigned to asymmetric stretching vibration
in theCH3 (as (CH3)) andCH2 radical (as(CH2)) and symmetric
stretching vibration in the CH2 radical (s(CH2)) [18], respectively.
The broad band at 3000cm−1 was observed, which was assigned to
as(NH3
+) and s(NH3
+) [18].
Fig. 6. Zeta-potential of quartz conditioned with DACL at constant ionic strength as
a function of concentration.
Author's personal copy
J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83 81
Fig. 7. Infrared transmission spectra of DACL.
Curves a and b in Fig. 8A show the absorption bands ranging
from 2750 to 3050cm−1 in the FTIR spectra of the species adsorbed
on quartz surface from the solutions of 4.5mMDACL at pH 6 and 9.5
and curves c and d show the same absorption bands with 0.45mM
DACL solution at pH 6 and 9.5, respectively. Curve e shows the
adsorption bands with quartz powder conditioned with ethanol
Fig. 8. Infrared transmission spectra of quartz after adsorption of DACL at different
concentrations and pHs for (A) alkyl chain region 2750–3050cm−1 and (B) (OH)
group and amine molecule region 3050–3500cm−1 (a) 4.5mM, pH 9.5; (b) 4.5mM,
pH 6; (c) 0.45mM, pH 9.5; (d) 0.45mM, pH 6; (e) quartz powder with ethanol.
Fig. 9. Infrared transmission spectra of quartz after adsorption of DACL at different
concentrations and pHs for the alkyl chain region 2750–3100cm−1 (a) 4.5mM, pH
9.5; (b) 4.5mM, pH 6; (c) 0.45mM, pH 9.5; (d) 0.45mM, pH 6; (e) quartz powder
with ethanol; (f) DACL.
solution as reference. Curves a–d show the absorption peaks at
2950, 2920, and 2850cm−1, which are undoubtedly caused by the
adsorption of alkyl chain on the quartz surface since they are consistent
with DACL spectra shown in Fig. 7. In comparison with curve
e, the intensities of these bands in Fig. 8 were found to increase with
increasing DACL concentration, which is in agreement with QCM-D
measurements and zeta-potential results.
According to Vidyadhar et al. [2], the amine molecule in the
bulk phase in neutral and protonated forms can be identified with
the appearance of the band at 3333cm−1 and the alcohol spectrum
shows a broad band around 3318cm−1 characteristic of (OH)
group. Fig. 8B (curves a–e) shows the spectra in the region from
3300 to 3350cm−1 of quartz treated with DACL solution under different
conditions and with ethanol solution. The broad band on the
spectra indicated the existence of H-bond of alcohol (OH) group
on quartz surface at all concentrations and pH’s as well as at the
higher concentration of amine.
The DACL spectrum in Fig. 7 also displays several bands in the
region of 1700–1000cm−1. The band at 1460–1480cm−1 is often
assigned to the methylene scissoring band ?(CH2) on the quartz surface.
But quartz also exhibits strong absorption in the same region
and thus it is hard to identify these bands undermono-or submonolayer
adsorption [18]. Fig. 9 shows the infrared transmission spectra
for the adsorption bands in the region from 1400 to 1800cm−1. It
was observed that the absorption of DACL gave rise to characteristic
absorption bands in the region of 1400–1700cm−1 on quartz
surface and the intensity of these absorption peaks increased with
increasing concentration.
The FTIR results did not show significant effects of pH on the
intensities of absorption peaks, which is inconsistent with flotation
performance of amine at different pH’s. In contrast the QCM-D
results clearly showed differences in adsorption behavior at different
pH’s. For example, theD−f plot shows K1 > K2 at high DACL
concentrations at pH 6 and K1 < K2 at pH 9.5. According to Sirkeci
[33] the coadsorption of ionic and molecular species is responsible
for quartz flotation with dodecylamine at high pHs. Sirkeci [33]
further claimed that pH 9.3 represents the equilibrium between
ionized and molecular amine species. Vieira and Peres [34] also
reported that the most favorable pH for quartz flotation was pH 9.0
Author's personal copy
82 J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83
for both ether monoamine and ether diamine. In good agreement
with previous studies QCM-D results have shown the significant
increase in adsorption density and adsorption layer stability as a
result of surface precipitation of molecular amine at pH 9.5 on the
quartz.
4. Discussion
The previous studies on the adsorption of dodecylamine on
quartz surface have revealed three adsorption mechanisms: (1)
single ion adsorption at lower concentration; (2) surface micelle
formation at intermediate concentration; and (3) multilayer
adsorption with molecular amine precipitation at higher concentration
[3]. The present QCM-D study provided for the first time the
real-time adsorption data and the kinetic D−f plots displayed
five distinguishable adsorption stages:
(1) At lower concentrations (0.07 and 0.11mM) at pH 6, a slight
decrease in f and a small D value of <1×10−6 with only
one slope K < 0.15 was observed, which indicates the adsorption
process formed a thin but rigid adsorption layer. The species
distribution diagram developed by Novich and Ring [3] for
0.04mM dodecylamine solution at different pH’s suggests that
cationic amine (RNH3
+) is the dominant species in pH 6 solution.
According to Vidyadhar and Rao [18], at pH from 2 to
7, where the surface potential of silicate is negative, ammonium
ions undergo physisorption and are electrostatically held
at the silicate–water interface much below the critical micelle
concentration (CMC, 10mM at pH 6), which is in agreement
with zeta-potential measurement results. Novich and Ring
[3] reported an adsorption density of 1.24×10−10 mol/cm2 at
first adsorption step, in agreement with QCM-D results which
showed 1.8×10−10 mol/cm2 adsorption density on quartz surface
at concentrations of 0.07 and 0.11mM.
(2) Three slopes (0.15 > K1 ≈K3 > 0.01 > K2) were observed at lower
DACL concentrations at pH 9.5, which revealed different
adsorption processes from at pH 6. According to the species distribution
diagram, the concentration of cationic amine starts
to decrease and a small amount of RNH2 is observed in the
0.04mM solution at pH 9.5. Yoon and Ravishankar [21] indicated
that neutral molecules of dodecylamine were coadsorbed
with dodecylammonium at low concentrations atpH9.5, which
may increase the packing density of hydrocarbon chains on the
surface and form a higher adsorption density but lower dissipative
adsorption layer. The existence of one small slope between
two higher slopes may indicate that the orientational changes
and layer compression step resulted in more rigid adsorption
layer.
(3) It should be mentioned that the QCM-D results show different
adsorption behavior at concentrations of 0.45mM at pH 6
and 1.13mM at pH 9.5. Although the D−f plot displayed
two slopes and higher adsorption density at 0.45mM and pH 6
and 1.13mM at pH 9.5, the dissipation shift D was <1×10−6.
These two adsorption conditions produced higher adsorption
density than at concentration <0.45mM, but lower dissipation
shift than at concentration >1.13mM. After the monolayer
coverage the formation of second layer with the hydrocarbon
tails oriented toward the surface and the charged amine heads
toward the solution accounted for the charge reversal of zetapotential
(0.36mM at pH 6 as shown in Fig. 6). According to
QCM-D results, the critical concentration for hemimicelle formation
increased with increasing pH from 6 to 9.5 and the
hemimicelle formed at this critical concentration was rigid and
well-ordered.
(4) At DACL concentrations higher than 0.45mM at pH 6, a significant
decrease in f happened simultaneously with a sharp
increase in D. Two slopes (K1 > 0.15 > K2) were observed on
D−f plots, suggesting the formation of a thick, dissipated
and less rigid structure followed by a compaction stage on the
quartz surface caused by surface precipitation. According to
Chernyshova et al. [8] the adsorbed film on quartz includes
both neutral and protonated amine H-bonded to the surface
silanols at higher concentration, which is in agreement with
FTIR results shown in Fig. 8 that also showed adsorption peak
intensity higher than at lower concentration and without DACL.
(5) At concentrations higher than 1.13mM and pH 9.5, a dramatic
decrease in f and a significant increase in D were observed
along with two slopes (K2 > 0.15 > K1) on D−f plots, which
suggests the formation of a rigid structure followed by a surface
precipitation that resulted in higher adsorption density. This
result is consistent with the theory proposed by Novich and
Ring [3] that there is a multilayer of adsorbates held together
by Van de Waals force associated with hydrocarbon chains.
McNamee et al. [35] also reported the Van der Waals attraction
between the cationic surfactant and negatively charged
silica surface based on AFM and zeta-potential measurements.
The poorly organized multilayer structure is reflected in high
dissipation shift D and low resonance frequency change f
observed in the present study. The DACL adsorption density in
2.25mM solution was estimated to be 47.3×10−10 mol/cm2 at
pH 9.5, which is very close to the value of 39.6×10−10 mol/cm2
Novich and Ring [3] reported as a result of multilayer adsorption.
5. Conclusions
In situ measurements of dodecylamine adsorption on quartz
surface at different concentrations and pH’s were conducted using
the QCM-D technique to investigate the adsorption density and
kinetics and the structural properties of the adsorption layer. The
following major conclusions were derived from this study:
(1) Three adsorption mechanisms, i.e., physisorption, formation
of hemimicelle and surface precipitation of neutral molecule,
were revealed by the QCM-D technique combined with Zetapotential
measurements and FTIR analysis for the adsorption
of dodecylamine on quartz surface. At low concentrations and
pH 6, physisorption of ammonium ion formed a rigid and thin
adsorption layer, resulting in a small change in both f and
D. The coadsorption of both dodecylammonium and dodecylamine
were observed in pH 9.5 solutions, which formed a
thicker but more rigid layer than at pH 6. Surface hemimicelle
with a higher density and a lower dissipation was formed
at intermediate concentrations, and the critical concentration
for hemimicelle formation increased from 0.45 to 1.13mM
with increasing pH from 6 to 9.5. Surface precipitation of neutral
amine molecules dominated at higher concentrations and
caused a rapid decrease in f and increase in D, which indicated
the formation of a thick but soft and dissipated adsorption
layer.
(2) Different adsorption behaviors can be readily derived from
the different slopes of D−f plots. At lower concentrations
the adsorption layer was thin and rigid with only
one small slope (K < 0.15) at pH 6 and three small slopes
(0.15 > K1 ≈K3 > 0.01 > K2) at pH 9.5. At higher concentrations
two slopes were observed with K1 > 0.15 > K2 at pH 6 and
K2 > 0.15 > K1 at pH 9.5, which indicated the formation of higher
adsorption density and more dissipated multilayer structure. A
critical concentration, i.e., 0.45mMat pH 6, and 1.13mMat pH
Author's personal copy
J. Kou et al. / Colloids and Surfaces A: Physicochem. Eng. Aspects 368 (2010) 75–83 83
9.5, was observed which indicated the occurrence of surface
hemimicelle with higher adsorption density and lower dissipation.
(3) Adsorption densities under different solution conditions were
calculated with Sauerbrey equation and Voigt model. At concentrations
≤0.45mM at pH 6 and ≤0.11mM at pH 9.5
the adsorption density was close to the theoretical monolayer
coverage of 6.64×10−10 mol/cm2; the adsorption density
increased significantly at concentrations higher than 0.45mM
at both pH’s.
(4) Zeta-potential and FTIR studies showed that the DACL adsorption
on quartz involved both physisorption as a result of the
electrostatic interaction at lower concentrations and surface
precipitation of amine molecules on quartz surface at higher
concentrations, which helped explain the different adsorption
processes observed in QCM-D studies.
References
[1] H. Sis, S. Chander, Reagents used in the flotation of phosphate ores: a critical
review, Miner. Eng. 16 (2003) 577–585.
[2] A. Vidyadhar, K.H. Rao, I.V. Chernyshova, Mechanisms of amine-feldspar interaction
in the absence and presence of alcohols studies by spectroscopic
methods, Colloids Surf. A: Physicochem. Eng. Aspects 214 (2003) 127–142.
[3] B.E. Novich, T.A. Ring, A predictive model for the alkylamine–quartz flotation
system, Langmuir 1 (1985) 701–708.
[4] H. Matsubara, T. Eguchi, H. Takumi, K. Tsuchiya, T. Takiue, M. Aratono, Surface
adsorption and aggregate formation of cationic gemini surfactant and longchain
alcohol mixtures, J. Phys. Chem. B 113 (2009) 8847–8853.
[5] Q. Zhou, P. Somasundaran, Synergistic adsorption of mixtures of cationic gemini
and nonionic sugar-based surfactant on silica, J. Colloid Interface Sci. 331 (2009)
288–294.
[6] J.H. Zhang, R.H. Yoon, J.C. Eriksson, AFM surface force measurements concluded
with silica in CnTACL solutions: effect of chain length on hydrophobic force,
Colloids Surf. A: Physicochem. Eng. Aspects 300 (2007) 335–345.
[7] D.W. Fuerstenau, H.M. Jang, On the nature of alkylsulfonate adsorption at the
rutile/water interface, Langmuir 7 (1991) 3138–3143.
[8] I.V. Chernyshova, K.H. Rao, A. Vidyadhar, Mechanism of adsorption of longchain
alkylamines on silicates. A spectroscopic study. 1. Quartz, Langmuir 16
(2000) 8071–8084.
[9] D.W. Fuerstenau, R. Jia, The adsorption of alkylpyridinium chlorides and their
effect on the interfacial behavior of quartz, Colloids surf. A: Physicochem. Eng.
Aspects 250 (2004) 223–231.
[10] I.V. Chernyshova, K.H. Rao, A. Vidyadhar, Mechanism of adsorption of longchain
alkylamines on silicates: a spectroscopic study. 2. Albite, Langmuir 17
(2001) 775–785.
[11] R.M.F. Lima, P.R.G. Brandao, A.E.C. Peres, The infrared spectra of amine collectors
used in the flotation of iron ores, Miner. Eng. 18 (2004) 267–273.
[12] I.V. Chernyshova, K.H. Rao, A new approach to the IR spectroscopic study of
molecular orientation and packing in adsorbed monolayers. Orientation and
packing of long-chain primary amines and alcohols on quartz, J. Phys. Chem. B
105 (2001) 810–820.
[13] I.V. Chernyshova, K.H. Rao, Mechanism of coadsorption of long-chain alkylamines
and alcohols on silicates: Fourier transform spectroscopy and X-ray
photoelectron spectroscopy studies, Longuire 17 (2001) 2711–2719.
[14] B.H. Bijsterbosch, Characteriztion of silica surfaces by adsorption from solution:
investigation into themechanismof adsorption of cationic surfactants, J. Colloid
Interface Sci. 47 (1973) 186–198.
[15] J.L. Menezes, J. Yan, M.M. Sharma, The mechanism of alteration of macroscopic
contact angles by the adsorption of surfactants, Colloids Surf. 38 (1989)
365–390.
[16] N.V. Churaev, I.P. Sergeeva, V.D. Sobolev, H.J. Jacobasch, P. Weidenhammer, F.J.
Schmitt, Modification of quartz surfaces using cationic surfactant solutions,
Colloids Surf. A: Physicochem. Eng. Aspects 164 (2000) 121–129.
[17] A.M. Gaudin, D.W. Fuerstenau, Quartz flotation with cationic collectors, Trans.
AIME 202 (1955) 958–962.
[18] A. Vidyadhar, K.H. Rao, Adsorption mechanism of mixed cationic/anionic collectors
in feldspar-quartz flotation system, J. Colloid Interface Sci. 305 (2007)
195–204.
[19] A. Vidyadhar, K.H. Rao, I.V. Chernyshova, K.S.E. Forssberg, Mechanisms of
amine-quartz interaction in the absence and presence of alcohols studied by
spectroscopic methods, J. Colloid Interface Sci. 256 (2002) 59–72.
[20] E. Chibowski, L. Holysz, Correlation of surface free energy changes and floatability
of quartz, J. Colloid Interface Sci. 112 (1985) 15–23.
[21] R.H. Yoon, S.A. Ravishankar, Long-range hydrophobic forces between mica surface
in alkaline dodecylammonium chloride solutions, J. Colloid Interface Sci.
179 (1996) 403–411.
[22] C.E. McNamee, H.J. Butt, K. Higashitani, I.U. Vakarelski, M. Kappl, Interaction of
cationic hydrophobic surfactants at negatively charged surfaces investigated
by atomic force microscopy, Longmuir 25 (2009) 11509–11515.
[23] M. Rodahl, F. Hook, C. Fredriksson, C.A. Keller, A. Krozer, P. Brzezinski, M.
Voinova, B. Kasemo, Simultaneous frequency and dissipation factor QCM measurements
of biomolecular adsorption and cell adhesion, Faraday Discuss. 107
(1997) 229–246.
[24] X. Turon, O.J. Rojas, R.S. Deinhammer, Enzymatic kinetics of cellulose hydrolysis:
a QCM-D study, Langmuir 24 (2008) 3880–3887.
[25] J.J.R. Stalgren, J. Eriksson, K. Boschkova, A comparative study of surfactant
adsorption on model surfaces using the quartz crystal microbalance and the
ellipsometer, J. Colloid Interface Sci. 253 (2002) 190–195.
[26] M.A. Plunkett, P.M. Claesson, M.W. Rutland, Adsorption of cationic polyelectrolyte
followed by surfactant-induced swelling, studied with a quartz crystal
microbalance, Langmuir 18 (2002) 1274–1280.
[27] B.D. Vogt, C.L. Soles, H.J. Lee, W.L. Wu, Moisture adsorption into ultrathin
hydrophilic polymer films on different substrate surface, Polymer 46 (2005)
1635–1642.
[28] M.V. Voinova, M. Rodahl, M. Jonson, B. Kasemo, Viscoelastic acoustic response
of layered polymer films at fluid–solid interfaces: continuum mechanics
approach, Phys. Scr. 59 (1999) 391–396.
[29] S. Paul, D. Paul, B. Tamara, A.K. Ray, Studies of adsorption and viscoelastic
properties of proteins onto liquid crystal phthalocyanine surface using quartz
crystal microbalance with dissipation technique, J. Phys. Chem. C 112 (2008)
11822–11830.
[30] P. Ekholm, E. Blomberg, P. Claesson, I.H. Auflem, J. Sjoblom, A. Kornfeldt, A
quartz crystal microbalance study of the adsorption of asphaltenes and resins
onto a hydrophilic surface, J. Colloid Interface Sci. 247 (2002) 342–350.
[31] F. Hook, M. Rodahl, P. Brzezinski, B. Kasemo, Energy dissipation kinetics for
protein and antibody-antigen adsorption under shear oscillation on a quartz
crystal microbalance, Langmuir 14 (1998) 729–734.
[32] J.A. Mielczarski, J.M. Cases, E. Bouquet, O. Barres, J.F. Delon, Nature and structure
of adsorption layer on apatite contacted with oleate solution. 1. Adsorption and
Fourier transform infrared reflection studies, Langmuir 9 (1993) 2370–2382.
[33] A.A. Sirkeci, Electrokinetic properties of pyrite, arsenopyrite and quartz in the
absence and presence of cationic collectors and their flotation behavior, Miner.
Eng. 13 (2000) 1037–1048.
[34] A.M. Vieira, A.E.C. Peres, The effect of amine type, pH, and size range in the
flotation of quartz, Miner. Eng. 20 (2007) 1008–1013.
[35] C.E. McNamee, M. Matsumoto, P.G. Hartley, P. Mulvaney, Y. Tsujii, M. Nakahara,
Interaction forces and zeta potentials of cationic polyelectrolyte coated silica
surfaces in water and in ethanol: effects of chain length and concentration of
perfluorinated anionic surfactants on their binding to the surface, Langmuir 17
(2001) 6220–6227.

400-833-6968
歡迎您的咨詢
我們將竭盡全力為您用心服務
關注公眾號
版權所有 ? 2025 佰奧林(上海)貿易有限公司  備案號:滬ICP備18035521號-3

TEL:18612271669

關注公眾號
性爱视频在线播放 | 色琪琪av男人的天堂 | 高清不卡伦理电影在线观看 | 亚洲一区二区观看播放 | 欧美精品黄页在线观看大全 | 亚洲人成77777A片张津瑜 | 风韵丰满熟妇啪啪区老熟熟女 | 色综合天天综合网站中国 | jizz曰本| 天天搞夜夜操 | 国产人妻人伦精品无码.麻豆 | 女性人体aa欣赏 | BT7086福利二区最新 | 日韩特黄特色大片免费视频 | 性一交一乱一伦一A片 | 91精品一区二区三区久久久久 | 亚洲精品偷拍的自拍的 | 久久九九久精品国产尤物 | 日日碰狠狠躁久久躁AV | 国产精品186在线观看在线播放 | 狠狠色丁香久久婷婷综合_中 | 不良女警免费看 | 天堂8资源在线官网资源 | 漂亮人妻被强中文字幕 | 日韩欧美一区二区三区视频 | 又爽又色又舒服A片免费 | 国产传媒在线观看视频免费观看 | 久久人妻AV一区二区软件 | 成人做爰69片免费看网站 | 亚洲一区在线观看视频 | 在线播放一区二区精品产 | 日本处888xxxx| 熟女人妻-蜜臀AV-首页 | 国产无套视频在线观看香蕉 | 有码中文 | 国产麻豆VIDEOXXXX实拍 | 四色777 | 黄色网页在线免费观看 | 69久久无码一区人妻A片 | 亚洲国色天香卡2卡3卡4 | 国产福利视频情侣视频 | 国产精品久久精品第一页不卡 | 亚洲高清在线不卡中文字幕网 | 护士97资源站共享总站 | 欧美精品做人一级爱免费 | 免费国产美女爽到喷出水来视频 | 深夜a级毛片免费视频 | 没有废话全色肉的黄文 | 欧美高清一区二区 | 久久婷婷激情综合中文字幕 | 欧美午夜一区 | 国产色情乱码久久久久一区二区 | 成人午夜A片产无码免费视频日本 | 女人18毛毛片兔费码A片 | 亚州老熟女A片AV色欲小说 | 国产1988精品A片 | 2018av在线 | 视频一区二区三区欧美日韩 | 最近中文字幕免费手机版 | 成人污污污WWW网站免费直播 | 天天爽夜爽免费精品视频 | 日本少妇裸体做爰高潮片 | 91精品手机国产在线观 | 中文有码视频 | videos欧美成人 | 黄色片之夜 | 成人美女免费网站视频 | 日本一道无马二区日本道专区 | 神马影院午夜理论二 | 无人在线观看免费高清直播视频 | 香港三级日本三级韩国三级 | 久久99精品国产99久久 | 免费aⅴ在线 | 色欲AV亚洲精品一区二区 | 毛片在线不卡 | 久久99热这里只频精品6 | 欧美精品无码久久久 | 欧美 国产日韩 综合在线 | 国内自拍小视频 | 免费视频片在线观看大片 | 欧洲精品不卡1卡2卡三卡 | 经典三级一区二区三区视频 | 少妇我被躁爽到高潮A片 | 国产99区 | 久久精品免费i 国产 | 凡人修仙传动漫53集免费 | 荡公乱妇翠翠 | 九九热国产精品视频 | 天天做天天爱天天综合网 | 欧美激情A片久久久久久 | 欧美成人精品一区二区综合A片 | 欧美日韩在线免费观看 | 国产综合有码无码中文字幕 | 日本无码精品一区二区三 | 午夜成人A片精品视频免费观看 | 泷川雅美的不雅视频 | 亚洲三级影院 | 日本少妇裸体做爰高潮片 | 色情图插插插 | 玩弄放荡人妇系列短篇下载 | 亚洲精品久久久久中文字幕二区 | 免费又粗又黄又爽又免费A片 | 狠狠久久免费视频在线 | 中文字幕一区2区 | 一个人在线视频免费观看www | 97碰在线视频 | 一级毛片美国 | 性色国产成人久久久精品一区二区 | 中文字幕一区二区三区在线观看 | 日本三级免费网站 | 国产偷窥熟妇高潮呻吟 | 91精品国产色综合久久不 | 四虎国产精品永久一区高清 | 91九色视频无限观看免费 | 国产一区二区亚洲精品 | 日本后进式猛烈xx00动态图 | 免费看污黄网站 | 成人A片动漫无码免费播放 成人A片免费看男人社区 | 99re在线精品 | 国精产品一区一区三区免费视频 | 国产福利美女福利视频免费看 | 亚洲一区二区三区四区五区六区 | 成人福利在线看 | 国产精品色无码AV在线观看 | WWW日韩AV免费高清看 | 草草视频在线播放 | 丰满少妇内射一区 | 精品国产成人亚洲午夜福利 | 免费观看亚洲视频 | 国产人妻人伦精品98 | 人人网站| 让人爽到湿的小黄书 | 插的痛的视频一卡二卡三卡 | 国产精品久久人妻互换毛片 | 最好看的最新的中文字幕 | 公车疯狂揉肉蒂高潮H失禁动态图 | 国产50岁熟妇露脸 | 玉蒲团之灯草和尚 | 欧美日韩第二页 | 久久不卡精品 | 亚洲一-卡2卡3卡乱码... | 污污内射久久一区二区欧美日韩 | 特级精品毛片免费观看 | 精品人妻无码一区二区三区婷婷 | 国产精品视频在线观看 | 国产熟睡乱子伦视频在线播放 | 少妇系列之白嫩人妻 | 国产a级一级久久毛片 | 丰满岳乱妇一区二区三区 | 国产蝌蚪| 亚洲AV成人影视在线观看 | 国产手机在线亚洲精品观看 | 国内色图| 久久精品欧美日韩精品 | 99视频30精品视频在线观看 | 欧美日韩欧美 | 性殴美高清视频 | 国产极品精频在线观看 | 日韩一区二区在线免费观看 | 黄页网站在线观看免费 | 国产毛片在线 | 99久久99久久精品免费看子 | 欧美日韩欧美日韩 | 四虎影视在线看免费 720p | 中字幕久久久人妻熟女 | 熟女人妻-蜜臀AV-首页 | 欧美黑人一区二区三区免费A片 | 日韩精美视频 | 男人的天堂中文字幕熟女人妻 | 高清中文字幕视频在线播 | 97在线视频免费观看97 | 国产美女裸露无遮挡双奶A片游戏 | 亚洲精品久久久久69影院 | 没有废话全色肉的黄文 | 欧美一区二区三区黄色 | 亚洲精品美女偷拍一区二 | 在线观看免费大片 | 亚洲日本一期二期三期精华液 | 2020精品极品国产色在线观看 | 午夜福利合集1000在线 | 99久久免费国产精品 | 日本一道人妻无码一区视频 | 日韩欧美国产免费看清风阁 | 菊色宫快播电影 | 欧美91精品久久久久网免费 | 国模少妇一区二区三区A片 国内精品A片XXX久久久 | 无尺码精品日本欧美 | 色爱区综合激情五月综合激情 | 星野亚希 快播 | 97久久精品国产成人影院 | 日本人丰满xxxxhd | 忘忧草在线影院日本图片 | 亚洲乱码AV中文一区二区 | 国产三级一区二区 | 国产亚洲欧美在线观看三区 | 国产乱码卡二卡三卡43 | 日韩A片无码毛片免费看久久 | 国产又粗又猛又爽又黄A片 国产又粗又猛又爽又黄A片漫 | 色欲天天天综合网免费 | 久久久国产精品免费A片分环卫 | 日本精高清区一 | 黄色一区二区三区 | 揄拍自拍日韩精品 | 日日噜噜夜夜躁躁狠狠 | 天天干天天操天天做 | 亚洲国产激情一区二区三区 | 久久久久久国产精品免费免费 | 蝌蚪窝在线观看高清 | 国产精品免费视频能看 | 日本真人做爰高潮全过程 | 韩国无码又爽又刺激的A片 韩日午夜在线资源一区二区 | 欧美性猛交一区二区三区精品 | 午夜精品福利在线导航小视频 | 国产精品久久久久国产精品三级 | 63jjj全国最大 | 天天操天天舔天天射 | 狠狠i撸| 玖玖精品视频在线 | 精品日韩卡1二2卡3卡4卡乱码 | chinese国产乱在线观看 | 夜夜操夜夜爽 | 波多野结衣中文在线观看 | 好爽又高潮了毛片免费下载 | 国产熟妇精品高潮一区二区三区 | 亚洲午夜精品久久久久久抢 | 久久久免费看少妇高潮A片特黄 | 人人看人人艹 | 色婷婷小说 | 白莲花乖腿打开h调教 | 波多野给衣一区二区三区 | 高辣H文短篇啪啪小说男男 高辣H文黄暴糙汉文H文 | 奇米777四色影视 | 西西人体全身祼体图片 | 亚洲精品成人AA片在线播 | 最新伦理电影_伦理电影 | 影音先锋影院中文无码 | 男人的天堂亚洲一线AV在线观看 | 中文字幕乱码熟妇五十中出 | 久久精品免费视频观看 | 日韩人妻少妇一区二区三区 | 亚洲 暴爽 AV人人爽日日碰 | 天天躁日日躁成人字幕aⅴ 天天躁日日躁aaaaxxxx | 五月丁香婷婷天堂 | 中文字幕5566看片资源 | 日韩亚洲国产高清免费视频 | 艳妇臀荡乳欲伦岳TXT免费下载 | 日韩精品三级 | 国产剧情一卡二卡麻豆 | 1区2区3区高清视频 2017国产小视频 | 日韩在线观看精品 | 色婷婷综合激情 | 亚洲黄页 | 一级做a爱过程免费视频超级 | 亚洲另类第一页 | 国产强伦姧人妻完整版 | 五色网| 国产亚洲欧美日韩v在线 | 亚洲性夜色噜噜噜在线观看不卡 | 秋霞在线观看视频一区二区三区 | 日韩艹 | 成年男人午夜片免费观看 | 久久国产大片 | 簧片在线免费观看 | 一本大道一卡二卡三卡四卡在线观 | 国产在线看不卡一区二区 | 国产一区二区精品在线观看 | 日产精品一二三四区气温 | 黄色网址www | 最新国产麻豆精品 | 国产精品无码人妻99999 | 久久丫精品忘忧草西安品 | 丁香花在线观看免费观看图片 | 日本波多野结衣在线观看 | 桃子视频在线观看高清免费视频 | 草莓ap| 久久精品免视看国产成人2021 | 97av视频| 中日文字字幕乱码视频 | 97久久综合九色综合 | 国产不卡在线看 | 影音先锋色情AV在线看片 | 日韩精品在线视频观看 | 国产一区免费在线观看 | 亚洲怡红院频在线视频 | 国产成人精品综合久久久软件 | 日本三级久久 | www.九色视频| 国产在线观看中文字幕 | 狠狠干狠狠操视频 | 国内精品久久久久久久小说 | 日本无码人妻精品一区二区蜜桃 | 久操手机在线视频 | 免费看黄色一级片 | 免费播放一区二区三区 | 色激情综合 | 黄色在线网站视频 | 综合色网站 | 黄色字幕网 | 亚洲情A成黄在线观看动漫软件 | 亚洲精品一区 | 女医学护士一级毛片 | 色丁香婷婷综合缴情综 | 波多在线 | 一本色道在线久88在线观看片 | 亚洲中文字幕无码一久久区 | 被伴郎的内捧猛烈进出H | 亚洲AV无码一区二区色情蜜芽 | 嗯好爽快点插我视频在线播放 | 婷婷四房综合激情五月性色 | 国产片AV片永久免费观看 | 诱人的女邻居在线观看 | 亚洲中文字幕日本在线观看 | 久久综合丁香 | 最好看的最新的中文字幕1 最好看的最新的中文字幕3 | 国产aaa免费视频国产 | 免费看片A级毛片免费看 | 97色婷婷 | 2024伊人查蕉在线观看 | 99久久综合精品国产 | 房奴试爱MBA | 国产乱子夫妻XX黑人XYX真爽 | 三级网站大全 | 韩国电影甜性涩爱 | 激情啪啪网站 | 全古装一级毛片大全 | 亚洲第一中文 | 波多野结衣免费线在线 | 亚洲日韩aⅴ在线视频 | 欧美阿v不卡资源在线 | 久久91精品国产91久久户 | 色情综合另类小说图片 | 日韩在线视频www色 日韩在线免费观看视频 | 日韩熟女精品一区二区三区 | 国产中文精品无码欧美综合小说 | 久久精品热2018 | 亚洲首页国产精品丝袜 | 少妇偷拍精品高潮少妇 | 色情无码永久免费网站WWW | 天堂 地址 在线 | 视频二区 中文字幕 欧美 | 久操视频在线免费观看 | 伊人大相蕉在线看青青 | 性生交大片免费看 | 国产成人精品亚洲2020 | 日本一卡二卡三卡四卡无卡免费播放 | 欧美性A片又大又长 | 一本三道a无线码一区v小说 | 99久久免费午夜国产精品 | 97神马| 6080新视觉理论在线观看 | 国产理论视频在线观看 | 国产黄的网站免费 | 欧美日韩中文国产一区二区三区 | 国产电影一曲二曲三曲图片 | 日韩免费一区 | 亚洲永久视频 | 亚洲精品乱码久久久久蜜桃 | 久草在在线免在线观看视频 | 亚洲欧美人成网站综合在线 | 91中文字幕在线 | 看黄免费在线 | 久久6699精品国产人妻 | 无码激情全黄做爰片 | 亚洲 视频 在线 国产 精品 | 韩国理伦三级做爰在线播放 | 少妇被猛烈进入A片 | 一级毛片美国 | 国产一级特黄高清免费大片dvd | 亚洲区激情区图片小说区 | 亚洲国产中文在线视频免费 | 韩国三级在线高速影院 | 肉肉描写很细致的黄文 | 欧美三级A做爰在线观看 | 夜色成人 | 亚洲精品久久99蜜芽尤物TV | 亚洲精品国产一区二区贰佰信息网 | 国产精品99精品无码视亚 | 黄页网站在线免费观看 | 午夜福利电影网站鲁片大全 | WWW国产亚洲精品久久久 | 水中色av综合 | 久久AV无码乱码A片无码软件 | 精品日韩免费视频在线观看 | 亚洲经典自拍 | 美女露出尿口让男生爽痛 | 免费看黄色一级 | 性一交一乱一乱A片AP88 | 亚洲欧美啪啪 | 波多野结衣 美乳人妻 | 女医学护士一级毛片 | 国产无遮挡A片又黄又爽小说 | 天天做天天爱天天干 | 欧美XXXXX高潮喷水麻豆 | 久久97久久99久久综合 | 亚洲成人福利网站 | 成人无码T髙潮喷水A片小说 | 老司机精品福利在线资源 | 最近最新中文字幕高清免费 | 蜜桃99影院 | 韩国污网站| 真人做爰片免费视频毛片中文 | 在线综合亚洲欧美网站 | 久久久亚洲天堂 | 一人看片WWW在线视频 | 浪妇叫床叫的很浪的小说 | 100国产精品人妻无码 | 小77论坛 唯美清纯 逍遥五月天 | 校花被房东C得合不拢腿H男男 | 五月婷婷丁香在线 | 爱你千万次韩剧在线观看 | 小视频免费观看在线 | 国产AV亚洲国产AV麻豆 | 四虎影在线 永久免费 | 最新日本中文字幕 | 中国老太婆bbwhd | 国产乱国产乱老熟300部视频 | 国产又色又爽在线观看 | 日日日日人人人夜夜夜2017 | xxxx性bbbb欧美野外 | 狠狠色丁香婷婷综合欧美 | 亚洲AV高清一区二区三区色欲 | 久草丁香 | 男女毛片免费视频看 | 欧美宗合网 | 亚洲日本在线观看 | 激情综合丁香婷婷色五月 | 色窝窝9797eecomehttp最新 | 91香蕉视频在线 | 咪咪久久| 久久三级网站 | 国产护士资源总站 | 丁香婷婷六月综合缴清 | 色综合天天综一个色天天综合网 | 国产精品 欧美在线 另类小说 | 妻子的背叛免费看 | 思思久久好好热精品国产 | 黄页网站视频免费 视频 | 好久被狂躁A片视频无码免费视频 | 欧美日韩精品一区二区三区视频播放 | 精品久久久久不卡无毒 | 成 人 网络| 欧美精品18videosex性欧 | 在线观看黄的网站 | 成 人 网 站免费观看 | tube jizz| 日韩AV爽爽爽久久久久久 | 伦理片无码电影在线看 | 色综合网 | 日本精品无码久久久久三级国产 | 性色AV一区二区三区V视界影院 | 国产精品乱码色情一区二区视频 | 性加拿大高清xxxxx | 影音先锋av看片资源库 | 日本中文字幕免费 | 亚洲欧洲日韩国产一区二区三区 | 亚洲色无色A片一区二区农夫山泉 | 青草草97久热精品视频 | 任你躁| 久久免费精彩视频 | 国产护士在病房AV做爰小说 | 麻豆一卡2卡三卡4卡网站 | 午夜精品福利在线导航小视频 | 久久精品国产免费中文 | 最近最新手机中文字幕在线看 | 亚洲高清视频一区 | 要看影院| 国产激情一级毛片久久久 | 中文国产乱码在线人妻一区二区 | 免费无码国产欧美久久18 | 四川丰满护士毛茸茸 | 五月丁香六月综合缴情基地 | 亚洲色欲色欲WWW在线成人网 | 久久国产亚洲电影天堂 | 色久久一个色综合在线 | 亚洲精品国产自在现线最新 | 美女免费视频一区二区三区 | 午夜影视在线视频观看免费 | 肥老熟妇伦子伦456视频 | 被黑人做的白浆直流 | 色网站在线看 | 丁香花在线观看视频在线 | 欧美又硬又粗进去好爽A片 欧美在线视频一区 | 日产精品卡二卡三卡四卡乱码视频 | 久久国产精品久久精 | 亚洲欧美中文字幕网站大全 | 免费看的一级毛片 | 韩国精品欧美一区二区三区 | 免费永久在线观看黄网站 | 男人使劲躁女人过程A片 | 满天星电影在线观看完整免费 | 精品一区二区三区影片 | 欧美日日夜夜 | 色欲天天婬色婬香视频综合网 | 日本亚洲免费无线码 | 福利区体验区120秒免费 | 欧美在线三级艳情网站 | 91在线一区二区 | 情色五月天首页 | 无码做爰全过程免费的床震 | 亚洲精品久久久久久中文 | 国产又黄又爽又刺激的免费网址 | 亚洲 无码 欧美 经典 | 黄视频网站免费看 | 国产在线观看免费观看不卡 | 影音先锋资源男人网 | 愉拍自拍视频在线播放 | 午夜电影在线观免费 | 网www天堂资源在线 王者荣耀最新更新公告 | 三级国产 | 一区二区乱子伦在线播放 | 免费精品国产人妻国语 | 少妇被躁爽到呻吟全过的小说 | 麻豆优品 | 亚洲色t图 | 精品黑人一区二区三区久久 | 亚洲精品久久久一区 | 欧美熟色妇 | 精品无码国产自产野外拍在线 | 日韩精品视频在线 | 久久综合结合久久很很很97色 | 色偷偷资源网 | 九九黄色网 | 免费女性裸身照无遮挡网站 | 国产女同视频 | 激烈娇喘叫床声床震爽文 | 欧洲裸妇图片大全 | 日本人×体xx艺术 | 日韩精品亚洲专在线电影 | 免费蜜芽官网网址永不失联 | 嗯宝贝叫的再浪一点 | 四虎2024最新免费观看 | 亚洲国产区男人本色在线观看 | 美国一级毛片片aa久久综合 | 一区二区三区精密机械 | 国产免费人成在线视频视频 | 91精品国产综合久久久久久 | 精品久久久影院 | 蜜桃在线线免费观看视频 | 国产黄色片在线播放 | 一二三区乱码不卡手机版 | 神兵小将第一季 | 亚洲AV国产精品无码市川京子 | 黑人狂躁日本妞免费视频 | 亚洲日本免费 | 刺激妇乱子伦短篇 | 欧美精品高清无码 | WWW婷婷AV久久久影片 | 亚洲国产天堂久久九九九 | 亚洲欧美一区二区三区久久 | 亚洲第一区第二区 | 欧美三级电影网 | 国产亚洲精品成人AA片 | 美女内射 | 青草青草久热精品视频在线网站 | 9999在线视频 | 日韩不卡视频在线 | 成人在线视频免费观看 | 久久精品国产亚洲麻豆小说 | 亚洲欧洲精品成人久久曰 | 亚洲成人99| 2021日日摸夜夜添夜夜添影院 | 一本色道久久综合亚洲精品加 | 日韩福利视频导航 | 黄色成年视频 | 亚洲乱色熟女一区二区三区蜜臀 | 国产精品高潮呻吟AV久久小说 | 欧美一区二区三区视频在线 | 中文字幕国产综合 | 野结衣在线 | 丁香婷婷亚洲六月综合色 | 成人福利免费观看体验区 | 国产又爽又黄无码无遮挡在线观看 | 俺也去网站 | 好硬啊进得太深了A片无码视频 | 奇米一区 | 日本免费网 | 国产福利美女福利视频免费看 | 亚洲国产精品久久又爽黄A片 | 国产一区二区三区四区五在线观看 | 免费看黄网站在线看 | 2022国产男人亚洲欧美天堂 | 99r在线 | 色社区| 亚洲成人一区 | 精品人妻无码一区二区三区4 | 婷婷网五月天天综合天天爱 | 日韩A片中文字幕视频免费 日韩MV欧美MV中文无码 | 色中色社区 | 亚洲精品偷拍影视在线观看 | 韩国最污禁网站免费观看 | 18禁无遮挡羞羞污污污污免费 | 国产福利小视频尤物98 | 日本爽爽爽爽爽爽在线观看免 | 文中字幕一区二区三区视频播放 | 最近免费更新中文在线观看 | 日本免费无码一区二区到五区 | 亚洲 日韩 在线 国产 精品 | 老司机免费视频福利0 | 无人在线观看高清视频单曲直播 | 9966久久精品免费看国产 | 久久久国产成人精品 | 精品人妻无码一区二区三区下载 | 精品欧美一区二区在线观看欧美熟 | 伦理片飘花手机在线 | 国产精品天干在线观看 | 日韩video| 欧美黑人添添高潮A片视频 欧美激情无码成人A片 | 欲女春潮bd | 久久精品视频9 | 欧美三级真做在线观看 | 国产在线无码不卡影视影院 | 日韩一区二区三区在线 | 亚洲国产天堂在线观看 | 纯肉高H肉辣浪荡NP论J | 3d肉蒲团观看地址 | 久久视频在线直播 | 午夜精品久久久久久毛片 | 精品亚洲国产成人A片在线观看 | 风平浪静 | 国产亚洲日韩欧美视频 | 国产在线视频自拍 | 国产在线不卡免费播放 | 三级久久| 亚洲欧美一区二区三区导航 | 丁香狠狠色婷婷久久综合 | 国产区视频在线 | 亚洲A片无码成人精品区 | 波多野一区二区 | 夜夜躁狠狠躁日日躁孕妇 | 黄色网址免费在线播放 | 在线观看成人网 | 日本性视频网站 | 九九国产视频 | 久久久久久久综合日本亚洲 | 午夜视频国产 | 国精产品999永久中国有限公司 | 欧美日韩一日韩一线不卡 | 国产亚洲精品在浅麻豆 | 精品一卡二卡三卡四卡视频区 | 蜜桃TV成人网站免费打开 | 日本VS中国VS亚洲看无码A | 国产aaaaaa | 在线伦理电影 | 交换邻居波多野结衣中文字幕 | 亚洲大尺度吃奶做爰电影 | 久久九九免费视频 | 久久精品免费电影 | 波多野结衣在线视频观看 | www.欧美色 | 日韩手机在线 | avtt2015天堂网 | 中文字幕理伦电影免费 | 蜜臀AV色欲A片精品一区 | 在线网站黄色 | 天天操天天摸天天曰天天干天天弄天天干 | 楚乔传第二部免费观看全集 | 国产午夜亚洲精品区 | 久久久久久国产精品免费 | 88av免费观看入口在线 | 色激情五月 | 2024亚洲电影最新地址 | 女人张开腿让男人桶爽免 | 国产真实乱人偷精品人妻69 | 97在线免费视频观看 | 国产麻豆精品乱码一区 | 深夜特黄a级毛片免费播放 深夜日韩 | 99精品久久久久久国产人妻 | 午夜激情影院 | 色妞网站 | 欧美激情中文字幕一区二区 | 免费的毛片视频 | 国语自产拍在线观看偷拍 | 中国xx爽69护士 | 91精品国产综合久久精品 | 亚洲情综合五月天 | 成人免费又大又爽A片视频 成人女人A级毛片免费软件 | 亚洲精品一区二区三区婷婷月色 | 91色老久久精品偷偷蜜臀 | 日韩亚洲欧美中文高清在线 | freev de0x x性欧美f| 天天操夜夜草 | 国产午夜男女爽爽爽爽爽 | 国产无遮挡A片无码免费 | 特级毛片绝黄A片免费播冫 特级毛片在线大全免费播放 | 色情成人吃奶激情视频在线播放 | 樱花草视频在线观看WWW在线观看 | 一道本视频一二三区 | 国产露脸无码A区久久 | 黄网在线视频 | 欧美色第一页 | 午夜日韩精品 | 97蜜桃网站 | 鲁丝一区二区三区不属 | 色婷婷色丁香 | 国产亚洲自拍一区 | 久久无码人妻AV精品一区 | 在线观看日本免费 | 免费 高清 日本视频 | 能播放女人另类ZOOZ0 | 欧美乱妇日本无乱码特黄大片 | 日本黄页88网成品网站 | 国自产拍偷拍福利精品啪啪 | 中文字幕在线看 | 亚洲精品免费在线 | 亚洲黄色网址 | 97色情在线观看免费高清 | 毛片大片免费看 | 国产精品久久精品 | 亚洲一级黄色毛片 | 免费观看欧美日韩亚洲 | 日韩精品综合 | 亚洲日本va中文字幕婷婷 | 宝贝乘把腿张开让我添动态图 | 97国产无遮挡A片又黄又爽小说 | 999国产精华是正规产品吗 | 国产一区二区不卡亚洲涩情 | 亚洲天堂精品视频 | 国产人妻久久精品一区 | 99精品噜噜噜成人AV | 最近最新中文字幕在线手机版 | 成人色网 | 国产日韩高清一区二区三区 | 无人高清视频免费观看在线下载 | 成全在线观看免费观看大全 | 口工绅士里番中文全彩 | 美妇吞吐粗长撞击迎合 | 奥様は元ヤリマン | 妺妺窝人体色WWW网站 | 2020色愉拍亚洲偷自拍 | 2024美女视频黄频大全视频 | 成人国产精品一级毛片了 | 日本韩国三级 | 欧美性猛交XXXX乱大交3 | 激情电影色影音先锋 | 国产中文字幕免费观看 | 国产一国产一区秋霞在线观看 | 国产精品人妻在线观看 | 无码人妻视频又大又粗欧美 | 黄色片网址在线观看 | 丁香综合网 | 偷偷色在线 男人天堂 | 国产精品人妻99一区二区 | 日本xxx在线观看免费播放 | 九九热在线免费观看 | 四虎影视214HU永久免费观看 | 国产三级电影网站 | 精品综合久久久久久蜜月 | 极品少妇粉嫩小泬啪啪小说 | 性一交一乱一美A片裸体 | 久久精品女人天堂 | 国产黄色片在线播放 | 欧美香蕉人人人人人人爱 | 老司机午夜剧场 | 亚洲成色A片202477在线小说 | 狠狠影院 | 国产 有码 无码 电影 | 天堂网2014a天堂网 | 特级A欧美做爰AAAAA片 | 波多野结衣中文字幕教师 | 欧美一区精品 | 最近最新免费中文字幕MV | 1769国产精品一区2区 | 久久aa毛片免费播放嗯啊 | 亚洲AV无码乱码A片无码18禁 | 天天草夜夜操 | 好男人视频社区精品免费 | 国产亚洲精品久久久久的角色 | 色中色 最新地址 | 国产97人妻人人做人碰人人爽 | 精品人妻伦一二三区久久AAA片 | 欧美亚洲国产激情一区二区 | 波多野结衣手机视频一区 | 激情五月黄色 | 亚洲激情中文字幕 | 无码人妻精品一区二区三区蜜臀 | 日本毛片高清免费视频 | 波多野结衣在线视频免费观看 | 粗大的内捧猛烈进出A片黄 粗大的内捧猛烈进出A片小说 | 韩剧 你是我的命运 | 欧洲一区 | 亚洲色一色噜一噜噜噜 | 狠狠干很很操 | 91精品人成在线观看 | 永久免费看黄A片无码软件 永久免费毛片 | 一女被多男灌满白浆受孕 | 中文一卡二卡三卡四卡免费 | 日韩精品资源 | 国产成人理在线观看视频 | 亚洲AV电影天堂男人的天堂 | 亚洲精品久久久无码大桥未久 | 韩国成人理伦片免费播放 | 四虎国产成人免费观看 | 日韩经典AV在线观看 | 超碰97人人无马 | 国产天天操 | 亚洲暴爽AV天天爽日日碰 | 狠狠色丁香久久综合五月 | 亚洲狼人综合网 | 翁公的粗大挺进我的密道 |