Mesoscale Particle-Based Model of Electrophoretic Deposition

20 Dec 2016 - We analyze particle configurations in order to observe how colloids accumulate at the ... We postulate ordering and stacking within the ...
1 downloads 0 Views 4MB Size
Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article

Mesoscale Particle-based Model of Electrophoretic Deposition Brian Giera, Luis A. Zepeda-Ruiz, Andrew J. Pascall, and Todd H. Weisgraber Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.6b04010 • Publication Date (Web): 20 Dec 2016 Downloaded from http://pubs.acs.org on December 24, 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 31

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

Langmuir

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

Mesoscale Particle-based Model of Electrophoretic Deposition

Brian Giera∗ , Luis A. Zepeda-Ruiz∗ , Andrew J. Pascall, and Todd H. Weisgraber Lawrence Livermore National Laboratory, Livermore, California 94550, USA Abstract We present and evaluate a semi-empirical particle-based model of electrophoretic deposition using extensive mesoscale simulations. We analyze particle configurations in order to observe how colloids accumulate at the electrode and arrange into deposits. In agreement with existing continuum models, the thickness of the deposit increases linearly in time during deposition. Resulting colloidal deposits exhibit a transition from highly ordered and bulk disordered region that can give rise to an appreciable density gradient under certain simulated conditions. The overall volume fraction increases and falls within a narrow range as the driving force due to the electric field increases and repulsive inter-colloidal interactions decrease. We postulate ordering and stacking within the initial layer(s) dramatically impacts the microstructure of the deposits. We find a combination of parameters, i.e. electric field and suspension properties, whose interplay enhances colloidal ordering beyond the commonly known approach of only reducing the driving force. ∗

To whom correspondence should be addressed. E-mail: [email protected]

1

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

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

1

Introduction

Electrophoretic deposition (EPD) occurs when electric field-driven colloidal particles suspended in a fluid migrate towards an electrode or liquid-liquid interface 1 where they assemble into a deposit. The deposition rate depends on many parameters, such as the applied field strength, volume fraction of colloids, and electrophoretic mobility. 2 The electrophoretic mobility is a function of the colloidal suspension properties, e.g. ionic strength, 3 permittivity, and viscosity of the suspension as well as the electric potential difference between the surface of the colloid and suspension (termed the zeta-potential). Due to the vast material set inherent to colloidal suspensions, 4 EPD has been demonstrated in systems of aqueous electrolytes, 5 non-polar solvents, 6,7 and ionic liquids 8 that stabilize colloids comprised of ceramics, 9 graphene-related materials, 10 nanometals, 11 polymers, 12 phosphors, 13 or biological materials. 14 EPD is central to numerous applications 15 that include industrial coating, 16 fabrication of supercapacitor electrodes 17 and quantum-dot light-emitting devices, 18 near-shape manufacturing of tubes or structured parts, 19 and controlled growth of artificial opals. 20 Recently, two additive manufacturing techniques were developed that relay on EPD with photoconductive 21 and coaxial cable 22 electrodes. Since its discovery more than two centuries ago, 5 EPD has maintained practical relevance in academia and industry because of its versatility and scalability. 23 Despite widespread usage, an understanding of the underlying mechanisms of deposit formation and microstructure evolution during EPD is presently lacking. 24 EPD kinetics are determined typically by weighing the deposited film at different time points, 25 which is a purely ex situ technique. While there are some emergent in situ probes that utilize a quartz crystal microbalance, 26 the only data collected is total mass of the deposit – no particle level details of deposition kinetics are obtained. Usually, the morphology of the deposit is determined by SEM or scattering measurements; however, this requires drying of the deposit that can introduce artifacts due to, e.g., capillary forces are well known to rearrange particles into crystals, cracking resulting from mismatched shrinkage of the film and substrate, 23 etc. 2

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

Furthermore, experimental measurements of deposits provide volumetric information and are always destructive. The limits of these experimental probes underscore the need for predictive models of EPD dynamics and the final structure of deposits. The majority of available models provide continuum level descriptions of EPD 27 that relate mass deposition rates to the electrophoretic velocity, deposition area, and a “sticking parameter” that accounts for the fraction of colloids that irreversibly incorporate into the deposit. 9,28–30 Film thickness models suggest nonuniform deposits occur when the colloidal particle permittivity exceeds that of the suspension 31 or near electrode edges where electric field singularities locally enhance deposition rates. 32 Continuum models inherently neglect inter-colloidal interactions within the bulk suspension, growth front, and deposit. Alternative particle-based modeling of EPD is still nascent, 33 but promises to offer more detailed predictive information about deposit formation and packing morphology than traditional continuum-level approaches. For instance, two-dimensional Discrete Element Method EPD simulations were used to construct “tentative phase diagrams” that imply particle area fraction of 2D deposits depends on electrolyte pH and concentration and electrode current density. 34 Brownian Dynamic simulations of “low dimensional” systems of colloids reveal grain boundary relaxation mechanisms in 2D colloidal crystal formation in an effort to tune voltages to avoid kinetic bottlenecks associated with slow grain boundary dynamics. 35 We recently presented a 3D mesoscale particle-based model of electrophoresis 36 that accurately captures the diffusive and electrophoretic motion of colloids over a comprehensive range of experimentally relevant conditions. Although our model successfully characterizes mass transport of bulk colloidal suspensions, it did not yet account for colloid-wall interactions, which are necessary to investigate particle rearrangement, defect formation, and the ultimate microstructure of the deposit. Here, we develop and evaluate a particle-based model of colloidal suspensions that undergo electrophoretic motion and deposition using an extensive set of mesoscale simulations representative of experimentally relevant colloidal suspensions. Since the model explicitly

3

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

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

computes inter-colloidal interactions, it is uniquely poised to elucidate how deposition conditions influence defect structures and particle rearrangement within electrophoretic deposits. We use the model to investigate how empirical parameters, such as electric field strength and the range of inter-colloidal repulsion between electric double layers, can be tuned in order to control the degree of colloidal ordering versus non-ordering that occurs during EPD. We observe spatial-temporal particle density profiles and examine colloidal ordering throughout the deposit. We find the deposit density increases for suspensions characterized by high ionic strength under strong driving forces. The degree of colloidal ordering increases for suspensions characterized by moderate ionic strength under weak driving forces.

2

Model Details

We perform mesoscale particle dynamics simulations of monodisperse colloidal suspensions in a uniform field, bounded by a flat electrode in order to model the electrophoretic motion and deposition of colloids. We use LAMMPS 37 to evaluate the internal (solvent-colloid and colloid-colloid) and external (electric field and wall) interactions that govern colloidal transport and deposit morphology within our simulations. Below, we describe these interaction potentials, the empirical properties they depend on, and how they are implemented.

2.1

Internal and External Interactions

In a system of Ncol colloids, the total energy acting on any particle is expressed as a sum of the internal Uinternal and external Uexternal energy components that depend on the pairwise separation distance between colloids r = |ri − rj | and each colloid’s distance from electrode surface z U (r

Ncol

)=

Ncol X

Uinternal (r) +

i 214, the resulting deposit is

20

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

Page 22 of 31

Page 23 of 31

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

indicating it is possible to optimize ordering within a range of Debye lengths. While for Pe=214 the fraction of ordered particles at 12.5 nm is nearly double that of all other λD , at Pe=26.5 the ordered fraction within 6.25 nm≤ λD ≤ 25 nm is three times as large as Nordered /N outside of that range. A surprising finding is the existence of such a window of Debye lengths that allows for a dramatic improvement in ordering as the driving force decreases. We believe that the range of optimal values in λD can be attributed to DLVO interactions. Since the DLVO potential contains appreciable attractive energy wells for λD ≤ 6.25 nm (see Fig. 1), there are many possible locations (at a lattice site or otherwise) for a colloid to incorporate into the deposit that results in sufficiently low energy configurations. Thus, colloids with attractive interactions are trapped where they land on the deposit and do not rearrange into more favorable sites. At the values of λD that favor ordering, repulsive forces are sufficient to drive depositing colloids towards locations within newly formed layers with lower energy. In the case of systems with λD = 50 nm, colloids in suspension experience mostly repulsive forces allowing them to incorporate into the deposit nearly everywhere with little difference in energy change. For λD = 12.5 nm, the extreme cases in ordering are shown for the lowest and highest Pe in Figures 7(b1-b3) and 7(c1-c3), respectively. The snapshot in Fig. 7b1 shows a higher number of ordered particles (blue) that extend further from the electrode than Fig. 7c1. The densification at the larger field correlates with a reduction in colloidal ordering. For clarity, we remove disordered particles (grey) to better view this difference in Figs. 7b2 and c2. In these panels, particles that appear isolated/floating are surrounded by a local environment characterized by short-ranged order. Finally, top views in Figs. 7b3 and c3 showing only ordered colloids away from electrode confirm domain sizes decrease with increasing Pe, which is consistent with observations of ordering at the electrode in Fig. 5.

22

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

Conclusions

We present a mesoscale model of EPD that explicitly accounts for inter-colloidal interactions to provide detailed predictions about the evolution and resulting microstructure of colloidal deposits. The model is formulated within a particle dynamics framework and assumes pairwise and external potentials that require only empirically measurable parameters as inputs. In its present form, the model is poised to explore systematically the vast interplay between the wide range of empirical parameters, e.g., electrophoretic mobility, particle radius, bulk volume fraction, suspension properties, among others, that govern EPD. We simulate a variety of monodisperse colloidal suspensions (listed in Table 1) that characterize suspensions with exclusively repulsive to partially attractive inter-colloidal interactions, under moderate to extreme uniform electric fields. We analyze in detail particle configurations in order to explore how colloids accumulate at the electrode. Simulated deposits increase in thickness linearly with time in accordance with traditional continuum-based theories. Fully deposited systems exhibit highest volume fractions for suspensions with thin double layers, e.g. concentrated electrolyte, under strong electric fields at large Pe. In the opposite limit, our simulations suggest a route to building deposits with graded density. Our model results indicate the volume fraction of the deposit can be tuned with the Péclet number and that the degree of tunability increases with Debye length. Unlike what is possible for standard (non-particle-based) modeling approaches, we use the model to investigate the degree of colloidal ordering versus non-ordering that occurs within the deposit during EPD. At lower fields and intermediate λD , domain sizes of ordered colloids at the electrodes is larger, which results in an increased overall fraction of ordered colloids within the deposit. We postulate that the degree of ordering at the electrode dictates the propensity of colloidal ordering in regions farther from the deposit. Enhanced control of the shape, composition, and performance of functional materials fabricated via EPD can benefit from computationally feasible models like the one presented here that can predict formation and morphology of colloidal depositions. The computational 23

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

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

expense of this model allows for EPD simulations across time and length scales that are amenable to direct experimental comparison. Here, we evaluate the model using a large set of simulations within a regime characterized by short times and thin deposits. For experimental validation schemes requiring thicker deposits, this model simply could be run for longer times or parameterized with larger particles, which permit longer time steps. Since the model predictions apply to wet deposits, comparisons against morphological features of dried deposits may not be reliable due to cracking that occurs from drying. Nevertheless, realtime particle tracking techniques 59 enable direct comparisons against simulated trajectories. Also, time series of weight measurements taken from wet or dry deposits offer a useful metric for assessing deposition rates. Colloidal ordering can be validated using scattering and/or tomography; however, sample preparation techniques such as gelling or freeze drying wet deposits could be used to preserve the microstructure. An immediate improvement to the model requires modification to account for changes in the resistance of and resulting potential drop across the deposit that strongly influences deposition kinetics, especially for thick deposits. Also, various treatments of the electrode/particle interface, e.g., a lattice of particles, an amorphous monolayer, surface charge heterogeneities, etc., along with spatial and temporal electric field variations can be investigated for their influence on defect structures, particle rearrangement, and deposit growth kinetics. Considering the variety and intricacies of actual colloidal suspensions, e.g. polydispersity and/or non-spherical colloids, it is fortunate that this model is extensible to simulate systems with size and/or surface potential distributions, ellipsoidal particles, or mixtures of different particle compositions. Furthermore, other non-DLVO systems, such as those with polarizable double layers, Janus or “patchy” particles, 60 or polymer coated colloids can be accounted for by including the appropriate particle interactions into the model.

Acknowledgement This work was performed under the auspices of the U.S. Department of Energy by Lawrence 24

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

Livermore National Laboratory under Contract DE-AC52-07NA27344.

References (1) Buie, C. Electrophoretic Deposition at the Interface of Immiscible Liquids. Presented at 5th International Conference on Electrophoretic Deposition: Fundamentals and Applications, Hernstein, Austria, October 5-10, 2014. (2) Smoluchowski, M. Contribution à la Théorie de l’Endosmose Électrique et de Quelques Phénomènes Corrélatifs. Pisma Mariana Smoluchowskiego 1924, 1, 403–420. (3) Hückel, E. Die Kataphorese der Kugel. Phys. Z 1924, 25, 204–210. (4) Russel, W. B.; Saville, D. A.; Schowalter, W. R. Colloidal Dispersions; Cambridge University Press, Cambridge, 1989. (5) Reuss, F. F. Sur un Nouvel Effet de l’Électricité Galvanique. Mem. Soc. Imp. Natur. Moscou 1809, 2, 327–337. (6) Corni, I.; Ryan, M. P.; Boccaccini, A. R. Electrophoretic Deposition: from Traditional Ceramics to Nanotechnology. J. Eur. Ceram. Soc. 2008, 28, 1353–1367. (7) Gonzalo-Juan, I.; Krejci, A. J.; Dickerson, J. H. Toward Dynamic Control over TiO2 Nanocrystal Monolayer-by-Monolayer Film Formation by Electrophoretic Deposition in Nonpolar Solvents. Langmuir 2012, 28, 5295–5301. (8) Große-Brauckmann, J.; Lair, V.; Argirusis, C. Application of Ionic Liquids in the Electrophoretic Deposition of Oxide Layers. Key Eng. Mat. 2009, 412, 27–32. (9) Sarkar, P.; Nicholson, P. S. Electrophoretic Deposition (EPD): Mechanisms, Kinetics, and Application to Ceramics. J. Am. Ceram. Soc. 1996, 79, 1987–2002.

25

ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31

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

(10) Diba, M.; Fam, D. W.; Boccaccini, A. R.; Shaffer, M. S. Electrophoretic Deposition of Graphene-Related Materials: A Review of the Fundamentals. Prog. Mater Sci. 2016, 82, 83–117. (11) Boccaccini, A. R.; Roether, J. A.; Thomas, B. J.; Shaffer, M. S.; Chavez, E.; Stoll, E.; Minay, E. J. The Electrophoretic Deposition of Inorganic Nanoscaled Materials– A Review. J. Ceram. Soc. Jpn. 2006, 114, 1–14. (12) Qian, F.; Pascall, A. J.; Bora, M.; Han, T. Y.-J.; Guo, S.; Ly, S. S.; Worsley, M. A.; Kuntz, J. D.; Olson, T. Y. On-Demand and Location Selective Particle Assembly via Electrophoretic Deposition for Fabricating Structures with Particle-to-Particle Precision. Langmuir 2015, 31, 3563–3568. (13) Talbot, J.; Sluzky, E.; Kurinec, S. Electrophoretic Deposition of Monochrome and Color Phosphor Screens for Information Displays. J. Mater. Sci. 2004, 39, 771–778. (14) Boccaccini, A.; Keim, S.; Ma, R.; Li, Y.; Zhitomirsky, I. Electrophoretic Deposition of Biomaterials. J. R. Soc. Interface 2010, 7, S581–S613. (15) Electrophoretic Deposition of Nanomaterials; Dickerson, J. H., Boccaccini, A. R., Eds.; Springer: New York, 2012; Vol. 50. (16) Davey, W. Process for Making and Applying Japan. US Patent 1,294,627, Feb 18, 1919. (17) Du, C.; Pan, N. High Power Density Supercapacitor Electrodes of Carbon Nanotube Films by Electrophoretic Deposition. Nanotechnology 2006, 17, 5314. (18) Song, K. W.; Costi, R.; Bulović, V. Electrophoretic Deposition of CdSe/ZnS Quantum Dots for Light-Emitting Devices. Adv. Mater. 2013, 25, 1420–1423. (19) Tabellion, J.; Clasen, R. Electrophoretic Deposition from Aqueous Suspensions for Near-shape Manufacturing of Advanced Ceramics and GlassesâĂŤApplications. J. Mater. Sci. 2004, 39, 803–811. 26

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

(20) Holgado, M.; García-Santamaría, F.; Blanco, A.; Ibisate, M.; Cintas, A.; Míguez, H.; Serna, C. J.; Molpeceres, C.; Requena, J.; Mifsud, A.; Meseguer, F.; López, C. Electrophoretic Deposition To Control Artificial Opal Growth. Langmuir 1999, 15, 4701– 4704. (21) Pascall, A. J.; Qian, F.; Wang, G.; Worsley, M. A.; Li, Y.; Kuntz, J. D. Light-Directed Electrophoretic Deposition: A New Additive Manufacturing Technique for Arbitrarily Patterned 3D Composites. Adv. Mater. 2014, 26, 2252–2256. (22) Nold, A.; Zeiner, J.; Assion, T.; Clasen, R. Electrophoretic Deposition as Rapid Prototyping Method. J. Eur. Ceram. Soc. 2010, 30, 1163–1170. (23) Bersa, L.; Liu, M. A Review on Fundamentals and Applications of Electrophoretic Deposition (EPD). Prog. Mater. Sci. 2007, 52, 1 – 61. (24) Van der Biest, O. O.; Vandeperre, L. J. Electrophoretic Deposition of Materials. Annu. Rev. Mater. Sci. 1999, 29, 327–352. (25) Sarkar, P.; De, D.; Rho, H. Synthesis and Microstructural Manipulation of Ceramics by Electrophoretic Deposition. J. Mater. Sci. 2004, 39, 819–823. (26) Fatisson, J.; Domingos, R. F.; Wilkinson, K. J.; Tufenkji, N. Deposition of TiO2 Nanoparticles onto Silica Measured Using a Quartz Crystal Microbalance with Dissipation Monitoring. Langmuir 2009, 25, 6062–6069. (27) Ferrari, B.; Moreno, R. EPD Kinetics: A Review. J. Eur. Ceram. Soc. 2010, 30, 1069 – 1078. (28) Hamaker, H. Formation of a Deposit by Electrophoresis. Trans. Faraday Soc. 1940, 35, 279–287. (29) Ferrari, B.; Moreno, R.; Cuesta, J. A. A Resistivity Model for Electrophoretic Deposition. Key Eng. Mat. 2006, 314, 175–180. 27

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

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

(30) Biesheuvel, P. M.; Verweij, H. Theory of Cast Formation in Electrophoretic Deposition. J. Am. Ceram. Soc. 1999, 82, 1451–1455. (31) Gonzalez-Cuenca, M.; Biesheuvel, P. M.; Verweij, H. Modeling Constant Voltage Electrophoretic Deposition from a Stirred Suspension. AIChE J. 2000, 46, 626–631. (32) Pascall, A. J.; Sullivan, K. T.; Kuntz, J. D. Morphology of Electrophoretically Deposited Films on Electrode Strips. J. Phys. Chem. B 2013, 117, 1702–1707. (33) Park, J. S.; Saintillan, D. Direct Numerical Simulations of Electrophoretic Deposition of Charged Colloidal Suspensions. Key Eng. Mat. 2012, 507, 47–51. (34) Cordelair, J.; Greil, P. Discrete Element Modeling of Solid Formation During Electrophoretic Deposition. J. Mater. Sci. 2004, 39, 1017–1021. (35) Edwards, T. D.; Yang, Y.; Beltran-Villegas, D. J.; Bevan, M. A. Colloidal crystal grain boundary formation and motion. Sci. Rep. 2014, 4 . (36) Giera, B.; Zepeda-Ruiz, L. A.; Pascall, A. J.; Kuntz, J. D.; Spadaccini, C. M.; Weisgraber, T. H. Mesoscale particle-based model of electrophoresis. JES 2015, 162, D3030– D3035. (37) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular-Dynamics. J. Comput. Phys. 1995, 117, 1–19. (38) Schunk, P. R.; Pierce, F.; Lechman, J. B.; Grillet, A. M.; in’t Veld, P. J.; Weiss, H.; Stoltz, C.; Heine, D. R. Performance of Mesoscale Modeling Methods for Predicting Rheological Properties of Charged Polystyrene/Water Suspensions. J. Rheol. 2012, 56, 353. (39) Ball, R. C.; Melrose, J. R. A Simulation Technique for Many Spheres in Quasi-static Motion under Frame-invariant Pair Drag and Brownian Forces. Physica A 1997, 247, 444–472. 28

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

(40) Kumar, A.; Higdon, J. J. L. Origins of the Anomalous Stress Behavior in Charged Colloidal Suspensions under Shear. Phys. Rev. E 2010, 82, 051401. (41) Ermak, D. L.; McCammon, J. A. Brownian Dynamics with Hydrodynamic Interactions. J. Chem. Phys. 1978, 69, 1352–10. (42) Kumar, A. Microscale Dynamics in Suspensions of Non-spherical Particles. Ph.D. thesis, University of Illinois at Urbana-Champaign, 2010. (43) Bybee, M. D. Hydrodynamic Simulations of Colloidal Gels: Microstructure, Dynamics, and Rheology; ProQuest, 2009. (44) Bolintineanu, D. S.; Grest, G. S.; Lechman, J. B. Particle Dynamics Modeling Methods for Colloid Suspensions. CPM 2014, 1, 321–356. (45) Derjaguin, B.; Landau, L. The Theory of Stability of Highly Charged Lyophobic Sols and Coalescence of Highly Charged Particles in Electrolyte Solutions. Acta Physicochim. URSS 1941, 14, 633–52. (46) Verwey, E.; Overbeek, G. Theory of the Stability of Lyophobic Colloids. Elsevier: Amsterdam 1948, 1, 205. (47) Safran, S. Statistical Thermodynamics of Surfaces, Interfaces, and Membranes; Addison-Wesley Pub., 1994. (48) Israelachvili, J. N. Intermolecular and Surface Forces: Revised Third Edition; Academic Press, 2011. (49) Everaers, R.; Ejtehadi, M. R. Interaction Potentials for Soft and Hard Ellipsoids. Phys. Rev. E 2003, 67, 041710. (50) Elimelech, M.; Gregory, J.; Jia, X. Particle Deposition and Aggregation: Measurement, Modeling and Simulation; Butterworth-Heinemann, 2013. 29

ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

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

(51) Probstein, R. F. Physicochemical Hydrodynamics: An Introduction; John Wiley & Sons, 2005. (52) Henry, D. C. The Cataphoresis of Suspended Particles Part I - The Equation of Cataphoresis. Proc. R. Soc. A 1931, 133, 106–129. (53) Biben, T.; Hansen, J.-P.; Barrat, J.-L. Density Profiles of Concentrated Colloidal Suspensions in Sedimentation Equilibrium. J. Chem. Phys. 1993, 98, 7330–7344. (54) Böhmer, M. In Situ Observation of 2-Dimensional Clustering during Electrophoretic Deposition. Langmuir 1996, 12, 5747–5750. (55) Trau, M.; Saville, D.; Aksay, I. Field-induced Layering of Colloidal Crystals. Science 1996, 272, 706. (56) Hayward, R.; Saville, D.; Aksay, I. Electrophoretic Assembly of Colloidal Crystals with Optically Tunable Micropatterns. Nature 2000, 404, 56–59. (57) Yeh, S.-R.; Seul, M.; Shraiman, B. I. Assembly of Ordered Colloidal Aggregrates by Electric-field-induced Fluid Flow. 1997, (58) Honeycutt, J. D.; Andersen, H. C. Molecular Dynamics Study of Melting and Freezing of Small Lennard-Jones Clusters. J. Phys. Chem. 1987, 91, 4950–4963. (59) Barnkob, R.; Kähler, C. J.; Rossi, M. General Defocusing Particle Tracking. Lab Chip 2015, 15, 3556–3560. (60) DeLaCruz-Araujo, R. A.;

Beltran-Villegas, D. J.;

Larson, R. G.;

Córdova-

Figueroa, U. M. Rich Janus Colloid Phase Behavior Under Steady Shear. Soft Matter 2016, 12, 4071–4081.

30

ACS Paragon Plus Environment