Simulation of Batch Nanoparticle Growth by the Generalized

diffusion, conservation of mass and Gibbs-Thomson equation.31 The generalized growth rate can. Page 3 of 44. ACS Paragon Plus Environment. The Journal...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF LOUISIANA

C: Physical Processes in Nanomaterials and Nanostructures

Simulation of Batch Nanoparticle Growth by the Generalized Diffusion Model Tianlong Wen, Xiaochen Zhang, Dainan Zhang, Chong Zhang, Qiye Wen, Huaiwu Zhang, and Zhiyong Zhong J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b07438 • Publication Date (Web): 31 Oct 2018 Downloaded from http://pubs.acs.org on November 5, 2018

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 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 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.

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 44 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

Simulation of Batch Nanoparticle Growth by the Generalized Diffusional Model Tianlong Wen*, Xiaochen Zhang, Dainan Zhang, Chong Zhang, Qiye Wen, Huaiwu Zhang, Zhiyong Zhong State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science and Technology of China, Chengdu, Sichuan, China 610054

KEYWORDS: Diffusional growth, monodisperse nanoparticles, batch nanoparticles, nanoparticle growth dynamics, Ostwald ripening, monomer concentration.

ACS Paragon Plus Environment

1

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 44

ABSTRACT. Diffusional growth in solution is the most efficient method to grow monodisperse colloidal nanoparticles. Instead of calculating a single nanoparticle, here the growth of an ensemble of nanoparticles has been numerically simulated by the generalized growth rate of nanoparticles, which is deduced in our previous work from the very fundamental laws and equations. By doing that, the dynamic interaction among nanoparticles through monomer attachment and detachment can be clearly revealed. The constant and diminishing monomer concentration have been simulated separately and compared together. The simulations give the growth dynamics of batch nanoparticles across a broad range of conditions, where diffusivity, adsorption and concentration of monomers play a significant role. The relations between average radius and time, variance and time, and variance and average radius have been obtained for different growth conditions. The average radius-time curve can be well fitted to 𝐸[𝑅] = 𝐸[𝑅0] +𝛼𝑡𝑛 andE[R] = E[𝑅0] +

(𝑅𝑠 ― 𝐸[𝑅0])𝐿𝑎𝑛𝑔𝑒𝑣𝑖𝑛(𝛽𝑡) for constant and diminishing monomer concentration respectively. By the simulation, many hidden details about diffusional growth have been unveiled, such as the details of Ostwald ripening and the role of monomer concentration.

ACS Paragon Plus Environment

2

Page 3 of 44 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

1. Introduction Monodisperse nanoparticles are essential to obtain nanoparticle crystals of long range order 1-7 and to achieve uniform and tunable physical properties of nanoparticle ensembles,8-12 which will deeply affect their many applications in such as display technology,13 lithography,14 flash memory,15 catalysis,16-17, terahertz modulators,18 magnetic nanoparticle imaging and hyperthermia,19-20 and magnetic nanocomposites.20-22 Solution-based synthesis of functional nanoparticles23-26 has many advantages over physical deposition in vacuum27-28 or air29 since the former method has better control over the size, composition, uniformity, stability and selfassembly of the nanoparticles. To foster a mature nanoparticle in solution, it involves two steps: (1) the birth of a nucleus by the aggregation of monomers (atom-surfactant complex) in solution and (2) the growth of the nucleus by consuming the supersaturated monomers in solution.30 To prevent the abortion of the nucleus, their size should be larger than the critical size for nucleation, which is dynamically determined by the monomer concentration.31 LaMer and Dinegar proposed that the nucleation of all particles should occur in a very short period so that they have nearly the same growth history to achieve monodispersity.30Under diffusion-layer assumption, Sugimoto later deduced the growth rate of nanoparticles under two extreme cases of diffusion-controlled and adsorptioncontrolled growth.32 The analytical expression of a single nanoparticle indicates that the size distribution of nanoparticle ensembles can be either focused in diffusion-controlled growth or broadened in adsorption-controlled growth. However, the growth rate under the intermediate conditions is not obtained in Sugimoto’s deduction. Without the diffusion-layer assumption, we deduced a generalized growth rate of spherical nanoparticles from the fundamental Fick’s law of diffusion, conservation of mass and Gibbs-Thomson equation.31 The generalized growth rate can

ACS Paragon Plus Environment

3

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 44

be automatically reduced to the Sugimoto’s expression under the extreme conditions.31 The generalized growth rate allows us to simulate the growth of batch nanoparticles under general growth conditions, which will provide us many clues for the fabrication of monodisperse nanoparticles. Equation 1 shows the generalized growth rate of nanoparticles, where R is the radius of nanoparticles at time t, Ω is the volume occupied by an atom, D is the diffusivity of the monomers, k is the adsorption constant, T is the absolute temperature, Cb is the monomer concentration in the bulk solution, Ceq is the equilibrium monomer concentration, and C∞ is the equilibrium monomer concentration at the flat surface. 31 𝑑𝑅 𝑑𝑡

where A =

=

𝐷𝐴Ω

(1)

𝑅2

𝑘𝑅2(𝐶𝑏 ― 𝐶𝑒𝑞)

(1a)

𝐷 + 𝑘𝑅

( )

and 𝐶𝑒𝑞 = 𝐶∞𝑒𝑥𝑝

2𝛾Ω 𝑅𝑘𝐵𝑇

(1b)

where equation (1a) is valid under quasi-static assumption, or if 4DΩkR(𝐶𝑏 ― 𝐶𝑒𝑞)/(𝐷 + 𝑘𝑅)2 ≪ 1, which is true if monomer concentration in solution is much lower than that in the solids.31 If solving the differential equation (1), the radius of a nanoparticle at any time t, R, can be expressed as 𝑅 = 𝑅(𝑡, 𝑅0,𝐷, 𝑘, 𝐶𝑏,𝑇, Ω, 𝛾, 𝐶∞)

(2)

where R0 is the initial radius of the nanoparticle right after the nucleation event. Namely, the growth of a given nanoparticle is determined by a group of parameters of t, R0, D, k, Cb, T, Ω, γ,

ACS Paragon Plus Environment

4

Page 5 of 44 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

and C∞. Among them, R0 is the nucleation parameter, D and k are the growth parameters, Cb is the precursor parameter, and T is the reaction condition parameters, and Ω, γ and C∞ are the material parameters. Consequently, the final product of nanoparticles is sensitive to the growth time, initial nucleation, solvent, type and concentration of precursors, reaction conditions and chosen materials. Variation of these parameters could lead to discrepancy of the final product. Apparently, it is of challenge to obtain the analytical expression of R in equation (2). In this work, instead of considering only a single nanoparticle, the growth of an ensemble nuclei is numerically simulated using the generalized growth model to unveil the many characteristics of the batch nanoparticle growth. Through the simulation, we can obtain the details of the growth dynamics of the batch nanoparticles across various growth conditions. By the simulation, we expect the growth of batch nanoparticles could be revealed clearly so that it can be finely tuned for better control over their size and uniformity. 2. Simulation Method In the simulation, we assume the synthesis follow the LaMer’s mechanism,30 namely, all nanoparticles nucleate simultaneously at the beginning, and the monomers are quickly consumed by nucleation to reach below the critical supersaturation to prevent secondary nucleation.30 As a result, the size of nanoparticles follows a unimodal distribution. Nanoparticles are subjected to a combination of finite diffusion and a reversible growth reaction with finite kinetics. The growth rate is given by equation (1), and in an ad hoc manner, the dissolution rate is assumed to follow equation (1) as well when Cb < Ceq. Further, during the growth, nanoparticles do not coalesce, aggregate, split, compete locally for monomers, or otherwise interact other than through their effect on the bulk monomer concentration.

ACS Paragon Plus Environment

5

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 44

After the one-time nucleation event, the radius of the nucleus ensemble is assumed to obey the unimodal log-normal distribution as follow,33

𝑓(𝑅0) =

(―

1 2𝜋𝜎𝑅0𝑒𝑥𝑝

(ln (𝑅0) ― 𝜇)2

)

2𝜎2

(

(3)

𝑉𝑎𝑟[𝑅0]

)

μ = 𝑙𝑛 (E[𝑅0]) ―0.5𝑙𝑛 1 + (𝐸[𝑅 ])2

(

𝑉𝑎𝑟[𝑅0]

0

)

𝜎2 = 𝑙𝑛 1 + (𝐸[𝑅 ])2 0

(3a)

(3b)

where E[𝑅0] and Var[𝑅0] are the expectation value and variance of the radius, which determine the parameters µ and σ in log-normal distribution as in equation (3a) and (3b). Among the ensemble, the nuclei of radius R0 𝑅𝐶). Let Nt and N be the total number of nuclei per unit volume after the one-time nucleation and the total number of nuclei per unit volume for growth respectively. Thus N = 𝑁𝑡(1 ― 𝐹(𝑅𝐶)), where F(R) is the cumulative distribution function of the corresponding log-normal distribution function. In this work, the software Matlab was used for numerical simulation of batch nanoparticles, and the numbers of the nucleus having radius R0i for growth in unit volume is digitized as dN(𝑅0𝑖) = 𝑁𝑡𝑓(𝑅0𝑖|𝑅0𝑖 > 𝑅𝐶)𝛿𝑅0

(4)

𝑅0(𝑖 + 1) = 𝑅0𝑖 +𝛿𝑅0

(4a)

ACS Paragon Plus Environment

6

Page 7 of 44 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

where i is positive integer, 𝑅01 = 𝑅𝐶 +𝑂(𝑅)≅𝑅𝐶, 𝛿𝑅0 = 0.02 𝑛𝑚. For a given material (Ω, γ and C∞) and reaction condition (T), Ceq is specified for each radius R0i according to equation 1(b). The size distribution of a specific nanoparticle ensemble at a given time t and constant temperature T is then determined by the growth parameters D, k and the precursor concentration Cb. In this work, two typical situations in the synthesis are considered, namely, (I) monomer is continuously supplied during the nanoparticle growth to keep the monomer concentration constant, and (II) monomer is gradually consumed by the nanoparticle growth without extra addition.

Case I: Monomer is continuously supplied during the growth, and Cb=constant For this case, equation (1b) is substituted into equation (1a), which is subsequently substituted into equation (1) to obtain the growth rate of nanoparticles at specific radius. The time is set at t = 0 second right after the nucleation event. The radius of nanoparticles at time t with initial radius of 𝑅0𝑖 is calculated accumulatively from 𝑅0𝑖 in a small time interval δt by

𝑅0𝑖(𝑡) = 𝑅0𝑖(𝑛𝛿𝑡) = 𝑅0𝑖((𝑛 ― 1)𝛿𝑡) +𝛿𝑡 ×

|

𝑑𝑅 𝑑𝑡 𝑅 = 𝑅 ((𝑛 ― 1)𝛿𝑡) 0𝑖

(5)

where 𝑅0𝑖(0) = 𝑅0𝑖, n is positive integers.

Case II: Monomer is gradually consumed, and Cb is gradually diminishing.

ACS Paragon Plus Environment

7

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 44

For this case, assume the total volume of the solution is Vb, and the total number of monomer in solution is equal to Cb(t)Vb at time t. The radius of a nanoparticle at time t is also calculated by equation (5). However, when calculate dR/dt at time t in equation (5), Cb is time dependent. And

[𝐶𝑏((𝑛 ― 1)𝛿𝑡) ― 𝐶𝑏(𝑛𝛿𝑡)]𝑉𝑏 1 = Ω



∑4𝜋3 [𝑅

3 0𝑖(𝑛𝛿𝑡)

― 𝑅30𝑖((𝑛 ― 1)𝛿𝑡)] ×

dN(𝑅0𝑖) × 𝑉𝑏 =

𝑖=1

1 Ω



∑4𝜋3 [𝑅

3 0𝑖(𝑛𝛿𝑡)

― 𝑅30𝑖((𝑛 ― 1)𝛿𝑡)]𝑁𝑡𝑓(𝑅0𝑖|𝑅0𝑖 > 𝑅𝐶)𝛿𝑅0𝑉𝑏

𝑖=1

=

𝑁𝑡𝑉𝑏𝛿𝑅04𝜋 Ω

3



∑ [𝑅

3 0𝑖(𝑛𝛿𝑡)

― 𝑅30𝑖((𝑛 ― 1)𝛿𝑡)]𝑓(𝑅0𝑖|𝑅0𝑖 > 𝑅𝐶)

𝑖=1

As a result

𝐶𝑏(𝑛𝛿𝑡) = 𝐶𝑏((𝑛 ― 1)𝛿𝑡) ―

4𝜋𝑁𝑡𝛿𝑅0 3Ω



∑𝑖 = 1[𝑅30𝑖(𝑛𝛿𝑡) ― 𝑅30𝑖((𝑛 ― 1)𝛿𝑡)]𝑓(𝑅0𝑖|𝑅

0𝑖

> 𝑅𝐶)

(6) And 𝐶𝑏(0) = 𝐶𝑏0 is the initial monomer concentration in bulk solution. As a result, the monomer concentration in this case is a function of nanoparticle concentration. As growth proceed, the monomer will be gradually consumed. For a given size, when 𝐶𝑏(𝑡) = 𝐶𝑒𝑞, the growth is ceased; when 𝐶𝑏(𝑡) < 𝐶𝑒𝑞, the nanoparticle will shrink according to equation 1. And the dissolved part of the shrinking nanoparticles will become monomers in the bulk solution again, where equation (6) is still valid. For shrinking nanoparticles, when R0i(t) is below the critical size for nucleation at the corresponding 𝐶𝑏(𝑡), the dN(R0i) is set to be zero, resulting in reduction of nanoparticle concentration.

ACS Paragon Plus Environment

8

Page 9 of 44 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

For both cases, the expectation value and variance of the nanoparticle ensemble at time t is given by equation (7a) and (7b) respectively. ∞

E[R(t)] = E[R(nδt)] = ∑𝑖 = 1𝑅0𝑖(𝑡) ×

{



Var[R(t)] = ∑𝑖 = 1𝑅20𝑖(𝑡) ×

𝑑𝑁(𝑅0𝑖) 𝑁

𝑑𝑁(𝑅0𝑖)

(7a)

𝑁

} ― {𝐸[𝑅(𝑡)}

2

(7b)

N and dN(R0i) are dynamically adjusted in the Matlab program if complete dissolving of nanoparticles occurs. To obtain the size distribution function, let F(R, t) be the cumulative function of the number of nanoparticles at size R and time t, where F(R, t) = ∑𝑅

0𝑖(𝑡)

≤𝑅

𝑑𝑁0𝑖

As a result, the size distribution function of the nanoparticle ensemble is given by

𝑓(R,t) =

∂𝐹(𝑅,𝑡) ∂𝑅

(8)

Integration of equation (8) from R=0 to R=∞ will give N, the total number of nanoparticles per unit volume at time t. 3. Results and Discussion 3.1 Simulation parameters The initial size distribution ((𝑓(𝑅0|𝑅0 > 𝑅𝐶)) and nucleus concentration (N) are determined by the operation of experimentalist as well as nucleation method, for example, injection method34 or rapid formation of monomers upon heating/reduction.35-36 In this work, we only concentrate on the growth rather than the nucleation. And for all simulations in this work, the expectation

ACS Paragon Plus Environment

9

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 10 of 44

value and variance of the nucleus ensemble are set at 2 nm and 1 nm respectively right after the one-time nucleation. And the nuclei of radius R Eaη ) monotonously with reaction temperature. As a result, solvent with viscosity of mild temperature dependence is preferred to obtain small D/k ratio and considerable growth rate according to equation (10)&(11) at elevated reaction temperature. Higher monomer concentration will effectively reduce the size distribution for large D/k ratio, whose effect is minute for small D/k ratio. To prevent broadening effect due to Ostwald ripening,

ACS Paragon Plus Environment

25

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 26 of 44

three strategies can be used: (1) collect nanoparticles well before the depletion of monomers for diminishing monomer concentration, (2) supply new monomers during synthesis by syringe pumping or producing monomers by slow chemical reaction in the solution, and (3) reduce the numbers of nuclei. More nuclei are often obtained for quick injection or large ramp rate of temperate for ‘hot-injection’ and ‘heating-up’ method respectively, and vice versa. Secondary nucleation event could occur during synthesis to yield a bimodal distribution of nanoparticle size. Under this circumstance, longer reaction time should be used to deplete the monomer concentration to the level such that the larger nanoparticle in the bimodal distribution still can grow while the smaller nanoparticles will dissolve. By doing this, the monodispersity could be achieved via Ostwald ripening.48 However, the process is usually slow compared to the focusing effect achieved by diffusional growth, and it can be expedited for large value of D and k (higher temperature). In this work, the temperature is maintained constant (T=573 K) during simulation. By substituting equation (12)-(14) into equation (1), the effect of reaction temperature for ‘hotinjection’ method and the temperature ramp rate for ‘heating-up’ method could be clearly revealed. The above bimodal distribution problem could also be simulated by the generalized growth rate of nanoparticles. Finally, if the generalized growth rate was extended to the diffusion and adsorption of two types of monomers, the growth of compound nanoparticle ensembles could also be simulated, which will have broader applicability. These works are in progress. 4. Summary The growth of batch nanoparticles, following unimodal log-normal size distribution, has been numerically simulated using the generalized diffusional model. The cases of constant and

ACS Paragon Plus Environment

26

Page 27 of 44 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

gradually reduced monomer concentration have been simulated. Our systematical simulations have unveiled the general rules of batch nanoparticle growth by diffusional growth. The E[R]-t, Var[R]-t, and Var[R]-E[R] relationships have been shown for a broad range of growth conditions. This work reveals many hidden details of batch nanoparticle growth, which cannot be easily revealed by experimental method, such as the two stages of Ostwald ripening and the role of monomer concentration at small D/k ratio. As a result, we expect the growth of batch nanoparticles by diffusional growth be well guided by our simulation. Future works include the simulation of temperature dependence, bimodal distribution and compound nanoparticles.

ACS Paragon Plus Environment

27

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 28 of 44

FIGURES.

Figure 1. The size distribution of the nanoparticle ensemble at different time under: (a) adsorption controlled growth (D=1X10-10 m2/s, k=5X10-7 m/s with D/k=2X105 nm), (b) diffusion controlled

ACS Paragon Plus Environment

28

Page 29 of 44 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

growth (D=1X10-15 m2/s, k=5X10-1 m/s with D/k=2X10-6 nm), and (c) the intermediate growth conditions (D=1.5X10-14 m2/s, k=1X10-6 m/s with D/k=15 nm). For the simulation, N=1X1015 /m3, C∞=3X1025 /m3, Cb=4.5X1025 /m3, γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

Figure 2. The variance and average radius of the nanoparticle ensembles as a function of time at different monomer concentration: (a)&(b) adsorption controlled growth (D=1X10-10 m2/s, k=5X10-7 m/s with D/k=2X105 nm), (c)&(d) diffusion controlled growth (D=1X10-15 m2/s, k=5X10-1 m/s with D/k=2X10-6 nm), and (e)&(f) the intermediate case (D=1.5X10-14 m2/s, k=1X10-6 m/s with D/k=15 nm). For this simulation, N=1X1015 /m3, C∞=3X1025 /m3, γ=1 J/m2, Ω =1X10-30 m3, T=573 K and δt=0.2 s. The following monomer concentration is simulated (otherwise specified): Cb=4.5X1025 /m3, 6X1025 /m3, 7.5X1025 /m3, 9X1025 /m3, 12X1025 /m3, and the orange arrows indicate the ascending order of the monomer concentration.

ACS Paragon Plus Environment

29

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 30 of 44

Figure 3 (a), (b) and (c) The variance of the nanoparticle ensembles as a function of radius at different D/k ratio and monomer concentration, where D and k were varied. The variance as a function of (d) growth time and (e) radius at D/k=15 nm and different D & k, the color dots in (e) indicate the size at 3600 seconds for corresponding D&k. (f) the Var[R]/E[R] ratio as a function of radius at D/k=15 nm. The variance as a function of radius at different monomer concentration for (g) D/k=0.01 nm (h) D/k=10 nm and (i) D/k=1000 nm respectively. For this simulation, C∞=3X1025 /m3, Cb=4.5X1025 /m3 (otherwise specified), γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

ACS Paragon Plus Environment

30

Page 31 of 44 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

Figure 4. The size distribution of the nanoparticle ensembles at 0, 1200, 2400, and 3600 s at different nanoparticle concentrations under different conditions: (a)(b)&(c) adsorption controlled growth (D=1X10-10 m2/s, k=5X10-7 m/s with D/k=2X105 nm), (d)(e)&(f) diffusion controlled growth (D=7X10-15 m2/s, k=1X10-1 m/s with D/k=7X10-5 nm), and (g)(h)&(i) the intermediate case (D=1.5X10-14 m2/s, k=1X10-6 m/s with D/k=15 nm). For the simulation, C∞=3X1025 /m3, Cb0=4.5X1025 /m3, γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

ACS Paragon Plus Environment

31

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 32 of 44

Figure 5. The average size and variance of the nanoparticle ensembles as a function of time at different nanoparticle concentration from 1015 to 1020 /m3 under: (a)(d) adsorption controlled growth (D=1X10-10 m2/s, k=5X10-7 m/s with D/k=2X105 nm), (b)(e) diffusion controlled growth (D=7X10-15 m2/s, k=1X10-1 m/s with D/k=7X10-5 nm), and (c)(f) the intermediate case (D=1.5X10-14 m2/s, k=1X10-6 m/s with D/k=15 nm). The orange arrows indicate the ascending order of nanoparticle concentration. The effect of monomer concentration on the E[R]-Var[R] curves is shown for (g) D/k=0.01 nm, (h) D/k=10 nm and (i) D/k=1000 nm at N=1018 /m3. For this simulation, C∞=3X1025 /m3, Cb0=4.5X1025 /m3 (otherwise specified), γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

ACS Paragon Plus Environment

32

Page 33 of 44 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

Figure 6. (a) The three-stage Var[R]-t curve obtained for intermediate D/k ratio (D=5X10-14 m2/s and k=5X10-6 m/s with D/k= 10 nm) and large nanoparticle concentration (N=1X1018 /m3). (b) shows details of the three stages. For the simulation, C∞=3X1025 /m3, Cb=4.5X1025 /m3, γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

ACS Paragon Plus Environment

33

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 34 of 44

Figure 7. The comparison of Var[R]-E[R] curves between the case I and II growth at (a) D/k=0.01, (b) D/k=10 nm and (c) D/k=1000 nm. (d) shows the two stages of Ostwald ripening in the Var[R]E[R] curve. For the simulation, N=1X1018 /m3, C∞=3X1025 /m3, Cb=4.5X1025 /m3 (otherwise specified), γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt=0.2 s.

ACS Paragon Plus Environment

34

Page 35 of 44 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

Figure 8. The 2D D-k map of variance at different size. (a), (b) an (c) show the variance of case I growth at 5 nm, 10 nm and 15 nm respectively, and (d), (e) an (f) show the variance of case II growth at 5 nm, 10 nm and 15 nm respectively. The lines of 10 seconds, 1 hour and 10 hours were indicated, on which the corresponding time should be taken to obtain the radius. For this simulation, N=1X1018 /m3, C∞=3X1025 /m3, Cb=4.5X1025 /m3, γ=1 J/m2, Ω=1X10-30 m3, T=573 K and δt is reduced for large D and k.

ACS Paragon Plus Environment

35

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 36 of 44

ASSOCIATED CONTENT Supporting Information. The average radius and variance obtained at different time step, the saturation value of variance as a function of monomer concentration for adsorption controlled growth in case I growth, the E[R]-t curve fitting, the Cb(t)-E[R] curve, and 2D D-k map of Var[R]/E[R] ratio. AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] (Prof. T. Wen) Notes The authors declare no competing financial interest ACKNOWLEDGMENT This work was financially supported by National Natural Science Foundation of China (Nos. 51772045, 51401046, 51572042, 61734002), National Key Research Development Program (No. 2016YFA0300801), Science Challenge Project (No. TZ2018003), and Fundamental Research Funds for the Central Universities (ZYGX2016J045).

ACS Paragon Plus Environment

36

Page 37 of 44 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

REFERENCES 1.

Murray, C. B.; Norris, D. J.; Bawendi, M. G., Synthesis and Characterization of Nearly

Monodisperse CdE (E = S, Se, Te) Semiconductor Nanocrystallites. J Am Chem Soc 1993, 115 (19), 8706-8715. 2.

Feng, J. G.; Song, Q.; Zhang, B.; Wu, Y. C.; Wang, T.; Jiang, L., Large-Scale, Long-

Range-Ordered Patterning of Nanocrystals via Capillary-Bridge Manipulation. Adv Mater 2017, 29 (46). 3.

Wen, T. L.; Majetich, S. A., Ultra-Large-Area Self-Assembled Mono layers of

Nanoparticles. ACS Nano 2011, 5 (11), 8868-8876. 4.

Ji, H. J.; Trevino, J.; Tu, R.; Knapp, E.; McQuade, J.; Yurkiv, V.; Mashayek, F.; Vuong,

L. T., Long-Range Self-Assembly via the Mutual Lorentz Force of Plasmon Radiation. Nano Lett 2018, 18 (4), 2564-2570. 5.

Majetich, S. A.; Wen, T.; Booth, R. A., Functional Magnetic Nanoparticle Assemblies:

Formation, Collective Behavior, and Future Directions. ACS Nano 2011, 5 (8), 6081-6084. 6.

Shevchenko, E. V.; Talapin, D. V.; Kotov, N. A.; O'Brien, S.; Murray, C. B., Structural

diversity in binary nanoparticle superlattices. Nature 2006, 439 (7072), 55-59. 7.

Wen, T. L.; Zhang, D. N.; Wen, Q. Y.; Zhang, H. W.; Liao, Y. L.; Li, Q.; Yang, Q. H.;

Bai, F. M.; Zhong, Z. Y., Magnetic Nanoparticle Assembly Arrays Prepared by Hierarchical Self-assembly on a Patterned Surface. Nanoscale 2015, 7 (11), 4906-4911.

ACS Paragon Plus Environment

37

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

8.

Page 38 of 44

Hines, M. A.; Scholes, G. D., Colloidal PbS Nanocrystals with Size-Tunable Near-

Infrared Emission: Observation of Post-Synthesis Self-Narrowing of the Particle Size Distribution. Adv Mater 2003, 15 (21), 1844-1849. 9.

Tong, S.; Quinto, C. A.; Zhang, L. L.; Mohindra, P.; Bao, G., Size-Dependent Heating of

Magnetic Iron Oxide Nanoparticles. ACS Nano 2017, 11 (7), 6808-6816. 10. Wen, T.; Li, Y.; Zhang, D.; et al., Manipulate the magnetic anisotropy of nanoparticle assemblies in arrays. J Colloid Interf Sci 2017, 497, 14-22. 11. Cademartiri, L.; Montanari, E.; Calestani, G.; Migliori, A.; Guagliardi, A.; Ozin, G. A., Size-Dependent Extinction Coefficients of PbS Quantum Dots. J Am Chem Soc 2006, 128 (31), 10337-10346. 12. Majetich, S. A.; Wen, T.; Mefford, O. T., Magnetic Nanoparticles. MRS Bull 2013, 38 (11), 899-903. 13. Kim, T. H.; Cho, K. S.; Lee, E. K.; et al., Full-Colour Quantum Dot Displays Fabricated by Transfer Printing. Nat Photonics 2011, 5 (3), 176-182. 14. Wen, T. L.; Booth, R. A.; Majetich, S. A., Ten-Nanometer Dense Hole Arrays Generated by Nanoparticle Lithography. Nano Lett 2012, 12 (11), 5873-5878. 15. Han, S. T.; Zhou, Y.; Xu, Z. X.; Huang, L. B.; Yang, X. B.; Roy, V. A. L., Microcontact Printing of Ultrahigh Density Gold Nanoparticle Monolayer for Flexible Flash Memories. Adv Mater 2012, 24 (26), 3556-3561.

ACS Paragon Plus Environment

38

Page 39 of 44 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. Kim, D.; Resasco, J.; Yu, Y.; Asiri, A. M.; Yang, P. D., Synergistic Geometric and Electronic Effects for Electrochemical Reduction of Carbon Dioxide Using Gold-Copper Bimetallic Nanoparticles. Nat. Commun. 2014, 5. 4948 17. Wu, L. H.; Li, Q.; Wu, C. H.; Zhu, H. Y.; Mendoza-Garcia, A.; Shen, B.; Guo, J. H.; Sun, S. H., Stable Cobalt Nanoparticles and Their Monolayer Array as an Efficient Electrocatalyst for Oxygen Evolution Reaction. J Am Chem Soc 2015. 137 (22), 7071-7074 18. Wen, T. L.; Zhang, D. N.; Wen, Q. Y.; et al., Enhanced Optical Modulation Depth of Terahertz Waves by Self-Assembled Monolayer of Plasmonic Gold Nanoparticles. Advanced Optical Materials 2016, 4 (12), 1974-1980. 19. Gleich, B.; Weizenecker, R., Tomographic Imaging Using the Nonlinear Response of Magnetic Particles. Nature 2005, 435 (7046), 1214-1217. 20. Bao, Y. P.; Wen, T. L.; Samia, A. C. S.; Khandhar, A.; Krishnan, K. M., Magnetic Nanoparticles: Material Engineering and Emerging Applications in Lithography and Biomedicine. J Mater Sci 2016, 51 (1), 513-553. 21. Wen, T. L.; Krishnan, K. M., Cobalt-Based Magnetic Nanocomposites: Fabrication, Fundamentals and Applications. J Phys D Appl Phys 2011, 44 (39), 393001. 22. Wen, T. L.; Liu, D.; Luscombe, C. K.; Krishnan, K. M., Granular Magnetoresistance in Cobalt/Poly (3-hexylthiophene, 2, 5-diyl) Hybrid Thin Films Prepared by a Wet Chemical Method. Appl Phys Lett 2009, 95 (8), 082509.

ACS Paragon Plus Environment

39

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 40 of 44

23. Yu, W. W.; Falkner, J. C.; Yavuz, C. T.; Colvin, V. L., Synthesis of Monodisperse Iron Oxide Nanocrystals by Thermal Decomposition of Iron Carboxylate Salts. Chem Commun 2004, (20), 2306-2307. 24. Hufschmid, R.; Arami, H.; Ferguson, R. M.; Gonzales, M.; Teeman, E.; Brush, L. N.; Browning, N. D.; Krishnan, K. M., Synthesis of Phase-Pure and Monodisperse Iron Oxide Nanoparticles by Thermal Decomposition. Nanoscale 2015, 7 (25), 11142-11154. 25. Paik, T.; Hong, S. H.; Gaulding, E. A.; Caglayan, H.; Gordon, T. R.; Engheta, N.; Kagan, C. R.; Murray, C. B., Solution-Processed Phase-Change VO2 Metamaterials from Colloidal Vanadium Oxide (VOx) Nanocrystals. ACS Nano 2014, 8 (1), 797-806. 26. Murray, C. B.; Kagan, C. R.; Bawendi, M. G., Synthesis and Characterization of Monodisperse Nanocrystals and Close-Packed Nanocrystal Assemblies. Annu. Rev. Mater. Sci. 2000, 30, 545-610. 27. Fang, X. Q.; Cheng, X. Q.; Zhang, Y. R.; Zhang, L. J. G.; Keidar, M., Single-Step Synthesis of Carbon Encapsulated Magnetic Nanoparticles in Arc Plasma and Potential Biomedical Applications. J Colloid Interf Sci 2018, 509, 414-421. 28. Liu, X.; Zhang, X. W.; Yin, Z. G.; Meng, J. H.; Gao, H. L.; Zhang, L. Q.; Zhao, Y. J.; Wang, H. L., Enhanced Efficiency of Graphene-Silicon Schottky Junction Solar Cells by Doping with Au Nanoparticles. Appl Phys Lett 2014, 105 (18), 183901. 29. Nedyalkov, N.; Nakajima, Y.; Terakawa, M., Magnetic Nanoparticle Composed Nanowires Fabricated by Ultrashort Laser Ablation in Air. Appl Phys Lett 2016, 108 (4), 043107.

ACS Paragon Plus Environment

40

Page 41 of 44 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

30. LaMer, V. K.; Dinegar, R. H., Theory, Production and Mechanism of Formation of Monodispersed Hydrosols. J Am Chem Soc 1950, 72 (11), 4847-4854. 31. Wen, T. L.; Brush, L. N.; Krishnan, K. M., A Generalized Diffusion Model for Growth of Nanoparticles Synthesized by Colloidal Methods. J Colloid Interf Sci 2014, 419, 79-85. 32. Sugimoto, T., Preparation of Monodispersed Colloidal Particles. Advances in Colloid and Interface Science 1987, 28, 65-108. 33. Bergmann, R. B.; Bill, A., On the Origin of Logarithmic-Normal Distributions: An Analytical Derivation, and Its Application to Nucleation and Growth Processes. J. Cryst. Growth 2008, 310 (13), 3135-3138. 34. Puntes, V. F.; Krishnan, K. M.; Alivisatos, A. P., Colloidal Nanocrystal Shape and Size Control: The case of cobalt. Science 2001, 291 (5511), 2115-2117. 35. Sun, S. H.; Zeng, H.; Robinson, D. B.; Raoux, S.; Rice, P. M.; Wang, S. X.; Li, G. X., Monodisperse MFe2O4 (M = Fe, Co, Mn) Nanoparticles. J Am Chem Soc 2004, 126 (1), 273279. 36. Newman, J. D. S.; Blanchard, G. J., Formation of Gold Nanoparticles Using Amine Reducing Agents. Langmuir 2006, 22 (13), 5882-5887. 37. He, M.; Protesescu, L.; Caputo, R.; Krumeich, F.; Kovalenko, M. V., A General Synthesis Strategy for Monodisperse Metallic and Metalloid Nanoparticles (In, Ga, Bi, Sb, Zn, Cu, Sn, and Their Alloys) via in Situ Formed Metal Long-Chain Amides. Chem Mater 2015, 27 (2), 635-647.

ACS Paragon Plus Environment

41

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 42 of 44

38. Park, J.; An, K. J.; Hwang, Y. S.; Park, J. G.; Noh, H. J.; Kim, J. Y.; Park, J. H.; Hwang, N. M.; Hyeon, T., Ultra-Large-Scale Syntheses of Monodisperse Nanocrystals. Nat Mater 2004, 3 (12), 891-895. 39. Niesner, R.; Heintz, A., Diffusion Coefficients of Aromatics in Aqueous Solution. J. Chem. Eng. Data 2000, 45 (6), 1121-1124. 40. Vreeland, E. C.; Watt, J.; Schober, G. B.; et al., Enhanced Nanoparticle Size Control by Extending LaMer's Mechanism. Chem Mater 2015, 27 (17), 6059-6066. 41. Jansons, A. W.; Hutchison, J. E., Continuous Growth of Metal Oxide Nanocrystals: Enhanced Control of Nanocrystal Size and Radial Dopant Distribution. ACS Nano 2016, 10 (7), 6942-6951. 42. Bastus, N. G.; Comenge, J.; Puntes, V., Kinetically Controlled Seeded Growth Synthesis of Citrate-Stabilized Gold Nanoparticles of up to 200 nm: Size Focusing versus Ostwald Ripening. Langmuir 2011, 27 (17), 11098-11105. 43. Bao, Y. P.; An, W.; Turner, C. H.; Krishnan, K. M., The Critical Role of Surfactants in the Growth of Cobalt Nanoparticles. Langmuir 2010, 26 (1), 478-483. 44. Dehsari, H. S.; Harris, R. A.; Ribeiro, A. H.; Tremel, W.; Asadi, K., Optimizing the Binding Energy of the Surfactant to Iron Oxide Yields Truly Monodisperse Nanoparticles. Langmuir 2018, 34 (22), 6582-6590. 45. Xia, Y. N.; Gates, B.; Yin, Y. D.; Lu, Y., Monodispersed Colloidal Spheres: Old Materials with New Applications. Adv Mater 2000, 12 (10), 693-713.

ACS Paragon Plus Environment

42

Page 43 of 44 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

46. Bian, X.; Kim, C.; Karniadakis, G. E., 111 Years of Brownian Motion. Soft Matter 2016, 12 (30), 6331-6346. 47. Bonhote, P.; Dias, A. P.; Papageorgiou, N.; Kalyanasundaram, K.; Gratzel, M., Hydrophobic, Highly Conductive Ambient-Temperature Molten Salts. Inorg. Chem. 1996, 35 (5), 1168-1178. 48. Zhang, Z. R.; Wang, Z. N.; He, S. N.; Wang, C. Q.; Jin, M. S.; Yin, Y. D., Redox Reaction Induced Ostwald Ripening for Size- and Shape-Focusing of Palladium Nanocrystals. Chem. Sci. 2015, 6 (9), 5197-5203.

ACS Paragon Plus Environment

43

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 44 of 44

Table of Contents Graphic

ACS Paragon Plus Environment

44