Understanding the Stabilization of Single-Walled Carbon Nanotubes

Dec 19, 2014 - Understanding Surfactant Stabilization of MoS2 Nanosheets in Aqueous ... Carbon Nanotubes for Biodelivery and Optical Sensing Applicati...
0 downloads 0 Views 2MB Size
Subscriber access provided by AKDENIZ UNIVERSITESI

Article

Understanding the Stabilization of Single-Walled Carbon Nanotubes and Graphene in Ionic Surfactant Aqueous Solutions: Large-Scale CoarseGrained Molecular Dynamics Simulations and Modified DLVO Theory Chih-Jen Shih, Shangchao Lin, Michael S. Strano, and Daniel Blankschtein J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/jp5093477 • Publication Date (Web): 19 Dec 2014 Downloaded from http://pubs.acs.org on December 21, 2014

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 49

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

Understanding the Stabilization of Single-Walled Carbon Nanotubes and Graphene in Ionic Surfactant Aqueous Solutions: Large-Scale CoarseGrained Molecular Dynamics Simulations-Assisted DLVO Theory Chih-Jen Shih,1,2 Shangchao Lin,1,3,4 Michael S. Strano,1and Daniel Blankschtein1* 1

Department of Chemical Engineering, MIT, Cambridge, MA 02139

2, Current Address 3

Department of Chemical Engineering, Stanford University, Stanford, CA 94305

Department of Civil and Environmental Engineering, MIT, Cambridge, MA 02139

4, Current Address

Department of Mechanical Engineering and Materials Science & Engineering Program, Florida State University, Tallahassee, Florida 32310

*Author to whom all correspondence should be addressed. Email: [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

ABSTRACT

Understanding the dispersion of single-walled carbon nanotubes (SWCNTs) and graphene in surfactant aqueous solutions is essential for processing of these materials. Herein, we develop the first theoretical framework, which combines large-scale coarse-grained (CG) molecular dynamics (MD) simulations with the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory and Langmuir isotherm, to inform the mechanisms of surfactant adsorption and induced colloidal stability. By carrying out large-scale CG-MD simulations (~1 million CPU hours), we successfully calculated the surface coverage of the widely-used ionic surfactant sodium cholate (SC) on different carbon nanomaterials at experimentally realistic SC concentrations. This was accomplished by simultaneously simulating SC micellization and adsorption in one simulation box. Our theoretical framework further allows us to quantify the surface electric potential and the potential energy barrier height that maintains colloidal stability of dispersed SWCNTs or graphene sheets as a function of SC concentration, C, and radius of the SWCNT, r. We found that for a specific carbon nanomaterial, there exists an optimal surfactant concentration, resulting from a trade-off between the increase in surfactant adsorption (i.e., surface charge density) and the increase in ionic strength (i.e., inverse Debye length) in the bulk aqueous phase. For the first time, we predict that small-radius SWCNTs promote higher surfactant adsorption than largeradius SWCNTs in the high SC concentration regime, and surprisingly, that the opposite behavior is exhibited in the low SC concentration regime. These findings suggest new experimental implications for the separation of SWCNTs based on their radius by tuning the SC concentration. In addition, we also predict that monolayer graphene exhibits better colloidal stability than multilayer graphene.

2 ACS Paragon Plus Environment

Page 2 of 49

Page 3 of 49

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 Carbon nanomaterials, such as single-walled carbon nanotubes (SWCNTs) and graphene,

have attracted considerable research interest due to their outstanding optoelectronic, thermal, and mechanical properties.1,2 These properties, however, are associated with single-chirality, isolated, SWCNTs or graphene, which requires advanced post-processes to debundle the SWCNTs, or to exfoliate graphene, for large-scale production. To this end, it is common to disperse these carbon nanomaterials in an aqueous surfactant solution,3-5 a method that has shown its superiority for further scaling-up, separation, and chemical functionalization.6-10 Indeed, this dispersion method represents the only viable route to sort SWCNTs by their chirality, following application of techniques like electrophoresis,11-14 density gradient ultracentrifugation (DGU),7,9 or gel chromatography6,15,16 to the solution. In addition, surfactant-induced stabilization and exfoliation of graphene becomes potentially important for the production of AB-stacked bilayer and trilayer graphene.8,17,18 This method can also be used as a high quality alternative to chemical vapor deposition.19 In combination with the appropriate surface chemistry on the target substrate, the interactions between the adsorbed surfactant molecules on the carbon nanomaterials and the functionalized substrate surface further enable the creation of device arrays for integrated circuits.20,21

The

key

challenge

for

the

implementation

of

this

method

involves

debundling/exfoliating the appropriate carbon nanomaterials over a specific range of surfactant concentrations, capable of dispersing the resulting SWCNTs/graphene sheets in a stable manner.3 The widely-recognized, ionic surfactant-aided dispersion process requires applying an external energy input (ultrasound) to separate the carbon nanomaterials at the bundled SWCNT ends (or at the graphite flake edges) by overcoming the strong van der Waals (vdW) attractions between them.3 Subsequently, the separated ends (edges) provide new adsorption sites for the 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

hydrophobic tails of the surfactant molecules.22 As a result, the repulsive potential energy resulting from the electrostatic interactions between the hydrophilic, ionic heads of the adsorbed surfactant molecules further enhances this separation process.22 Eventually, individually-isolated, surfactant-coated carbon nanomaterials are released to and dispersed in the bulk aqueous phase. The ionic surfactant molecules adsorbed on the nanomaterial surfaces generate a diffuse layer of counter-ions present in the bulk solution phase to form an electrical double layer (EDL).23 The diffusive nature of the counter-ions results in an effective surface charge, which can be quantified in terms of the  potential, the surface electric potential under an electrophoretic field.23 The interactions between nearby EDLs induce a repulsive potential energy, which offsets the attractive vdW potential energy, thereby stabilizing the suspended carbon nanomaterials. The competitive repulsive and attractive potential energies can be quantified using the DerjaguinLandau-Verwey-Overbeek (DLVO) theory.23 In recent years, significant research effort has been devoted to: (i) screen through many types of surfactants for the preparation of high weight fraction, aqueous solutions of SWCNTs, including ranking of surfactants in terms of their ability to solubilize suspended SWCNTs and graphene,18,24-29 and (ii) understand the self-assembled structures of the adsorbed surfactant molecules around SWCNTs and graphene,29-34 including rationalizing the surfactant-induced colloidal stability as a function of the surfactant concentration, C,18,35-42 and the SWCNT radius, r (note that, for graphene, r   ).36,43-45 For example, it has been experimentally shown that the widely-used anionic surfactant SC preferentially stabilizes SWCNTs possessing smaller radii,44 and is therefore particularly suitable for SWCNT separation via DGU.7 In addition, the stability of carbon nanomaterials can be optimized when the SC concentration, C, is in the range 10 × CMC > C > CMC, where CMC is the critical micelle concentration of SC,18,27 although micellization is not a prerequisite to obtain

4 ACS Paragon Plus Environment

Page 4 of 49

Page 5 of 49

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

stable colloidal dispersions.37 Indeed, these experimental findings clearly suggest that the adsorption of surfactant molecules highly depends on the geometry of the carbon nanomaterials and the surfactant concentration. In this respect, however, very little is known about the molecular details of the interactions between carbon nanomaterials and surfactant molecules, including correlating these interactions with the surfactant-assisted colloidal stability. In particular, there is a lack of fundamental understanding about the adsorption isotherms that describe the behavior at high surfactant concentrations (above the CMC), where typical experiments of carbon nanomaterials dispersions are conducted. In order to bridge macroscopic variables (e.g., C) and the nanoscale self-assembly, molecular simulations have been used to elucidate the interactions between carbon nanomaterials and surfactants. For example, molecular dynamics (MD) simulations have been carried out to resolve the self-assembled structures of surfactants, including sodium dodecyl sulfate (SDS)43,46,47, single-tailed and double-tailed phosphatidylcholine,48-50 and SC,51,52 on SWCNTs and graphene, by using all-atomistic (AA) or coarse-grained (CG) force fields. In addition, interactions between two surfactant-coated SWCNTs / graphene sheets were studied using MD simulations to elucidate the contributions from electrostatics and vdW interactions to the simulated potential of mean force (PMF) between two parallel SWCNTs or graphene sheets.43,51,52 The density differences of SC-SWCNT assemblies having different SWCNT radii have been also discussed.53,54 Despite significant progress on modeling/rationalizing local structural features of self-assembled surfactant molecules on carbon nanomaterials, the surface coverage of surfactant molecules on a carbon nanomaterial as a function of C, as well as the effect of the carbon nanomaterial geometry, have not been investigated computationally. The main reason is that modeling experimentally-relevant surfactant concentrations, typically in the

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

CMC range, requires simultaneous simulation of surfactant micellization and adsorption in one simulation box, which is extremely computationally expensive. Indeed, using AA-MD simulations, one can only simulate a length scale of ~10 nm, which is too small to demonstrate micelle formation, requiring length scales of the order of ~50 nm. In addition, the time scale for capturing the dynamics of micelle formation is also very demanding (~ 100 ns).55 As a result, most MD simulations can only be used to model/rationalize local structural features of selfassembled surfactant molecules on carbon nanomaterials under a “known” surfactant surface coverage (e.g., Ref. 43), while the dependence of the surface coverage on r and C, as well as the effect of micellization, remain undetermined. Consequently, the mechanism underlying the surfactant-aided stabilization of aqueous dispersions of carbon nanomaterials is still not sufficiently understood to permit a rational design of suitable surfactants for different carbon nanomaterials at the molecular level. With the above background in mind, in this report, we develop the first theoretical framework, which combines CG-MD simulations and the DLVO theory and Langmuir isotherm, to inform the mechanisms of surfactant adsorption and the induced colloidal stability, as well as their dependence on surfactant concentration, C, and radius of the SWCNT, r. In order to reduce computational expense, in addition to using CG simulations, the planar anionic surfactant – sodium cholate (SC), which has a small micelle aggregation number (n = 2 ~ 6),56 is considered. As a result, SC micellization and SC adsorption can be simulated simultaneously in one simulation box, thereby enabling simulation of graphene/SWCNT aqueous solutions at experimentally realistic concentrations (C = 1 – 200 mM). Specifically, we have carried out large-scale CG-MD simulations (~ 1 million CPU hours), and have successfully calculated the SC surface coverage of different carbon nanomaterials, including SWCNTs (6,6), (12,12), and

6 ACS Paragon Plus Environment

Page 6 of 49

Page 7 of 49

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

(20,20), as well as monolayer, bilayer, and trilayer graphene (MLG, BLG, and TLG), as a function of SC concentration. The simulated surface coverage allows us to generate a continuum spectrum of the surface electric potential, 0, as a function of r and C. Combining the simulation results with the DLVO theory and Langmuir isotherm, we have further quantified the potential energy barrier height required to maintain colloidal stability between two parallel identical SWCNTs or graphene sheets in aqueous SC solutions as a function of r and C. The theoretical methodology presented here allows us to rationalize many experimental findings, as well as to provide practical implications for separating different carbon nanomaterials from their isomers in a solution phase. Our investigation also provides fundamental insight on the molecular-level design of efficient surfactant or co-surfactant systems7,57 to enhance the dispersion and separation of SWCNTs and graphene in the solution phase.

2.

Simulation Method We carried out MD simulations under the NPT (constant number of atoms, constant

pressure of 1.0 bar, and constant temperature of 298.15 K) ensemble using the GROMACS 4.558 software package and the MARTINI-CG force field.59 The SWCNT and graphene species, including SWCNTs (6,6), (12,12), (20,20), MLG, BLG, and TLG, were modeled using the 3:1 mapping scheme proposed by Wallace and Sansom,49 and described by the uncharged SC359 Lennard-Jones (LJ) spheres in the MARTINI force field. The target carbon nanomaterial was fixed in the middle of the simulation box with 3D periodic boundary conditions. For the SWCNT simulations, the length of the simulation box along the x-axis, Lx, was kept constant in order to mimic an infinitely-long nanotube along this axis, and similarly, for the graphene simulations, the lengths of the simulation box along the x-axis and y-axis, Lx and Ly, respectively, were kept

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

constant to simulate infinitely-large graphene sheets in the xy plane. In terms of the water model, it is noteworthy that we observed a serious freezing of CG water when the original MARTINI water beads or even antifreeze beads59 are used, particularly for graphene simulations due to flatgeometry-enforced ordering effects. Therefore, to avoid this artifact, the polarizable MARTINI water model is used,60 which models four water molecules with three CG beads. Sodium cholate, unlike conventional linear surfactants like SDS, is a rigid facial amphiphile, often referred to as a “two-faced detergent”. Figures 1a and 1b show its chemical structure and all-atomistic molecular structure, respectively. The SC surfactant possesses a quasi-planar, slightly bent but rigid steroid ring with a hydrophilic face (the hydroxyl groups and the charged carboxylate group), and a hydrophobic face (the methyl groups and the tetracyclic carbon backbone) residing back-to-back. The CG model of SC was adopted from Marrink,61 and the resulting CG structure is shown in Fig. 1c. All the non-bonded interactions (LJ and Coulombic) were treated using the default values in the MARTINI library.59 The velocity-rescaled Berendsen thermostat was implemented to maintain the system at a constant temperature.62 The pressure was coupled to a semiisotropic Berendsen barostat (Py = Pz = 1 bar for the SWCNT case, and Pz = 1 bar for the graphene case).63 The LJ and Coulombic interactions were treated with a cutoff distance of 1.2 nm. The equation of motion was integrated at a time step of 0.02 ps for the self-assembly process. For each simulation with a different carbon nanomaterial, we varied the number of surfactant molecules in the simulation box by distributing cholate ion monomers and sodium counter-ions randomly in the bulk phase. Subsequently, depending on the surfactant concentration considered, SC monomers can stay in the bulk, form micelles, or adsorb on the carbon nanomaterial surface freely. For each simulation, an ~ 8 × 8 × 24 nm3 simulation box, which contains ~12,500 coarse-

8 ACS Paragon Plus Environment

Page 8 of 49

Page 9 of 49

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

grained water molecules, is used. For example, Fig. 1d (top) shows a post-equilibrium snapshot of a MLG-SC simulation, which is comprised of SC monomers, micelles, and a sheet of SCadsorbed MLG. In the system considered, the density profile of cholate ions as a function of z was calculated, as shown in Fig. 1d (bottom), and z = 0 corresponds to the location of the MLG. One can observe a peak on each side of the two graphene surfaces, which is associated with monolayer adsorption of the cholate molecules. The density profile of the cholate ions indicates that they prefer to distribute in a very compact manner within 1 nm on the graphene surface (where the density profile decays to the bulk concentration beyond z = ±1 nm). Note that throughout the concentration range considered in this report, we did not observe multilayer adsorption. Using the obtained concentration profile, the surfactant concentration in the bulk region, C (in units of mM), was therefore quantified by averaging the concentration of cholate ions in the bulk region (z > 1 or < ‒1 nm), while the surface coverage of surfactants on carbon nanomaterials,  (in units of nm-2), was obtained by integrating the density of cholate ions along the axis normal to the nanomaterial surface (the z-axis in the case of graphene and the radial axis in the case of SWCNTs) in the adsorption region (‒1 < z < 1 nm). It is noteworthy that in the low SC concentration regime, particularly in the case of graphene, almost all the inserted SC molecules are adsorbed, leaving only a few SC monomers in the simulation box. Consequently, the assumptions of constant concentration in the bulk phase, as well as of thermodynamic equilibrium of the system, may be less accurate. The accuracy of the simulations in the low SC concentration regime can be improved by considering a larger simulation box, which is more computationally expensive and beyond the scope of this report. The simulations considered in this report typically took at least 2 - 5 s to reach equilibrium, by monitoring the calculated  as a function of time. Therefore, each simulated system was equilibrated for a total of 20 s, and

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

Page 10 of 49

only the last 5 s of simulation were used for data analysis. Considering the time scale and the MARTINI mapping schemes used in this manuscript, we estimate that the computational load is only about 1% to 5% of that using AA simulation methods. In order to carry out large-scale calculations (note that ~1 million CPU hours were used), we believe that the CG approach is still necessary. For the statistical analysis of micelle formation in the bulk region, we assume that surfactant molecules belong to the same micelle if the distance between any hydrophobic beads (SC1 and SC3 beads) within their structures is smaller than the cutoff distance of 7.5 Å.64 The large-scale CG-MD simulations allow us to investigate the effects of surfactant concentration and carbon nanomaterial geometry for the first time.

3.

Results and Discussion

3.1

Large-Scale CG-MD Simulations

3.1.1

Surfactant Adsorption at Low Surfactant Concentration Representative post-equilibrium simulation snapshots of SWCNTs (6,6), (12,12), and

(20,20) covered with cholate ions at two specific concentrations (C = 10 and 60 mM) are shown in Figs. 2a, 2b, and 2c, respectively. Similar to the AA MD-simulated self-assembly structure of cholate ions on SWCNTs and graphene,51,52 the cholate molecules form an adsorbed monolayer on the surface. An important feature in the low-concentration (10 mM) simulations is that, instead of distributing uniformly on the SWCNT surfaces, the adsorbed cholate molecules prefer to self-assemble as clusters, exposing a fraction of the SWCNT surface to the aqueous phase, particularly for SWCNT (6,6). Interestingly, in the three configurations associated with C = 10 mM (see Figs. 2a, 2b, and 2c. all left), it appears that the extent of surface coverage for SWCNT (6,6) is much lower than those for SWCNTs (12,12) and (20,20). Since more surfactant

10 ACS Paragon Plus Environment

Page 11 of 49

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

adsorption generally implies a lower interfacial tension and stronger repulsive electrostatic interactions between the dispersed nanomaterials (for additional discussions, see Section 3.2),3 this finding predicts that at low SC concentration, large-radius SWCNTs and graphene (because r   for graphene) should exhibit higher colloidal stability. It further suggests that dispersing

SWCNTs at low SC concentration potentially allows separation of large-radius SWCNTs from small-radius SWCNTs. To our knowledge, this method has not been proposed before, and we look forward to its experimental confirmation in the future. At low SC concentration, we believe that the preferential adsorption of cholate ions onto large-radius SWCNTs results from the rigid facial structure of the cholate ion. In general, adsorption of surfactant monomers from the bulk phase to the hydrophobic carbon nanomaterial surface is determined by the competition between enthalpic and entropic contributions. For a surfactant monomer, in order to minimize the contact area between its hydrophobic tail and water, there is an enthalpic driving force for adsorption onto the hydrophobic carbon nanomaterial surface. However, since SC is a planar surfactant, in order to cover the surface, it has to bend to some extent, depending on the local curvature of the surface, with a higher extent of bending resulting in a reduction in the enthalpic contribution, thereby making the SC adsorption less enthalpically favorable. On the other hand, when the surfactant monomer is adsorbed, it is confined on a 2-dimensional surface, leading to a reduction in its entropy. In the case of adsorption onto large-radius SWCNTs (see Figs. 2b and 2c left), the local curvature is small, such that the adsorption is more enthalpically favorable, and the enthalpic contribution dominates the entropic one, thereby driving more adsorption of SC molecules onto the SWCNT’s. On the other hand, in the case of adsorption onto small-radius SWCNTs at C = 10 mM (see Fig. 2a left), the enthalpic contribution is diminished because the cholate molecules

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 49

need to bend to a higher extent in order to cover the more curved surface of the small-radius SWCNT. As a result, the SC monomers prefer to stay in the bulk phase, and for the adsorbed SC molecules, they tend to form aggregates to diminish the enthalpic bending penalty. The same rationale can also be used to explain the SC adsorption at high SC concentration (C = 60 mM, see Figs. 2a, 2b, and 2c right). The SC concentration increase in the bulk phase diminishes the entropic contribution, such that the surface is covered to a greater extent by the cholate molecules. On a small-radius SWCNT, there is an enthalpic bending penalty that hinders extensive SC adsorption, while the formation of SC aggregates on the surface, on the other hand, can stabilizes the adsorbed SC molecules and reduces the enthalpic bending penalty. On the other hand, on a large-radius SWCNT, the enthalpic bending penalty is small, such that the adsorbed SC molecules tend not to aggregate on the surface. Additional discussion is presented in Section 3.2.

3.1.2

Micellization in the Bulk Phase As discussed in Section 1, one of the most difficult challenges associated with simulating

surfactant adsorption onto the surfaces of nanomaterials is to accommodate micellization and adsorption in one simulation box. In order to monitor micellization in the bulk phase, we statistically analyzed the population of SC micelles as a function of the total SC concentration, C. Figure 3a shows the simulated concentration of free SC monomers (not part of a micelle) at the CMC in the bulk phase, C1, as a function of C. The dashed vertical line corresponds to the experimental CMC of SC (~ 10 mM)65, the solid line corresponds to C1 = C, and the blue circles correspond to the simulated C1. Recall that for an ideal micellization process, when the surfactant

12 ACS Paragon Plus Environment

Page 13 of 49

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

concentration exceeds the CMC, it becomes thermodynamically unfavorable for the added surfactant molecules to remain dispersed as free monomers. Accordingly, C1 should deviate significantly from the line C1 = C, reflecting micelle formation. In other words, when C > CMC, C1 should reach a plateau and remain constant at the CMC value.23 Figure 3b shows the simulated populations of micelles as a function of micelle aggregation number for the two SC concentrations considered (10 mM and 200 mM). For C = 200mM, more than 50% of the population consists of micelles (n > 2), while for C = 10mM, monomers are clearly the dominant species (> 95%). However, the CG-MD-simulated C1 (the blue circles in Figure 3a) only exhibit a gradual deviation from the C1 = C line. As a result, it becomes challenging to infer the CMC from Fig. 3a, since the monomer-to-micelle transition is too gradual. Two observations can be made to rationalize the behavior shown in Fig. 3b: (i) The size of the SC micelles is inherently small (the reported experimental aggregation number, n = 2 ~ 6, which is consistent with the simulated range of small n values in Fig. 3b),56 because its facial structure does not allow the formation of large aggregates. Consequently, a gradual monomer-to-micelle transition is expected,23 and (ii) The CG force-field parameters of water and SC do not allow a very accurate CMC prediction, based on the free monomer concentration.64 Several recent reports have suggested that in the case of linear surfactants like SDS, which form micelles having larger aggregation numbers (n ~ 50),23 the CG-MD-simulated free monomer concentration strongly depends on the total surfactant concentration considered, and the extracted CMC values are often systematically underestimated.64,66 Observations (i) and (ii) are also reflected in the simulated aggregation number distributions (see Fig. 3b), which do not show a peak aggregation number, with larger micelles becoming more favorable at higher surfactant concentrations.64 Nevertheless, to our knowledge, the results presented here represent the first CG-MD study to directly

13 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 14 of 49

investigate the monomer-to-micelle transition near C = CMC. This was possible because the small aggregation number of the SC micelles significantly reduces the computational expense incurred at low SC concentrations. We hope that advances in the development of graphics processing units (GPU)-based algorithms,64 as well as further refinements in the development of CG force fields, will facilitate simulating the micellization behavior of linear surfactants like SDS, which form micelles having large aggregation numbers. This will allow developing a better understanding of the mechanisms underlying micelle structural transitions at the molecular level.

3.1.3

Modeling SC Adsorption and Micellization in One Simulation Box The various symbols in Figure 4a correspond to the simulated SC surface coverage,  (in

units of nm-2), as a function of C for the three SWCNTs considered. First, as discussed in Section 3.1.1, when C < 10 mM, the simulated  vs. C for SWCNT (6,6) is only about half of those for SWCNTs (12,12) and (20,20) due to the large local curvature of SWCNT (6,6), which significantly reduces the enthalpy change associated with SC adsorption. Moreover, the simulated s for SWCNTs (12,12) and (20,20) are very close, suggesting that the curvature effect becomes negligible when r exceeds some threshold value. As discussed in Section 3.1.1, we believe that understanding the SWCNT radius-dependence of the SC surface coverage at low SC concentration may be useful for separating large-radius SWCNTs from small-radius ones. Note that for C < 10 mM, almost all the SC molecules in the bulk phase are free SC monomers (see Fig. 3a). Therefore, the increase in  with SC concentration shown in Fig. 4a results from the equilibrium between free SC monomers and adsorbed SC molecules.

14 ACS Paragon Plus Environment

Page 15 of 49

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

Next, by further increasing C (C > 20mM), the population of micelles increases gradually (see Fig. 3a). The CG-MD simulated  for SWCNT (6,6) increases rapidly and exceeds those for SWCNTs (12,12) and (20,20) at C ~ 20 mM, corresponding to the vertical dashed line in Fig. 4a. In addition, the surface coverage for SWCNT (20,20) appears to increase more slowly, and begins deviating from that for SWCNT (12,12). In other words, relative to the simulations at low SC concentrations, the trend of radius dependence is reversed, that is, the small-radius SWCNTs exhibit higher SC surface coverage. More quantitative discussions will be presented in Section 3.2. Our simulated SC surface coverage is consistent with the recent experimental finding that SC preferentially stabilizes SWCNTs with smaller radius at a high SC concentration (specifically, a ~ 30 mM SC aqueous solution was used).44 As suggested in Fig. 4a, we further predict that a higher SC concentration would be even more favorable to differentiate between SWCNTs with different radii. Considering the fact that a sufficiently high SC concentration can reduce the solubility of SWCNTs27 (additional mechanistic discussions are provided in Section 3.2), based on our simulation results, we propose that a SC concentration of ~ 40 – 80 mM (or 1.7 – 3.5 wt%) would be desirable to separate small-radius SWCNTs from large-radius ones. In combination with post-separation techniques (e.g., DGU7,9 and electrophoresis11-14), we believe that surfactant concentration is an useful parameter to tune the solubility of SWCNTs based on their radius, at least for the case of SC. In order to model the dependence of the simulated  on SC concentration and SWCNT radius at high SC concentration, we utilize the Langmuir isotherm to fit the simulated  values for C > 20 mM. Specifically,23

(r , X )  max (r )

K ads (r )X 1  K ads (r )X

15 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 16 of 49

where r is the radius of the SWCNT, max (r) is the radius-dependent saturated SC surface coverage, X is the molar fraction of cholate in the bulk phase (i.e., X ~ 18C/1000, where C is in the unit of M), and Kads (r) is the radius-dependent equilibrium constant, which can be calculated using the free-energy change associated with the adsorption of a SC molecule from the bulk phase to the surface of a SWCNT of radius r, Gads (r), as follows:   G ads (r )   K ads (r )  exp   kBT 

(2)

where kB is the Boltzmann constant, and T is the absolute temperature. Note that at constant temperature, the two parameters, max (r) and Gads (r), depend solely on r. Note also that Eqs. (1) and (2) are based on the assumptions of: (i) cholate monolayer adsorption, and (ii) no interactions between the adsorbed cholate molecules.23 Assumption (i) is generally valid (see Fig. 1d), while assumption (ii) is approximate, particularly at low SC concentrations, where as discussed above, the adsorbed cholate molecules tend to self-assemble as clusters (see Fig. 2a). Consequently, to carry out the required fitting, we only considered our simulation results at high SC concentration (C > 20 mM). Specifically, by least-square fitting Eqs. (1) and (2) to the MD simulation results for C > 20 mM, we extracted values for max (r) and Gads (r) for each SWCNT considered, as shown in Table 1. The solid lines in Fig. 4a correspond to the fitted Langmuir isotherms for each SWCNT considered, where the color of each line and symbols match. Figure 4a shows that each Langmuir isotherm can describe the SC adsorption behavior for C > 20 mM reasonably well, suggesting that: (i) SC adsorption reaches equilibrium with the SC micelles and the SC monomers in the bulk phase, and (ii) interactions between the adsorbed cholate molecules are either negligible or weak relative to the SWCNT-SC interactions. Note, however, that the isotherms in Fig. 4a, which are obtained by fitting the simulated  for C > 20 mM, do not

16 ACS Paragon Plus Environment

Page 17 of 49

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

accurately describe the SC adsorption at low SC concentration (C < 10 mM). In particular, they do not allow prediction of the preferential adsorption of SC on large-radius SWCNTs (see symbols in Fig. 4a and Section 3.1.1). Moreover, the adsorption isotherms fitted at high SC concentration tend to underestimate the SC surface coverage at low SC concentration, except for the small-radius SWCNTs (e.g., SWCNT (6,6) in Fig. 4a). As discussed above, we believe that SC micellization in the bulk phase is playing an important role. Nevertheless, in order to attain a better solubility of carbon nanotubes, typically, it is experimentally desirable to disperse SWCNTs and graphene above the CMC,3 as we described in Section 1. Hereafter, we will utilize the fitted Langmuir isotherms to model the adsorption behavior. The various symbols in Fig. 4b correspond to the CG-MD simulated  for MLG, BLG, and TLG. The Langmuir isotherms (the various lines having colors corresponding to those of the symbols) were obtained by fitting the simulated  for C > 20 mM. The extracted values of max and Gads are shown in Table 1. Note that because graphene allows SC adsorption on both sides, the simulated surface coverages for graphene are inherently higher. The simulated SC surface coverages for MLG, BLG, and TLG do not change significantly with the SC concentration. Interestingly, as Fig. 4b shows, the simulated SC surface coverages on graphene at high SC concentrations are significantly limited, as will be discussed later. In addition, the simulated s for MLG, BLG, and TLG are very close, indicating that SC adsorption is not strongly influenced by the substances adsorbed on the opposite side of the first graphene layer, including the adsorbed cholate molecules or the additional graphitic layers. This finding confirms that MLG is a “weakly-transparent” material for the transmission of vdW interactions through it.67,68 As indicated above, Table 1 summarizes the extracted values of max and Gads for the three SWCNTs and three graphene species considered. The extracted values of Gads are all very 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 49

close, suggesting that the adsorption free-energy change is mainly determined by the interactions between the SC molecules and the carbon atoms on the nanomaterial surface, irrespective of the local curvature. Therefore, we can further simplify Eqs. (1) and (2) by assuming that Gads is a constant, which is independent of the radius r of the SWCNT (or graphene, corresponding to an infinite r value). Accordingly, the average value of Gads = -8.5 kBT (calculated by averaging the six Gads values reported in Table 1) will be used in all the calculations which follow. On the other hand, Table 1 shows that max is r-dependent (for graphene, r   ), and more importantly, a small-radius SWCNT tends to have higher max (the surface coverage when C   ), as stressed earlier. In order to understand the resulting radius dependence, in Fig 5b,

the extracted max values are plotted as a function of 1/r. Note that for graphene, 1/r is equal to 0, and in order to make a proper comparison, the calculated max for graphene should be divided by 2, since graphene potentially allows SC adsorption on both sides. Interestingly, we find that the extracted max is linearly proportional to 1/r, suggesting a semi-empirical relation for max(r), which is given by:

max (r )  a 

b r

(3)

where a and b are fitting parameters ( a = 0.947 nm-2 and b = 0.725 nm-1 are obtained for the various nanomaterials considered in Fig. 5b). Along with Eqs. (1) and (2), the extracted values of a, b, and Gads allow us to generate a continuous description of  as a function of r and C, based on the CG-MD simulated results. Note that as discussed above, Eqs. (1) - (3) should provide more accurate predictions at high SC concentration (C > CMC). The mechanism that is responsible for the radius-dependence in the high SC concentration regime is schematically illustrated in Fig. 5a. Specifically, by increasing the SC

18 ACS Paragon Plus Environment

Page 19 of 49

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

concentration, the SC surface coverage increases to a great extent and the interactions between the charged hydrophilic heads of the adsorbed SC molecules become more significant. These interactions are dominated by electrostatic repulsions, which thermodynamically limits further adsorption of SC molecules. Figure 5a depicts schematically the adsorption of cholate molecules on small-radius (top) and large-radius (bottom) SWCNTs having the same surface coverage. Note that although the structure of the cholate molecule is relatively facial, following adsorption, the negatively-charged carboxyl group (the cyan bead in Fig. 1c) does not completely attach to the carbon nanomaterial surfaces (see Figs. 2a – 2c), due to steric effects. Therefore, in Fig. 5a, the extension attached to each hydrophilic SC head group only represents part of its hydrophobic tail. As shown in Fig. 5a, the small-radius SWCNT allows a larger distance between the charged head groups, thereby enabling a higher saturated SC surface coverage. The reverse is observed for the large-radius SWCNTs. In other words, the electrostatic repulsions between the adsorbed surfactant clusters increase with increasing r. Consequently, small-radius SWCNTs allow more SC adsorption at high SC concentration, where max is inversely proportional to the radius of the SWCNT, as discussed above. Interestingly, recall that at low SC concentration, there is an enthalpic bending penalty for the adsorption on small-radius SWCNTs (see Section 3.1.1). However, it appears that at high SC concentration, the electrostatic repulsion between the adsorbed SC head groups becomes more dominant. One can observe that the adsorbed SC molecules on a small-radius SWCNT tend to aggregate (see Fig. 2), in order to gain more enthalpic reduction. By increasing the SC concentration in the bulk phase, more SC aggregates are formed, thereby stabilizing SC adsorption on small-radius SWCNTs. Consequently, we expect that the effect of the enthapic

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 49

bending penalty is diminished and overwhelmed by the electrostatic repulsions between the adsorbed SC head groups. We can therefore conclude that the radii of SWCNTs and graphene determine the SC surface coverage through different mechanisms. In the low-SC concentration regime, the SC surface coverage for small-radius SWCNTs is lower, since the cholate molecules need to bend to a greater extent in order to cover the increasingly curved surface. On the other hand, in the highSC concentration regime, the electrostatic repulsions between the adsorbed SC clusters increase with increasing r, resulting in a higher saturated SC surface coverage for small-radius SWCNTs. Because, experimentally, it is generally desirable to disperse SWCNTs and graphene above the CMC of SC, in the next section, we will utilize a simplified adsorption model (Eqs. (1) – (3)) to describe (r,C), which actually provides a better prediction of SC adsorption in the high-SC concentration regime.

3.2

Simulation-Assisted DLVO Theory

3.2.1

Modeling the Surface Electric Potential

The theoretical analysis presented so far provides a thermodynamic description of surfactant (SC) adsorption on an isolated SWCNT or graphene sheet. However, for practical purposes, in order to correlate the calculated (r,C) to the resulting colloidal stability observed in actual experiments, additional modeling which accounts for the interactions between two surfactant-coated SWCNTs or graphene sheets is required. With this in mind, we next present a theoretical analysis that combines the (r,C) results from our CG-MD simulations with a simulation-assisted DLVO theory23 of colloidal stability. This analysis will allow an overall

20 ACS Paragon Plus Environment

Page 21 of 49

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

understanding of the surfactant-assisted stabilization of SWCNT and graphene aqueous dispersions. In order to elucidate the colloidal stability of SWCNTs/graphene sheets with different chiralities, in this section, we model the interactions between two parallel, identical, SC-coated SWCNTs or graphene sheets in water. The following assumptions were made to simplify the model: (i) The thickness of the adsorbed cholate molecules is negligible, implying that the adsorbed SC molecules are treated as effective surface charges uniformly distributed on the surface, (ii) The SC molecules in micellar form contribute to the solution ionic strength as the SC monomers do. This assumption is reasonable because: (1) the simulated SC micelle size is generally small (n ~ 2 – 6). As a result, the electric field generated by the negatively-charged SC micelle is expected to be weak, resulting in small binding of positively-charged sodium counterions, (2) the SC micelles are dynamic entities, and (3) According to our CG-MD simulations, even at a very high SC concentration (C = 200 mM), more than 50% of the cholate molecules exist as monomers (see Fig. 3b). We anticipate that assumption (ii) may overestimate the solution ionic strength at a very high SC concentration, (iii) The dominant interactions between two carbon nanomaterials (SWCNTs or graphene sheets) include: (a) vdW attractions, and (b) electrostatic repulsions resulting from the adsorbed SC molecules, consistent with the description underlying the DLVO theory. Assumption (i) above implies that the charge density on the carbon nanomaterial surface as a function of r and C, (r,C), is given by:

 (r ,C )  (r ,C )  ( e)

(4)

where e is the elementary charge. The exact solution of the Poisson-Boltzmann equation for a cylinder of radius r in a 1:1 electrolyte (where the electrolyte refers to SC) aqueous solution

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 49

provides an analytical description to relate the surface electrostatic potential, 0(r,C), to (r,C) as follows:69  e (r ,C )   2 0 k BT  2  1  (r ,C )  sinh  0 1      2 e cosh (e0 / 4k BT )   2k BT  

1/2

(5)

where ε0 is the dielectric permittivity of vacuum, ε is the relative dielectric permittivity of water, and  is a shape (cylindrical) parameter given by:

 

K 0 ( r ) K 1 ( r )

(6)

where Ki (i=0 and 1) is the i-th-order modified Bessel function of the second kind,

 

e 2  is the inverse of the Debye-Hückel screening length, and   is the SC number  0k BT

concentration in the bulk. Assumption (ii) above implies that    2CN A , where NA is Avogadro’s number. It is noteworthy that Eq. (5) reduces to the well-known Goüy-Chapman solution23 of the Poisson-Boltzmann equation in the case of graphene ( r   ). Indeed, lim   1 , which reduces the last term in Eq. (5) to unity. In combination with the proposed

r 

adsorption model (Eqs. (1) – (3)) and the three parameters a, b, and Gads deduced in Section 3.1.3 (Gads = -8.5 kBT, a = 0.947 nm-2, and b = 0.725 nm-1), the surface electrostatic potential as a function of r and C, 0(r,C), can be calculated by solving Eq. (5) numerically. A continuous spectrum of 0, spanning a concentration and reciprocal radius space, can therefore be generated, as shown in Fig. 6. The scale of 0 ranges from -80 mV (blue) to -160 mV (red) (see the color scale bar in Fig. 6). The reciprocal radii corresponding to SWCNTs (6,6), (12,12), and (20,20), and to graphene are shown as dashed vertical lines for reference (note that for graphene, 1/r = 0, corresponding to the y-axis). For each carbon nanomaterial considered (along the dashed vertical

22 ACS Paragon Plus Environment

Page 23 of 49

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

lines for the three SWCNTs and the y-axis for graphene), there exists an optimal surfactant concentration which depends on 1/r, Copt(1/r), corresponding to the maximum value of 0, as shown by the solid white line in Fig. 6. Interestingly, the optimal concentration for graphene, Copt(0), corresponds to C ~ 10 mM, consistent with the experimental CMC of SC.65 This finding is actually coincidental, since as we discussed above, our CG-MD simulations do not accurately describe the micellization of SC in the bulk aqueous phase. In fact, we found that our model (Eqs. (1) – (6)) mathematically demands that the value of Copt(0) be equal to the molar fraction of ln(G ads ) (i.e., X opt (0)  ln(G ads ) ), by analytically determining the maximum of lim 0 (r ,C ) . r 

The value of Copt(1/r) slightly increases with 1/r due to the shape (cylindrical) factor in Eq. (5) (the term associated with ). The optimal concentration results from the competition between the SC adsorption and the ionic strength in the bulk aqueous phase. Qualitatively speaking, one can simplify Eq. (5) by considering the limits of graphene (   1 ) and of the Debye-Hückel approximation ( e0  k BT ), which yield    00 .23 Combining this result with Eq. (4) and the definitions of  and   given above, it follows that: 0   / C

(7)

Indeed, along the dashed vertical line for graphene (1/r = 0) in Fig. 6, one can clearly see the competition between  and

C . Specifically, at low SC concentration, as described by the

Langmuir isotherm (see Eq. (1)),  increases rapidly with C and dominates the C term in the denominator of Eq. (7), resulting in the increase of 0 (more negative) with SC concentration. On the other hand, at high SC concentration (C > 10 mM), SC adsorption approaches its saturated surface coverage. Accordingly, the increase of  with SC concentration becomes much 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

slower, and the

Page 24 of 49

C term becomes dominant, thereby reducing 0 (less negative) with

concentration. Although mathematically more complex in the case of SWCNTs, because of the finite curvature of the SWCNT surface and the nonlinearity of the Poisson-Boltzmann equation when  is not unity (see Eq. (5)), a similar competition still exists between the SC adsorption (the  in the numerator of Eq. (7)) and the ionic strength in the bulk aqueous phase (the C in the denominator of Eq. (7)). Theoretically, the predicted 0 is closely related to the experimentally-measured zeta potential, , which corresponds to the electric potential in the EDL at the location of the slipping plane, residing a few -1 away from the charged surface.23 As a result, it is possible to relate  to the stability of dispersions of carbon nanomaterials.27,29 For example, a highly-negative zeta potential typically suggests better colloidal stability, due to a stronger electrostatic repulsion between the dispersed nanomaterials.23 Therefore, if we assume that   0 , Fig. 6 clearly rationalizes the experimental findings that: (i) at any given r, the colloidal stability of SWCNTs and graphene is significantly reduced when C >> CMC (0 becomes more bluish or less negative), and (ii) there is an optimal SC concentration range, which maximizes colloidal stability, at around 1 to 10 times the CMC,18,27,37 corresponding to the reddish area in Fig. 6. Furthermore, the theoretical predictions presented in Fig. 6 also have important implications for the separation of SWCNTs based on their radius. Specifically, in addition to the separation of small-radius SWCNTs from large-radius SWCNT’s at high C, as we suggested earlier, it is also possible to separate SWCNTs based on their radius using electrophoresis11-14 by tuning the surfactant concentration. To this end, one needs to experimentally investigate the zeta potential as a function of r and C, which to our knowledge, has not yet been pursued extensively. In the future, in addition to developing new techniques to prepare large-scale, single-chirality SWCNT 24 ACS Paragon Plus Environment

Page 25 of 49

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

or graphene dispersions,6 we hope that a better understanding can be advanced to correlate the predicted 0 to the experimentally-measured .

3.2.2

Modeling the Energy Barrier Height

In order to further understand the surfactant-assisted stabilization of dispersions of carbon nanomaterials, we next quantify the interactions between two identical, parallel SWCNTs or graphene sheets, as a function of the intertube or intersheet separation, d. In combination with the calculated 0(r,C), the repulsive electrostatic potential energy between two identical, parallel SWCNTs per unit length, rep, as a function of r, C, and d, can be obtained by considering the interactions between two EDLs,23 which is given by:23 ( r )1 / 2  64 0 (k BT / e)2 tanh 2 e0 (r , C ) / 4k BT  exp(  d ) 2 

 rep (r , C , d ) 

(8)

As assumed above, the attractive vdW potential energy between two identical, parallel SWCNTs per unit length, att, is not influenced by surfactant adsorption. Therefore, att is independent of C, and can be obtained by integrating the Lennard-Jones (LJ) potential over the surfaces of two SWCNTs. This yields:70  att (r , d ) 

3 2 21B ( A I A  IB ) 3 8r 32r 6

(9)

where A and B are the attractive and repulsive Lennard-Jones constants (A = 15.2 eV×Å6, and B = 24.1×103 eV×Å12),70   4 /

3a 2 is the mean surface density of carbon atoms,70 and IA and IB

are integrals given by:70 IA 

  (cos 

IB 

  (cos 

2 2

0

0

2

 cos 1 )2  (sin  2  sin 1  d / r )2



5 / 2

d 1d  2

(10) 2 2

0

0

2

 cos 1 )2  (sin  2  sin 1  d / r )2

25 ACS Paragon Plus Environment



 11 / 2

d 1d  2

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 49

where 1 and 2 are angular coordinates, ranging from 0 to 2. Therefore, the total potential energy between two identical, parallel SWCNTs per unit length, , is given by:23 (r , C , d )   att (r , d )   rep (r ,C , d )

(11)

Note that the unit of  in Eq. (11) is energy per unit length, which is inherently higher for largeradius SWCNTs and diverges in the graphene limit ( r   ). Therefore, graphene is studied separately and will be discussed below. Using Eqs. (8) - (10), with known 0(r,C) (see Fig. 6), one can calculate  rep (r ,C , d ) ,  att (r , d ) , and (r ,C , d ) . For each SWCNT species considered (known r), the continuum

spectrum of as a function of C and d was generated. As an illustration, the calculated total potential energy, , between two parallel (6,6) SWCNTs per unit length as a function of the intertube separation, d, and the surfactant concentration, C, is shown in Fig. 7a. The energy scale of  ranges from zero (blue) to 27 meV/nm (red) (see the color scale bar in Fig. 7a). Note that in order to focus on the energy barrier height, the energy scale for  < 0 is not presented here (i.e., the color scale for  < 0 remains blue). At a fixed surfactant concentration (constant y, e.g., C = 10 mM), the total energy is zero when the two SWCNTs are far apart ( d   ). By reducing d, the competition between att and rep leads to an increase of  (repulsive) first, followed by a rapid decrease to negative (attractive), when the two SWCNTs are too close, as described by the classical DLVO theory.3,23 This is not surprising, since the vdW attractions are relatively shortranged, and therefore, when the two SWCNTs are far apart, the repulsive electrostatic potential dominates, thereby kinetically stabilizing the SWCNTs. In other words, one can view the reddish region for d > 10Å as the “boundary layer”, responsible for protecting the SWCNT from re-

26 ACS Paragon Plus Environment

Page 27 of 49

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

aggregation, at a specific surfactant concentration. Since att is independent of C (see Eq. (9)), the thickness of the boundary layer directly relates to the Debye-Hückel screening length, -1, which characterizes the length scale of electrostatic influence in an electrolyte solution, i.e., a larger -1 results in a wider boundary layer thickness. Consequently, as we discussed earlier, because -1 decreases with C, Fig. 7a clearly shows that the boundary layer thickness of  decreases with C. At a fixed surfactant concentration, we can also calculate the height of the energy barrier hindering the re-aggregation of the SWCNTs, Barrier, as the maximum potential energy attainable, corresponding to d = dmax. Accordingly, the following two equations are utilized to determine Barrier:       d r ,C

 0 ;  Barrier (r , C )  (C , r , d max )

(12)

d dmax

The calculated height of the energy barrier, as a function of SWCNT radius, r, and SC concentration, C, is shown in Fig. 7b. The energy barrier height scales from zero (blue) to 50 meV/nm (red) (see the color bar scale in Fig. 7b). Typically, a higher value of Barrier implies better colloidal stability.23 However, again, the unit used for Barrier (energy per unit length) magnifies the values for the large-radius SWCNTs. To interpret Fig. 7b properly, one should compare the relative values of Barrier for each SWCNT considered (along the vertical dashed lines in Fig. 7b). Accordingly, for each SWCNT, we are able to clearly define an optimal concentration, corresponding to the highest energy barrier height, located at C ~ 10 mM, which is close to that described by the Copt curve in Fig. 6. Another important finding is that the energy barrier height of the SWCNTs at a very high SC concentration (C > 200 mM) is even lower (more bluish) than that at a low SC concentration (C ~ 1 mM). This clearly suggests that 27 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 28 of 49

stabilization of SWCNTs at a low SC concentration should be feasible, while a high SC concentration could lead to an even lower colloidal stability of SWCNTs, consistent with most experimental observations.37,38 Following a modeling procedure similar to the one discussed above in the SWCNT case, a simulation-assisted DLVO theory was developed to model the stabilization of graphene sheets in a SC aqueous solution. First, the repulsive electrostatic potential energy between two parallel graphene sheets per unit area, rep, as a function of SC concentration, C, and intersheet separation, d, is given by:23  Rep (C , d )  64k BT   1 tanh 2 e0 (, C ) / k BT  exp( d )

(13)

where 0(∞,C) was obtained by solving Eq. (5) in the limiting case  = 1. Note that in deriving Eq. (13), we have assumed that graphene is a continuous conductor that can perfectly screen the electric field exerted by the EDL on the opposite side of graphene. In other words, we have assumed that graphene is “opaque” to the electric field, and Eq. (13) is therefore independent of the layer number of the graphene sheet considered. On the other hand, the attractive vdW interactions can be partially transmitted through the graphene layers,67,68 suggesting that the layer number of graphene needs to be considered to formulate the attractive vdW potential energy per unit area between two parallel graphene sheets. Based on a model which accounts for the vdW interactions between two parallel MLG sheets,70 we have derived the expression of the vdW interactions between two parallel, N-layer and M-Layer graphene sheets per unit area as a function of d, att,NM(d), as follows: N 1 M 1  i    B A   att, NM (d )   2      4 10    5(d  jd 0 )    i  0  j  i  2(d  jd 0 )

28 ACS Paragon Plus Environment

(14)

Page 29 of 49

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 d0 = 3.35 Å is the interlayer distance between two graphitic layers in an N-layer ( N  2 ) graphene sheet. In the subsequent calculations, we consider the cases of N = M = 1 (MLG), N = M = 2 (BLG), or N = M= 3 (TLG). Finally, similar to the SWCNT formulation, the total potential energy between two identical, parallel, N-layer graphene sheets per unit area, NN, is

given by:23  NN (C , d )   att, NN (d )   rep (C , d )

(15)

Following a procedure similar to the one discussed above in the SWCNT case, rep, att , and  are obtained by solving Eqs. (13), (14), and (15), respectively. In addition, the height of the energy barrier hindering the re-aggregation of two identical, N-layer graphene sheets, Barrier,NN, can be determined using the following equations:   NN     d C

 0 ;  Barrier, NN (C )  (C , dmax )

(16)

d dmax

The calculated total potential energy between two identical, parallel, MLG sheets per unit area, 11, as a function of C and d is shown in Fig. 8a, which exhibits characteristics which are very similar to those shown in Fig. 7a. The energy scale of 11 ranges from zero (blue) to 27 meV/nm (red). In particular, at a specific SC concentration, there exists a boundary layer that can hinder the re-aggregation of MLG sheets, where the boundary-layer thickness decreases significantly with C, due to the increased ionic strength in the bulk phase. Figure 8b shows the calculated Barrier,NN for MLG (N = 1), BLG (N = 2), and TLG (N = 3). Again, similar to the SWCNT case, we find that the colloidal stability of graphene sheets becomes very low when C > 200 mM, while at a relatively low SC concentration (C ~ 1 – 10 mM), an adequate colloidal stability can still be imparted. Most importantly, in the graphene case, we predict that the colloidal stability significantly decreases with the layer number, a very 29 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 30 of 49

interesting finding which has never been reported experimentally. According to our predictions, monolayer graphene can provide about a 40% higher energy barrier than multilayer graphene, thereby exhibiting better colloidal stability. It is noteworthy that our CG-MD simulation results suggest that SC adsorption is only weakly affected by the layer number (see Fig. 4b). Therefore, it would appear that the reduced colloidal stability as the layer number increases results from an increase in the vdW interactions with layer number. Since research on the stability of graphene aqueous dispersions is still at an early stage, we hope that our results will stimulate experimental work to confirm their validity.

4.

Conclusions

The mechanism of surfactant-assisted stabilization of SWCNTs and graphene sheets in an aqueous solution of the surfactant SC was investigated using large-scale CG-MD simulations. In addition, a simulation-assisted DLVO theory which quantifies the energy barrier hindering the re-aggregation of the carbon nanomaterials was proposed to elucidate their radius and surfactant concentration dependences. By carrying out large-scale CG-MD simulations, we successfully calculated the SC surface coverage on SWCNTs (6,6), (12,12), (20,20), and MLG, BLG, and TLG in SC aqueous solutions at experimentally realistic SC concentrations (C = 1 – 200 mM). Most importantly, the phenomena of SC adsorption and SC micellization were considered in one simulation box when C > CMC. According to our simulation results, for the first time, we predict that in a SC aqueous solution in the low-SC concentration regime (C < 10 mM), small-radius SWCNTs have lower SC surface coverages than large-radius SWCNTs, because their large local curvature can reduce the adsorption enthalpy change. On the other hand, in the high-SC

30 ACS Paragon Plus Environment

Page 31 of 49

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

concentration regime (C > 20 mM), small-radius SWCNTs have higher SC surface coverages than large-radius SWCNTs due to the reduced electrostatic repulsions between the negativelycharged cholate head groups of the adsorbed SC molecules. These findings provide new experimental implications to separate SWCNTs based on their radius by tuning surfactant (SC) concentration. In order to describe the CG-MD simulated surface coverages, we developed a new adsorption model which explicitly accounts for the nanomaterial radius, r, and the surfactant concentration, C. In combination with analytical solutions of the Poisson-Boltzmann equation, a continuum spectrum of the surface electric potential, 0, as a function of r and C, was generated, suggesting that separating SWCNTs based on their radius using electrophoresis may be feasible. Finally, we developed a simulation-assisted DLVO theory to quantify the potential energy barrier height required to impart colloidal stability to two identical SWCNTs or graphene sheets in a SC aqueous solution as a function of C and r. We found that, for a given chirality of SWCNT or graphene, there exists an optimal SC concentration, corresponding to the highest energy barrier height. This results from a competition between the increase in SC adsorption and the increase in the ionic strength in the bulk phase, consistent with previous experimental findings. We also found that monolayer graphene provides about a 40% higher energy barrier than multilayer graphene, thereby suggesting better colloidal stability for monolayer graphene dispersions in a SC aqueous solution. The mechanism and theoretical methodology presented here are expected to be useful in understanding the stabilization of SWCNTs and graphene sheets in surfactant aqueous solutions, as well as in allowing development of new separation methodologies aimed at preparing singlechirality carbon nanomaterials on a large scale. In addition, we hope that the molecular-level design of better surfactant or co-surfactant7,57 systems to disperse and sort SWCNT and graphene

31 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 32 of 49

dispersions will be facilitated by the fundamental principles and theoretical methods presented here.

Acknowledgements

The authors are grateful for the financial support from the National Science Foundation CBET1133813. M.S. Strano acknowledges funding from the 2009 U.S. Office of Naval Research Multi University Research Initiative (ONR-MURI) on Graphene Advanced Terahertz Engineering (GATE) at MIT, Harvard and Boston University. C.J. Shih is grateful for partial financial support from the Chyn Duog Shiah memorial Fellowship awarded by Massachusetts Institute of Technology. S. Lin is grateful for the startup funding from the Energy Materials Initiative at Florida State University.

References

(1)

Baughman, R. H.; Zakhidov, A. A.; de Heer, W. A. Carbon Nanotubes--the Route toward Applications. Science 2002, 297, 787-792. 32 ACS Paragon Plus Environment

Page 33 of 49

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

(2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18)

(19)

Geim, A. K.; Novoselov, K. S. The Rise of Graphene. Nat Mater 2007, 6, 183-191. Coleman, J. N. Liquid-Phase Exfoliation of Nanotubes and Graphene. Advanced Functional Materials 2009, 19, 3680-3695. Vaisman, L.; Wagner, H. D.; Marom, G. The Role of Surfactants in Dispersion of Carbon Nanotubes. Advances in Colloid and Interface Science 2006, 128(30), 37-46. Wang, H. Dispersing Carbon Nanotubes Using Surfactants. Current Opinion in Colloid & Interface Science 2009, 14, 364-371. Liu, H.; Nishide, D.; Tanaka, T.; Kataura, H. Large-Scale Single-Chirality Separation of Single-Wall Carbon Nanotubes by Simple Gel Chromatography. Nat Commun 2011, 2, 309. Arnold, M. S.; Green, A. A.; Hulvat, J. F.; Stupp, S. I.; Hersam, M. C. Sorting Carbon Nanotubes by Electronic Structure Using Density Differentiation. Nat Nano 2006, 1, 6065. Shih, C.-J.; Vijayaraghavan, A.; Krishnan, R.; Sharma, R.; Han, J.-H.; Ham, M.-H.; Jin, Z.; Lin, S.; Paulus, G. L. C.; Reuel, N. F.; Wang, Q. H.; Blankschtein, D.; Strano, M. S. Bi- and Trilayer Graphene Solutions. Nat Nano 2011, 6, 439-445. Ghosh, S.; Bachilo, S. M.; Weisman, R. B. Advanced Sorting of Single-Walled Carbon Nanotubes by Nonlinear Density-Gradient Ultracentrifugation. Nat Nano 2010, 5, 443450. Usrey, M. L.; Strano, M. S. Controlling Single-Walled Carbon Nanotube Surface Adsorption with Covalent and Noncovalent Functionalization. The Journal of Physical Chemistry C 2009, 113, 12443-12453. Krupke, R.; Hennrich, F.; Kappes, M. M.; v. Lohneysen, H. Surface Conductance Induced Dielectrophoresis of Semiconducting Single-Walled Carbon Nanotubes. Nano Letters 2004, 4, 1395-1399. Krupke, R.; Hennrich, F.; Lohneysen, H. v.; Kappes, M. M. Separation of Metallic from Semiconducting Single-Walled Carbon Nanotubes. Science 2003, 301, 344-347. Lustig, S. R.; Jagota, A.; Khripin, C.; Zheng, M. Theory of Structure-Based Carbon Nanotube Separations by Ion-Exchange Chromatography of DNA/Cnt Hybrids. The Journal of Physical Chemistry B 2005, 109, 2559-2566. Tu, X.; Manohar, S.; Jagota, A.; Zheng, M. DNA Sequence Motifs for Structure-Specific Recognition and Separation of Carbon Nanotubes. Nature 2009, 460, 250-253. Tvrdy, K.; Jain, R. M.; Han, R.; Hilmer, A. J.; McNicholas, T. P.; Strano, M. S. A Kinetic Model for the Deterministic Prediction of Gel-Based Single-Chirality SingleWalled Carbon Nanotube Separation. ACS Nano 2013, 7, 1779-1789. Flavel, B. S.; Kappes, M. M.; Krupke, R.; Hennrich, F. Separation of Single-Walled Carbon Nanotubes by 1-Dodecanol-Mediated Size-Exclusion Chromatography. ACS Nano 2013, 7, 3557-3564. Green, A. A.; Hersam, M. C. Solution Phase Production of Graphene with Controlled Thickness Via Density Differentiation. Nano Letters 2009, 9, 4031-4036. Lotya, M.; Hernandez, Y.; King, P. J.; Smith, R. J.; Nicolosi, V.; Karlsson, L. S.; Blighe, F. M.; De, S.; Wang, Z.; McGovern, I. T.; Duesberg, G. S.; Coleman, J. N. Liquid Phase Production of Graphene by Exfoliation of Graphite in Surfactant/Water Solutions. Journal of the American Chemical Society 2009, 131, 3611-3620. Mattevi, C.; Kim, H.; Chhowalla, M. A Review of Chemical Vapour Deposition of Graphene on Copper. Journal of Materials Chemistry 2011, 21, 3324-3334. 33 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

(20) (21) (22)

(23) (24) (25) (26) (27)

(28)

(29) (30)

(31) (32) (33)

Page 34 of 49

Park, H.; Afzali, A.; Han, S.-J.; Tulevski, G. S.; Franklin, A. D.; Tersoff, J.; Hannon, J. B.; Haensch, W. High-Density Integration of Carbon Nanotubes Via Chemical SelfAssembly. Nat Nano 2012, 7, 787-791. Vijayaraghavan, A.; Sciascia, C.; Dehm, S.; Lombardo, A.; Bonetti, A.; Ferrari, A. C.; Krupke, R. Dielectrophoretic Assembly of High-Density Arrays of Individual Graphene Devices for Rapid Screening. ACS Nano 2009, 3, 1729-1734. Strano, M. S.; Moore, V. C.; Miller, M. K.; Allen, M. J.; Haroz, E. H.; Kittrell, C.; Hauge, R. H.; Smalley, R. E. The Role of Surfactant Adsorption During Ultrasonication in the Dispersion of Single-Walled Carbon Nanotubes. Journal of Nanoscience and Nanotechnology 2003, 3, 81-86. Israelachvili, J. N. Intermolecular and Surface Forces; 3rd ed.; Academic Press: Burlington, MA, 2011. White, B.; Banerjee, S.; O'Brien, S.; Turro, N. J.; Herman, I. P. Zeta-Potential Measurements of Surfactant-Wrapped Individual Single-Walled Carbon Nanotubes. The Journal of Physical Chemistry C 2007, 111, 13684-13690. Hertel, T.; Hagen, A.; Talalaev, V.; Arnold, K.; Hennrich, F.; Kappes, M.; Rosenthal, S.; McBride, J.; Ulbricht, H.; Flahaut, E. Spectroscopy of Single- and Double-Wall Carbon Nanotubes in Different Environments. Nano Letters 2005, 5, 511-514. Moore, V. C.; Strano, M. S.; Haroz, E. H.; Hauge, R. H.; Smalley, R. E.; Schmidt, J.; Talmon, Y. Individually Suspended Single-Walled Carbon Nanotubes in Various Surfactants. Nano Letters 2003, 3, 1379-1382. Sun, Z.; Nicolosi, V.; Rickard, D.; Bergin, S. D.; Aherne, D.; Coleman, J. N. Quantitative Evaluation of Surfactant-Stabilized Single-Walled Carbon Nanotubes: Dispersion Quality and Its Correlation with Zeta Potential. The Journal of Physical Chemistry C 2008, 112, 10692-10699. Gubitosi, M.; Trillo, J. V.; Alfaro Vargas, A.; Pavel, N. V.; Gazzoli, D.; Sennato, S.; Jover, A.; Meijide, F.; Galantini, L. Characterization of Carbon Nanotube Dispersions in Solutions of Bile Salts and Derivatives Containing Aromatic Substituents. The Journal of Physical Chemistry B 2014, 118, 1012-1021. Smith, R. J.; Lotya, M.; Coleman, J. N. The Importance of Repulsive Potential Barriers for the Dispersion of Graphene Using Surfactants. New Journal of Physics 2010, 12, 125008. O'Connell, M. J.; Bachilo, S. M.; Huffman, C. B.; Moore, V. C.; Strano, M. S.; Haroz, E. H.; Rialon, K. L.; Boul, P. J.; Noon, W. H.; Kittrell, C.; Ma, J.; Hauge, R. H.; Weisman, R. B.; Smalley, R. E. Band Gap Fluorescence from Individual Single-Walled Carbon Nanotubes. Science 2002, 297, 593-596. Matarredona, O.; Rhoads, H.; Li, Z.; Harwell, J. H.; Balzano, L.; Resasco, D. E. Dispersion of Single-Walled Carbon Nanotubes in Aqueous Solutions of the Anionic Surfactant Naddbs. The Journal of Physical Chemistry B 2003, 107, 13357-13367. Islam, M. F.; Rojas, E.; Bergey, D. M.; Johnson, A. T.; Yodh, A. G. High Weight Fraction Surfactant Solubilization of Single-Wall Carbon Nanotubes in Water. Nano Letters 2003, 3, 269-273. Richard, C.; Balavoine, F.; Schultz, P.; Ebbesen, T. W.; Mioskowski, C. Supramolecular Self-Assembly of Lipid Derivatives on Carbon Nanotubes. Science 2003, 300, 775-778.

34 ACS Paragon Plus Environment

Page 35 of 49

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

(34) (35) (36) (37) (38) (39)

(40) (41) (42) (43) (44) (45)

(46) (47) (48) (49)

Yurekli, K.; Mitchell, C. A.; Krishnamoorti, R. Small-Angle Neutron Scattering from Surfactant-Assisted Aqueous Dispersions of Carbon Nanotubes. Journal of the American Chemical Society 2004, 126, 9902-9903. Angelikopoulos, P.; Bock, H. Directed Self-Assembly of Surfactants in Carbon Nanotube Materials. The Journal of Physical Chemistry B 2008, 112, 13793-13801. Calvaresi, M.; Dallavalle, M.; Zerbetto, F. Wrapping Nanotubes with Micelles, Hemimicelles, and Cylindrical Micelles. Small 2009, 5, 2191-2198. Blanch, A. J.; Lenehan, C. E.; Quinton, J. S. Optimizing Surfactant Concentrations for Dispersion of Single-Walled Carbon Nanotubes in Aqueous Solution. The Journal of Physical Chemistry B 2010, 114, 9805-9811. Angelikopoulos, P.; Gromov, A.; Leen, A.; Nerushev, O.; Bock, H.; Campbell, E. E. B. Dispersing Individual Single-Wall Carbon Nanotubes in Aqueous Surfactant Solutions Below the Cmc. The Journal of Physical Chemistry C 2009, 114, 2-9. Utsumi, S.; Kanamaru, M.; Honda, H.; Kanoh, H.; Tanaka, H.; Ohkubo, T.; Sakai, H.; Abe, M.; Kaneko, K. Rbm Band Shift-Evidenced Dispersion Mechanism of Single-Wall Carbon Nanotube Bundles with Naddbs. Journal of Colloid and Interface Science 2007, 308, 276-284. Duque, J. G.; Densmore, C. G.; Doorn, S. K. Saturation of Surfactant Structure at the Single-Walled Carbon Nanotube Surface. Journal of the American Chemical Society 2010, 132, 16165-16175. Bouchard, D.; Zhang, W.; Powell, T.; Rattanaudompol, U. s. Aggregation Kinetics and Transport of Single-Walled Carbon Nanotubes at Low Surfactant Concentrations. Environmental Science & Technology 2012, 46, 4458-4465. Niyogi, S.; Densmore, C. G.; Doorn, S. K. Electrolyte Tuning of Surfactant Interfacial Behavior for Enhanced Density-Based Separations of Single-Walled Carbon Nanotubes. Journal of the American Chemical Society 2008, 131, 1144-1153. Xu, Z.; Yang, X.; Yang, Z. A Molecular Simulation Probing of Structure and Interaction for Supramolecular Sodium Dodecyl Sulfate/Single-Wall Carbon Nanotube Assemblies. Nano Letters 2010, 10, 985-991. McDonald, T. J.; Engtrakul, C.; Jones, M.; Rumbles, G.; Heben, M. J. Kinetics of Pl Quenching During Single-Walled Carbon Nanotube Rebundling and Diameter-Dependent Surfactant Interactions. The Journal of Physical Chemistry B 2006, 110, 25339-25346. Duque, J. G.; Parra-Vasquez, A. N. G.; Behabtu, N.; Green, M. J.; Higginbotham, A. L.; Price, B. K.; Leonard, A. D.; Schmidt, H. K.; Lounis, B.; Tour, J. M.; Doorn, S. K.; Cognet, L.; Pasquali, M. Diameter-Dependent Solubility of Single-Walled Carbon Nanotubes. ACS Nano 2010, 4, 3063-3072. Tummala, N. R.; Striolo, A. Sds Surfactants on Carbon Nanotubes: Aggregate Morphology. ACS Nano 2009, 3, 595-602. Tummala, N. R.; Striolo, A. Curvature Effects on the Adsorption of Aqueous SodiumDodecyl-Sulfate Surfactants on Carbonaceous Substrates: Structural Features and Counterion Dynamics. Physical Review E 2009, 80, 021408. Qiao, R.; Ke, P. C. Lipid-Carbon Nanotube Self-Assembly in Aqueous Solution. Journal of the American Chemical Society 2006, 128, 13656-13657. Wallace, E. J.; Sansom, M. S. P. Carbon Nanotube/Detergent Interactions Via CoarseGrained Molecular Dynamics. Nano Letters 2007, 7, 1923-1928.

35 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

(50) (51)

(52)

(53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66)

Page 36 of 49

Wallace, E. J.; Mark, S. P. S. Carbon Nanotube Self-Assembly with Lipids and Detergent: A Molecular Dynamics Study. Nanotechnology 2009, 20, 045101. Lin, S.; Blankschtein, D. Role of the Bile Salt Surfactant Sodium Cholate in Enhancing the Aqueous Dispersion Stability of Single-Walled Carbon Nanotubes: A Molecular Dynamics Simulation Study. The Journal of Physical Chemistry B 2010, 114, 1561615625. Lin, S.; Shih, C.-J.; Strano, M. S.; Blankschtein, D. Molecular Insights into the Surface Morphology, Layering Structure, and Aggregation Kinetics of Surfactant-Stabilized Graphene Dispersions. Journal of the American Chemical Society 2011, 133, 1281012823. Carvalho, E. J. F.; dos Santos, M. C. Role of Surfactants in Carbon Nanotubes Density Gradient Separation. ACS Nano 2010, 4, 765-770. Quintilla, A.; Hennrich, F.; Lebedkin, S.; Kappes, M. M.; Wenzel, W. Influence of Endohedral Water on Diameter Sorting of Single-Walled Carbon Nanotubes by Density Gradient Centrifugation. Physical Chemistry Chemical Physics 2010, 12, 902-908. Sanders, S. A.; Sammalkorpi, M.; Panagiotopoulos, A. Z. Atomistic Simulations of Micellization of Sodium Hexyl, Heptyl, Octyl, and Nonyl Sulfates. The Journal of Physical Chemistry B 2012, 116, 2430-2437. Garidel, P.; Hildebrand, A.; Neubert, R.; Blume, A. Thermodynamic Characterization of Bile Salt Aggregation as a Function of Temperature and Ionic Strength Using Isothermal Titration Calorimetry. Langmuir 2000, 16, 5267-5275. Hersam, M. C. Progress Towards Monodisperse Single-Walled Carbon Nanotubes. Nat Nano 2008, 3, 387-394. Van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. C. Gromacs: Fast, Flexible, and Free. Journal of Computational Chemistry 2005, 26, 17011718. Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; de Vries, A. H. The Martini Force Field:Coarse Grained Model for Biomolecular Simulations. The Journal of Physical Chemistry B 2007, 111, 7812-7824. Yesylevskyy, S. O.; Schäfer, L. V.; Sengupta, D.; Marrink, S. J. Polarizable Water Model for the Coarse-Grained Martini Force Field. PLoS Comput Biol 2010, 6, e1000810. Marrink, S. J. Molecular Dynamics Simulations of Cholesterol Nucleation in Mixed Micelles Modelling Human Bile. Proceedings of Falk Symposium 2004, 139, 98 - 115. Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. The Journal of Chemical Physics 2007, 126, 014101. Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular Dynamics with Coupling to an External Bath. The Journal of Chemical Physics 1984, 81, 3684-3690. LeBard, D. N.; Levine, B. G.; Mertmann, P.; Barr, S. A.; Jusufi, A.; Sanders, S.; Klein, M. L.; Panagiotopoulos, A. Z. Self-Assembly of Coarse-Grained Ionic Surfactants Accelerated by Graphics Processing Units. Soft Matter 2012, 8, 2385-2397. Reis, S.; Moutinho, C. G. e.; Matos, C.; de Castro, B.; Gameiro, P.; Lima, J. L. F. C. Noninvasive Methods to Determine the Critical Micelle Concentration of Some Bile Acid Salts. Analytical Biochemistry 2004, 334, 117-126. Shinoda, W.; DeVane, R.; Klein, M. L. Coarse-Grained Molecular Modeling of Non-Ionic Surfactant Self-Assembly. Soft Matter 2008, 4, 2454-2462. 36 ACS Paragon Plus Environment

Page 37 of 49

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

(67) (68) (69)

(70)

Shih, C.-J.; Wang, Q. H.; Lin, S.; Park, K.-C.; Jin, Z.; Strano, M. S.; Blankschtein, D. Breakdown in the Wetting Transparency of Graphene. Physical Review Letters 2012, 109, 176101. Shih, C.-J.; Strano, M. S.; Blankschtein, D. Wetting Translucency of Graphene. Nat Mater 2013, 12, 866-869. Ohshima, H.; Healy, T. W.; White, L. R. Accurate Analytic Expressions for the Surface Charge Density/Surface Potential Relationship and Double-Layer Potential Distribution for a Spherical Colloidal Particle. Journal of Colloid and Interface Science 1982, 90, 1726. Girifalco, L. A.; Hodak, M.; Lee, R. S. Carbon Nanotubes, Buckyballs, Ropes, and a Universal Graphitic Potential. Physical Review B 2000, 62, 13104-13110.

Table

37 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

Table 1.

Page 38 of 49

Extracted values of max and Gads for the SWCNTs and graphene sheets

considered in Fig. 4. SWNT (6,6)

SWNT (12,12)

SWNT (20,20)

MLG

BLG

TLG

max [nm-2]

2.75

1.77

1.47

0.95

0.98

0.99

Gads [kBT]

-8.59

-8.46

-8.39

-8.98

-8.37

-8.19

Figure Captions Figure 1.

Coarse-grained MD simulations of the stabilization of a carbon nanomaterial in a

surfactant (sodium cholate - SC) aqueous solution accounting simultaneously for the phenomena of SC micellization and SC adsorption in one simulation box. (a) Chemical structure of SC showing the rigid, steroid-ring backbone of this facial surfactant, (b) AA molecular model of SC showing its hydrophobic and hydrophilic faces , (c) CG molecular model of SC showing the grouped beads used in our simulations. (d) Representative post-equilibrium simulation snapshot (top) and the simulated concentration profile (bottom) of a SC-coated MLG sheet. Color code for AA model: red – oxygen, cyan – carbon, green – sodium, and white – hydrogen. Color code for the CG model: purple – PW.

Figure 2.

Representative post-equilibrium simulation snapshots (side and front views) for

SWCNTs (a) (6,6), (b) (12,12), and (c) (20,20) at C = 10 mM (left) and 60 mM (right).

Figure 3.

Micellization of sodium cholate molecules in the bulk phase. (a) Simulated free

monomer concentration, C1, as a function of C (the blue circles). The solid and dashed lines

38 ACS Paragon Plus Environment

Page 39 of 49

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

correspond to C1 = C and to the experimental CMC, respectively. (b) Simulated aggregation number distributions for C = 10 mM and 200 mM.

Figure 4.

Adsorption of sodium cholate molecules on SWCNTs and graphene sheets. The

various symbols denote the simulated SC surface coverages on single, isolated (a) SWCNT (6,6), SWCNT (12,12), and SWCNT (20,20), and (b) MLG, BLG, and TLG sheets as a function of C, and the various lines are the fitted adsorption isotherms using Eq. (1), where the colors used for each symbol and line match.

Figure 5.

Preferential adsorption of cholate molecules on small-radius SWCNTs in the

high-SC concentration regime (C > 20 mM). (a) Schematic illustration of the reduction of electrostatic repulsions between the adsorbed SC molecules on a small-radius SWCNT. (b) Fitted values of saturated surface coverage, max, from the CG-MD simulations as a function of 1/r for the carbon nanomaterials considered (symbols) and their linear fit (line).

Figure 6.

Calculated surface electrical potential, 0, as a function of 1/r and C. The solid

white line is the calculated optimal concentration maximizing 0, Copt, as a function of 1/r, the dashed black lines correspond to the different carbon nanomaterials considered, and the white dashed line denotes the experimental CMC.

Figure 7.

Simulation-Assisted DLVO theory accounting for the interactions between two

identical, parallel SWCNTs. (a) Calculated total potential energy between two SWCNTs (6,6) per unit length as a function of intertube separation, d, and C. (b) Calculated energy barrier

39 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 40 of 49

height hindering the re-aggregation of SWCNTs, Barrier, as a function of r and C. The solid white line corresponds to the calculated Copt as a function of r.

Figure 8.

Simulation-Assisted DLVO theory accounting for the interactions between two

identical, parallel graphene sheets. (a) Calculated total potential energy between two MLG sheets per unit area as a function of intersheet separation, d, and C. (b) Calculated energy barrier height hindering the re-aggregation of the graphene sheets, Barrier, for MLG, BLG, and TLG as a function of C.

40 ACS Paragon Plus Environment

Page 41 of 49

a

b

c Hydrophobic Side

Hydrophilic Side

d

z

20

C [nm-3]

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

Bulk Region

Adsorption Region

Bulk Region

10

0 -10

-5

0

5

10

z [nm]

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

a

b SWCNT (6,6) 10 mM

60 mM

Page 42 of 49

c SWCNT (12,12) 10 mM

60 mM

SWCNT (20,20) 10 mM

60 mM

Fig. 2 ACS Paragon Plus Environment

Page 43 of 49

a Exp. CMC*

C1 [mM]

100

CG MD

10

1 1

10 C [mM]

b

100

100

Percentage [%]

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

10 C = 200 mM

1 C = 10 mM

0.1 0.01

1

2

3

4

n

5

6

7

8

Fig. 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

a

3

 [nm−2]

(6,6)

2

(12,12)

1

(20,20)

0 1

b

 [nm−2]

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 49

10 C [mM]

100

3

2

MLG BLG TLG

1

0 1

10 C [mM]

100

Fig. 4

ACS Paragon Plus Environment

Page 45 of 49

a Small r

r

Large r

b 3

max [nm-2]

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

(6,6) (12,12)

2 (20,20) TLG

1

BLG MLG Linear Fit

0 0

1 1/r [nm-1]

2

Fig. 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

-80

0 [mV]

C [mM]

(12,12) SWCNT

(20,20) SWCNT

100

-160 (6,6) SWCNT

1000

Graphene

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 46 of 49

Copt

10

Exp. CMC

1 0

1

2

1/r

3

4

[nm-1]

Fig. 6 ACS Paragon Plus Environment

Page 47 of 49

a 0

 [meV/nm]

27

1000 r

r d SWCNT (6,6)

C [mM]

100

Exp. CMC

10

1

3

20

40

60

80

100

d [Å]

b 0

Barrier [meV/nm]

50

1000

C [mM] 10

(20,20)

(12,12)

100

(6,6)

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

Exp. CMC

1 5

10

15

20

r [Å]

Fig. 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

a

0

 [meV/nm2]

10

1000

MLG

C [mM]

100

Exp. CMC

10

1 4

20

40

60

80

100

d [Å]

b MLG

barrier,NN [meV/nm2]

10

BLG TLG 5 Exp. CMC

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 48 of 49

0 1

10

C [mM]

100

1000

Fig. 8 ACS Paragon Plus Environment

Page 49 of 49

TOC Figure

-80

0 [mV]

C [mM]

(20,20) SWCNT

(12,12) SWCNT

100

-160

(6,6) SWCNT

1000

Graphene

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

Copt

10

Exp. CMC

1 0

1

2

3

4

1/r [nm-1]

ACS Paragon Plus Environment