Modeling the Morphology and Phase Behavior of One-Component

Nov 17, 2017 - Polymer-grafted or “hairy” nanoparticles (HNP) represent an important and relatively new class of materials. Traditionally, polymer...
6 downloads 8 Views 4MB Size
Article Cite This: Macromolecules XXXX, XXX, XXX-XXX

pubs.acs.org/Macromolecules

Modeling the Morphology and Phase Behavior of One-Component Polymer-Grafted Nanoparticle Systems Valeriy V. Ginzburg* Materials Science and Engineering, The Dow Chemical Company, Building 1702, Midland, Michigan 48674, United States S Supporting Information *

ABSTRACT: Polymer-grafted or “hairy” nanoparticles (HNP) represent an important and relatively new class of materials. Traditionally, polymers or oligomers were grafted onto the particle (silica, metal, or semiconductor) surface to improve the dispersion of particles in a polymer matrix. Recently, the scope of research has broadened substantially, as it was shown that such nanoparticles can form anisotropic structures ranging from self-assembled wires to sheets to networks. Furthermore, it has been demonstrated that one could make hybrid polymer−inorganic materials with HNPs alone, without using a separate matrix polymer. Such one-component hybrid materials are not prone to macroscopic phase separation and can, in principle, have a variety of interesting microphase-separated, anisotropic morphologies, similar to surfactants or block copolymers. Here, we develop a new selfconsistent field theory describing the behavior of such one-component HNP systems and apply it to predict morphology as a function of the ligand molecular weight and grafting density. As in the case of block copolymers, we observe lamellar, cylindrical, and spherical morphologies and elucidate phase boundaries as a function of the core (nanoparticle) volume fraction and the ratio of the particle radius to the ligand radius of gyration. We also observe the formation of a novel phase, labeled as “sheets”, where the lamellar-like ordering of particle-rich and ligand-rich layers is additionally characterized by the hexagonal ordering of the particles within the particle-rich layer. Our theoretical approach can be easily extended to other HNPs, including those with mixed ligands and block copolymer ligands.



INTRODUCTION In recent years, polymer-grafted nanoparticles (PGNP), alternatively described as polymer-brushed nanoparticles (PBNP) or hairy nanoparticles (HNP), have been the subject of increased interest in the polymer science field.1−5 Initially, grafted ligands were primarily used as tools to improve nanoparticle dispersion in polymer matrices.6−11 Thus, grafting organic ligands on the particles became necessary to form functional nanocomposites in various applications, from quantum dot (QD) containing resins to polymer−clay thermoplastic composites. However, several years ago, it was shown that the behavior of HNPs in polymers was much richer than originally thought, with anisotropic self-assembly caused by a complicated interplay between the enthalpy of particle− polymer repulsion and the entropy of ligand and matrix chain conformational arrangements.4,12−16 The ligand “cloud” around the nanoparticleeven if originally isotropiccan “polarize” in the presence of other HNPs, forming structures such as “wires” and “sheets”, in addition to more familiar morphologies such as uniformly dispersed particles or, at the other extreme, large aggregates completely phase-separated from the matrix. In principle, it is not even necessary to have a separate polymer matrix. With the versatility of chemical “grafting to” and “grafting from” methods, one can create HNPs with any polymer-to-filler volume ratio, from very high (with higher© XXXX American Chemical Society

molecular-weight ligands) to very low (with low-molecularweight ligand and/or very low grafting density). By making HNPs in solution and then casting films out of those solutions, one can create “one-component nanocomposites”, also referred to as “nanoparticle organic hybrid materials” (NOHM).5,17−23 Since these materials are one-component, theylike block copolymersshould generally exhibit microphase separation, as opposed to macrophase separation. Changing the ligand characteristics (e.g., using block copolymer ligands or mixture of ligands of various molecular weights and chemical compositions) could expand the complexity of possible morphologies, while changing the particle nature (metal, semiconductor, or polymeric) could determine which application the composite would be best suited for. Thus, HNPs could become ideal “building blocks” for versatile nanostructure engineering, as envisioned, e.g., by Cheng et al.24−27 and Glotzer et al.28−33 Experimentally, researchers have observed three-dimensional crystals such as face-centered cubic (FCC)19 as well as local anisotropic ordering5 and disordered liquid.23 Other morphologies likely can be found as well; to understand them, we need the right theoretical description. Received: September 5, 2017 Revised: October 31, 2017

A

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules Theory and modeling of polymer/nanoparticle composites has been a very active field for several decades.4,7,34−40 Theoretical and computational methods used include strong segregation theory,16 Monte Carlo simulations,16 molecular dynamics, 15 polymer reference interaction site model (PRISM),35,41−45 and self-consistent field theory (SCFT) based approaches.46−49 For the pure HNP systems and their solutions (“giant surfactants”), particle-based simulations of Glotzer and co-workers offered many interesting insights into their morphologies and phase behavior.11,31,32,50 The transition between FCC crystal and disordered liquid in one-component HNP systems was investigated by Choi et al.23 using scaling theory and by Koch and co-workers51−55 by means of density functional theory (DFT). Chremos and co-workers modeled the HNP behavior using coarse-grained molecular dynamics (CG-MD).56−59 They also compared CG-MD with DFT54 and found that in DFT (and similar field-based theories) the grafted brush stretching is somewhat underestimated. Jayaraman and Schweizer used PRISM approach to compute pair correlation functions and determine microphase separation boundaries for one- and multitethered nanoparticles.43−45 Lee et al.60 utilized the self-consistent-field/density functional theory (SCF-DFT) approach61,62 to predict microphase-separated lamellar and hexagonal morphologies of “tadpole” (coil−sphere) diblock copolymers (one-tethered nanoparticles). A slightly different version of the SCF-DFT approach was developed by Hur et al.63 for the case of multiple ligands. Here, we adapt the SCFDFT approach of refs 60−62 to model the phase behavior of one-component HNP assemblies where the nanoparticles can have arbitrary number of grafted ligands. We use this model to investigate how the HNP morphologies depend on the ligand grafting density and chain length and sketch a simple HNP phase map.

Q HNP =

1 V

(2a)

∫ dr′ ρP (r′)Θ(2RP − |r − r′|)

(2b)

∫ dr ρP (r)

(2c)

4x − 3x 2 (1 − x)2

(2d)

∫ dr (ZL(r))n exp[−wP(r)]

(3)

where ZL(r) =

∫ dr′ q(r′, 1)δ(|r − r′| − RP) (4)

A

Here A is the surface area of the particle, and δ(x) is the standard (one-dimensional) Dirac delta function. The propagator q(r,s) is computed by solving the so-called modified diffusion equation (MDE) ∂q(r, s) = (R g 2∇2 − wL(r))q(r, s) ∂s

(5) Na 2

with the “initial condition” q(r,0) = 1, and R g 2 = 6 , where a is the Kuhn length. The index s runs from 0 to 1, with 0 corresponding to the free end and 1 corresponding to the bound end; thus, q(r,1) is related to the probability of finding the bound end of the ligand at point r. The integral in eq 4 is the sum over all configurations of a ligand grafted on a surface of a particle whose center is located at point r. Equations 15 describe the HNP system. We note here that it is straightforward to extend this formalism to other HNPbased systems (e.g., mixtures of HNPs with homopolymers, mixtures of HNPs with block copolymers, HNPs with block copolymer ligands, HNPs with mixed homopolymer ligands resulting in patchy particles, etc.) For the time being, we concentrate on the HNPs with homopolymer ligands. To calculate the morphologies of these systems, we need to find the saddle point of the free energy functional with respect to the densities φL(r) and ρP(r), as well as the conjugate fields wL(r) and wP(r), and the Lagrange multiplier (“pressure”) field ξ(r). The self-consistency equations, again, are very similar to those obtained by Thompson et al.61,62

⎛Q ⎞ FNv = −nHNP ln⎜ HNP ⎟ ⎝ V ⎠ kBTV 1 + dr[(χN ){φL(r) − ϕL}{φP (r) − ϕP } V







1 2 Nv 3

∫ dr′ ρP (r′)Θ(RP − |r − r′|)

Here Θ(x) is the Heaviside step function: Θ(x) = 0, if x < 0, Θ(0) = 0.5, and Θ(x) = 1, if x > 0. The particle center probability, ρP(r), and the overall number of particles per unit volume, nHNP, are multiplied by the ligand volume, Nv, following the notation of the earlier papers. The Flory− Huggins parameter χ describes the local interaction between the ligand and the particle species. Finally, we need to express the partition function of the hairy nanoparticles, QHNP, as functional of the fields wL and wP. For that we assume that the n ligands grafted on the particle surface are mobile in a sense that their first segments are able to travel freely on the nanoparticle surface but unable to leave it. In that case, the partition function can be written as

)

− ξ(r)(1 − φL(r) − φP (r))] 1 − dr[wP(r)ρP (r) + wL(r)φL(r)] V 1 + dr[ρP (r)ΨHS(φP̿ (r))] V

φP̿ (r) =

ΨHS(x) =

METHOD AND IMPLEMENTATION We start from a generic version of the SCF-DFT mean-field free energy similar to that of Thompson et al.60−62 The system we consider is a one-component “melt” of polymer-grafted or “hairy” nanoparticles (HNP), occupying total volume V. Each HNP consists of the “core” (spherical nanoparticle with radius 4 RP and volume VP = 3 πRP 3) and the “shell” (n homopolymer ligands with chain length N). The particle species is denoted as P and ligand as L. The overall volume of each HNP is given by VHNP = VP + nNv, where v is the monomer (repeat unit) volume. The overall particle volume fraction, 4 4 ϕP = 3 πRP 3 / 3 πRP 3 + nNv , and the overall ligand volume fraction ϕL = 1 − ϕP. The free energy is then written as46,61,62

)(

1 Nv

nHNP =



(

φP (r) =

(1)

wL(r) = (χN ){φP (r) − ϕP } + ξ(r)

Here, the following definitions are used: B

(6a)

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules = (r)) + wP(r) = ΨCS(φ P

1 N

⎧ V [φL(r)]new = [φL(r)]old + λ⎨ϕL ⎩ QL

∫ dr′[(χN ){φL(r′) − ϕL}





1 dr′ + ξ(r′)]Θ(RP − |r − r′|)+ 3 2 Nv = (r′))]Θ(2R − |r − r′|) × [ρP (r′)Ψ′CS(φ P P



ρP (r) = nHNP

φL(r) = ϕL

V Q HNP

V QL

∫0

1

× (6b)

(6d) (6e)

The local particle density φP(r) is computed by integrating particle center probability ρP(r) according to eq 2a. To complete the set of equations, we need to determine the ligand partition function QL and the conjugate propagator q†L(r,s). The single-ligand partition function is given by

∫ qL†(r,0) dr

The propagator equation ∂qL†(r,s) ∂s

(7)

is described by the modified diffusion

= −(R g 2∇2 − wL(r))qL†(r,s)

(8a)

with the “initial condition” reflecting the fact that the “bound” end of each ligand must be on the particle surface which it shares with (n − 1) other ligands

∫ dr′(ZL(r′))n−1 exp[−wP(r′)]δ(|r − r′| − RP)

qL†(r,1) = A−1

(8b)

In the limit of an infinitely small particle (RP → 0), eq 8b would reduce to the familiar expression for the propagator of an n-arm star.64 To find numerical solutions of the self-consistency equations, we discretize them on a cubic grid, with mesh size dx. The Θ and δ functions are approximated as follows: Θ(RP − |r|) ≈

δ(|r| − RP) ≈

⎛ |r| − RP ⎞⎤ 1⎡ ⎟⎥ ⎢1 − tanh⎜K ⎝ 2⎣ d x ⎠⎦

(9a)

⎧ ⎛ |r| − R ⎞2 ⎫ K P ⎟⎬ exp⎨−K ⎜ π dx ⎩ ⎝ dx ⎠ ⎭ ⎪







χ=v

(9b)

⎧ V [ρP (r)]new = [ρP (r)]old + λP ⎨nHNP Q ⎩ HNP ⎪











(10b) (10c)

( eP − eL )2 (δP − δL)2 =v RT RT

(11)

where eP and eL denote cohesive energy densities (CED, usually expressed in units of J/cm3) of particle and ligand, respectively, and δL , P = eL , P are the respective solubility parameters. Cohesive energy densities of inorganic particles are much greater than those for typical polymers, and thus the polymer− particle interaction parameters are expected to be large (see Supporting Information). As χ is increased, the segregation between each particle and its surroundings becomes stronger, even though at high particle volume fractions, particles can still “fuse” into wires or sheets.46 Here, we used χ chosen in a way that for all compositions χN is much greater than the corresponding χNODT for flexible block copolymers (see

where K is sufficiently large to ensure that the transition from the “inside” to the “outside” occurs over a single layer of grid points (in our calculations, we used K = 1.5). We solve eqs 6a−6e in an iterative manner, using Picard algorithm to update the densities, as follows:

⎫ × exp[−wP(r)](ZL(r)n ) − [ρP (r)]old ⎬ ⎭

⎫ qL(r, s)qL†(r,s) ds − [φL(r)]old ⎬ ⎭

After updating the densities and pressures according to eqs 10a−10c, we then compute the new fields wP and wL according to eqs 6a and 6b and repeat the process until convergence. The Rasmussen−Kalosakas65 pseudospectral approach is employed to speed up the calculations. The use of pseudospectral method and utilizing Fourier transforms for the computation of convolution integrals makes it possible to solve the ligand propagator MDE (eqs 5, 8a, and 8b) relatively fast, so that all computations can be performed within hours (exact time depends on the specific hardware, size of the simulation box, and the number of discretization points used to describe the ligand chains). We use convergence criterion ⟨|φL(r) + φP(r) − 1|⟩ < 10−3. Other simulation parameters are as follows: λ = 0.025, λP = 0.007, and ε = 30.0. Although the approach to finding the locally equilibrium morphology in SCF-DFT is similar to that in a pure SCFT, there are some notable challenges specific to hybrid nanoparticle−polymer systems. In SCF-DFT, particles are not introduced as solid objects whose interior is impenetrable to the polymers; rather, regions with high φP (“particles”) and low φP (“polymers”) emerge from the calculations. The resulting morphology is usually more or less smeared, representing ensemble-averaged particle density profiles, rather than instantaneous configurations. While so-called “hybrid” models66−70 preserve the particle shape and disallow particle− particle or particle−polymer overlap at any given moment, they often suffer from significant perturbations in the polymer field when particles are moved and require lengthy equilibration process. Here, we attempt to utilize SCF-DFT while minimizing the particle deformation and overlap. We follow the ideas of Chen et al.,71 who proposed to add a “particle cohesive energy” term in the free energy. The addition of this term results in sharpening the peaks of ρP(r) and preserving the spherical shapes of individual particles. For the case of a twospecies system, increasing the cohesive energy density of the particle relative to the polymer is equivalent to increasing the Flory−Huggins parameter. Indeed, the Flory−Huggins parameter χ is related to the cohesive energy densities of the polymer and the particle via72

(6c)

qL(r, s)qL†(r,s) ds

q†L(r,s)

1

[ξ(r)]new = [ξ(r)]old + λε{φL(r) + φP (r) − 1}

exp[−wP(r)](ZL(r))n

φL(r) + φP (r) = 1

QL =

∫0

(10a) C

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules

Figure 1. Various HNPs having the same core volume fraction (φP) of 0.5. Left to right: one grafted ligand with N = 120; two grafted ligands with N = 60; four grafted ligands with N = 30; and eight grafted ligands with N = 15.

Figure 2. (a) Three-dimensional view of lamellar morphology, with red representing particle-rich layers. The surfaces correspond to ϕ = 0.5. The blue arrow goes from the center of the particle-rich layer to the center of the ligand layer. (b) Density profiles of particle (red) and ligand (blue) species along the direction of the blue arrow. Here, n = 8 and N = 15.

further analysis, we also define three important parameters: the n grafting density σ = 4πR 2 , the ratio of the particle radius to the

Supporting Information). Thus, we expect that all calculations below correspond to strong segregation behavior and the morphology is then χN-independent. The calculation of the phase diagram is done using a realspace approach (see, e.g., Drolet and Fredrickson73), with the Rasmussen−Kalosakas65 pseudospectral approach to solving the modified diffusion equation (MDE). In this initial study, we considered three classical morphologies ubiquitous in block copolymers, i.e., spherical (cubic), cylindrical (hexagonal), and lamellar. For each composition, we computed the free energy of the face-centered cubic (FCC) structure (if stable) and also found the free energy minima of the cylindrical and lamellar structures with respect to the unit cell dimensions (for more details, see the Supporting Information). Interestingly, we found that in addition to the classical phases there is also a “layered sheet” phase, similar to lamellar but with 2D-hexagonal ordering of the particles within the particle-rich layer. The detailed description of those morphologies is given below.



P

(unperturbed) radius of gyration of a single ligand chain, RP/Rg, and the volume fraction of the nanoparticle core, ϕP. The first two parameters are similar to those used in studies of HNPs in polymer matrices and solvents,15,16 while the third one is more akin to the composition parameter f in the theory of block copolymers.74 In choosing the parameter space for this investigation, we concentrated on the region of intermediate particle volume fractions, 0.2 < ϕP < 0.5. It is known that in colloid−polymer mixtures, at low ϕP, particles are dispersed uniformly, while at high ϕP, they form either gel or crystal, depending on the details of the polymer−particle interactions.75−78 In the case where there is a strong depletion attraction between the particles, there is also a broad two-phase coexistence region between particle-rich and particle-poor phases; however, such a macrophase separation is not possible for the single-component HNP system, as discussed above. Thus, we expect that the HNP material is in a homogeneous liquid state at low ϕP and in a gel or glassy state at high ϕP. At intermediate ϕP, formation of periodic soft-crystalline anisotropic structures is possible, and it is this region we would like to probe with SCF-DFT calculations. This expectation is reinforced by recent theoretical

RESULTS AND DISCUSSION

In this study, we fix the nanoparticle radius RP = 3.05 nm, so the nanoparticle volume VP = 120 nm3. The monomer volume is set to v = a3 = 1 nm3. The number of ligands per particle, n, and the ligand chain length, N, are varied as shown in Figure 1 and tabulated in Table S4 (Supporting Information). For D

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules

Figure 3. (a) Three-dimensional view of the “sheet” structure. The surfaces correspond to ϕ = 0.5. (b) Cross section through the center of the particle-rich domainthe color map representing the particle density in the XY-plane. Here, n = 10 and N = 30.

Figure 4. (a) Three-dimensional view of cylindrical morphology, with red representing particle-rich regions. The surfaces correspond to ϕ = 0.75. The blue arrow goes from the center of a cylinder to the point halfway between two adjacent cylinders. (b) Density profiles of particle (red) and ligand (blue) species along the direction of the blue arrow. Here, n = 11 and N = 15.

At lower grafting densities (higher particle volume fraction), nanoparticles attempt to reduce the overall interfacial energy of the system by organizing themselves into two-dimensional sheets or lamellae. The ligands are expelled from those monolayers and form separate particle-free domains, giving rise to lamellar morphology. For example, the system shown in Figure 2 has n = 8 and N = 15, corresponding to σ = 0.068 chains/nm2, RP/Rg = 1.93, and φP = 0.5. In rendering 3D images (Figure 2a and others below) the box size was doubled in each of the three directions. The lamellar period is approximately 2.4RP. The degree of interpenetration between the two species is substantially higher than it is for the case of flexible block copolymer with f = φP = 0.5 and similar degree of

studies showing formation of anisotropic structures in onecomponent hard core/soft shell sphere dispersions.79,80 The morphologies of the HNP melts obtained in the SCFDFT calculations are reminiscent of the theoretical and computational results of Pryamitsyn et al.16 (even though that work was devoted to solutions of HNPs, not melts, the main factors determining anisotropic assembly are very similar). The critical competition is between the interfacial energy of particle/polymer contacts (favoring the nanoparticle clustering and aggregation) and the conformational entropy of the ligands (favoring a more uniform dispersion). The result of this competition is different at different ligand lengths and grafting densities. We will now describe the various morphologies in more detail. E

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules

Figure 5. (a) Three-dimensional view of the FCC morphology, with red representing particle-rich regions. The surfaces correspond to ϕ = 0.5. The blue arrow goes from the center of a sphere to the point halfway between two spheres in the same layer. (b) Density profiles of particle (red) and ligand (blue) species along the direction of the blue arrow. Here, n = 24 and N = 15.

incompatibility χN, in excellent agreement with earlier study of Lee and co-workers.60 If one increases the grafting density of the ligands, the overall particle volume fraction is reduced, and one would expect a transition to a morphology where the particles are packed less tightly. While in flexible-chain block copolymers, such a morphology usually is a cylindrical one, here we also find a structure that can be labeled as “sheets”: the particle-rich and particle-free layer alternate similar to lamellar structure, but within the particle-rich layer, two-dimensional crystallization of spheres is taking place (Figure 3). The structure of this layer is similar to sheets observed in simulations of Akcora et al.15 (HNPs in polymer melts), Glotzer and co-workers28,50 (tethered nanospheres), and experiments of Jaeger and coworkers81 (nanoparticle membranes assembled at oil−water interfaces). Figure 3 shows the particle arrangement for the system with n = 10 and N = 30, corresponding to σ = 0.086 chains/nm2, RP/Rg = 1.36, and φP = 0.29. The arrangement of the particles within the layer is consistent with two-dimensional (monolayer) hexagonal crystal. The interparticle distance within the layer depends on the particle volume fraction and is significantly smaller than the distance between neighboring particles in adjacent layers. Most of the ligands are expelled from the particle-rich layers and form a particle-free layer between them, similar to the lamellar case. At lower particle volume fractions (higher grafting densities), a cylindrical morphology is observed. The particles now are arranged into strings, and the strings, in turn, are assembled into a 2D hexagonal crystalline structure. The system shown in Figure 4 has n = 11 and N = 15, corresponding to σ = 0.094 chains/nm2, RP/Rg = 1.93, and φP = 0.42. The intercylinder spacing is approximately 2.8RP. As in the lamellar case, there is substantial interpenetration between the particle core and the polymer matrix, once again in agreement with the results of Lee and co-workers.60 Further increase in the grafting density (decrease in the particle volume fraction) leads to the formation of threedimensional crystalline structures of spheres. Before discussing this morphology in more detail, we caution that the SCF-DFT free energy functional is not accurate enough to distinguish

between crystalline structures like FCC, BCC, and HCP. In this study, the FCC structure was used as a representative one. It is possible that the choice of FCC or BCC would give slightly different results, but the main featuretransition from cylinders to dispersed sphereswould remain. The topic of determining the exact crystalline morphologies of sphereforming block copolymers (including the so-called Frank− Kasper phases82) recently emerged as a new frontier in polymer science, and similar challenges would apply to HNP systems; it is, however, out of scope for this study. The morphology of an FCC crystal corresponding to the system with n = 24 and N = 15, corresponding to σ = 0.205 chains/nm2, RP/Rg = 1.93, and φP = 0.25, is shown in Figure 5. The distance between two neighboring spheres is approximately 3.0RP. The ligands are now distributed in a nearly isotropic fashion, so that each HNP can be modeled as a core− shell sphere, with hard core and soft shell. To our knowledge, three-dimensional assemblies of HNPs have not yet been probed with SCF-DFT but have been extensively studied using particle-based28,30−32,50 and integral equation-based41,42,45 approaches. To plot the phase map, we need to compare the free energies of all morphologies for each value of N and σ. The resulting equilibrium (lowest-free energy) morphologies are shown in Figure 6. From Figure 6, one can see that if the ligands are sufficiently long (large N), they cause segregation between particle-rich and particle-free regions. Then, at relatively low grafting density, the HNP system is arranged into lamellae; further increase in the grafting density leads to the formation of sheets with twodimensional crystallization within the particle-rich layers. As the grafting density is increased further, one can expect those layers to melt and change into isotropic, disordered liquid. In the case of shorter ligands, the sheet morphology is pre-empted by the cylindrical and FCC structures. We can gain further insights by replotting Figure 6 using (φP, RP/Rg) coordinates instead of (N, σ) coordinates (Figure 7). Such a representation would facilitate comparison between the behavior of tethered nanoparticles and conventional diblock copolymers and also allow comparison with experimental data. F

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules

system prefers to segregate the particle-rich and particle-free regions as either lamellar or sheet morphology. The transition from sheet to lamellar shifts to higher φP as RP/Rg is increased. At higher RP/Rg, the one-dimensional segregation between particle-rich and particle-free regions becomes less favorable, and the sheet morphology is pre-empted by morphologies with higher average interfacial curvature, hexagonal cylindrical, or FCC. It is important to caution that the current phase map is based on the specific choice of the particle size and the specific sphere free energy functional; thus, one can expect that the phase boundaries could be refined by future studies. We can also compare our phase map to some recent experimental studies (see Table 1 for detailed description of the nanocomposite particles). Koerner et al.5 observed the formation of anisotropic structures for their HNP system with the “intermediate” grafting density (φP = 0.18, RP/Rg = 0.82). According to SCF-DFT predictions, this HNP system is indeed expected to exhibit anisotropic ordering (“sheet” morphology to be more precise), as marked by the gray star labeled “1”. Goel et al.19 prepared a number of HNP systems forming a weak FCC crystalall those systems would be located in the upper left corner of our phase map (large RP/Rg and/or small φP), as schematically depicted by the blue star labeled “2” and the two arrows pointing leftward and upward. As was mentioned above, several earlier theoretical models analyzed in detail the transition from crystal (FCC) to isotropic liquid.23,51,53 Within those theories, it was predicted that the FCC crystal would “melt” upon increase in the degree of polymerization N and grafting density σ. As the “reduced” grafting density, σ* = σNa2, is increased, the polymer layer undergoes transitions from “mushroom” to “semidilute polymer brush” (SDPB) to “concentrated polymer brush” (CPB). Upon further increase in σ*, the polymer shell becomes layered, with the CPB layer surrounding the core and the SDPB layer facing the neighboring particles. As proposed e.g., by Choi et al.,23 the melting of the FCC crystals takes place where the volume fraction of SDPB becomes greater than the “free volume” of the FCC lattice (∼0.26). However, exact calculation of the CPB and SDPB volumes is challenging even in the limit of isolated HNP in a good solvent or polymer melt. Previous theoretical models were based on the scaling theory of Daoud and Cotton,83 originally developed for solutions of star polymers, and appropriately adapted to the HNP case.84 Other authors recently demonstrated that the use of selfconsistent field theory (SCFT) following the approach of Wijmans and Zhulina85 can result in a reasonably good agreement with experiment.86,87 Our prediction for the FCC melting (dashed thick blue line in Figure 7) can be ⎛ R ⎞* approximated by the following expression, ⎜ RP ⎟ ≈ 3.0ϕP1/2 , ⎝ g⎠

Figure 6. Phase map for the HNP melt. Here, N is ligand length and σ is grafting density. Symbols represent simulation results, with equilibrium (minimum free energy) morphologies shown as follows: yellow = lamellar, gray = sheets, orange = cylindrical, and blue = FCC spherical. Dashed lines are boundaries of the region studied in these simulations, while solid lines are approximate locations of transitions between various phases.

Figure 7. Phase map of the HNP one-component systems expressed in terms of φP and RP/Rg. Symbols are labeled as in Figure 6. Solid lines are guide to the eye depicting approximate phase boundaries. Stars with numbers represent the phase map locations of HNP systems from refs 5 and 19 (see text for more details). The dashed blue line is the predicted FCC melting transition.

Within the SCF-DFT theory in the strong segregation regime, HNP systems having the same scaled parameters (φP, RP/Rg) are expected to have the same morphologies; thus, we can compare results for systems with different core radii and superimpose them on the same phase map. Examination of the phase map of Figure 7 shows several interesting features. For small RP/Rg (longer ligands), the

or, alternatively, N * =

2 3

2

( Ra ) ϕP−1 P

for the transition chain

length assuming a given core volume fraction and nanoparticle radius. This scaling is qualitatively similar to the analysis of

Table 1. Experimental HNP Systems Referenced in This Study HNP system, as labeled in the original paper

core material

core radius (nm)

polymer shell

polymer MW (kg/mol)

grafting density (chains/nm2)

estimated φP

estimated RP/Rg

morphology

ref

PBA-4K PBA-25K SiPS0p05

SiO2 SiO2 SiO2

8.25 8.25 8.0

PBA PBA PS

3.9 24.7 120

0.8 0.8 0.05

0.33 0.08 0.18

4.65 1.85 0.82

FCC FCC aniso

19 19 5

G

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules Choi et al.,23 although detailed comparison is beyond the scope of this paper. One crucial question is whether the predicted anisotropic microphase-separated phases could be observed experimentally. Unlike in the case of block copolymers, microphase-separated ordered phases like hexagonal cylinders or lamellae have yet to be seen, even though anisotropic aggregates with short-range order are now known. It is possible that the HNP assemblies are not able to fully equilibrate and the morphologies observed in experiments are metastable or kinetically trapped. Alternatively, it is possible that the inability of the HNPs to organize into long-range ordered anisotropic structures is due to polydispersity (with regard to both particle size and ligand length). A recent study by Bachhar et al.88 shows that in polydisperse HNP systems (at least in solution or in a polymer matrix) there is a tendency to fractionate, rather than organize based on the “average” particle size and ligand length. On the basis of the modeling, one can hypothesize that “single-component nanocomposites” consisting of nanoparticles with one or more grafted polymer ligands offer wide variety of complex and potentially unexplored structures. In the simplest case where the ligands are pure monodisperse oligomers or homopolymers (and, in principle, we can graft linear block copolymers, star-polymers, blends of ligands, etc.), the morphology map is somewhat similar to that of diblock copolymers. Furthermore, the characteristic dimensions of those structures should be smaller than those of conventional block copolymers. For example, a lamellar period of a 104 kg/ mol PS-b-PMMA is approximately 48 nm,89 and even for higher-χ polymers such as PS-b-PDMS, structures with characteristic dimensions below 10 nm are difficult to prepare. For the HNP materials, characteristic dimensions are about 2− 3 times the particle radius, so sub-10 nm lamellar domains could be obtained in a fairly straightforward way. Furthermore, the two domain types (particle-rich and particle-free) will have extremely different properties (mechanical, thermal, physical, electrical, and optical). And finally, one could potentially find ways of making different types of porous materials by removing either the organic portion (perhaps by high-temperature treatment that should also ensure that particles have sufficiently sintered) or the particles themselves. The SCF-DFT approach, while fast and simple, has some limitations. As already mentioned, the free energy functional is based on the ensemble-averaged description of the particle density, rather than instantaneous particle positions. As a result, the detailed investigation of particle arrangements, e.g., within the monolayers, is not straightforward. Moreover, within this approach, one is likely limited to exploration of “intermediate” (φP ∼ 0.2−0.5) core volume fractions. At small φP, the particles become more delocalized, and density fluctuations (not captured in SCF-DFT) would become the dominant contribution to thermodynamics. At large φP, while the overall approach could still work, one would need a more elaborate density functional and substantially higher spatial resolution. Likewise, the model is also somewhat limited in terms of RP/Rg range. For very large or very small RP/Rg, the need to increase the spatial resolution so as to accurately describe either polymers or the particle would make computations expensive and impractical. Even so, the SCF-DFT approach can provide initial screening, to be followed by more extensive particlebased simulations. Furthermore, a full exploration of a block copolymer phase diagram using SCFT (or, in this case, SCFDFT) requires extensive optimization of the unit cells for

various candidate morphologies, followed by very accurate free energy calculation and selection of the lowest-free energy morphology for each point in the phase space.90,91 In particular, it is likely that gyroid, double-diamond, or other bicontinuous structures could be found at least for some parts of the phase map. Such a program is beyond the scope of the current paper and will be explored in the future.



CONCLUSIONS We developed a new modeling framework aimed at predicting the morphology of polymer-grafted or “hairy” nanoparticles. The model is based on the SCF-DFT approach originally developed by Thompson et al.61 and expanded to account explicitly for the ligand configurations. The simulations demonstrated substantial similarity between the HNP onecomponent systems and the block copolymer melts, with spherical, cylindrical, and lamellar morphologies present in both material types. The simulation results suggest that the arrangement of nanoparticles in single-component HNP nanocomposites can be isotropic or nearly isotropic (FCC), quasi-one-dimensional (cylindrical), or quasi-two-dimensional (sheets/lamellar). The anisotropic morphologies are formed as the ligands rearrange along the surfaces of the particles, thus allowing more intimate particle−particle contacts. Note that we assumed a weak bonding between the ligand and the particle surface; it is possible that the phase diagram would change if the ligands are “locked” in their positions (although one would assume that the importance of the ligand binding strength should diminish as the ligand chain length is increased). The method proposed here can be extended to a number of other HNP systems (for example, considering the cases of block copolymer ligands or mixed homopolymer ligands). This will be the subject of future studies.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.macromol.7b01922. Selection of the Flory−Huggins parameters, morphology optimization procedure (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Valeriy V. Ginzburg: 0000-0002-2775-5492 Notes

The author declares no competing financial interest.



ACKNOWLEDGMENTS This work was supported by The Dow Chemical Company. The author thanks Prof. Arthi Jayaraman (University of Delaware) and Drs. Jeffrey Weinhold, YuanQiao Rao, and Peter Trefonas (Dow) for stimulating discussions. I am also thankful to anonymous reviewers for many helpful suggestions.



REFERENCES

(1) Chakkalakal, G. L.; Ramakrishnan, S.; Bockstaller, M. R. PolymerTethered Nanoparticle MaterialsAn Emerging Platform for Multi-

H

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules functional Hybrid Materials. In Hybrid and Hierarchical Composite Materials; Springer: 2015; pp 65−94. (2) Ruan, Y.; Gao, L.; Yao, D.; Zhang, K.; Zhang, B.; Chen, Y.; Liu, C.-Y. Polymer-Grafted Nanoparticles with Precisely Controlled Structures. ACS Macro Lett. 2015, 4 (10), 1067−1071. (3) Li, Y.; Wang, L.; Natarajan, B.; Tao, P.; Benicewicz, B. C.; Ullal, C.; Schadler, L. S. Bimodal “matrix-free” polymer nanocomposites. RSC Adv. 2015, 5 (19), 14788−14795. (4) Kumar, S. K.; Jouault, N.; Benicewicz, B.; Neely, T. Nanocomposites with Polymer Grafted Nanoparticles. Macromolecules 2013, 46 (9), 3199−3214. (5) Koerner, H.; Drummy, L.; Benicewicz, B.; Li, Y.; Vaia, R. Nonisotropic Self-Organization of Single-Component Hairy Nanoparticle Assemblies. ACS Macro Lett. 2013, 2, 670−676. (6) Chevigny, C.; Dalmas, F.; Di Cola, E.; Gigmes, D.; Bertin, D.; Boué, F. o.; Jestin, J. Polymer-Grafted-Nanoparticles Nanocomposites: Dispersion, Grafted Chain Conformation, and Rheological Behavior. Macromolecules 2011, 44 (1), 122−133. (7) Ganesan, V.; Jayaraman, A. Theory and simulation studies of effective interactions, phase behavior and morphology in polymer nanocomposites. Soft Matter 2014, 10 (1), 13−38. (8) Martin, T. B.; Jayaraman, A. Identifying the Ideal Characteristics of the Grafted Polymer Chain Length Distribution for Maximizing Dispersion of Polymer Grafted Nanoparticles in a Polymer Matrix. Macromolecules 2013, 46 (22), 9144−9150. (9) Ginzburg, V. Self-Consistent Field Theory Modeling of Polymer Nanocomposites. In Modeling and Prediction of Polymer Nanocomposite Properties; Mittal, V., Ed.; Wiley-VCH: Weinheim, Germany, 2013; pp 11−38. (10) Jog, P.; Ginzburg, V.; Srivastava, R.; Weinhold, J.; Jain, S.; Chapman, W. G. Application of Mesoscale Field-Based Models to Predict Stability of Particle Dispersions in Polymer Melts. In Advances in Chemical Engineering; West, D. H., Ed.; Elsevier: Amsterdam, 2010; Vol. 39, pp 131−164. (11) Balazs, A. C.; Bicerano, J.; Ginzburg, V. V. Polyolefin/Clay Nanocomposites: Theory and Simulation. In Polyolefin Composites. John Wiley & Sons, Inc. 2007, 415−448. (12) Sunday, D.; Ilavsky, J.; Green, D. L. A Phase Diagram for Polymer-Grafted Nanoparticles in Homopolymer Matrices. Macromolecules 2012, 45 (9), 4007−4011. (13) Jiao, Y.; Akcora, P. Assembly of Polymer-Grafted Magnetic Nanoparticles in Polymer Melts. Macromolecules 2012, 45 (8), 3463− 3470. (14) Wang, X.; Foltz, V. J.; Rackaitis, M.; Böhm, G. G. A. Dispersing hairy nanoparticles in polymer melts. Polymer 2008, 49 (26), 5683− 5691. (15) Akcora, P.; Liu, H.; Kumar, S. K.; Moll, J.; Li, Y.; Benicewicz, B. C.; Schadler, L. S.; Acehan, D.; Panagiotopoulos, A. Z.; Pryamitsyn, V.; Ganesan, V.; Ilavsky, J.; Thiyagarajan, P.; Colby, R. H.; Douglas, J. F. Anisotropic self-assembly of spherical polymer-grafted nanoparticles. Nat. Mater. 2009, 8 (4), 354−U121. (16) Pryamtisyn, V.; Ganesan, V.; Panagiotopoulos, A. Z.; Liu, H. J.; Kumar, S. K., Modeling the anisotropic self-assembly of spherical polymer-grafted nanoparticles. J. Chem. Phys. 2009, 131 (22).22110210.1063/1.3267729. (17) Koerner, H.; Opsitnick, E.; Grabowski, C. A.; Drummy, L. F.; Hsiao, M.-S.; Che, J.; Pike, M.; Person, V.; Bockstaller, M. R.; Meth, J. S. Physical aging and glass transition of hairy nanoparticle assemblies. J. Polym. Sci., Part B: Polym. Phys. 2016, 54 (2), 319−330. (18) Fernandes, N. J.; Wallin, T. J.; Vaia, R. A.; Koerner, H.; Giannelis, E. P. Nanoscale Ionic Materials. Chem. Mater. 2014, 26 (1), 84−96. (19) Goel, V.; Pietrasik, J.; Dong, H. C.; Sharma, J.; Matyjaszewski, K.; Krishnamoorti, R. Structure of Polymer Tethered Highly Grafted Nanoparticles. Macromolecules 2011, 44 (20), 8129−8135. (20) Bourlinos, A. B.; Chowdhury, S. R.; Herrera, R.; Jiang, D. D.; Zhang, Q.; Archer, L. A.; Giannelis, E. P. Functionalized Nanostructures with Liquid-Like Behavior: Expanding the Gallery of Available Nanostructures. Adv. Funct. Mater. 2005, 15 (8), 1285−1290.

(21) Bourlinos, A. B.; Herrera, R.; Chalkias, N.; Jiang, D. D.; Zhang, Q.; Archer, L. A.; Giannelis, E. P. Surface-Functionalized Nanoparticles with Liquid-Like Behavior. Adv. Mater. 2005, 17 (2), 234− 237. (22) Fernandes, N. J.; Koerner, H.; Giannelis, E. P.; Vaia, R. A. Hairy nanoparticle assemblies as one-component functional polymer nanocomposites: opportunities and challenges. MRS Commun. 2013, 3 (1), 13−29. (23) Choi, J.; Hui, C. M.; Schmitt, M.; Pietrasik, J.; Margel, S.; Matyjazsewski, K.; Bockstaller, M. R. Effect of polymer-graft modification on the order formation in particle assembly structures. Langmuir 2013, 29 (21), 6452−6459. (24) Zhang, W.-B.; Yu, X.; Wang, C.-L.; Sun, H.-J.; Hsieh, I.-F.; Li, Y.; Dong, X.-H.; Yue, K.; Van Horn, R.; Cheng, S. Z. Molecular nanoparticles are unique elements for macromolecular science: From “nanoatoms” to giant molecules. Macromolecules 2014, 47 (4), 1221− 1239. (25) Yu, X.; Yue, K.; Hsieh, I.-F.; Li, Y.; Dong, X.-H.; Liu, C.; Xin, Y.; Wang, H.-F.; Shi, A.-C.; Newkome, G. R.; Ho, R.-M.; Chen, E.-Q.; Zhang, W.-B.; Cheng, S. Z. D. Giant surfactants provide a versatile platform for sub-10-nm nanostructure engineering. Proc. Natl. Acad. Sci. U. S. A. 2013, 110 (25), 10078−10083. (26) Su, H.; Zheng, J.; Wang, Z.; Lin, F.; Feng, X.; Dong, X.-H.; Becker, M. L.; Cheng, S. Z.; Zhang, W.-B.; Li, Y. Sequential Triple “Click” Approach toward Polyhedral Oligomeric Silsesquioxane-Based Multiheaded and Multitailed Giant Surfactants. ACS Macro Lett. 2013, 2 (8), 645−650. (27) Yue, K.; Liu, C.; Guo, K.; Yu, X.; Huang, M.; Li, Y.; Wesdemiotis, C.; Cheng, S. Z.; Zhang, W.-B. Sequential “click” approach to polyhedral oligomeric silsesquioxane-based shape amphiphiles. Macromolecules 2012, 45 (20), 8126−8134. (28) Zhang, Z. L.; Horsch, M. A.; Lamm, M. H.; Glotzer, S. C. Tethered nano building blocks: Toward a conceptual framework for nanoparticle self-assembly. Nano Lett. 2003, 3 (10), 1341−1346. (29) Phillips, C. L.; Iacovella, C. R.; Glotzer, S. C. Stability of the double gyroid phase to nanoparticle polydispersity in polymertethered nanosphere systems. Soft Matter 2010, 6 (8), 1693−1703. (30) Glotzer, S. C.; Anderson, J. A. Nanoparticle assembly: made to order. Nat. Mater. 2010, 9 (11), 885−887. (31) Glotzer, S. C.; Solomon, M. J. Anisotropy of building blocks and their assembly into complex structures. Nat. Mater. 2007, 6 (8), 557− 562. (32) Chan, E. R.; Ho, L. C.; Glotzer, S. C. Computer simulations of block copolymer tethered nanoparticle self-assembly. J. Chem. Phys. 2006, 125 (6), 064905. (33) Glotzer, S. C.; Horsch, M. A.; Iacovella, C. R.; Zhang, Z. L.; Chan, E. R.; Zhang, X. Self-assembly of anisotropic tethered nanoparticle shape amphiphiles. Curr. Opin. Colloid Interface Sci. 2005, 10 (5−6), 287−295. (34) Jayaraman, A.; Schweizer, K. S. Liquid state theory of the structure and phase behaviour of polymer-tethered nanoparticles in dense suspensions, melts and nanocomposites. Mol. Simul. 2009, 35 (10−11), 835−848. (35) Hall, L. M.; Jayaraman, A.; Schweizer, K. S. Molecular theories of polymer nanocomposites. Curr. Opin. Solid State Mater. Sci. 2010, 14 (2), 38−48. (36) Ginzburg, V. Nanoparticle/Polymer Blends: Theory and Modeling. In Encyclopedia of Polymer Blends; Isayev, A. I., Ed.; Wiley-VCH: Weinheim, 2011; Vol. 1, pp 233−267. (37) Estridge, C. E.; Jayaraman, A. Effect of homopolymer matrix on diblock copolymer grafted nanoparticle conformation and potential of mean force: a molecular simulation study. J. Polym. Sci., Part B: Polym. Phys. 2015, 53 (1), 76−88. (38) Martin, T. B.; Jayaraman, A. Effect of matrix bidispersity on the morphology of polymer-grafted nanoparticle-filled polymer nanocomposites. J. Polym. Sci., Part B: Polym. Phys. 2014, 52 (24), 1661− 1668. (39) Martin, T. B.; Dodd, P. M.; Jayaraman, A. Polydispersity for tuning the potential of mean force between polymer grafted I

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules nanoparticles in a polymer matrix. Phys. Rev. Lett. 2013, 110 (1), 018301. (40) Jayaraman, A. Polymer grafted nanoparticles: Effect of chemical and physical heterogeneity in polymer grafts on particle assembly and dispersion. J. Polym. Sci., Part B: Polym. Phys. 2013, 51 (7), 524−534. (41) Jayaraman, A.; Nair, N. Integrating PRISM theory and Monte Carlo simulation to study polymer-functionalised particles and polymer nanocomposites. Mol. Simul. 2012, 38 (8−9), 751−761. (42) Jayaraman, A.; Schweizer, K. S. Effective interactions and selfassembly of hybrid polymer grafted nanoparticles in a homopolymer matrix. Macromolecules 2009, 42 (21), 8423−8434. (43) Jayaraman, A.; Schweizer, K. S. Effective Interactions, Structure, and Phase Behavior of Lightly Tethered Nanoparticles in Polymer Melts. Macromolecules 2008, 41 (23), 9430−9438. (44) Jayaraman, A.; Schweizer, K. S. Structure and assembly of dense solutions and melts of single tethered nanoparticles. J. Chem. Phys. 2008, 128 (16), 164904. (45) Jayaraman, A.; Schweizer, K. S. Effect of the number and placement of polymer tethers on the structure of concentrated solutions and melts of hybrid nanoparticles. Langmuir 2008, 24 (19), 11119−11130. (46) Ginzburg, V. V. Polymer-Grafted Nanoparticles in Polymer Melts: Modeling Using the Combined SCFT−DFT Approach. Macromolecules 2013, 46 (24), 9798−9805. (47) Chao, H.; Hagberg, B.; Riggleman, R. The distribution of homogeneously grafted nanoparticles in polymer thin films and blends. Soft Matter 2014, 10 (40), 8083−8094. (48) Koski, J.; Chao, H.; Riggleman, R. A. Predicting the structure and interfacial activity of diblock brush, mixed brush, and Janus-grafted nanoparticles. Chem. Commun. 2015, 51 (25), 5440−5443. (49) Koski, J.; Chao, H.; Riggleman, R. A. Field theoretic simulations of polymer nanocomposites. J. Chem. Phys. 2013, 139 (24), 244911. (50) Marson, R. L.; Nguyen, T. D.; Glotzer, S. C. Rational design of nanomaterials from assembly and reconfigurability of polymertethered nanoparticles. MRS Commun. 2015, 5 (03), 397−406. (51) Yu, H.-Y.; Koch, D. L. Structure of Solvent-Free Nanoparticle− Organic Hybrid Materials. Langmuir 2010, 26 (22), 16801−16811. (52) Yu, H.-Y.; Srivastava, S.; Archer, L. A.; Koch, D. L. Structure factor of blends of solvent-free nanoparticle−organic hybrid materials: density-functional theory and small angle X-ray scattering. Soft Matter 2014, 10 (45), 9120−9135. (53) Yu, H.-Y.; Koch, D. L. Predicting the Disorder−Order Transition of Solvent-Free Nanoparticle−Organic Hybrid Materials. Langmuir 2013, 29 (26), 8197−8202. (54) Chremos, A.; Panagiotopoulos, A. Z.; Yu, H.-Y.; Koch, D. L. Structure of solvent-free grafted nanoparticles: Molecular dynamics and density-functional theory. J. Chem. Phys. 2011, 135 (11), 114901. (55) Yu, H.-Y.; Koch, D. L. Structure of solvent-free nanoparticleorganic hybrid materials, 2010. (56) Chremos, A.; Douglas, J. F. Self-assembly of polymer-grafted nanoparticles in solvent-free conditions. Soft Matter 2016, 12 (47), 9527−9537. (57) Hong, B. B.; Chremos, A.; Panagiotopoulos, A. Z. Dynamics in coarse-grained models for oligomer-grafted silica nanoparticles. J. Chem. Phys. 2012, 136, 204904. (58) Chremos, A.; Panagiotopoulos, A. Z.; Yu, H. Y.; Koch, D. L. Structure of solvent-free grafted nanoparticles: Molecular dynamics and density-functional theory. J. Chem. Phys. 2011, 135, 114901. (59) Chremos, A.; Panagiotopoulos, A. Z. Structural Transitions of Solvent-Free Oligomer-Grafted Nanoparticles. Phys. Rev. Lett. 2011, 107, 105503. (60) Lee, J. Y.; Balazs, A. C.; Thompson, R. B.; Hill, R. M. Selfassembly of amphiphilic nanoparticle-coil “tadpole” macromolecules. Macromolecules 2004, 37 (10), 3536−3539. (61) Thompson, R. B.; Ginzburg, V. V.; Matsen, M. W.; Balazs, A. C. Predicting the mesophases of copolymer-nanoparticle composites. Science 2001, 292 (5526), 2469−2472. (62) Thompson, R. B.; Ginzburg, V. V.; Matsen, M. W.; Balazs, A. C. Block copolymer-directed assembly of nanoparticles: Forming

mesoscopically ordered hybrid materials. Macromolecules 2002, 35 (3), 1060−1071. (63) Hur, K.; Hennig, R. G.; Escobedo, F. A.; Wiesner, U. Mesoscopic structure prediction of nanoparticle assembly and coassembly: theoretical foundation. J. Chem. Phys. 2010, 133 (19), 194108. (64) Matsen, M. W.; Schick, M. Microphase Separation in Starblock Copolymer Melts. Macromolecules 1994, 27 (23), 6761−6767. (65) Rasmussen, K. O.; Kalosakas, G. Improved numerical algorithm for exploring block copolymer mesophases. J. Polym. Sci., Part B: Polym. Phys. 2002, 40 (16), 1777−1783. (66) Reister, E.; Fredrickson, G. H. Nanoparticles in a diblock copolymer background: The potential of mean force. Macromolecules 2004, 37 (12), 4718−4730. (67) Reister, E.; Fredrickson, G. H. Phase behavior of a blend of polymer-tethered nanoparticles with diblock copolymers. J. Chem. Phys. 2005, 123 (21), 214903. (68) Sides, S. W.; Kim, B. J.; Kramer, E. J.; Fredrickson, G. H. Hybrid particle-field simulations of polymer nanocomposites. Phys. Rev. Lett. 2006, 96, 250601. (69) Balazs, A. C.; Ginzburg, V. V.; Qiu, F.; Peng, G. W.; Jasnow, D. Multi-scale model for binary mixtures containing nanoscopic particles. J. Phys. Chem. B 2000, 104 (15), 3411−3422. (70) Ginzburg, V. V.; Gibbons, C.; Qiu, F.; Peng, G. W.; Balazs, A. C. Modeling the dynamic behavior of diblock copolymer/particle composites. Macromolecules 2000, 33 (16), 6140−6147. (71) Chen, K.; Li, H.-s.; Zhang, B.-k.; Li, J.; Tian, W.-d. A New SelfConsistent Field Model of Polymer/Nanoparticle Mixture. Sci. Rep. 2016, 6, 20355. (72) Bicerano, J. Prediction of Polymer Properties, 3rd ed.; Marcel Dekker: New York, 2002. (73) Drolet, F.; Fredrickson, G. H. Combinatorial screening of complex block copolymer assembly with self-consistent field theory. Phys. Rev. Lett. 1999, 83 (21), 4317−4320. (74) Bates, F. S.; Fredrickson, G. H. Block copolymersdesigner soft materials. Phys. Today 1999, 52 (2), 32−38. (75) Zaccarelli, E. Colloidal gels: equilibrium and non-equilibrium routes. J. Phys.: Condens. Matter 2007, 19 (32), 323101. (76) Poon, W. The physics of a model colloid−polymer mixture. J. Phys.: Condens. Matter 2002, 14 (33), R859. (77) Fuchs, M.; Schweizer, K. S. Structure of colloid-polymer suspensions. J. Phys.: Condens. Matter 2002, 14 (12), R239. (78) Löwen, H. Melting, freezing and colloidal suspensions. Phys. Rep. 1994, 237 (5), 249−324. (79) Pauschenwein, G. J.; Kahl, G. Zero temperature phase diagram of the square-shoulder system. J. Chem. Phys. 2008, 129 (17), 174107. (80) Glaser, M. A.; Grason, G. M.; Kamien, R. D.; Košmrlj, A.; Santangelo, C. D.; Ziherl, P. Soft spheres make more mesophases. EPL (Europhysics Letters) 2007, 78 (4), 46004. (81) Jiang, Z.; He, J.; Deshmukh, S. A.; Kanjanaboos, P.; Kamath, G.; Wang, Y.; Sankaranarayanan, S. K.; Wang, J.; Jaeger, H. M.; Lin, X.-M. Subnanometre ligand-shell asymmetry leads to Janus-like nanoparticle membranes. Nat. Mater. 2015, 14 (9), 912. (82) Lee, S.; Bluemle, M. J.; Bates, F. S. Discovery of a Frank-Kasper σ phase in sphere-forming block copolymer melts. Science 2010, 330 (6002), 349−353. (83) Daoud, M.; Cotton, J. Star shaped polymers: a model for the conformation and its concentration dependence. J. Phys. 1982, 43 (3), 531−538. (84) Ohno, K.; Morinaga, T.; Takeno, S.; Tsujii, Y.; Fukuda, T. Suspensions of silica particles grafted with concentrated polymer brush: Effects of graft chain length on brush layer thickness and colloidal crystallization. Macromolecules 2007, 40 (25), 9143−9150. (85) Wijmans, C. M.; Zhulina, E. B. Polymer Brushes at Curved Surfaces. Macromolecules 1993, 26 (26), 7214−7224. (86) Dukes, D.; Li, Y.; Lewis, S.; Benicewicz, B.; Schadler, L.; Kumar, S. K. Conformational Transitions of Spherical Polymer Brushes: Synthesis, Characterization, and Theory. Macromolecules 2010, 43 (3), 1564−1570. J

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX

Article

Macromolecules (87) Hore, M. J. A.; Ford, J.; Ohno, K.; Composto, R. J.; Hammouda, B. Direct Measurements of Polymer Brush Conformation Using SmallAngle Neutron Scattering (SANS) from Highly Grafted Iron Oxide Nanoparticles in Homopolymer Melts. Macromolecules 2013, 46 (23), 9341−9348. (88) Bachhar, N.; Jiao, Y.; Asai, M.; Akcora, P.; Bandyopadhyaya, R.; Kumar, S. K. Impact of the Distributions of Core Size and Grafting Density on the Self-Assembly of Polymer Grafted Nanoparticles. Macromolecules 2017, 50 (19), 7730−7738. (89) Han, E.; Kang, H. M.; Liu, C. C.; Nealey, P. F.; Gopalan, P. Graphoepitaxial Assembly of Symmetric Block Copolymers on Weakly Preferential Substrates. Adv. Mater. 2010, 22 (38), 4325−4329. (90) Arora, A.; Morse, D. C.; Bates, F. S.; Dorfman, K. D. Accelerating self-consistent field theory of block polymers in a variable unit cell. J. Chem. Phys. 2017, 146 (24), 244902. (91) Arora, A.; Qin, J.; Morse, D. C.; Delaney, K. T.; Fredrickson, G. H.; Bates, F. S.; Dorfman, K. D. Broadly accessible self-consistent field theory for block polymer materials discovery. Macromolecules 2016, 49 (13), 4675−4690.

K

DOI: 10.1021/acs.macromol.7b01922 Macromolecules XXXX, XXX, XXX−XXX