Number Density Distribution of Small Particles around a Large Particle

Sep 29, 2016 - colloidal suspension where small and large colloidal particles exist and we experimentally obtain the number density distribution of sm...
0 downloads 0 Views 2MB Size
Subscriber access provided by Northern Illinois University

Article

Number Density Distribution of Small Particles around a Large Particle: Structural Analysis of a Colloidal Suspension Ken-ichi Amano, Mitsuhiro Iwaki, Kota Hashimoto, Kazuhiro Fukami, Naoya Nishi, Ohgi Takahashi, and Tetsuo Sakka Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.6b02628 • Publication Date (Web): 29 Sep 2016 Downloaded from http://pubs.acs.org on October 4, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Langmuir is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Number Density Distribution of Small Particles around a Large Particle: Structural Analysis of a Colloidal Suspension

Ken-ichi Amano, *a Mitsuhiro Iwaki, bc Kota Hashimoto,a Kazuhiro Fukami,d Naoya Nishi, a Ohgi Takahashi,e and Tetsuo Sakkaa

a

Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan. E-mail: [email protected] b Quantitative Biology Center, RIKEN, Suita, Osaka 565-0874, Japan. c Graduate School of Frontier Biosciences, Osaka University, Suita, Osaka, 565-0874, Japan. d Department of Materials Science and Engineering, Graduate School of Engineering, Kyoto University, Kyoto 606-8501, Japan. e Faculty of Pharmaceutical Sciences, Tohoku Medical and Pharmaceutical University, Sendai 981-8558, Japan.

ABSTRACT Some colloidal suspensions contain two types of particles—small and large particles—to improve the lubricating ability, light absorptivity, etc. Structural and chemical analyses of such colloidal suspensions are often performed to understand their properties. In a structural analysis study, the observation of the number density distribution of small particles around a large particle (g LS ) is difficult because these particles are randomly moving within the colloidal suspension by Brownian motion. We obtained g LS using the data from a line optical tweezer (LOT) which can measure the potential of mean force between two large colloidal particles (Φ LL). We propose a theory that transforms Φ LL into g LS . The transform theory is explained in detail and tested. We demonstrate for the first time that LOT can be used for the structural analysis of a colloidal suspension. LOT combined with the transform theory will facilitate structural analyses of the colloidal suspensions, which is important for both understanding colloidal properties and developing colloidal products.

1

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1 Introduction Colloidal suspensions are used in various industrial products, such as coating materials, cosmetics, pharmaceuticals, and food products. Optimization of the physical and chemical properties of colloidal suspensions is indispensable to improve the quality of the products, and information about the microscopic structure of the colloidal particles is important to understand the various properties of the colloidal suspensions. The structural information for mixtures of two or more different colloidal particles is also important, and such mixtures are often used to improve the lubricating ability, light absorptivity, heat capacity, durability, etc., of the colloidal suspension. In the present paper, we focus on the structure of a binary colloidal suspension where small and large colloidal particles exist, and we experimentally obtain the number density distribution of small colloidal particles around a large colloidal particle (g LS ). However, the observation of g LS is difficult because the particles are randomly moving by Brownian motion. This difficulty is overcome by using a line optical tweezer (LOT)1 and proposing a transform theory. LOT can measure the potential of mean force between two large colloidal particles in a sea of small colloidal particles, and the transform theory can provide g LS from the potential of mean force. The combination of LOT and the transform theory is useful for the observation of g LS . The neutron contrast variation method of small-angle neutron scattering (SANS)2 is known as an experimental method for measuring g LS . SANS is helpful for measuring the structural properties of colloidal suspensions; however, it cannot be performed in a general lab space. Hence, LOT could provide a convenient method for the measurement of g LS. LOT can measure the potential of mean force between two large colloidal particles (Φ LL), but until now, g LS has not been obtained from Φ LL. For that reason, we propose a theory that can transform Φ LL into g LS. The transform theory is derived from a statistical mechanics a simple liquid. 3 The small and large colloidal particles are modeled as small and large spheres, respectively. In the transform theory, the dispersing media (liquid) is treated as a continuous fluid. The combination of LOT and the transform theory enables us to measure g LS in a general lab space without deuteration. This method will facilitate structural analyses of the colloidal suspensions, which is important for both understanding and developing colloidal products. Several apparatuses can measure subtle interactions between two surfaces. For example, LOT, one of the aforementioned optical tweezers, can measure the potential of mean force between two large colloidal particles. 1 Stepwise movements of kinesins4 and myosins5,6 have been measured by optical tweezers (kinesin and myosin are linear motor proteins). The forces during the stretching of a DNA chain 7 and a protein8,9 have also been measured by optical tweezers. Another technique for measuring this type of subtle interaction is atomic force microscopy (AFM). AFM can measure the oscillatory force 10–14 originating from the solvation layers on the surface of the substrate, which is 2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

the so-called solvation force. Furthermore, with a colloidal probe, AFM can measure the force between two colloidal particles. 15 Surface force apparatus (SFA) can also measure the oscillatory forces in liquids 16–19 and colloidal suspensions. 20 LOT, AFM, and SFA can measure the subtle interactions between two surfaces; however, the observation of the number density distribution of small spheres on a sample surface from the measured subtle interaction is one of the difficult tasks in this area of study. The theories that transform from the subtle interactions to the number density distributions in the research fields of AFM 21–24 and SFA 25,26 have already been studied. However, such a study has not been performed for LOT. The layer structures of colloidal particles on some surfaces have been studied from theoretical and experimental points of view. For example, fundamental measure theory (FMT)27–29 and integral equation theory (IET)30–32 are well known within the field. The simplest and most constitutive model of a colloidal suspension is an ensemble of small hard spheres, and FMT and IET can very precisely calculate the layer structure using the model. From these theories, it has been found that the granularity (i.e., the volumes of the colloidal particles themselves) is essential for the formation of the layer structure. Experimentally, the layer structure is studied by observing the oscillatory force curve between a colloidal prove and a sample surface using colloidal probe AFM. 33–35 The oscillatory force curve represents the existence of the layer structure. The influences of the volume fraction, surface roughness, etc., on the oscillatory force curve are also studied using colloidal probe AFM. The observation of the oscillatory force is a remarkable achievement; however, it would be valuable if we could also evaluate the number density distribution of the colloidal particles on the sample surface from the force data. We believe that the proposition of a theory that can transform from the measured force curve into the number density distribution would enable us to analyze the layer structure in more detail. Moreover, the transform theory would bridge the gap between theory and experiment. In this paper, we propose a transform theory to calculate g LS (the number density distribution of small colloidal particles around a large colloidal particle) from the input data of Φ LL (the potential of mean force between the two large colloidal particles). This paper is composed of five chapters. In Chapter 2, the transform theory is explained in detail. We propose two processes for the calculation of g LS . In one process, we focus on the contact number density of the small colloidal particles around the large colloidal particle. Thus, it is called the contact density (CD) process. In the other process, we focus on the pressure on a surface element of the large colloidal particle’s excluded volume. Hence, it is called the surface element (SE) process. Eventually, both processes are proven to be the same. Hence, we collectively call them particle pair (PP) transform. In our opinion, both processes are important because the CD process justifies the unclear Derjaguin’s assumption 36–38 used in the SE process. In Chapter 3, the computational details for the verification tests of the PP transform are described. In 3

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chapter 4, the results of the theoretical (computational) tests of the PP transform are shown. In addition, the experimentally obtained Φ LL1 is transformed into g LS. Finally, we present the conclusions in Chapter 5.

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

2 Theory 2.1 Calculation of the density distribution using the CD process We illustrate the model system for the CD (contact density) process in Fig. 1, where large spheres 1 and 2 with radii, r L, are immersed in an ensemble of small spheres. The radius of the small sphere is rS . The bulk number density of the small spheres is ρ 0 . The separation between the centers of the two large spheres is expressed as s. The distance between the center of large sphere 1 and the center of the small sphere in contact with large sphere 2 is expressed as l. The meaning of θ is illustrated in Fig. 1.

Figure 1. System configuration of the CD process.

The large sphere has a radial density distribution of small spheres around itself, and the radial density distributions of the two large spheres overlap when they approach each other. The overlap creates a complex distribution of small spheres within the confined space. In a simple form, the overlap can be expressed by the Kirkwood superposition approximation: 39–41   ;         , 1

where g W, g LS1 , and g LS2 represent the normalized number density distribution of the small spheres in the whole system and the normalized number density distributions around large spheres 1 and 2, respectively. We notice that g LS1 and g LS2 are the same (g LS1 = g LS2 = g LS ), and they are the distributions in the bulk (in the isolated condition). Hereafter, we simply mention the normalized number density distribution as the density distribution. The three-dimensional (3D) vector, ɤ, starts from the center of large sphere 1, and it indicates an arbitrary position from the center. The 3D vector, s, 5

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

represents the separation between the centers of large spheres 1 and 2 and the direction. The Kirkwood superposition approximation is exact when large spheres 1 and 2 are sufficiently separated and/or when there are only three bodies (i.e., one small sphere and two large spheres) in the system. However, as the bulk number density of the small spheres increases, the approximation becomes invalid. Next, we employ an approximation that the force (f) acting on large sphere 2 depends only on the contact density of the small spheres. 25,29,42 That is, the two-body potential between the large and small spheres is approximated as a rigid potential. In this case, f is expressed as follows: 20,25,29,42

 = 2      

(/



   ,  !sin cos '   , 2

where f simply means:

 = 

' *   + 

. 3

'

The functions Φ LL and u LL represent the potential of mean force and the (bare) two-body potential between the two large spheres, respectively. The subtracted product (Φ LL − u LL) represents the potential induced by the small spheres. Φ LL is exactly related to the normalized number density between the two spheres (g LL) as follows:

 = exp1* /2  3, 4

where k B and T are the Boltzmann constant and absolute temperature, respectively. The symbol g LSC represents the number density at the contact point; π is the circular constant and r ≡ r L + rS . The symbol, f B (constant), is the force acting on the back (right side) of large sphere 2, which is expressed as:

 = 2      

(/



sin cos ' =       . 5

As shown in Fig. 1, l can be expressed as:

 = 6  +    2 sin . 6

Differentiation of l with respect to θ is written as:

6

ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

64       +       ' = . 7

' 2 Furthermore, sinθ and cosθ are written as:

  +     , 8

2

sin

=

cos

=;

4       +       . 9

4   

Consequently, Eq. (2) can be rewritten as: √>  +  =     >BA

? @A ?

  +     

   ' . 10



Eq. (10) can be expressed in the form of a matrix operation, as follows:

D∗ = FG , 11

where F * corresponds to the left-hand side of Eq. (10). G1 and H correspond to g LS1 (l) and all the remaining parts, respectively. We use sufficiently large matrices of F *, G 1 , and H . The boundary conditions of F * are (f(2r) + f B )/(ρ 0 k B Tπg LSC ) at the upper side and f B /(ρ 0 k B Tπg LSC ) at the lower side. The latter boundary condition originates from f(s∞ ) = 0, where s ∞ represents a sufficiently long separation. The boundary conditions of G 1 are g LS1 (r) at the upper side and 1 at the lower side. g LS1 (r) is the contact value of the density distribution that is the same as g LSC , and it is an unknown value at this stage. To solve Eq. (11), we set H as a square matrix whose variables are l and s. A cell of matrix H is composed of l(s2 + r 2 + l 2 )dl/s 2 inside the integral range and 0 outside the integral range. At the lower right area of H , where correspondent g LS1 (l) converges to 1, the integral range protrudes from the matrix area. For such an area, the block matrix is a square unit matrix multiplied by f B /(ρ 0 k B Tπg LSC ). To obtain G1 , an arbitrary initial value is substituted into g LSC , and a trial F * is prepared. G1 is numerically predicted by multiplying the inverse matrix of H by F *. Because large spheres 1 and 2 are equivalent, when the output g LS1 (r) at the contact is equal to g LSC of the input value, the numerical calculation converges. The convergence is, for instance, operated by a bisection method.

7

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

2.2 Calculation of the density distribution using the SE process The outline of the SE (surface element) process is as follows. First, the force curve between the two large spheres is converted to pressure on a projected surface element of the large sphere’s excluded volume (see Fig. 2). The excluded volume of the large sphere is defined as the volume that the center of the small sphere cannot enter. We call this conversion the FPSE conversion, where FPSE means force to pressure on a surface element. Second, the pressure on the surface element is transformed into a density distribution of small spheres around the large sphere by applying the transform theory recently proposed by Amano. 25

Figure 2 . System configuration of the SE process.

The FPSE conversion is an exact version of the Derjaguin approximation. 36–38 The Derjaguin approximation has been applied in many studies due to its universality and validity. However, its applicability is limited to large particles, and it is restricted to very short surface–surface separations. Bhattacharjee and Elimelech also proposed an exact version of the Derjaguin approximation called the surface element integration (SEI). 43,44 The difference between the SEI and FPSE conversion is that SEI can convert the interaction between two parallel flat substrates into the interaction between two spheres, while the FPSE conversion can change its inverse transformation. The applicability of the FPSE conversion is not limited to large particles and is not restricted to very short surface–surface separations because the FPSE conversion is an exact version of the Derjaguin approximation. (We note that the three methods, the Derjaguin approximation, SEI conversion, and FPSE conversion, each contain Derjaguin’s assumption, which is explained in the next paragraph and Eq. (12).) The system configuration for the SE process is shown in Fig. 2. There are many small spheres with a bulk density, ρ 0 ; large spheres 1 and 2 are immersed in the ensemble of small spheres. The space where the center of the small sphere cannot enter 8

ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

is the excluded volume of the large sphere, and r is the radius of the excluded volume. rS and r L are the radii of the small and large spheres, respectively. The meaning of θ is illustrated in Fig. 2. The separation between the centers of large spheres 1 and 2 is represented as s. The length of the horizontal line between a pair of surface elements of the excluded volumes is represented as lʹ. If the force between the two large spheres, f, is expressed by a summation of pressures on the surface elements. Thus, the force can be written as:

 = H I J KL ′;  , 12

NO

where P is the pressure along the x axis and A 2 x is the projected flat area of the surface element of the excluded volume of large sphere 2. We note that Eq. (12) represents Derjaguin’s assumption. 36–38 In the present case, Eq. (12) can be rewritten as: P 

 = 2  I′ sin cos ' . 13





Here, P is equal to the pressure from the left side plus that from the right side of large sphere 2. In accordance with Fig. 2, l′ can be written as:

 J =   2sin . 14

Hence, the following two expressions are obtained:

cos ' = 1/2 '′, 15

sin =   ′ /2 . 16

Then, Eq. (13) is rewritten as:

2  / = 

>

>BA

I J    J '′. 17

If we differentiate both sides of Eq. (17) with respect to s, we obtain: > 2 Q 

= 2I  2 +  I J ' J . 18  Q >BA

9

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Moreover, if we eliminate the integral term of Eq. (18), the Derjaguin approximation for the granular system 37 can be obtained. The Derjaguin approximation for the granular system indicates that (i) the medium around the large sphere is not a continuum fluid but an ensemble of small spheres and (ii) the surface element of the excluded volume is considered as the surface element. However, in the SE process, we would like to obtain P from f without eliminating the integral term. Then, we make use of the matrix operation below:

D# = ST, 19

where F # corresponds to the left-hand side of Eq. (17). P and J correspond to P(lʹ) and all the remaining parts, respectively. We use sufficiently large matrixes of F #, P, and J . The boundary conditions of F # are 2f(2r)/π at the upper side and 0 at the lower side because f(s) = 0 at a sufficiently large separation. The boundary conditions of P are P(0) at the upper side and 0 at the lower side. To obtain P from Eq. (19), we set J as a square matrix whose variables are lʹ and s. A cell of matrix J is composed of (s – lʹ)dlʹ inside the integral range and 0 outside the integral range. At the lower right of J , the integral range protrudes from the matrix area. For such an area, the block matrix is a square unit matrix multiplied by 2r 2. P is numerically obtained by multiplying the inverse matrix of J by F # . This is an example of the solving process of Eq. (19), and the FPSE conversion is finished here. We found that the result of the FPSE conversion is similar to the result of the Derjaguin approximation for the granular system; when r L >> rS (r ≈ r L), both results are almost the same. Next, we transform P (i.e., P(lʹ)) into g LS1 . In this step, we approximate that the force acting on large sphere 2 depends only on the contact density of the small spheres. 25,29,42 By applying this approximation, g LS1 is readily calculated as follows:25

  |V + WJ | =

I ′

+ 1, 20

   

where

  =

1 + 61 + 4I0 / 

. 21

2

The bold characters r and l′ represent the 3D vectors of r and l′, respectively. The start point of r is the center of large sphere 1 and its end point is connected to the start point of l′.

10

ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

2.3 Relation between the CD and SE processes Here, we demonstrate that the CD and SE processes are identical, where we suppose that the force acting on the large sphere depends only on the contact density of the small spheres. We show that Eq. (17) in the SE process leads to Eq. (10) in the CD process. First, the pressure on the surface element of large sphere 2 is given by: 25

I′ =  X ′ /KL ′ =       |V + WJ |      . 22

Note that Eq. (22) indicates that the force acting on the surface element (f SE) is the sum of the forces from the right and left sides. Substituting Eq. (22) into Eq. (17), we obtain:

  =

>       |V + WJ |  1    J '′. 23

2    >BA 

Next, the relation between l and lʹ can be obtained by comparing Figs. 1 and 2:

J =      /. 24

Thus, Eq. (23) is rewritten as:

 =     

√>? @A ?

>BA

  |V + WJ |  1   +     

'. 25



Since the value of g LS1(|r + l ʹ|) is equal to that of g LS1 (l), Eq. (25) is rewritten as:

 =     

√>? @A ?

>BA

   +     

   '   . 26 

Eq. (26) is the same as Eq. (10). Therefore, the values of g LS1 calculated by the CD and SE processes are the same. In our opinion, the CD process is important for corroborating the validity of Derjaguin’s assumption in the SE process.

11

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3 Computational details Because the CD (contact density) and SE (surface element) processes produce the same results, we collectively call these processes the PP (particle pair) transform. In this chapter, we explain the computational details of the tests for the PP transform. In the tests, we employ rigid systems and a one-dimensional Ornstein–Zernike equation coupled with a hypernetted-chain closure (1D-OZ-HNC). 3,45–47 A bridge function proposed by Verlet, Choudhury, and Ghosh (VCG)30 is incorporated in 1D-OZ-HNC to precisely perform the verification test. The incorporated VCG bridge function is useful because the function becomes a bridge function of the Verlet type 31,32 (a highly sophisticated bridge function for the rigid system) when the system is rigid. The 1D-OZ-HNC with the VCG bridge function generates reliable results for the tests. The verification test is performed in a computationally closed cycle. The outline of the verification test is as follows. In step (A), an ensemble of small spheres is prepared in a radially symmetric space, and a pair correlation function between the two small spheres is computed. In step (B), a single large sphere is immersed in the ensemble of small spheres, and the density distribution of the small spheres around the large sphere (g O ) is calculated. We mention g O as the original density distribution. In step (C), one additional large sphere is immersed in the system, and the potential of mean force between the two large spheres (Φ LL) is computed. In step (D), Φ LL is changed to the mean force through Eq. (3). Then, in step (E), the density distribution of the small spheres around the large sphere (g LS ) is calculated from the force curve using the PP transform. The reproducibility of the PP transform is validated when the g LS calculated using the PP transform agrees with the original one (g O). The model systems for the verification tests are as follows. The diameters of the large and small spheres are d L and d S, respectively. Here, d S is an arbitrary short length, which is a unit length for the calculation. Several values of d L were tested, but in this paper, only the results obtained in the case of d L = 10d S are shown. (Almost the same reproducibilities were also obtained when d L = d S, 2d S , 3d S, 4d S , 5d S , 15d S , 20dS , etc.) The grid spacing is set at dS /50, and the number of grid points is 2048. The number of grid points (2048) is sufficiently large because the center of the large sphere is placed on the boundary of the grid line. The bulk number density of the small spheres is ρ 0 . The tested volume fractions (φ S) are 0.1, 0.2, 0.3, and 0.4 (φ S = ρ 0 d S 3 π/6). The tested results obtained for the rigid systems are shown in Section 4.1, and those for the non-rigid systems are shown in the Supporting Information. In the derivation process for the PP transform (CD and SE processes), the two-body potential between the large and small spheres is supposed to be rigid, but the other two-body potentials were arbitrary. Hence, the PP transform can be used in an ensemble of soft small spheres. We found that the PP transform can also work well in the ensemble of soft small spheres. Moreover, the PP transform can be used even when the two-body 12

ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

potential between the large and small spheres is soft. However, we found that when the depth of the potential minimum is much lower than –k B T, the reproducibility of the PP transform is not very high. As a result, we conclude that the PP transform works well when the two-body potential between the large and small spheres is close to the rigid potential.

13

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4 Results and discussion 4.1 Tests for the PP transform We show the results of the tests conducted for the four volume fractions. Fig. 3 shows the potential of mean force between the large spheres (Φ LL) as a function of the distance between the centers of the large spheres calculated by 1D-OZ-HNC-VCG. As observed in Fig. 3, the amplitude of the curve increases as the volume fraction (φ S ) increases. By using the PP transform, the curves are transformed into the number density distributions of the small spheres around the large sphere (see Fig. 4). The bold solid line represents the original number density, gO, which is the benchmark structure. The thin solid line represents g LS obtained from Φ LL. As shown in Fig. 4, the reproducibility of the number density is high when the volume fraction is low. This trend originates from the Kirkwood superposition approximation because the approximation becomes precise as the volume fraction decreases.

Figure 3 . Potential of mean force between the two large spheres calculated by 1D-OZ-HNC-VCG.

14

ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 4 . Comparison among g O, g LS , and ğ LS . φS = 0.1, 0.2, 0.3, and 0.4 in (A), (B), (C), and (D), respectively.

The thin solid line in Fig. 4(D) shows the wrongly underestimated minimum around r LS /d S = 5.8, the shape of which is unusual compared to the benchmark structure. To improve the g LS , it is approximated as follows:

  = exp1*  /2  3 1  *  /2  , 27

where Φ LS is the potential of mean force between the large and small spheres. Then, the modified density distribution, ğ LS , is empirically expressed as:

15

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ğ  =   when   ≥ 1, 28a ğ  = exp1   13 when   < 1. 28b The introduction of the method corrects the wrongly underestimated g LS. In all of the tested cases, we found that ğ LS is better than g LS. The layer spacings of g LS (ğ LS ) are shorter than that of gO. The reason for the decrease is as follows. When the small spheres are confined between the two large spheres, the layer structure is compressed. However, the PP transform cannot correct for the compression (shortening) of the layer spacings because it uses the Kirkwood superposition approximation to treat the density distribution of the small spheres in the confined space. To obtain more precise results, this problem should be solved in the future.

16

ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

4.2 Practical use of the PP transform In this section, the PP transform is applied to LOT (line optical tweezer) to obtain the density distribution of small colloidal particles around a large colloidal particle from the experimental data. Crocker et al. 1 measured the potential of mean force between poly(methyl methacrylate) (PMMA) particles in a colloidal suspension of polystyrene (PS) particles. The diameters of the PS and PMMA particles are 83 nm (Seradyn, Inc.) and 1100 ± 15 nm (Bangs Labs, Inc.) in products, respectively. The colloidal suspension contains 5 mM NaCl and 5 mM sodium dodecyl sulfate, which prevents aggregation. 1 The screening lengths of the repulsive interactions between the colloid particles are about 3 nm. 1 Since the screening length is so small compared to the diameters of the colloid particles, they are treated as hard spheres. 1 However, the repulsive interaction between the PS particles causes its effective diameter to be slightly larger than the actual diameters. Crocker et al. applied the Asakura-Oosawa theory, 48,49 and estimated the effective diameter of a PS particle as 97 ± 6 nm. 1 However, it is likely that the effective diameter of a PS particle is overestimated because the diameter was determined in the low volume fractions. Therefore, we assume the effective diameter of a PS particle as 91 nm (97 nm – 6 nm, taken from “± 6”). However, we do not need to determine the effective diameter of a PMMA particle here because the PS-induced potential between the PMMA particles (Fig. 5) almost automatically yields the value. The PS-induced potential between the PMMA particles in φ PS = 0.34 (φPS : volume fraction of PS particles) is calculated by subtracting the (bare) two-body potential between the PMMA particles from the potential of mean force between the PMMA particles. The (bare) two-body potential is obtained from the potential of mean force measured as φ PS = 0.00. The volume fraction of PMMA particles is less than 10 –7 in both cases (φ PS = 0.00 and 0.34). 1 The fitting is performed using a polynomial function of sexp , where sexp is the separation between the centers of the PMMA particles.

17

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5 . Potential between PMMA particles in the colloidal suspension. The thin and bold solid lines are PS-induced potentials (the same thing). 1 The bold solid line is used for the input of the PP transform. The dotted line is an arbitrarily drawn potential curve that contains PS-induced and steric repulsive potentials.

The PP transform is performed using the bold solid line in Fig. 5. At first, the PS-induced potential (bold solid line) is transformed into the pressure (P(l′)) on the projected surface element of PMMA’s excluded volume using the FPSE conversion. It corresponds to the pressure between the two flat surfaces. Secondly, the point of the first maximum peak of P(l′) is determined, which corresponds to P(0). When the separation between two flat surfaces is just equal to the diameter of a small colloidal particle, the force between the two flat surfaces takes a maximum value. 50,51 Thirdly, the contact density of PS particles around a PMMA particle is calculated using Eq. (21). Finally, the density distribution of PS particles around a PMMA particle is obtained using Eq. (20). Fig. 6 is the density distribution of PS particles around a PMMA particle (ğ exp ), where the modification (Eq. (28)) is used. The ratios of the amplitudes of the potential and the density are not the same (by comparing Figs. 5 and 6). There are three observable layers in ğexp , which is very similar to the liquid structure around a solute. 52–58 Since three PS layers are semi-adsorbed around a PMMA particle, the 18

ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

layers may move with the PMMA particle. That is, when the dynamics of the PMMA particle are studied in the PS colloidal suspension, the layer structure around PMMA should be considered.

Figure 6. Normalized number density of PS particles around a PMMA particle.

Next, we use Eq. (1) to calculate the density distribution of the PS particles around a pair of PMMA particles. Fig. 7(A) shows the density distribution at sexp = 1.3 µm. The repulsive force arises between the PMMA particles when s exp = 1.3 µm (from the slope of Fig. 5). In this situation, the contact density of PS particles on the inside of PMMA particles is expected to be high compared to that on the outside of the PMMA particles. Actually, as shown in Fig. 7(B), the inside contact density is relatively high compared to the outside contact density. This density distribution generates the repulsive force between the PMMA particles.

19

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Normalized number density of PS particles around a pair of PMMA particles (A), and its magnified view (B). The separation between the centers of PMMA particles (s exp) is 1.3 µm. The white areas represent the excluded volumes of the PMMA particles, in which the centers of the PS particles cannot enter. (C) Color scale of the normalized number density.

Crocker et al. predicted “the formation of shells around the large sphere” 1 from their experimental results. Consequently, we were able to visualize their prediction. Furthermore, they discovered a dramatic slowing of Brownian motion for the trapped PMMA particles for a high volume fraction of PS particles. In our opinion, the slowing is closely related to the layer structure of PS particles around a PMMA particle.

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

5 Conclusions We have proposed a theory (PP transform) that transforms from the potential of mean force between two large colloidal particles to the density distribution of small colloidal particles around the large colloidal particles. The PP transform was derived through two processes: the CD (contact density) and SE (surface element) processes. We demonstrated that the CD and SE processes are mathematically identical. Verification tests of the PP transform were conducted in a computationally closed cycle. The PP transform reproduced the density distribution of small spheres around a large sphere. The potential of mean force between the PMMA particles measured by LOT (line optical tweezer) was transformed into the density distribution of PS particles around a PMMA particle for the first attempt using the PP transform for practical use. We demonstrated that LOT combined with the PP transform can be used for the structural analysis of colloidal suspensions. We believe that LOT combined with the PP transform will help studies of shear force, lubricating ability, etc., because the density distribution of the colloidal particles is closely related to the mechanisms. In addition, we expect that the combination of fluorescence microscopy and the PP transform can also be used for the structural analyses of colloidal suspensions because fluorescence microscopy can also detect the potential of mean force between two labeled colloidal particles. We believe that this study makes the structural analyses of the colloidal suspensions easier and facilitates developments of colloidal products (e.g., coating materials, cosmetics, pharmaceuticals, and food products) because LOT combined with the PP transform can measure the density distribution without SANS and deuteration.

21

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Associated Content Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: . Results of the computational test of PP transform in non-rigid systems (PDF)

Author Information Corresponding author *E-mail: [email protected] Notes The authors declare no competing financial interest

Acknowledgements We appreciate support from Masahiro Kinoshita, Kozue Maruta, and Ryosuke Sawazumi. This work was supported by a Grant-in-Aid for Young Scientists (B) from Japan Society for the Promotion of Science (15K21100) and partly supported by a Grant-in-Aid for Scientific Research (B) from Japan Society for the Promotion of Science (15H03877).

References (1) Crocker, J. C.; Matteo, J. A.; Dinsmore, A. D.; Yodh, A. G. Entropic Attraction and Repulsion in Binary Colloids Probed with a Line Optical Tweezer. Phys. Rev. Lett. 1999 , 82 (21), 4352–4355. (2) Lutterbach, N.; Versmold, H.; Reus, V.; Belloni, L.; Zemb, T.; Linder, P. Charge-Stabilized Liquidlike Ordered Binary Colloidal Suspensions. 2. Partial Structure Factors Determined by Small-Angle Neutron Scattering. Langmuir 1999 , 15 (6), 345–352. (3) Hansen, J.-P.; McDonald, I. R. Theory of Simple Liquids; Academic Press: London, 1986 . 22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(4) Taniguchi, Y.; Nishiyama, M.; Ishii, Y.; Yanagida, T. Entropy Rectifies the Brownian Steps of Kinesin. Nat. Chem. Biol. 2005 , 1 (6), 342–347. (5) Oguchi, Y.; Mikhailenko, S. V; Ohki, T.; Olivares, A. O.; De La Cruz, E. M.; Ishiwata, S. Load-Dependent ADP Binding to Myosins V and VI: Implications for Subunit Coordination and Function. Proc. Natl. Acad. Sci. U. S. A. 2008 , 105 (22), 7714–7719. (6) Yanagida, T.; Iwaki, M.; Ishii, Y. Single Molecule Measurements and Molecular Motors. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 2008, 363 (1500), 2123–2134. (7) Williams, M. C.; Wenner, J. R.; Rouzina, I.; Bloomfield, V. A. Entropy and Heat Capacity of DNA Melting from Temperature Dependence of Single Molecule Stretching. Biophys. J. 2001 , 80 (4), 1932–1939. (8) Elms, P. J.; Chodera, J. D.; Bustamante, C.; Marqusee, S. The Molten Globule State Is Unusually Deformable under Mechanical Force. Proc. Natl. Acad. Sci. 2012 , 109 (10), 3796–3801. (9) Gebhardt, J. C. M.; Bornschlögl, T.; Rief, M. Full Distance-Resolved Folding Energy Landscape of One Single Protein Molecule. Proc. Natl. Acad. Sci. U. S. A. 2010 , 107 (5), 2013–2018. (10) Fukuma, T. Subnanometer-Resolution Frequency Modulation Atomic Force Microscopy in Liquid for Biological Applications. Jpn. J. Appl. Phys. 2009 , 48, 08JA01. (11) Hiasa, T.; Kimura, K.; Onishi, H.; Ohta, M.; Watanabe, K.; Kokawa, R.; Oyabu, N.; Kobayashi, K.; Yamada, H. Aqueous Solution Structure over α-Al2O3(0112) Probed by Frequency-Modulated Atomic Force Microscopy. J. Phys. Chem. C 2010 , 114, 21423–21426. (12) Fukuma, T.; Onishi, K.; Kobayashi, N.; Matsuki, A.; Asakawa, H. Atomic-Resolution Imaging in Liquid by Frequency Modulation Atomic Force Microscopy Using Small Cantilevers with Megahertz-Order Resonance Frequencies. Nanotechnology 2012 , 23 (13), 135706. (13) Hiasa, T.; Kimura, K.; Onishi, H. Minitips in Frequency-Modulation Atomic Force Microscopy at Liquid-Solid Interfaces. Jpn. J. Appl. Phys. 2012 , 51, 025703. (14) Nishioka, R.; Hiasa, T.; Kimura, K.; Onishi, H. Specific Hydration on P -Nitroaniline Crystal Studied by Atomic Force Microscopy. J. Phys. Chem. C 2013 , 117 (6), 2939–2943. (15) Borkovec, M.; Szilagyi, I.; Popa, I.; Finessi, M.; Sinha, P.; Maroni, P.; Papastavrou, G. Investigating Forces between Charged Particles in the Presence of Oppositely Charged Polyelectrolytes with the Multi-Particle Colloidal Probe Technique. Adv. Colloid Interface Sci. 2012 , 179-182, 85–98. (16) Horn, R. G.; Israelachvili, J. N. Direct Measurement of Structural Forces between Two Surfaces in a Nonpolar Liquid. J. Chem. Phys. 1981 , 75 (3), 1400.

23

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(17) Christenson, H. K. Forces between Solid Surfaces in a Binary Mixture of Non-Polar Liquids. Chem. Phys. Lett. 1985 , 118 (5), 455–458. (18) Robert, H.; Gregory, G.; Warr, G.; Atkin, G. Physical Chemistry of Ionic Liquids. Phys. Chem. Chem. Phys. 2010, 12 (8), 1709–1723. (19) Perkin, S. Ionic Liquids in Confined Geometries. Phys. Chem. Chem. Phys. 2012 , 14 (15), 5052. (20) Richetti, P.; Kekicheff, P. Direct Measurement of Depletion and Structural Force in a Micellar System. Phys. Rev. Lett. 1992 , 68 (12), 1951–1954. (21) Kimura, K.; Ido, S.; Oyabu, N.; Kobayashi, K.; Hirata, Y.; Imai, T.; Yamada, H. Visualizing Water Molecule Distribution by Atomic Force Microscopy. J. Chem. Phys. 2010 , 132 (19), 194705. (22) Watkins, M.; Reischl, B. A Simple Approximation for Forces Exerted on an AFM Tip in Liquid. J. Chem. Phys. 2013, 138 (15), 154703. (23) Amano, K.; Suzuki, K.; Fukuma, T.; Takahashi, O.; Onishi, H. The Relationship between Local Liquid Density and Force Applied on a Tip of Atomic Force Microscope: A Theoretical Analysis for Simple Liquids. J. Chem. Phys. 2013 , 139 (22), 224710. (24) Amano, K.; Liang, Y.; Miyazawa, K.; Kobayashi, K.; Hashimoto, K.; Fukami, K.; Nishi, N.; Sakka, T.; Onishi, H.; Fukuma, T. Number Density Distribution of Solvent on a Substrate: A Transform Theory for Atomic Force Microscopy. Phys. Chem. Chem. Phys. 2016 , 18 (23), 15534–15544. (25) Amano, K.; Takahashi, O. Simple Transform Methods of a Force Curve Obtained by Surface Force Apparatus to the Density Distribution of a Liquid near a Surface. Physica A 2015 , 425, 79–89. (26) Amano, K.; Tanaka, E.; Kobayashi, K.; Onishi, H.; Nishi, N.; Sakka, T. Force Measurement Reveals Structure of a Confined Liquid: Observation of the Impenetrable Space. Surf. Sci. 2015, 641, 242–246. (27) Roth, R.; Evans, R.; Lang, A.; Kahl, G. Fundamental Measure Theory for Hard-Sphere Mixtures Revisited: The White Bear Version. J. Phys. Condens. Matter 2002 , 14 (46), 12063. (28) Yu, Y. X.; Wu, J. Structures of Hard-Sphere Fluids from a Modified Fundamental-Measure Theory. J. Chem. Phys. 2002 , 117 (22), 10156–10164. (29) Roth, R. Fundamental Measure Theory for Hard-Sphere Mixtures: A Review. J. Phys. Condens. Matter 2010 , 22 (6), 063102. (30) Choudhury, N.; Ghosh, S. K. Integral Equation Theory of Lennard-Jones Fluids: A Modified Verlet Bridge Function Approach. J. Chem. Phys. 2002 , 116 (19), 8517– 8522. (31) Verlet, L. Integral Equations for Classical Fluids: I. The Hard Sphere Case. Mol. Phys. 1980 , 41, 183–190.

24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(32) Duh, D.-M.; Henderson, D. Integral Equation Theory for Lennard-Jones Fluids: The Bridge Function and Applications to Pure Fluids and Mixtures. J. Chem. Phys. 1996 , 104 (17), 6742. (33) Piech, M.; Walz, J. Y. Direct Measurement of Depletion and Structural Forces in Polydisperse, Charged Systems. J. Colloid Interface Sci. 2002 , 253 (1), 117–129. (34) Zeng, Y.; Von Klitzing, R. Oscillatory Forces of Nanoparticle Suspensions Confined between Rough Surfaces Modified with Polyelectrolytes via the Layer-by-Layer Technique. Langmuir 2012 , 28 (15), 6313–6321. (35) Zeng, Y.; Schön, S.; von Klitzing, R. Silica Nanoparticle Suspensions under Confinement of Thin Liquid Films. J. Colloid Interface Sci. 2015, 449, 522–529. (36) Derjaguin, B. V.; Rabinovich, Y. I.; Churaev, N. V. Measurement of Forces of Molecular Attraction of Crossed Fibres as a Function of Width of Air Gap. Nature 1977, 265 (5594), 520–521. (37) Oettel, M. Depletion Force between Two Large Spheres Suspended in a Bath of Small Spheres: Onset of the Derjaguin Limit. Phys. Rev. E 2004 , 69, 041404. (38) Rentsch, S.; Pericet-Camara, R.; Papastavrou, G.; Borkovec, M. Probing the Validity of the Derjaguin Approximation for Heterogeneous Colloidal Particles. Phys. Chem. Chem. Phys. 2006, 8 (21), 2531–2538. (39) Kirkwood, J. G.; Boggs, E. M. The Radial Distribution Function in Liquids. J. Chem. Phys. 1942 , 10 (6), 394–402. (40) Karino, Y.; Akiyama, R.; Kinoshita, M. Three-Dimensional Density Profiles of Small and Medium Spheres near a Pair of Large Spheres: Relevance to Entropie Interaction Induced between Large Spheres. J. Phys. Soc. Japan 2009 , 78 (4), 044801. (41) Kubota, Y.; Akiyama, R. Assessment of Three-Dimensional Distribution Functions and Triplet Distribution Functions for Hard Spheres Calculated Using a Three-Dimensional Ornstein-Zernike Equation with Hypernetted-Chain Closure. J. Phys. Soc. Japan 2012 , 81, SA017. (42) Attard, P. Spherically Inhomogeneous Fluids. II. Hard-Sphere Solute in a Hard-Sphere Solvent. J. Chem. Phys. 1989, 91 (5), 3083. (43) Bhattacharjee, S.; Elimelech, M. Surface Element Integration: A Novel Technique for Evaluation of DLVO Interaction between a Particle and a Flat Plate. J. Colloid Interface Sci. 1997 , 193 (2), 273–285. (44) Bhattacharjee, S.; Elimelech, M.; Borkovec, M. DLVO Interaction between Colloidal Particles: Beyond Derjaguin’s Approximation. Croat. Chem. Acta 1998 , 71 (4), 883–903. (45) Amano, K.; Kinoshita, M. Entropic Insertion of a Big Sphere into a Cylindrical Vessel. Chem. Phys. Lett. 2010 , 488 (1-3), 1–6. (46) Amano, K.; Yoshidome, T.; Iwaki, M.; Suzuki, M.; Kinoshita, M. Entropic Potential Field Formed for a Linear-Motor Protein near a Filament:

25

ACS Paragon Plus Environment

Langmuir

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Statistical-Mechanical Analyses Using Simple Models. J. Chem. Phys. 2010 , 133 (4), 045103. (47) Amano, K.; Oshima, H.; Kinoshita, M. Potential of Mean Force between a Large Solute and a Biomolecular Complex: A Model Analysis on Protein Flux through Chaperonin System. J. Chem. Phys. 2011 , 135, 185101. (48) Asakura, S.; Oosawa, F. On Interaction Between Two Bodies Immersed In a Solution of Macromolecules. J. Chem. Phys. 1954 , 22 (1954), 1255–1256. (49) Asakura, S.; Oosawa, F. Interaction between Particles Suspended in Solutions of Macromolecules. J. Polym. Sci. 1958 , 33 (126), 183–192. (50) Trokhymchuk, A.; Henderson, D.; Nikolov, A.; Wasan, D. T. Simple Calculation of Structural and Depletion Forces for Fluids/Suspensions Confined in a Film. Langmuir 2001 , 17, 4940–4947. (51) Hayes, R.; Warr, G. G.; Atkin, R. Physical Chemistry of Ionic Liquids. Phys. Chem. Chem. Phys. 2010, 12, 1709–1723. (52) Fries, P. H.; Patey, G. N. The Solution of the Hypernettedchain Approximation for Fluids of Nonspherical Particles. A General Method with Application to Dipolar Hard Spheres. J. Chem. Phys. 1985, 82, 429–440. (53) Kusalik, P. G.; Patey, G. N. The Solution of the Reference Hypernetted-Chain Approximation for Water-like Models. Mol. Phys. 1988 , 65 (5), 1105–1119. (54) Kovalenko, A.; Hirata, F. Three-Dimensional Density Profiles of Water in Contact with a Solute of Arbitrary Shape: A RISM Approach. Chem. Phys. Lett. 1998, 290 (1-3), 237–244. (55) Matubayasi, N.; Nakahara, M. Super- and Subcritical Hydration of Nonpolar Solutes. I. Thermodynamics of Hydration. J. Chem. Phys. 2000 , 112 (18), 8089. (56) Sato, H.; Hirata, F.; Sakaki, S. Distortion of Electronic Structure in Solvated Molecules: Tautomeric Equilibrium of 2-Pyridone and 2-Hydroxypridine in Water Studied by the RISM-SCF Method. J. Phys. Chem. A 2004, 108 (11), 2097–2102. (57) Sato, H.; Sakaki, S. Comparison of Electronic Structure Theories for Solvated Molecules: RISM-SCF versus PCM. J. Phys. Chem. A 2004 , 108 (9), 1629–1634. (58) Kinoshita, M.; Yoshidome, T. Molecular Origin of the Negative Heat Capacity of Hydrophilic Hydration. J. Chem. Phys. 2009 , 130, 144705.

26

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Graphical abstract (TOC graphic)

27

ACS Paragon Plus Environment