Characterization of the Stratified Morphology of Nanoparticle

Aug 22, 2016 - Nanoparticles surrounded by gas agglomerate in a hierarchical fashion. From production until powder processing in the gas phase, ...
0 downloads 0 Views 2MB Size
Subscriber access provided by Northern Illinois University

Article

Characterization of the Stratified Morphology of Nanoparticle Agglomerates Andrea Fabre, Teun Steur, Wim G. Bouwman, Michiel Thijs Kreutzer, and J. Ruud van Ommen J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b07437 • Publication Date (Web): 22 Aug 2016 Downloaded from http://pubs.acs.org on August 28, 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.

The Journal of Physical Chemistry C 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 41

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

The Journal of Physical Chemistry

Characterization of the Stratified Morphology of Nanoparticle Agglomerates Andrea Fabrea , Teun Steura , Wim G. Bouwmanb , Michiel T. Kreutzera , J. Ruud van Ommena,∗ a

Delft University of Technology, Department of Chemical Engineering, Van der Maasweg 9, 2629 HZ Delft, the Netherlands

b

Delft University of Technology, Department of Radiation, Science and Technology, Mekelweg 15, 2629 JB Delft, The Netherlands Phone: +31 15 278 2133. E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 2 of 41

Abstract Nanoparticles surrounded by gas agglomerate in a hierarchical fashion. From production until powder processing in the gas phase, nanoparticles go from individual particles to aggregates, simple agglomerates, and complex agglomerates. Even though the structures at each level have unique properties, they are commonly assessed as a whole. Additionally, the effect of external factors during gas processing on the morphology of these structures is not well understood, and challenging to study due to the limited techniques for in situ analysis of the dynamic phenomenon. Here, we study three materials in their hydrophobic and hydrophilic version. We describe the structural characteristics of each hierarchical level of complex agglomerate formation obtained from two in situ techniques. The first scale, namely aggregates, are open structures with a fractal dimension of about 1.5, which then form simple agglomerates with a fractal dimension close to 3, that later cluster into complex agglomerates that present a fractal dimension of about 2. Furthermore, gas dynamics were found to densify the simple agglomerates increasing their fractal dimension by more than 0.1.

2

ACS Paragon Plus Environment

Page 3 of 41

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

The Journal of Physical Chemistry

Introduction Nanoparticle agglomeration is of tremendous importance in medicine, energy storage, catalysis, consumer products, etc 1 . This process develops in a stepwise fashion, where the structures formed at each step have unique properties and thus, react differently to the mechanical interactions they are exposed to 2–4 . This paper presents underlying morphological differences among the nanoparticle clusters generated throughout the agglomeration process, and the effect of external factors during powder processing; more specifically, during nanopowder suspension in an upward gas flow (i.e. fluidization). Nanoparticles are of great interest in academia and industry due to their unique properties which make them suitable for applications such as in textiles, electronics, medicine, and renewable energy 3,5–11 . Some of these applications require the adaptation of the nanoparticles to the desired function 12 . Due to its potential scalability and large gas-solid contact area, fluidization is commonly used for nanoparticle processing 2,8,13–17 . Nonetheless, particles at this scale fluidize as clusters due to the relatively strong adhesion forces 3,18–25 , a behavior that could hamper characteristics of the individual particles. Nanoparticle agglomeration is a hierarchical process 2,3 . During particle production at high temperatures, nanoparticles come into contact and form chain-like structures of a few hundreds of nanometers called aggregates, connected by solid necks and assumed to be very hard, if not impossible to break 4,26,27 . While stored, these aggregates form highly porous clusters named simple agglomerates, which can reach tens of micrometers in size. The stored nanopowder is then

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 4 of 41

processed in a fluidized bed, where it creates complex agglomerates from the simple agglomerates already present. As a stepwise growth, each of the three structures have unique properties. The nanoparticle cluster is commonly characterized by its fractal-like structure 3,18,28–33 . The Np number of monodisperse particles with size dp forming a mass fractal cluster of size da scales to the power of the ratio da /dp according to

da dp

Df

da dp

Df −3

 Np = kn

(1)

which in terms of density translates to

 ρa = ρp kn

(2)

where ρa is the cluster density,ρp is the particle density, kn is a prefactor, and Df is the mass fractal dimension 34,35 . The prefactor is used to give physical meaning to the expression at the limits of the single nanoparticle and extremely large clusters. The mass fractal dimension gives information on the nanoparticle distribution within the structure, which is a key feature to consider when studying cluster morphology. The agglomerate morphology can be affected by external factors such as the fluidizing gas or collisions with other solids 7,18,36 . Studying the structure of the fluidized agglomerates in situ is quite challenging. Visualization of fluidized agglomerates has been achieved before 7,8,13,37,38 . Nonetheless, it should be noted that visual analysis is always difficult since the data consists of 2D images of 3D objects. Moreover, all camera-based techniques are limited by the maximum resolution of the image. Other techniques such as γ-rays adsorbtion densitometry have been used 4

ACS Paragon Plus Environment

Page 5 of 41

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

The Journal of Physical Chemistry

to study the fluidized bed 14 ; however, the detailed morphology of the fluidized structures could not be determined from this method. Previously, de Martin et al. 30 described the morphology of the multiple agglomerate scales of static SiO2 (R812) and TiO2 (P25), and fluidized TiO2 (P25). Titania P25 consists of hydrophilic 21 nm primary particles that fluidize with channels and bubbles. Since only one nanopowder was analyzed in the fluidized state, the applicability of the results to a wider range of powders was open to debate. Since the inherent morphology of the agglomerate depends on the properties of the primary particles, which include material and surface characteristics such as the presence of hydroxide groups 3,25 , data analysis of a single powder sample suggested interpretations that required further research. Here, an in situ study of the morphology of the structures at the different agglomeration levels of six metal oxide nanopowders is presented. We compare the static and fluidized states to find the scale at which the fluidization dynamics modify the morphology of the agglomerates, which presumably should not reach the strong aggregates. Moreover, a preliminary comparison between hydrophobic and hydrophilic powders is made since the presence of hydroxide groups can strengthen the clusters via hydrogen bonds; hence, enhancing their resistance to morphological changes. This study determines the static and dynamic fractal dimension of nanoparticle aggregates, simple agglomerates, and complex agglomerates using the spin-echo small-angle neutron scattering (SESANS) 30,39,40 technique for the first two scales, and the settling tube technique 35,41 for the third one.

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 6 of 41

Methods Materials Hydrophilic (P) and hydrophobic (A) versions of TiO2 , Al2 O3 , SiO2 commercial nanopwders from Evonik are used in this work; the same set used in previous studies 25,42 . The specifications given by the provider are shown in Table 1. All six powders were fluidized and studied with the settling tube technique. However, fluidization experiments with SESANS were only performed for titania and alumina since silica adhered to the walls, giving unreliable measurements. Static SESANS measurements were done with all three hydrophilic powders, plus hydrophobic titania and alumina. Table 1: Properties of the fluidized nanopowders. The particle density and size are provided by the manufacturer.

Powder Commercial name TiO2 -P P25 TiO2 -A T805 Al2 O3 -P AluC Al2 O3 -A AluC805 SiO2 -P A130 SiO2 -A R972

Surface ρp (kg/m3 ) dp (nm) hydrophilic 4000 21 hydrophobic 4000 21 hydrophilic 3600 13 hydrophobic 3600 13 hydrophilic 2200 16 hydrophobic 2200 16

Previous to fluidization, all powders were sieved to remove agglomerates larger than 400 µm in size to improve fluidization 8 . The powders were fluidized within the next 24 hours after sieving to prevent further formation of large clusters. The measurements with SESANS and the settling tube were made once the fluidized bed had reached a pseudo-steady state with minimal height fluctuations.

6

ACS Paragon Plus Environment

Page 7 of 41

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

The Journal of Physical Chemistry

SESANS The properties of the nanoparticle aggregates and simple agglomerates were studied in situ by the spin-echo small-angle neutron scattering (SESANS) technique which covers the size range from 5nm to 20µm 39 . The time required for each experiment varied between 1 and 5 hours. The measurements are obtained from the polarization of a neutron beam which passes through a precession device with tilted interface before crossing the sample, and through a second precession device with opposite magnetic field direction after the sample, before reaching the detector 30,39 . The setup for SESANS consists of a fluidization column made by a single piece of quartz with a 30 × 25 mm rectangular cross-section. The rectangular cross-section was chosen to ensure that the whole neutron beam would travel through the sample exactly the same distance (Fig.1). The area of the beam is 16 × 10 mm with a wavelength of 0.203±0.001 nm. The size of the column was determined from preliminary measurements on neutron depolarization of the different nanopowders using cuvettes. The distance traveled by the beam through the sample affects the level of neutron depolarization, and the optimum sample depth depends on the material properties; thus, varies among the powders used. From the cuvette experiments it was determined that all powders required a sample thickness of 30 mm, with the exception of hydrophobic TiO2 -A, which needed 25 mm.

Figure 1: Schematic of the SESANS setup. Main parts are labeled with the neutron beam going from left to right, passing through the sample placed at the middle of the beam trajectory.

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 8 of 41

The column had distributor plates at the top and bottom; the bottom one to spread the fluidizing gas, and the top one to prevent particle loss. After sieving, the powders were placed inside the column to a height of 6±1 cm, enough to guarantee full sample coverage of the neutron beam. The powders were fluidized at 4.5 cm/s-8.3 cm/s with (99.999% pure) N2 gas for about 5 hours. Measurements were taken before, during, and in some cases, after fluidization for TiO2 -P, TiO2 -A, Al2 O3 -P, and Al2 O3 -A powders. Settling Tube The settling tube setup comprises a 45×45 mm cross-sectional area quartz column, 15 cm high, with a distributor plate at the bottom, and a bubbler and HEPA filter attached to the top to trap any entrained nanopowder. One of the column’s side walls has two openings for the connections of the settling tube. This settling tube is a black box with dimensions 10×6×30 mm, an opening at the top to catch the falling agglomerates, and two side openings for the borescope-high-speed camera system and the flushing gas that empties the box before each experiment (Fig.2). A rigid Olympus (R040-021-000-60 S5) borescope and a high-speed Phantom v9.1 camera are used for in situ fluidized agglomerate visualization 35,41 in the range of 30 µm to 1 mm.

Figure 2: Settling tube with the agglomerate inlet at the top, and the borescope and flushing gas inlet on the side. A sample image of the recorded fluidized agglomerates is shown.

8

ACS Paragon Plus Environment

Page 9 of 41

All six nanopowders were analyzed with the settling tube technique. The nanopowders were fluidized with N2 at 0.13 m/s after being sieved to remove structures larger than 400 µm. The recordings were taken 15 minutes after fluidization, once the bed height was steady.

Data Analysis The data obtained from SESANS is the polarization (P ) as a function of the spin-echo length (z) (Fig.3). The spin-echo length covers the range of 5 nm-20 µm. The relation between polarization and spin-echo length is characterized by the correlation function, which describes the density of the sample. 1

1

Al2O3-P static

0.9

0.8 P/P0

0.8 0.7

0.7

0.6

0.6

Df1=1.1 ± 0.7 Df2=2.55 ± 0.01

0.5 0.4

TiO2-P fluidized

0.9

P/P0

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

The Journal of Physical Chemistry

rc1=0.02 ± 0.01 [µm]

0

2

4

6

Df1=1.7 ± 0.3

0.5 8 10 z [µm]

12

14

16

0.4

18

Df2=2.87 ± 0.02 rc1=0.13 ± 0.08 [µm]

0

2

4

(a)

6

8 10 z [µm]

12

14

16

18

(b)

Figure 3: Sample plots of static Al2 O3 -P (a) and fluidized TiO2 -P (b) obtained from SESANS. Polarization is normalized by that of an empty container (P0 ). The best fit from equation 8 is shown as a red line. The first and second fractal dimensions (Df 1 and Df 2 , respectively), and transition size (rc1 ) obtained from each plot are presented.

Since SESANS does not give straightforward information about the sample properties, thorough data processing and qualitative understanding of the sample morphology are necessary.The solid fraction of the fluidized samples was determined from the voidage of the static powder, and transmission values from the empty column and fluidized bed from the following equation: 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

T = AeB(1−φf )

Page 10 of 41

(3)

where T is the SESANS transmission of the sample, A is the transmission of the empty column, B is the voidage of the static samples estimated as B = 1 − ρb /ρp , and φf is he solid fraction of the fluidized powder. To relate the experimental data to the sample properties, a theoretical polarization based on known fluidization conditions and powder properties is evaluated, compared, and reiterated until agreeing with the experimental polarization. The estimation of the theoretical polarization requires knowledge on the density-density correlation function for the sample, the projection of the densitydensity correlation function in the form of an Abel transform, and the solid fraction of the sample. A detailed description of the SESANS data analysis can be found in the publications by de Martin et al. 30,43 and Andersson et al. 40,44 . The density-density correlation function for a bifractal system is described as:

 γ∝

γ∝

r a

r +1 rp

+1

Df,1 −3 r < rc1

Df,2 −3

r

e− ξ

r > rc1

(4)

(5)

Where rc1 is a characteristic length scale indicating the transition from one density-density correlation function to the other, and not a radius. The correlation length of the density distribution ξ is equal to the width of the distribution, which means that: 10

ACS Paragon Plus Environment

Page 11 of 41

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

The Journal of Physical Chemistry



Z

γ(r)dr

ξ=2

(6)

z

By an Abel transform, the projection of the density-density correlation function reads as follows:

2 G(z) = ξ



Z z

γ(r)r √ dr r2 − z 2

(7)

To analyze the data extracted from SESANS, a theoretical polarization of the neutrons is estimated and compared to the experimental values for optimization of the initial parameters (Df 1 , Df 2 , and rc1 ). The expression for the theoretical polarization is:

P

P (z) = e

t (G(z)−1)

(8)

where Σt depends on the sample thickness t, neutron wavelength λ, neutron scattering length density difference in the sample ∆ρ0 , correlation length of the density distribution ξ, solid fraction φ1 , and void fraction φ2 as shown below:

Σt = tλ2 ∆ρ20 φ1 φ2 ξ

(9)

Once all the values are calculated, the error difference between the theoretical and experimental polarizations is evaluated. If the difference surpasses the threshold determined by the experimental uncertainties, the loop of calculation is restarted with new values. In this manner, the fractal

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 12 of 41

dimension and transition size for the first and second scales can be obtained from the SESANS measurements. Video analysis of the settling tube experiments is done using in-house MATLAB scripts 41 . The videos have a resolution of 5.6µm per pixel, giving a minimum recognizable size of ≈12µm. A background frame per video is created after MATLAB evaluates all the recorded frames, and then subtracts it to remove the noise. The group of background-less frames are looped to recognize all the agglomerates, track them, and filter the ones that attach to the window or get out of focus. Once a list of labeled agglomerates is made, with their respective frame location, size, and sphericity, the settling velocity Ut and gyration diameter da are estimated assuming that the agglomerates follow Stokes Law 35,41 . The gyration diameter is estimated as 2/3 the outer diameter measured from the video which holds for fractal aggregates in the Stokes regime 34 . Since the average sphericity δ of the agglomerates is about 0.9, they are assumed to be spheres, thus, the gyration diameter da , area Aa , and volume Va are directly obtained from the videos. With this data, the Reynold’s number Re is evaluated, which then is used to calculate the drag coefficient CD , completing all the parameters needed to determine the agglomerate density ρa as follows:

ρa =

Ut2 Aa CD ρf + ρf g2Va

(10)

The logarithmic form of Eq. 2 gives a line which depends on the size and density of the nanoparticle and agglomerate. The agglomerate size and density are obtained from the videos, while those of the primary particles are given by the provider. Then, the prefactor kˆn3 and fractal dimension Df 3 of the recorded fluidized agglomerates is calculated. 12

ACS Paragon Plus Environment

Page 13 of 41

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

The Journal of Physical Chemistry

Results and Discussion SESANS SESANS was used to study the structure of the first two formation levels of complex fluidized agglomerates, which we call aggregates and simple agglomerates. The varying parameters were the material, hydrophobicity, and dynamic state of the powder. Aggregates are known to be connected by solid necks that can withstand any external forces 2 , which means that the transition length scale rc1 obtained for the static powder is fixed for the fluidized samples. The stable aggregates are characterized by Df 1 (Fig.4(a)), which should then be the same for the static and fluidized samples. Hence, Df 1 was left as a fitting parameter for the fluidized system to confirm the value obtained for the static sample; indeed, we found the same Df 1 with experimental accuracy for the static and fluidized systems. To verify the experimental reproducibility, in addition to Df 1 , rc1 was also kept as a fitting parameter and compared for the static and fluidized samples (Supporting Information). It was found that rc1 and Df 1 were consistent between the static and fluidized samples, resulting in a first fractal dimension below 2 and an aggregate size below 500 nm for all the materials, hydrophobicity, and dynamic state. Hence, fluidization dynamics have no effect at the aggregate scale. On the other hand, the values obtained for Df 2 allow for the comparison among materials, properties, and dynamic state (Fig.4(b)). There is no significant difference among the materials; however, as seen in Figure 4(b), the hydrophobic version of the nanopowders has a higher fractal dimension than the hydrophilic one, showing a difference slightly higher than the confidence interval. The divergence between the hydrophobic and hydrophilic materials can be attributed 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

to the influence of the OH groups on the interactions among the particles 25,42 . The presence of OH groups affects the formation of H-bonds, which increase the force required to rearrange the particles within the agglomerate, thus making it harder to modify.

Static Fluidized

2.5 2

1.5

2.5 2

1.5

1

0.5 0

3

Second Fractal Dimension (Df2)

3

First Fractal Dimension (Df1)

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

Page 14 of 41

1

0.5

TiO2 – P

TiO2 – A

Al2O3 – P

Al2O3 – A

0

SiO2 – P

TiO2 – P

(a)

TiO2 – A

Al2O3 – P

Al2O3 – A

SiO2 – P

(b)

Figure 4: First (a) and second (b) fractal dimensions of static (solid) and fluidized (light patterned) nanopowders. Values for TiO2 -P, TiO2 -A, Al2 O3 -P, Al2 O3 -A, and SiO2 -P are presented. The error bars are the best fit confidence interval of equation A.5.

It is also shown that, in contrast with the TiO2 -P results of de Martin el al. 30 , the fluidized nanopowders (light bars) have a higher value for Df 2 than the static sample (solid bars) (Fig.4(b)). During fluidization, agglomerates collide, break, coalesce, and deform; a dynamic that has been confirmed by fluidizing nanopowders of three different colors (red, green, and white) and finding agglomerates with a combination of all the colors 13 . The impact of this unsteadiness on the morphology of the agglomerate was unknown until now. The evident divergence on the second fractal dimension of the static and fluidized samples suggests that fluidization densifies the structures in the 1-20 µm size range. As the simple agglomerates collide in the fluidized bed, numerous outcomes are possible depending on the relative collision velocity and angle, and the properties of the colliding structures. Some collisions have enough energy to compress the agglomerates by rearranging the aggregates within them, and leading to a bounce or coalesce of the colliding clus14

ACS Paragon Plus Environment

Page 15 of 41

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

The Journal of Physical Chemistry

ters. Collisions with less energy most likely result in the coalescence of the simple agglomerates, loosely joined by interpaticle forces among protruding aggregates, thus creating the complex agglomerates with very low density known to be in the fluidized bed and observed with the settling tube.

Settling Tube From the recordings, the size and settling velocity of the agglomerates can be extracted, and used to estimate the fractal dimension and prefactor when combined with the primary particle properties (Supporting Information). The size range of the recorded agglomerates is 50-980 µm, which agrees with the size found in literature of 70-900 µm 7,8,13,32,37 . The prefactor and fractal dimension are estimated from a linear fit using the logarithm of Eq. 2 to the experimental data. In all cases, except for SiO2 -P, the coefficient of determination was above 0.8, which indicates a considerably good fit to the experimental data. For SiO2 -P, the coefficient was 0.53 (Fig.5 (c)). The sparse cloud of data points below the denser cloud of SiO2 -P could be the agglomerates falling close to the window with an inclined trajectory, at a seemingly lower velocity than their corresponding settling velocity. The prefactor and fractal dimension calculated from the best fit line are shown in Table 2. The error on the prefactor can be quite large since it is obtained from the line intercept on the log(ρa /ρp ) axis and the error propagation from the uncertainty on the agglomerate density results in the given values.

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 16 of 41

(a)

(b)

(c) Figure 5: The plots show the data points of the complex agglomerates as recorded by the settling tube, with the best fit line from which the prefactor and fractal dimension was extracted. Data from TiO2 -P (left) and TiO2 -A (right) (a); Al2 O3 -P (left) and Al2 O3 -A (right) (b), and SiO2 -P (left) and SiO2 -A (right) (c).

16

ACS Paragon Plus Environment

Page 17 of 41

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

The Journal of Physical Chemistry

Table 2: Prefactor and fractal dimension of complex nanoparticle agglomerates estimated from the number of data points indicated on the last column.

Powder TiO2 -P TiO2 -A Al2 O3 -P Al2 O3 -A SiO2 -P SiO2 -A

kˆn3 (-) Df 3 (-) Data Points 51 ± 3 2.12 ± 0.01 3276 150 ± 6 2.04 ± 0.01 8880 367 ± 76 1.89 ± 0.02 295 584 ± 51 1.92 ± 0.01 3483 38 ± 11 2.19 ± 0.03 683 747 ± 117 1.90 ± 0.02 819

There is no trend among the hydrophobicity and material of the different powders. TiO2 -P and SiO2 -P have a higher fractal dimension than TiO2 -A and SiO2 -A; whereas for Al2 O3 , the opposite seems to be true. Additionally, the pattern on the values according to the material is not clear given the large discrepancy between the SiO2 powders. Nevertheless, it is clear that the structure of the agglomerates is rather dilute at this scale.

Combining all three scales The results from both techniques, SESANS and settling tube, were merged to give an overall picture of the fluidized agglomerate structure (Fig.6). An estimate of the second scale prefactor was obtained by equating Eq. 2 for the first and second scales with the fractal dimensions and transition size determined from SESANS. Likewise, the size at the intersection between the second and third scales was estimated from Eq. 2 using the previously determined fractal dimensions and prefactors (Supporting Information). A compilation of all three fractal dimensions for each powder is shown in Table 3.

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 18 of 41

Table 3: Fractal dimension of all three hierarchical levels of the agglomeration process during fluidization.

Powder TiO2 -P TiO2 -A Al2 O3 -P Al2 O3 -A SiO2 -P SiO2 -A

Df 1 (-) 1.9 ± 0.2 1.7 ± 0.2 1.1 ± 0.3 1.6 ± 0.2 1.6 ± 0.2

Df 2 (-) 2.76 ± 0.04 2.98 ± 0.02 2.65 ± 0.03 2.74 ± 0.04 2.64 ± 0.09

Df 3 (-) 2.12 ± 0.01 2.04 ± 0.01 1.89 ± 0.02 1.92 ± 0.01 2.19 ± 0.03 1.90 ± 0.02

The superposition of the data points from the settling tube to the continuous line estimated from the values mentioned before shows an extension of the settling tube data points over the predicted line for the simple agglomerates (Fig.6(a)(c)(d)). This misfit could arise from the large uncertainty of the first fractal dimension and the transition size obtained from SESANS (Supporting Information). A smoother continuous prediction of the size to density correlation could be made with a smaller transition size, or a higher first fractal dimension. Therefore, further studies on the first scale will confirm the adequate adjustment for the precise values to predict the full agglomeration range.

18

ACS Paragon Plus Environment

Page 19 of 41

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

The Journal of Physical Chemistry

(a)

(b)

(c)

(d)

Figure 6: Relation between density and size for all three agglomeration scales. a) TiO2 -P, b) TiO2 -A, c) Al2 O3 -P, d) Al2 O3 -A

As mentioned before, Df 1 and rc1 should be the same for static and fluidized powders since aggregates are connected by solid necks, resistant to any external force. On the other hand, the increase of Df 2 on the fluidized samples fits the results from the settling tube. If the second scale was kept the same for the static and fluidized powders, the connection between the simple and complex fluidized agglomerates would not be possible (Fig.7, dotted blue line). The static powder does not make individual complex agglomerates since the whole batch of powder is a network without a third agglomeration scale and an overall powder density. The tapped density of TiO2 -P given by the supplier is 100-180 kg/m3 corresponding to the size range 8-30 µm. Therefore, this is the size expected for nanoparticle agglomerates formed during storage.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 20 of 41

Figure 7: Plot of hydrophilic TiO2 -P comparing the static (dotted blue line) and fluidized (continuous red line) predicted size-to-density correlation.

Conclusions In this paper, we employ in situ techniques to analyze the morphology of the three structural levels, namely aggregates, simple agglomerates, and complex agglomerates, for commercial metal-oxide nanopowders. The first two scales are studied by spin-echo small-angle neutron scattering (SESANS), while the third one involves a settling tube with a borescope-high-speed camera system. The results show that the first fractal dimension for the aggregate structure is less than 2 for all the powders. The second fractal dimension lays above 2.6 with a maximum at 3, and the third one decreases back to values between 1.8 and 2.2. The data obtained from both methods was combined to understand the correlation between size and density throughout the nanopowder agglomeration size scales in a fluidized bed. Hence, the conclusions on the fractal dimension of the specific samples studied by de Martin et al. 30 are rather general to fluidization of commercial nanopowder. It was found that fluidization affects the second scale by densifying the structure. The fractal dimension calculated from SESANS gave higher values for the fluidized samples, which according to Eq. 2, results in a denser structure. This outcome was supported by the experiments done 20

ACS Paragon Plus Environment

Page 21 of 41

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

The Journal of Physical Chemistry

with the settling tube, which required an increase on the second fractal dimension to have a continuous correlation among all three scales in the fluidized state. Therefore, it can be concluded that fluidization dynamics affect the second scale, creating denser structures that later form large complex agglomerates.

Supporting Information Confirmation of the aggregate stability by setting the transition size (first-second scales) and first fractal dimension as fitting parameters. Equations used to evaluate the second prefactor and transition size between simple and complex agglomerates. This material is available free of charge via the Internet at http://pubs.acs.org/.

Acknowledgments The research leading to this paper has received funding from the European Research Council under the European Unions Seventh Framework Programme (FP/2007-2013)/ERC Grant, agreement no. 279632. We also want to acknowledge Chris Duif for the help with the SESANS measurements. [1] Sokolov, S.V.; K¨atelh¨on, E.; Compton, R.G. A thermodynamic view of agglomeration. J. Phys. Chem. C . 2015. 119, 25093–25099. [2] Yao, W.; Guangsheng, G.; Fei, W.; Jun, W. Fluidization and agglomerate structures of SiO2 nanoparticles. Powder Technol . 2002. 124, 152–159. [3] van Ommen, J.R.: Manufacturing core-shell nanoparticles by atomic layer deposition in a fluidized bed. International Symposium on Chemical Reaction Engineering. 2012 .

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 22 of 41

[4] Teleki, A.; Wengeler, R.; Wengeler, L.; Nirschl, H.; Pratsinis, S. Distinguishing between aggregates and agglomerates of flame-made TiO2 by high-pressure dispersion. Powder Technol . 2008. 181, 292 – 300. [5] Heggen, M.; Oezaslan, M.; Houben, L.; Strasser, P. Formation and analysis of core-shell fine structures in Pt bimetallic nanoparticle fuel cell electrocatalysts. J. Phys.Chem. C . 2012. 116, 19073–19083. [6] Biswas, N.; Bhattacharya, R.; Saha, A.; Jana, N.R.; Basu, J.K. Interplay of electrostatics and lipid packing determines the binding of charged polymer coated nanoparticles to model membranes. Phys. Chem. Chem. Phys. . 2015. 17, 24238–24247. [7] Valverde, J. M.and Quintanilla, M.A.S.; Castellanos, A.; Lepek, D.; Quevedo, J.; Dave, R.N.; Pfeffer, R. Fluidization of fine and ultrafine particles using nitrogen and neon as fluidizing gases. AIChE Journal . 2008. 54, 86–103. [8] Zhu, C.; Yu, Q.; Dave, R.N.; Pfeffer, R. Gas fluidization characteristics of nanoparticle agglomerates. AIChE Journal . 2005. 51, 426–439. [9] Heikkil¨a, E.; Martinez-Seara, H.; Gurtovenko, A.A.; Javanainen, M.; H¨akkinen, H.; Vattulainen, I.; Akola, J. Cationic Au nanoparticle binding with plasma membrane-like lipid bilayers: Potential mechanism for spontaneous permeation to cells revealed by atomistic simulations. J. Phys. Chem. C . 2014. 118, 11131–11141. [10] Vaidya, T.S. Fluidization behavior of alumina nano-particles. App. Mech. and Mat. . 2011. 110-116, 1833–1840. [11] Hu, P.; Zhang, X.; Zhang, C.; Chen, Z. Molecular interactions between gold nanoparticles and model cell membranes. Phys. Chem. Chem. Phys. . 2015. 17, 9873–9884. [12] Suhendi, A.; Nandiyanto, A.B.D.; Ogi, T.; Okuyama, K. Agglomeration-free core-shell polystyrene/silica particles preparation using an electrospray method and additive-free cationic polystyrene core. Mater. Lett. . 2013. 91, 161–164. [13] Hakim, L.F.; Portman, J.L.; Casper, M.D.; Weimer, A.W. Aggregation behavior of nanoparticles in fluidized beds. Powder Technol . 2005. 160, 149–160. [14] Jung, J.; Gidaspow, D. Fluidization of nano-size particles. J. Nanopart. Res . 2002. 4, 483–497. [15] Dadkhah, M.; Peglow, M.; Tsotsas, E. Characterization of the internal morphology of agglomerates produced in a spray fluidized bed by x-ray tomography. Powder Technol . 2012. 228, 349–358.

22

ACS Paragon Plus Environment

Page 23 of 41

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

The Journal of Physical Chemistry

[16] Valverde, J.M.; Pontiga, F.; Soria-Hoyo, C.; Quintanilla, M.A.S.; Moreno, H.; Duran, F.J.; Espin, M.J. Improving the gas-solids contact efficiency in a fluidized bed of CO2 adsorbent fine particles. Phys. Chem. Chem. Phys. . 2011. 13, 14906–14909. [17] Quevedo, J.A.; Omosebi, A.; Pfeffer, R. Fluidization enhancement of agglomerates of metal oxide nanopowders by microjets. AIChE Journal . 2010. 56, 1456–1468. [18] Shabanian, J.; Jafari, R.; Chaouki, J. Fluidization of ultrafine powders. Inter. Review of Chem. Eng. . 2012. 4, 16–50. [19] Parveen, F.; Berruti, F.; Briens, C.; McMillan, J. Effect of fluidized bed particle properties and agglomerate shape on the stability of agglomerates in a fluidized bed. Powder Technol . 2013. 237, 46–52. [20] Zhou, T.; Li, H. Estimation of agglomerate size for cohesive particles during fluidization. Powder Technol . 1999. 101, 57–62. [21] Khadilkar, A.; Rozelle, P.L.; Pisupati, S.V. Models of agglomerate growth in fluidized bed reactors: Critical review, status and applications. Powder Technol . 2014. 264, 216 – 228. [22] Liang, Y.; Ozawa, M.; Krueger, A. A general procedure to functionalize agglomerating nanoparticles demonstrated on nanodiamond. ACS Nano . 2009. 3, 2288–2296. [23] Salameh, S.; Schneider, J.; Laube, J. Adhesion mechanisms of the contact interface of TiO2 nanoparticles in films and aggregates. Langmuir . 2012. 28, 11457–11464. [24] Laube, J.; Salameh, S.; Kappl, M.; Madler, L.; Ciacchi, L.C. Contact forces between TiO2 nanoparticles governed by an interplay of adsorbed water layers and roughness. Langmuir . 2015. 31, 11288–11295. [25] Tahmasebpoor, M.; de Martin, L.; Talebi, M.; Mostoufi, N.; van Ommen, J.R. The role of the hydrogen bond in dense nanoparticle-gas suspensions. Phys. Chem. Chem. Phys. . 2013. 15, 5788–5793. [26] Seipenbusch, M.; Rothenbacher, S.; Kirchhoff, M.; Schmid, H.J.; Kasper, G.; Weber, A.P. Interparticle forces in silica nanoparticle agglomerates. J. Nanopart. Res . 2010. 12, 2037–2044. [27] Yan, W.; Li, S.; Zhang, Y.; Yao, Q.; Tse, S.D. Effects of dipole moment and temperature on the interaction dynamics of titania nanoparticles during agglomeration. J. Phys. Chem. C . 2010. 114, 10755–10760. [28] Vicsek, T.: Fractal Growth Phenomena. World Scientific Pub Co Inc. 1992.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

Page 24 of 41

[29] Friedlander, S.K.: Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics. Oxford University Press. 2000. [30] de Martin, L.; Bouwman, W.; van Ommen, J.R. Multidimensional nature of fluidized nanoparticle agglomerates. Langmuir . 2014. 30, 12696–12702. [31] Hu, D.; Zhuang, J.; Ding, M. A review of studies on the granular agglomeration mechanisms and antiagglomeration methods. Key Eng. Mater. . 2012. 501, 515–519. [32] Nam, C.H.; Pfeffer, R.; Dave, R.N.; Sundaresan, S. Aerated vibrofluidization of silica nanoparticles. AIChE Journal . 2004. 50, 1776–1785. [33] Valverde, J.M.; Castellanos, A. Types of gas fluidization of cohesive granular materials. Phys. Rev. E . 2007. 75. [34] Bushell, G.; Yan, Y.; Woodfield, D.; Raper, J.; Amal, R. On techniques for the measurement of the mass fractal dimension of aggregates. Adv. Colloid Interface Sci. . 2002. 95, 1 – 50. [35] de Martin, L.; Fabre, A.; van Ommen, J.R. The fractal scaling of fluidized nanoparticle agglomerates. Chem. Eng. Sci. . 2014. 112, 79 – 86. [36] Morooka, S.; Kusakabe, K.; Kobata, A.; Kato, Y. Fluidization state of ultrafine powders. J. Chem. Eng. Japan . 1988. 21, 41–46. [37] Wang, X.S.; Palero, V.; Soria, J.; Rhodes, M.J. Laser-based planar imaging of nano-particle fluidization: Part I-determination of aggregate size and shape. Chem. Eng. Sci. . 2006. 61, 5476–5486. [38] Quintanilla, M.A.S.; Valverde, J.M.; Espin, M.J.; Castellanos, A. Electrofluidization of silica nanoparticle agglomerates. Ind. & Eng. Chem. Res. . 2012. 51, 531–538. [39] Rekveldt, M.T.; Plomp, J.; Bouwman, W.G.; Kraan, W.H.; Grigoriev, S.; Blaauw, M. Spin-echo small angle neutron scattering in delft. Revi. Sci. Instrum . 2005. 76, 033901. [40] Andersson, R.; Bouwman, W.G.; Luding, S.; Schepper, I.M. Structure in cohesive powders studied with spinecho small angle neutron scattering. Granul Matter . 2008. 10, 407–414. [41] de Martin, L.; Sanchez-Prieto, J.; Hernandez-Jimenez, F.; van Ommen, J. A settling tube to determine the terminal velocity and size distribution of fluidized nanoparticle agglomerates. J. Nanopart. Res . 2013. 16, 2183.

24

ACS Paragon Plus Environment

Page 25 of 41

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

The Journal of Physical Chemistry

[42] Xanthakis, E.; van Ommen, J.R.; Ahrn´e, L. Flowability characterization of nanopowders. Powder Technol . 2015. 286, 156 – 163. [43] de Mart´ın, L.; Bouwman, W.G.; van Ommen, J.R.: Multidimensionality in fluidized nanopowder agglomerates. Powders and Grains, AIP Conf. Proc. 2013 . [44] Andersson, R.; Bouwman, W.; Plomp, J.; Mulder, F.; Schimmel, H.; Schepper, I.D. Structure, anisotropy and fractals in compressed cohesive powders. Powder Technol . 2009. 189, 6 – 13.

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

26

ACS Paragon Plus Environment

Page 26 of 41

Page 27 of 41

1 2 3 4 5 6 7

The Journal of Physical Chemistry

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 ACS Paragon Plus Environment

Page 28 of 41

Page 29 of 41

The Journal of Physical Chemistry

1

Al2O3-P static

0.9 0.8

P/P0

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

0.7 0.6

Df1=1.1 ± 0.7 Df2=2.55 ± 0.01

0.5 0.4

rc1=0.02 ± 0.01 [µm]

0

2

4

6

8 10 z [µm]

12

14

ACS Paragon Plus Environment

16

18

The Journal of Physical Chemistry

1

TiO2-P fluidized

0.9 0.8

P/P0

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

Page 30 of 41

0.7 0.6

Df1=1.7 ± 0.3

0.5 0.4

Df2=2.87 ± 0.02 rc1=0.13 ± 0.08 [µm]

0

2

4

6

8 10 z [µm]

12

14

ACS Paragon Plus Environment

16

18

Page 31 of 41

3

First Fractal Dimension (Df1)

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

The Journal of Physical Chemistry

Static Fluidized

2.5 2

1.5 1

0.5 0

TiO2 – P

TiO2 – A

Al2O3 – P

ACS Paragon Plus Environment

Al2O3 – A

SiO2 – P

The Journal of Physical Chemistry

3

Second Fractal Dimension (Df2)

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

Page 32 of 41

2.5 2

1.5 1

0.5 0

TiO2 – P

TiO2 – A

Al2O3 – P

ACS Paragon Plus Environment

Al2O3 – A

SiO2 – P

Page 33 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

The Journal of Physical Chemistry

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

ACS Paragon Plus Environment

Page 34 of 41

Page 35 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

The Journal of Physical Chemistry

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14

ACS Paragon Plus Environment

Page 36 of 41

Page 37 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

The Journal of Physical Chemistry

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14

ACS Paragon Plus Environment

Page 38 of 41

Page 39 of 41

1 2 3 4 5 6 7 8 9 10 11 12 13 14

The Journal of Physical Chemistry

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14

ACS Paragon Plus Environment

Page 40 of 41

Page 41 of 41

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

The Journal of Physical Chemistry

ACS Paragon Plus Environment