Theoretical Perspectives on the Influence of Solution-Phase Additives

Feature Article pubs.acs.org/JPCC

Cite This: J. Phys. Chem. C XXXX, XXX, XXX−XXX

Theoretical Perspectives on the Influence of Solution-Phase Additives in Shape-Controlled Nanocrystal Synthesis Xin Qi,† Tonnam Balankura,† and Kristen A. Fichthorn*,†,‡ †

Department of Chemical Engineering and ‡Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802, United States ABSTRACT: Shape-selective solution-phase nanocrystal growth is facilitated by capping agents, or structure-directing agents (SDAs), which guide shape evolution. It is often stated that these chemical additives impart shape selectivity by promoting nanocrystals with a majority of facets to which they bind most strongly. However, little is known on the mechanisms through which they impart shape selectivity. In this Feature Article, we highlight our recent studies aimed at understanding the thermodynamic and kinetic influences of SDAs using theory and computational tools. We review our studies of the poly(vinylpyrrolidone) (PVP)-mediated synthesis of {100}faceted Ag nanocubes in ethylene glycol solution, which has been studied experimentally. Our studies of the interfacial free energy of Ag−PVP solution interfaces show that while solution-phase PVP does bind more strongly to Ag(100) than to Ag(111), this selectivity is not sufficient to thermodynamically change the Wulff shape of a PVP-covered Ag nanocrystal in solution from that of the bare metal in vacuum. These studies indicate that a strong facet binding selectivity is needed for a SDA to thermodynamically alter the solution-phase crystal shape from that of the bare metal. Interestingly, the binding selectivity of PVP for Ag(100) is sufficient to regulate the atom deposition fluxes to Ag(100) and Ag(111), so that cubic Ag(100) nanocrystals form kinetically. Altogether, our studies indicate that kinetic control of metal nanocrystal shapes is likely more prevalent than thermodynamic control. We outline some current challenges in understanding shape-selective solution-phase nanocrystal syntheses.

■

INTRODUCTION Nanocrystals with well-defined shapes have the capability to advance many modern technologies, including selective catalysis,1−6 electronic devices,7 energy storage,8 and medicine.9,10 Noble metals have attracted attention because they uniquely enhance light scattering and absorption within the visible light range via localized surface-plasmon resonance (LSPR).11−13 This fascinating property allows their application in areas such as effective photovoltaics14−18 and high-precision imaging and sensing.11,19,20 Because these applications are highly sensitive to the sizes, shapes, and symmetry of the nanocrystals, numerous shape-controlled synthesis techniques have been developed in the last decade, among which solutionphase syntheses appear to be especially successful for a variety of metals, such as Ag,21−23 Au,24,25 Cu,26−28 Pt,29−31 and Pd.32,33 Despite many reports of successful syntheses and some quantitative efforts to understand them, it is still challenging to predict the outcome of a chemical synthesis of metal nanocrystals. Major challenges include understanding the complex chemistry involved when a metal salt is reduced, how metal atoms and ions nucleate seeds, and how seeds grow into well-defined shapes, possibly with the help of solutionphase additives called capping agents, or structure-directing agents (SDAs). © XXXX American Chemical Society

Theory and simulation can be beneficial in efforts to understand these syntheses. However, a complete theoretical description of an entire synthesis is hampered by numerous complexities. First, the chemical complexities in these systems practically mandate that any theoretical description should be based on quantum mechanics. Though quantum mechanical studies, particularly those based on density-functional theory (DFT), have been useful in understanding how solution-phase species interact with metal crystals as they grow,34−42 such calculations are typically done in a zero-temperature and vacuum environment, which is far removed from the solutionphase environment of experimental syntheses. DFT-based ab initio molecular dynamics (AIMD) studies can and have been used to study small, capped nanoparticles at finite temperatures.43 However, AIMD simulations are still limited to ∼1000 atoms and picosecond time scales. These difficulties make it unlikely that quantum calculations alone will lead to significant advances in understanding syntheses. Classical molecular-dynamics (MD) simulations can extend the range of first-principles DFT calculations, to further advance our understanding of solution-phase syntheses. The accuracy of these simulations is governed by the underlying Received: January 17, 2018 Revised: April 5, 2018 Published: April 17, 2018 A

DOI: 10.1021/acs.jpcc.8b00562 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Feature Article

The Journal of Physical Chemistry C force field, which is typically fit to results from quantum calculations or experiment. Typical length and time scales that can be reached in MD simulations are currently in the nanoscale, and it is not uncommon to use such simulations to probe microscale lengths and times. With MD simulations, we can rigorously compute free energies, which is difficult in DFT without resorting to approximations for the liquid phase. Though MD simulations may still not reach experimental scales, many different theoretical techniques have been developed to allow for enhanced MD sampling of thermodynamic and kinetic quantities, such as interfacial free energies and reaction rates. These quantities can be incorporated into continuum theories and coarse-grained simulations to extend the range of lengths and times that can be probed. Thus, MD simulations can be useful in evaluating various aspects of solution-phase metal crystal growth. In this Feature Article, we review our MD studies44−48 aimed at understanding the origins of {100}-faceted Ag nanostructures that are observed experimentally21,22,49−51 in the solutionphase polyol process. In particular, we investigated the role of poly(vinylpyrrolidone) (PVP), a commonly used SDA in solution-phase syntheses, in directing the growth of monocrystalline Ag seeds into {100}-faceted Ag cubes. Our studies are based on the metal−organic many-body (MOMB) force field that we developed to bridge DFT and MD.52−54 This force field encompasses both metal−organic direct-bonding interactions through a Morse potential and a many-body term based on the embedded-atom method potential, as well as pairwise van der Waals interactions. By including a many-body term for direct chemical bonding of certain organic atoms to metal surfaces, it captures the facet-selective binding of PVP on Ag(100), as well as facet-selective binding of alkylamines on Cu(100)55 that are seen in DFT calculations. We examine the commonly held idea that SDAs promote shapes with a predominance of a particular crystal facet by binding most strongly to that facet. We note that our firstprinciples studies with DFT indeed confirm this idea, as these studies indicate that PVP does bind more strongly to Ag(100) than to Ag(111).36,37 However, DFT studies do not indicate the thermodynamic and kinetic ramifications that this selective binding may have for crystal shapes. We note that other MD studies have proven useful in resolving the anisotropic growth of Au nanorods,56−59 peptide-mediated growth of Pt nanocubes,60 etc.61 A hallmark of our work is that we rigorously evaluate free energies associated with the kinetics and thermodynamics of nanocrystal growth at liquid−solid interfaces. Below, we review our studies.

In the Wulff construction, facets in the equilibrium crystal shape follow the relation hi =λ γsl, i

(1)

where hi is the orthogonal distance from the crystal center to facet i, γsl,i is the solid−liquid interfacial free energy of facet i, and λ is a constant for the structure. Thus, given γsl,i for each relevant facet i, we can predict the equilibrium shape of a crystal in liquid solution. Experimentally, γsl can be estimated from the contact angle of a liquid drop on a solid surface65,66 or backcalculated from nucleation or crystal growth rates.67−70 Interfacial energies obtained via such means tend not to be facet-specific,67−70 so they would be of little use in the Wulff construction. Theoretical and computational tools can be an excellent way to evaluate γ and its facet dependence.71−75 First, we should clarify that although the interfacial tension σ is occasionally used to represent the interfacial free energy γ, this mathematical equivalence only holds true for fluid systems and not for fluid− solid interfaces. The interfacial free energy is related to the interfacial tension via σ = γ + A(dγ /dA)

(2) 76,77

where A is the interfacial area. This equation indicates that for stress-sensitive surfaces, such as the surfaces of solids, γ varies with A such that dγ/dA is nonzero. It is evident that γsl cannot be approximated using σsl for such surfaces and that it is the interfacial free energy that is relevant in determining crystal shape. We recently introduced an MD-based multi-scheme thermodynamic integration (TI) method44 to estimate γsl. Our method is inspired by prior “cleaving-wall” methods,72−74 and it is targeted to resolve previous difficulties in applying cleaving-wall methods to all-atom simulations in multicomponent systems.44 This is a six-step method in which we reversibly transform a solid−liquid interfacial system to individual solid and liquid bulk phases, as is outlined in Figure 1. In each step i, the free-energy change ΔFi is calculated from TI, in which ΔFi is given by ΔFi =

∫

∂U(λ) dλ ∂λ

(3)

Here λ is a scaling parameter that varies between 0 and 1. λ is used to manipulate the interaction(s) that accomplish the needed transition. U(λ) is the potential energy of the system at an intermediate state, whose energy is influenced by λ, and the angled brackets denote an ensemble average. We briefly summarize our method below and note that a detailed description can be found in ref 44. The multi-scheme TI method begins with a solid−liquid interface, in which the solid is represented using a slab and the rest of the simulation cell is filled with liquid/solution molecules at the bulk liquid density. In step 1, we gradually switch on two superimposed, short-range repulsive walls at each of the two solid−liquid interfaces that occur in a simulation cell with periodic boundary conditions (PBCs). One of the superimposed walls interacts only with the solid, and the other interacts only with the liquid. The positions of these walls are fixed throughout all steps of the simulation. The solid− liquid interaction is gradually turned off in step 2. Since the solid and liquid are confined by repulsive walls, they behave as

■

THERMODYNAMIC INFLUENCE OF CAPPING MOLECULES Under thermodynamic control, the solid−liquid interfacial free energy, γsl, is a key parameter in determining the equilibrium nanocrystal shape. When crystal growth reaches equilibrium, the total interfacial free energy is minimized, and the Wulff construction yields the equilibrium shape.62,63 For bare, facecentered cubic (fcc) metals, we have γ{111} < γ{100} < γ{110}, and the Wulff shapes primarily express {111} facets.64 Various adsorbed species, including solvent and capping molecules, can alter γsl of crystal facets in the solution phase and potentially influence the equilibrium crystal shape. Consequently, a prevailing hypothesis is that preferential binding of a capping molecule to a particular crystal facet can lead to a nanocrystal shape that predominantly expresses that facet. B


Feature Article

The Journal of Physical Chemistry C

where Lx and Ly are the x- and y-dimensions of the simulation box, the product of which gives the interfacial area, and the factor of 2 in the denominator indicates there are two interfaces in a simulation cell with PBCs. We applied the multi-scheme TI method to understand the role of PVP in the solution-phase growth of {100}-faceted Ag nanocrystals45 by calculating γsl for Ag(100) and Ag(111), the most relevant facets in this synthesis. We calculated the surface free energy γ of the bare metal facets (i.e., Ag in contact with vacuum), as well as γsl of the surfaces in pure solvent (ethylene glycol, EG) and in PVP−EG solution with various PVP concentrations, where we use PVP icosamer (PVP20) to represent long-chain PVP in experiment. The results for vacuum−Ag, EG−Ag, and PVP−EG−Ag interfaces are summarized in Figure 2. Our results for Ag surfaces in contact

Figure 2. Interfacial/surface free energy of the studied interfaces. The predicted thermodynamic equilibrium shape for each system is nearly identical and is shown as the structure in the figure. Reproduced from ref 45 with permission from the Royal Society of Chemistry. Copyright 2017 Royal Society of Chemistry.

with vacuum confirm the expectation that fcc metals have a lower γ for {111} facets than {100} facets. Ag facets in contact with EG solvent have significantly lower values of γsl. Although the adsorption of PVP can further lower γsl from the puresolvent values, the ratio of interfacial free energies between Ag(100) and Ag(111) remains relatively unaffected. Using the Wulff construction (eq 1), we predict that nanocrystals in EG solvent or PVP−EG solution will have essentially the same shape as crystals in vacuum−a truncated octahedron or cuboctahedron. We find that the stronger binding of PVP to Ag(100) facets only lowers γsl,Ag(100) by 5.3 mJ/m2 more than γsl,Ag(111), whereas γAg(100) is 72.4 ± 1.1 mJ/m2 higher than γAg(111). Therefore, it is unlikely that PVP induces thermodynamic Ag nanocubes. Although it is evident that PVP is incapable of altering the equilibrium shape of an Ag crystal in solution, this result raises the question of how facet-selective the binding of a capping molecule needs to be to effectively induce an equilibrium crystal-shape change. To probe this question, we performed an analysis using generic capping molecules, represented as Lennard-Jones (LJ) particles with tunable affinities for {100} and {111} facets.45 By varying the LJ interaction-strength parameters for these model capping agents and calculating γsl for Ag(100) and Ag(111) facets, we can predict equilibrium crystal shapes for a range of binding-affinity combinations. From the results of these calculations, we constructed a crystalshape map (Figure 3) to link equilibrium, solution-phase crystal shapes to zero-temperature SDA binding energies. We see in Figure 3 that a relatively large difference in the binding energies

Figure 1. Schematic illustration of the multi-scheme TI method. Details on steps 1 through 6 are given in the text and in ref 44. Reprinted from ref 44, with the permission of AIP Publishing. Copyright 2016 American Institute of Physics.

non-interacting solid and liquid slabs shown in Figure 1 (2b) at the end of step 2. In steps 3 and 4, we recover a bulk liquid by changing the PBC of the liquid cell in the direction normal to the confining walls, then we remove the walls. In steps 5 and 6 we recover a bulk solid. We note that we employ different schemes to recover the bulk liquid and solid and that these differences are discussed in ref 44. The interfacial free energy per unit area is then given by ⎛ ΔF + ΔF + ΔF + ΔF + ΔF + ΔF ⎞ 2 3 4 5 6⎟ γsl = −⎜⎜ 1 ⎟ 2LxLy ⎠ ⎝

(4) C


Feature Article


{100}-faceted cubes (lower right corner of Figure 3). Moderate or weak differences in facet binding energies lead to truncated octahedra, cuboctahedra, or truncated cubes, which fall on the diagonal that runs between the cubic and octahedral regions of Figure 3. We note that Figure 3 was constructed for SDA binding to Ag, so the surface energies of facets in a solution environment reflect the bare-surface energies of Ag (cf., Figure 1). Since Ag(111) has a lower energy than Ag(100), the map in Figure 3 is somewhat skewed toward {111}-faceted crystals, as discussed in ref 45. The crystal-shape map in Figure 3 can be helpful in predicting equilibrium shapes of solution-phase crystals when the solvent−surface and solvent−SDA interactions are similar to those in our studies and when SDA binding energies are known from DFT calculations. To test the veracity of our predictions, we examine two chemical species whose DFT binding energies and the resulting crystal morphologies in solution-phase (EG) syntheses are known. One is 1-methyl-2pyrrolidone (1M2P),52 an analog of PVP monomer. The DFT binding energies for this monomer, along with its predicted surface energies and morphology, are indicated in Figure 3. These predictions fall in line with experimental observations that 1-methyl-2-pyrrolidone is not an effective structuredirecting agent in the shape-selective synthesis of Ag

Figure 3. A crystal-shape map that relates the thermodynamic equilibrium shape of a nanocrystal to the SDA binding energies on Ag(100) and Ag(111). Shapes on the transition curve indicate the predicted shape in each color region. Reproduced from ref 45 with permission from the Royal Society of Chemistry. Copyright 2017 Royal Society of Chemistry.

of SDAs to {100} and {111} surfaces is needed to promote {111}-faceted octahedra (upper left corner of Figure 3) or

Figure 4. Snapshots from MD simulations of Ag deposition on Ag(100) depicting the key steps in the deposition mechanism. For clarity, EG solvent molecules (pink) are downsized in panel a and removed in panel b. (a) An example initial configuration. (b) Key steps in deposition. A solutionphase Ag atom (yellow) arrives at the PVP−EG interface (A). The Ag atom further reaches the PVP layer through interaction with extended segments, especially through Ag−O (red) interaction (B). The immediate deposition is blocked by PVP and the Ag atom can diffuse inside the PVP layer (C). The deposition is accomplished when Ag can channel through a “hole” with suitable size (D). Reprinted with permission from ref 46. Copyright 2015 American Chemical Society. D


Feature Article


Figure 5. (a) PMF for a Ag atom to approach PVP10-covered Ag(100) and Ag(111) surfaces and (b) the density profile of O atoms in the corresponding PVP10 film. The inset shows a zoom-in feature at the “tail” part of the PVP layer. Reprinted with permission from ref 46. Copyright 2015 American Chemical Society.

nanostructures78 and are also in agreement with our full calculation with PVP20, indicated in Figure 2. The second species tested is citric acid (CA),34 the acid form of sodium citrate, which helps to make Ag octahedra in experiments.79 As we can see in Figure 3, the strong binding of CA to Ag(111) relative to Ag(100) places it in the range where an octahedron with slight corner truncation is expected. This consistency with experiment suggests that CA induces the formation of octahedra via thermodynamic control.

It is helpful to qualitatively characterize the deposition of a Ag atom onto PVP-covered Ag(100) and Ag(111) facets in EG solution using direct MD simulations.46 We initiated deposition in these simulations by inserting a free Ag atom with a randomly chosen velocity from the Maxwell−Boltzmann distribution near the PVP−EG interface at z [cf., Figure 4a] and following its trajectory. We find that deposition occurs in four steps. First the Ag atom diffuses from the EG solution to the PVP film. When the atom reaches the film, it is attracted to the PVP segments, and it diffuses into the PVP layer. The atom is held in the PVP layer for some time, during which it can diffuse parallel to the surface, until a hole opens in the PVP layer that allows the atom to access the Ag surface below. Details of each of these stages are illustrated in Figure 4b. The potential of mean force (PMF), or the free-energy profile that characterizes the process by which a solution-phase atom accesses the metal substrate, provides quantitative insight into the observed mechanisms. We obtained the PMF using umbrella sampling81 and umbrella integration82 by sampling the deposition path normal to the surface. The PMFs are shown in Figure 5a, where we notice three main features. First, as the Ag atom approaches the PVP layer, the free energy begins to decrease at a certain height above the surface, corresponding to Ag-atom trapping in the PVP layer. We note that the decrease in the free energy commences further from the surface for Ag(111) than for Ag(100). The free energy continues to decrease until it reaches a local minimum, which corresponds to the Ag atom being trapped in the PVP layer. The atom must surmount a free-energy barrier to access the Ag surface. We note that the barrier to access the surface is significantly lower than the free-energy difference for moving back to the solution phase. Thus, we expect Ag atom adsorption to be virtually irreversible. We also note that the energy barrier to reach the surface is lower for Ag(111) than Ag(100). Differences between the free-energy profiles for Ag(100) and Ag(111) can be attributed to the stronger binding of PVP to Ag(100) and than to Ag(111). These differences are also reflected in the PVP segment density profiles near the two surfaces. We show density profiles of O atoms in PVP (which represents the segment density) on the two surfaces in Figure 5b. At the interface between the PVP layer and the bulk solution, PVP has a greater segment density far from the surface on Ag(111) than on Ag(100) because PVP is less strongly bound on Ag(111) and it forms a looser layer there. These extended segments can more easily attract solution-phase

■

KINETIC SHAPE CONTROL VIA SDA REGULATION OF DEPOSITION FLUX The thermodynamic equilibrium of a nanocrystal is achieved via intra- and interfacet atom migration. If atom diffusion is slow on the time scale of deposition, we expect to observe steadystate kinetic nanocrystal shapes. These kinetic shapes can be predicted by the kinetic Wulff construction, in which the relative linear growth rates, Gi, of various facets i determine the shape.80 The linear facet growth rate of facet i can occur via several different processes: the nucleation of facet i on a facet of a different type, diffusion onto facet i from another facet, or deposition onto facet i from solution. In our studies to date,46−48 we considered facet growth due to deposition. Such a scenario would occur if interfacet diffusion is slow, as would be the case for the relatively large nanocrystals observed in relevant growth experiments,22 and if the nucleation and growth of new facets on top of existing ones is slow or does not occur. If Gi is determined by the solution-phase atom deposition flux to crystal facet i, the relative G̅ i normalized to a reference facet is equivalent to the normalized deposition flux F̅i. Then, the normalized distance xi̅ along a vector, which originates at the center of the crystal and is normal to facet i, is given by xi̅ = Gi = Fi ,

i = 1, ..., N − 1

(5)

where N is the number of facets in the crystal. For an fcc crystal with two possible facets, {100} and {111}, an octahedron is G predicted when G111 ≤ 1/ 3 , and a cube is predicted when 100

G111 G100

≥

3 .46,47 Below, we discuss our investigations with MD

simulations on how the adsorption of capping molecules can influence kinetic nanocrystal shapes by regulating relative deposition fluxes. E


Feature Article


contribute significantly to the MFPT; therefore, Ag atom deposition is neither diffusion-limited nor reaction-limited. In Figure 6, we see that the MFPTs for the systems with solvent only (i.e., no PVP) are significantly larger than those for systems with PVP. This is because PVP attracts incoming Ag atoms and, in the absence of PVP, EG solvent molecules form a tight network on the surface that results in a free-energy barrier for solution-phase Ag atoms to reach the surface.47 The MFPTs are longer for Ag(100) than for Ag(111) for the PVP10 and PVP20 systems, which mostly results from the higher energy barrier in the PMF for Ag(100) surfaces [cf., Figure 5a]. We note that the results in Figure 6 are within a factor of 2 of values obtained from direct MD simulations,46 which can be regarded as being exact for the system simulated. Thus, eq 6 provides a satisfactory description of the MFPT for these systems. The kinetic Wulff shapes that correspond to computed growth-rate ratios for each PVP system are shown in Figure 7.

atoms, so that the free energy begins its decrease further from the surface on Ag(111) than on Ag(100). Similarly, due to the higher binding affinity of PVP for Ag(100), PVP segments are denser and closer to the surface on Ag(100), which leads to a higher free-energy barrier for Ag to access Ag(100). These differences in the PMFs and PVP segment (O atom) density profiles imply that we will have different deposition rates on the two facets, which could lead to steady-state kinetic shapes that are different from thermodynamic Wulff shapes. We obtain the ratio of deposition rates, FAg(111)/FAg(100), that will enable us to calculate kinetic Wulff shapes via eq 5. Here, we utilize the reciprocal relation between the deposition rate and the mean first-passage time (MFPT), tm, for an atom to reach the surface from the bulk solution phase, such that FAg(111)/FAg(100) = tm,Ag(100)/tm,Ag(111). For a Ag atom starting at a height of z from the surface [Figure 4a] and ending at the absorbing boundary at zf at the Ag surface, tm is given by solution of the one-dimensional Smoluchowski equation83 tm(z , z f ) =

zf

∫z ∫z

z

′

exp[−W (z″)/(kT )] dz″

0

exp[W (z′)/(kT )] d z′ D(z′)

(6)

where W(z) is the PMF, D(z) is the z-dependent diffusion coefficient of the Ag atom, and z0 represents the reflecting boundary that accounts for the bulk-solution region. We obtained z-dependent diffusion coefficients from the umbrellasampling trajectories using the simplified Woolf−Roux equation.84,85 The MFPTs obtained from eq 6 for deposition on Ag(100) and Ag(111) are shown in Figure 6. To provide insight into the

Figure 7. Kinetic Wulff shapes of Ag nanocrystals in different solution environments predicted from MFPT calculations. The {100} facets are green, and the {111} facets are orange. Reprinted from ref 47 with the permission of AIP Publishing. Copyright 2016 American Institute of Physics.

For Ag surfaces with a full monolayer coverage of PVP10 and PVP20, {100}-faceted cubes with no {111} corner facets are predicted, while we predict truncated octahedra for PVP5 films. Our results suggest the formation of truncated cubes in EG solvent only; however, the kinetic Wulff construction does not take into account nanocrystal aggregation. We would expect irregular nanostructures to be formed experimentally in the absence of PVP because EG solvent molecules do not adsorb strongly enough to Ag surfaces to protect the nanocrystals from aggregating. These trends are also observed from the direct MD approach.46 Further, in our direct MD study, we observed that cubes were not formed with less-than-monolayer PVP coverage.46 Experimentally, Ag nanocubes are obtained with sufficient PVP concentration,49 so our predictions are in agreement with experiment. We note that Meena and Sulpizi have also proposed for the growth of Au nanorods that the more dense packing of CTAB molecules on the (110) and (100) surfaces with respect to the (111) surface could result in greater growth of the nanorod in the (111) direction.56

Figure 6. Mean first passage times obtained from eq 6 for various adsorbing films on Ag(100) and Ag(111) surfaces.47 For the PVP systems, the MFPTs are partitioned into diffusion and reaction times. Reprinted from ref 47 with the permission of AIP Publishing. Copyright 2016 American Institute of Physics.

deposition mechanism, we divided the MFPTs for the PVPcovered surfaces into diffusion and reaction times, as seen in Figure 6. The diffusion time corresponds to the average time required for the atom to diffuse from bulk solution into the PVP film, and the reaction time characterizes the average time required for the atom to reach the Ag surface once it has entered the PVP film. We note that both diffusion and reaction F


Feature Article


■

CONCLUSIONS AND OUTLOOK In this Feature Article, we summarized our recent theoretical studies aimed at understanding the role of solution-phase additives in shape-controlled metal nanocrystal synthesis, with a focus on the role of PVP in colloidal Ag nanocube formation in EG solution. We examined both the thermodynamic and kinetic influences of PVP on nanocrystal shape by using appropriate computational tools. To examine the thermodynamic influence of SDAs in nanocrystal synthesis, we developed a method to calculate the interfacial free energy for multicomponent molecular solid−liquid systems.44 We applied this method to resolve the effect of PVP in altering the interfacial free energies of solution-phase, PVP-covered Ag surfaces.45 Using our calculated interfacial free energies to predict Wulff shapes of solution-phase Ag nanocrystals, we found that PVP does not significantly influence the equilibrium shapes of these nanocrystals, so they essentially assume the Wulff shape of bare Ag in vacuum.45 We then considered the question of how facet-selective SDA binding needs to be to influence equilibrium nanocrystal shapes, by employing a general model. We found that only a strong binding preference to a particular facet can induce an equilibrium shape that is significantly different from that of the bare metal. This is because the formation of the bare metal surface has the largest contribution to the metal−solvent−SDA interfacial free energy, in which the contribution of the SDA− metal interaction is small in comparison.45 We also investigated the kinetic influence of PVP on Ag nanocrystal shapes, by assuming that the linear facet growth rates of Ag nanocrystals are determined by the Ag atom deposition fluxes to the various facets.46,47 We obtained MFPTs for Ag atoms to reach Ag(100) and Ag(111) surfaces with several PVP coverages. By using MFPT ratios on the two surfaces to predict kinetic Wulff shapes, we found that {100}faceted PVP-covered Ag nanocubes are the preferred kinetic shapes when the PVP coverage is sufficiently high and the PVP chains are sufficiently long. This result is in agreement with experiment.49 Thus, we find that SDA binding must exhibit considerable selectivity to a certain facet to thermodynamically alter the shape of a nanocrystal from that of the bare metal. However, the binding selectivity needed to kinetically influence nanocrystal shape is much weaker. There has been much discussion regarding possible thermodynamic nanocrystal shapes, as has been evidenced in two recent reviews.86,87 Our results indicate that kinetic control of metal nanocrystal shapes is likely more prevalent than thermodynamic control and provide a rational basis for understanding the requirements for thermodynamic control. We note that there is much more to be done in the quest to understand shape-selective nanocrystal growth. Our work has shown that SDA-induced differences in facet deposition rates can alter kinetic nanocrystal shapes, and this idea has been embraced by others.56 However, we note that the differences in atom deposition rates that we observe here cannot explain the shapes of 5-fold-twinned Ag nanowires that grow from decahedral seeds in similar synthesis conditions.88 Preliminary calculations indicate that we would need flux ratios (cf., Figure 7) of more than 1000 to achieve experimentally observed nanowire aspect ratios. There are other possible kinetic steps that can determine nanocrystal shapes in certain synthesis protocols, such as the surface diffusion of adsorbed atoms on

buried interfaces and nanocrystal aggregation. More research is needed to understand the role of these phenomena in promoting nanocrystal shapes. It is still not fully understood which attributes of SDAs or solvent molecules lead to certain nanocrystal shapes for a given initial seed. Recent work from Wall and co-workers59 has shown that surfactants may not be needed at all in the growth of Au nanocrystals. They developed a theoretical framework for predicting linear facet growth rates based on atom diffusion and surface nucleation kinetics, which may be a promising approach for future studies. Since the complexity of the solution-phase chemistry that accompanies nanocrystal nucleation and growth rivals that of gas-phase combustion reaction networks, it is still difficult to ascertain key species or intermediates in certain synthesis protocols. It is a current theoretical challenge to describe the solution-phase chemistry associated with these syntheses in an accurate and efficient way, so understanding how seed crystals nucleate and grow into selective shapesan enormously important problem as the shape of the final crystal is often determined by the shape of the seedis still beyond our current theoretical grasp. Continued progress along these lines, with the help of detailed experimental investigations, will be fruitful in efforts to grow nanocrystal shapes by design.

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: fi[email protected]. ORCID

Kristen A. Fichthorn: 0000-0002-4256-714X Notes

The authors declare no competing financial interest. Biographies

Xin Qi graduated from the University of Iowa with a B.S.E. in Chemical and Biochemical Engineering with distinction and is currently pursuing a Ph.D. in Chemical Engineering at the Pennsylvania State University. Her Ph.D. work under Dr. Kristen Fichthorn focuses on investigating the growth mechanisms of monocrystalline and multiply twinned nanocrystals using computational methods. She is interested in using various theoretical approaches to unravel complex growth phenomena from multiple perspectives, as well as method development. G



■

Feature Article

REFERENCES

(1) Bratlie, K. M.; Lee, H.; Komvopoulos, K.; Yang, P.; Somorjai, G. A. Platinum Nanoparticle Shape Effects on Benzene Hydrogenation Selectivity. Nano Lett. 2007, 7, 3097−3101. (2) Lim, B.; Jiang, M.; Camargo, P. H. C.; Cho, E. C.; Tao, J.; Lu, X.; Zhu, Y.; Xia, Y. Pd-Pt Bimetallic Nanodendrites with High Activity for Oxygen Reduction. Science 2009, 324, 1302−1305. (3) Prieto, G.; Zečević, J.; Friedrich, H.; de Jong, K. P.; de Jongh, P. E. Towards Stable Catalysts by Controlling Collective Properties of Supported Metal Nanoparticles. Nat. Mater. 2013, 12, 34−39. (4) Michalak, W. D.; Krier, J. M.; Komvopoulos, K.; Somorjai, G. A. Structure Sensitivity in Pt Nanoparticle Catalysts for Hydrogenation of 1,3-Butadiene: in situ Study of Reaction Intermediates Using SFG Vibrational Spectroscopy. J. Phys. Chem. C 2013, 117, 1809−1817. (5) Kim, M.; Kim, Y.; Hong, J. W.; Ahn, S.; Kim, W. Y.; Han, S. W. The Facet-Dependent Enhanced Catalytic Activity of Pd Nanocrystals. Chem. Commun. 2014, 50, 9454−9457. (6) Niu, Z.; Becknell, N.; Yu, Y.; Kim, D.; Chen, C.; Kornienko, N.; Somorjai, G. A.; Yang, P. Anisotropic Phase Segregation and Migration of Pt in Nanocrystals en route to Nanoframe Catalysts. Nat. Mater. 2016, 15, 1188−1194. (7) Choi, J.-H.; Wang, H.; Oh, S. J.; Paik, T.; Sung, P.; Sung, J.; Ye, X.; Zhao, T.; Diroll, B. T.; Murray, C. B.; et al. Exploiting the Colloidal Nanocrystal Library to Construct Electronic Devices. Science 2016, 352, 205−208. (8) Jena, P. Materials for Hydrogen Storage: Past, Present, and Future. J. Phys. Chem. Lett. 2011, 2, 206−211. (9) Liu, L.; Ni, F.; Zhang, J.; Jiang, X.; Lu, X.; Guo, Z.; Xu, R. Silver Nanocrystals Sensitize Magnetic-Nanoparticle-Mediated ThermoInduced Killing of Cancer Cells. Acta Biochim. Biophys. Sin. 2011, 43, 316−323. (10) Howes, P. D.; Chandrawati, R.; Stevens, M. M. Colloidal Nanoparticles as Advanced Biological Sensors. Science 2014, 346, 1247390. (11) Lee, K.-S.; El-Sayed, M. A. Gold and Silver Nanoparticles in Sensing and Imaging: Sensitivity of Plasmon Response to Size, Shape, and Metal Composition. J. Phys. Chem. B 2006, 110, 19220−19225. (12) McLellan, J. M.; Siekkinen, A.; Chen, J.; Xia, Y. Comparison of the Surface-Enhanced Raman Scattering on Sharp and Truncated Silver Nanocubes. Chem. Phys. Lett. 2006, 427, 122−126. (13) Mayer, K. M.; Hafner, J. H. Localized Surface Plasmon Resonance Sensors. Chem. Rev. 2011, 111, 3828−3857. (14) Kamat, P. V. Meeting the Clean Energy Demand:? Nanostructure Architectures for Solar Energy Conversion. J. Phys. Chem. C 2007, 111, 2834−2860. (15) Beck, F. J.; Polman, A.; Catchpole, K. R. Tunable Light Trapping for Solar Cells using Localized Surface Plasmons. J. Appl. Phys. 2009, 105, 114310. (16) Atwater, H. A.; Polman, A. Plasmonics for Improved Photovoltaic Devices. Nat. Mater. 2010, 9, 205−213. (17) Linic, S.; Christopher, P.; Ingram, D. B. Plasmonic-Metal Nanostructures for Efficient Conversion of Solar to Chemical Energy. Nat. Mater. 2011, 10, 911−921. (18) Kholmicheva, N.; Moroz, P.; Rijal, U.; Bastola, E.; Uprety, P.; Liyanage, G.; Razgoniaev, A.; Ostrowski, A. D.; Zamkov, M. Plasmonic Nanocrystal Solar Cells Utilizing Strongly Confined Radiation. ACS Nano 2014, 8, 12549−12559. (19) Anker, J. N.; Hall, W. P.; Lyandres, O.; Shah, N. C.; Zhao, J.; Van Duyne, R. P. Biosensing with Plasmonic Nanosensors. Nat. Mater. 2008, 7, 442−453. (20) Abalde-Cela, S.; Carregal-Romero, S.; Coelho, J. P.; GuerreroMartínez, A. Recent Progress on Colloidal Metal Nanoparticles as Signal Enhancers in Nanosensing. Adv. Colloid Interface Sci. 2016, 233, 255−270. (21) Zeng, J.; Xia, X.; Rycenga, M.; Henneghan, P.; Li, Q.; Xia, Y. Successive Deposition of Silver on Silver Nanoplates: Lateral versus Vertical Growth. Angew. Chem., Int. Ed. 2011, 50, 244−249.

Tonnam Balankura received a B.Eng. from Sirindhorn International Institute of Technology, Thailand, in 2013 and a Ph.D. with Prof. Kristen A. Fichthorn at the Pennsylvania State University in 2018, both in Chemical Engineering. His area of research during his Ph.D. studies included shape-control phenomena in nanocrystal synthesis, interfacial phenomena, rare-event sampling, and molecular simulation.

Kristen Fichthorn is the Merrell Fenske Professor of Chemical Engineering and a Professor of Physics at the Pennsylvania State University. She received a B.S. in Chemical Engineering from the University of Pennsylvania in 1985 and a Ph.D. (with Robert M. Ziff and Erdogan Gulari) in Chemical Engineering from the University of Michigan in 1989. She spent one year as an IBM Postdoctoral Fellow (with W. Henry Weinberg) in the Department of Chemical Engineering at the University of California at Santa Barbara before joining the Department of Chemical Engineering at Penn State as an Assistant Professor in 1990. She was promoted to Associate Professor and Professor in 1996 and 1999, respectively. Professor Fichthorn’s research is primarily in multiscale simulation of fluid−solid interfaces, in which she applies theoretical techniques ranging from quantum density functional theory to molecular dynamics, Monte Carlo methods, and continuum theories to a diverse array of applications in nanoscale materials, thin-film and crystal growth, colloidal assembly, and wetting.

■

ACKNOWLEDGMENTS This work is funded by the Department of Energy, Office of Basic Energy Sciences, Materials Science Division, Grant Number DE-FG02-07ER46414. T.B. acknowledges training provided by the Computational Materials Education and Training (CoMET) NSF Research Traineeship (Grant Number DGE-1449785). This work used the Extreme Science and Engineering Discovery Environment (XSEDE) supported by NSF/OCI-1053575. H


Feature Article


Platinum Nanoparticle Shape Control by Polyvinylpyrrolidone. J. Phys. Chem. C 2016, 120, 7532−7542. (42) Almora-Barrios, N.; Novell-Leruth, G.; Whiting, P.; Liz-Marzán, L. M.; López, N. Theoretical Description of the Role of Halides, Silver, and Surfactants on the Structure of Gold Nanorods. Nano Lett. 2014, 14, 871−875. (43) Greenwood, A. R.; Vörös, M.; Giberti, F.; Galli, G. Emergent Electronic and Dielectric Properties of Interacting Nanoparticles at Finite Temperature. Nano Lett. 2018, 18, 255−261. (44) Qi, X.; Zhou, Y.; Fichthorn, K. A. Obtaining the Solid-Liquid Interfacial Free Energy via Multi-Scheme Thermodynamic Integration: Ag-Ethylene Glycol Interfaces. J. Chem. Phys. 2016, 145, 194108. (45) Qi, X.; Fichthorn, K. A. Theory of the Thermodynamic Influence of Solution-Phase Additives in Shape-Controlled Nanocrystal Synthesis. Nanoscale 2017, 9, 15635−15642. (46) Qi, X.; Balankura, T.; Zhou, Y.; Fichthorn, K. A. How StructureDirecting Agents Control Nanocrystal Shape: PolyvinylpyrrolidoneMediated Growth of Ag Nanocubes. Nano Lett. 2015, 15, 7711−7717. (47) Balankura, T.; Qi, X.; Zhou, Y.; Fichthorn, K. A. Predicting Kinetic Nanocrystal Shapes through Multi-Scale Theory and Simulation: Polyvinylpyrrolidone-Mediated Growth of Ag Nanocrystals. J. Chem. Phys. 2016, 145, 144106. (48) Fichthorn, K. A.; Balankura, T.; Qi, X. Multi-Scale Theory and Simulation of Shape-Selective Nanocrystal Growth. CrystEngComm 2016, 18, 5410−5417. (49) Xia, X.; Zeng, J.; Oetjen, L. K.; Li, Q.; Xia, Y. Quantitative Analysis of the Role Played by Poly(vinylpyrrolidone) in SeedMediated Growth of Ag Nanocrystals. J. Am. Chem. Soc. 2012, 134, 1793−1801. (50) Sun, Y. G.; Yin, Y. D.; Mayers, B. T.; Herricks, T.; Xia, Y. Uniform Silver Nanowires Synthesis by Reducing AgNO3 with Ethylene Glycol in the Presence of Seeds and Poly(vinyl pyrrolidone). Chem. Mater. 2002, 14, 4736−4745. (51) Xia, Y.; Xiong, Y.; Lim, B.; Skrabalak, S. E. Shape-Controlled Synthesis of Metal Nanocrystals: Simple Chemistry Meets Complex Physics? Angew. Chem., Int. Ed. 2009, 48, 60−103. (52) Zhou, Y.; Saidi, W. A.; Fichthorn, K. A. A Force Field for Describing the Polyvinylpyrrolidone-Mediated Solution-Phase Synthesis of Shape-Selective Ag Nanoparticles. J. Phys. Chem. C 2014, 118, 3366−3374. (53) Williams, P. L.; Mishin, Y.; Hamilton, J. C. An Embedded-Atom Potential for the Cu-Ag System. Modell. Simul. Mater. Sci. Eng. 2006, 14, 817−833. (54) Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; Mackerell, A. D. CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields. J. Comput. Chem. 2010, 31, 671−690. (55) Liu, S.-H.; Fichthorn, K. Interaction of Alkylamines with Cu Surfaces: A Metal-Organic Many-Body Force Field. J. Phys. Chem. C 2017, 121, 22531. (56) Meena, S. K.; Sulpizi, M. Understanding the Microscopic Origin of Gold Nanoparticle Anisotropic Growth from Molecular Dynamics Simulations. Langmuir 2013, 29, 14954−14961. (57) Meena, S. K.; Celiksoy, S.; Schäfer, P.; Henkel, A.; Sönnichsen, C.; Sulpizi, M. The Role of Halide Ions in the Anisotropic Growth of Gold Nanoparticles: a Microscopic, Atomistic Perspective. Phys. Chem. Chem. Phys. 2016, 18, 13246−13254. (58) Meena, S. K.; Sulpizi, M. From Gold Nanoseeds to Nanorods: The Microscopic Origin of the Anisotropic Growth. Angew. Chem., Int. Ed. 2016, 55, 11960−11964. (59) Wall, M. A.; Harmsen, S.; Pal, S.; Zhang, L.; Arianna, G.; Lombardi, J. R.; Drain, C. M.; Kircher, M. F. Surfactant-Free Shape Control of Gold Nanoparticles Enabled by Unified Theoretical Framework of Nanocrystal Synthesis. Adv. Mater. 2017, 29, 1605622. (60) Ramezani-Dakhel, H.; Ruan, L.; Huang, Y.; Heinz, H. Molecular Mechanism of Specific Recognition of Cubic Pt Nanocrystals by Peptides and of the Concentration-Dependent Formation from Seed Crystals. Adv. Funct. Mater. 2015, 25, 1374−1384.

(22) Xia, X.; Zeng, J.; Zhang, Q.; Moran, C. H.; Xia, Y. Recent Developments in Shape-Controlled Synthesis of Silver Nanocrystals. J. Phys. Chem. C 2012, 116, 21647−21656. (23) da Silva, R. R.; Yang, M.; Choi, S.-I.; Chi, M.; Luo, M.; Zhang, C.; Li, Z. Y.; Camargo, P. H. C.; Ribeiro, S. J. L.; Xia, Y. Facile Synthesis of Sub-20 nm Silver Nanowires through a BromideMediated Polyol Method. ACS Nano 2016, 10, 7892−7900. (24) Zhang, J.; Langille, M. R.; Personick, M. L.; Zhang, K.; Li, S.; Mirkin, C. A. Concave Cubic Gold Nanocrystals with High-Index Facets. J. Am. Chem. Soc. 2010, 132, 14012−14014. (25) Casillas, G.; Velázquez-Salazar, J. J.; Jose-Yacaman, M. A New Mechanism of Stabilization of Large Decahedral Nanoparticles. J. Phys. Chem. C 2012, 116, 8844−8848. (26) Rathmell, A. R.; Wiley, B. J. The Synthesis and Coating of Long, Thin Copper Nanowires to Make Flexible, Transparent Conducting Films on Plastic Substrates. Adv. Mater. 2011, 23, 4798−4803. (27) Jin, M.; He, G.; Zhang, H.; Zeng, J.; Xie, Z.; Xia, Y. ShapeControlled Synthesis of Copper Nanocrystals in an Aqueous Solution with Glucose as a Reducing Agent and Hexadecylamine as a Capping Agent. Angew. Chem., Int. Ed. 2011, 50, 10560−10564. (28) Luo, M.; Zhou, M.; Rosa da Silva, R.; Tao, J.; Figueroa Cosme, L.; Gilroy, K. D.; Peng, H.-C.; He, Z.; Xia, Y. Pentatwinned Cu Nanowires with Ultrathin Diameters below 20 nm and Their Use as Templates for the Synthesis of Au-Based Nanotubes. ChemNanoMat 2017, 3, 190−195. (29) Kang, Y.; Pyo, J. B.; Ye, X.; Diaz, R. E.; Gordon, T. R.; Stach, E. A.; Murray, C. B. Shape-Controlled Synthesis of Pt Nanocrystals: The Role of Metal Carbonyls. ACS Nano 2013, 7, 645−653. (30) Liao, H.-G.; Zherebetskyy, D.; Xin, H.; Czarnik, C.; Ercius, P.; Elmlund, H.; Pan, M.; Wang, L.-W.; Zheng, H. Facet Development during Platinum Nanocube Growth. Science 2014, 345, 916−919. (31) Ye, H.; Liu, Y.; Chhabra, A.; Lilla, E.; Xia, X. Polyvinylpyrrolidone (PVP)-Capped Pt Nanocubes with Superior Peroxidase-Like Activity. ChemNanoMat 2017, 3, 33−38. (32) Xiong, Y.; Xia, Y. Shape-Controlled Synthesis of Metal Nanostructures: The Case of Palladium. Adv. Mater. 2007, 19, 3385−3391. (33) Lim, B.; Wang, J.; Camargo, P. H. C.; Cobley, C. M.; Kim, M. J.; Xia, Y. Twin-Induced Growth of Palladium−Platinum Alloy Nanocrystals. Angew. Chem., Int. Ed. 2009, 48, 6304−6308. (34) Kilin, D. S.; Prezhdo, O. V.; Xia, Y. Shape-Controlled Synthesis of Silver Nanoparticles: Ab initio Study of Preferential Surface Coordination with Citric Acid. Chem. Phys. Lett. 2008, 458, 113−116. (35) Mdluli, P. S.; Sosibo, N. M.; Revaprasadu, N.; Karamanis, P.; Leszczynski, J. Surface Enhanced Raman Spectroscopy (SERS) and Density Functional Theory (DFT) Study for Understanding the Regioselective Adsorption of Pyrrolidinone on the Surface of Silver and Gold Colloids. J. Mol. Struct. 2009, 935, 32−38. (36) Al-Saidi, W. A.; Feng, H.; Fichthorn, K. A. Adsorption of Polyvinylpyrrolidone on Ag Surfaces: Insight into a StructureDirecting Agent. Nano Lett. 2012, 12, 997−1001. (37) Saidi, W. A.; Feng, H.; Fichthorn, K. A. Binding of Polyvinylpyrrolidone to Ag Surfaces: Insight into a Structure-Directing Agent from Dispersion-Corrected Density Functional Theory. J. Phys. Chem. C 2013, 117, 1163−1171. (38) Zhou, Y.; Saidi, W. A.; Fichthorn, K. A. Comparison of the Binding of Polyvinylpyrrolidone and Polyethylene Oxide to Ag Surfaces: Elements of a Successful Structure-Directing Agent. J. Phys. Chem. C 2013, 117, 11444−11448. (39) Liu, S.-H.; Saidi, W. A.; Zhou, Y.; Fichthorn, K. A. Synthesis of 111-Faceted Au Nanocrystals Mediated by Polyvinylpyrrolidone: Insights from Density-Functional Theory and Molecular Dynamics. J. Phys. Chem. C 2015, 119, 11982−11990. (40) Liu, S.-H.; Balankura, T.; Fichthorn, K. A. Self-Assembled Monolayer Structures of Hexadecylamine on Cu Surfaces: DensityFunctional Theory. Phys. Chem. Chem. Phys. 2016, 18, 32753−32761. (41) Ye, J.-Y.; Attard, G. A.; Brew, A.; Zhou, Z.-Y.; Sun, S.-G.; Morgan, D. J.; Willock, D. J. Explicit Detection of the Mechanism of I


Feature Article


(83) Hänggi, P.; Talkner, P.; Borkovec, M. Reaction-Rate Theory: Fifty Years after Kramers. Rev. Mod. Phys. 1990, 62, 251. (84) Woolf, T. B.; Roux, B. Conformational Flexibility of oPhosphorylcholine and o-Phosphorylethanolamine: A Molecular Dynamics Study of Solvation Effects. J. Am. Chem. Soc. 1994, 116, 5916−5926. (85) Hummer, G. Position-Dependent Diffusion Coefficients and Free Energies from Bayesian Analysis of Equilibrium and Replica Molecular Dynamics Simulations. New J. Phys. 2005, 7, 34. (86) Xia, Y.; Xia, X.; Peng, H. C. Shape-Controlled Synthesis of Colloidal Metal Nanocrystals: Thermodynamic versus Kinetic Products. J. Am. Chem. Soc. 2015, 137, 7947−7966. (87) Marks, L. D.; Peng, L. Nanoparticle Shape, Thermodynamics and Kinetics. J. Phys.: Condens. Matter 2016, 28, 053001. (88) Wiley, B.; Sun, Y.; Mayers, B.; Xia, Y. Shape-Controlled Synthesis of Metal Nanostructures: The Case of Silver. Chem. - Eur. J. 2005, 11, 454−463.

(61) Heinz, H.; Pramanik, C.; Heinz, O.; Ding, Y.; Mishra, R. K.; Marchon, D.; Flatt, R. J.; Estrela-Lopis, I.; Llop, J.; Moya, S.; et al. Nanoparticle Decoration with Surfactants: Molecular Interactions, Assembly, and Applications. Surf. Sci. Rep. 2017, 72, 1−58. (62) von Laue, M. Der Wulffsche Satz fur die Gleichgewichtsform von Kristallen. Z. Kristallogr. - Cryst. Mater. 1943, 105, 124−133. (63) Dinghas, A. Uber Einen Geometrischen sat von Wulff fur die Gleichgewichtsform von Kristallen. Z. Kristallogr. - Cryst. Mater. 1943, 105, 304−314. (64) Tran, R.; Xu, Z.; Radhakrishnan, B.; Sun, W.; Persson, K. A.; Ong, S. P.; Winston, D. Surface Energies of Elemental Crystals. Sci. Data 2016, 3, 160080. (65) Good, R. J.; Girifalco, L. A. A Theory for Estimation of Surface and Interfacial Energies. III. Estimation of Surface Energies of Solids from Contact Angle Data. J. Phys. Chem. 1960, 64, 561. (66) Owens, D. K.; Wendt, R. C. Estimation of the Surface Free Energy of Polymers. J. Appl. Polym. Sci. 1969, 13, 1741. (67) Louhichi, A.; Tamborini, E.; Ghofraniha, N.; Caton, F.; Roux, D.; Oberdisse, J.; Cipelletti, L.; Ramos, L. Nucleation and Growth of Micellar Polycrystals under Time-Dependent Volume Fraction Conditions. Phys. Rev. E 2013, 87, 032306. (68) Bokeloh, J.; Rozas, R. E.; Horbach, J.; Wilde, G. Nucleation Barriers for the Liquid-to-Crystal Transition in Ni: Experiment and Simulation. Phys. Rev. Lett. 2011, 107, 145701. (69) Bokeloh, J.; Wilde, G.; Rozas, R. E.; Benjamin, R.; Horbach, J. Nucleation Barriers for the Liquid-to-Crystal Transition in Simple Metals: Experiment vs. Simulation. Eur. Phys. J.: Spec. Top. 2014, 223, 511. (70) Nguyen, V. D.; Hu, Z.; Schall, P. Single cCrystal Growth and Anisotropic Crystal-Fluid Interfacial Free Energy in Soft Colloidal Systems. Phys. Rev. E 2011, 84, 011607. (71) Leroy, F.; Müller-Plathe, F. Solid-Liquid Surface Free Energy of Lennard-Jones Liquid on Smooth and Rough Surfaces Computed by Molecular Dynamics using the Phantom-Wall Method. J. Chem. Phys. 2010, 133, 044110. (72) Broughton, J. Q.; Gilmer, G. H. Molecular Dynamics Investigation of the Crystal−Fluid Interface. VI. Excess Surface Free Energies of Crystal−Liquid Systems. J. Chem. Phys. 1986, 84, 5759− 5768. (73) Davidchack, R. L.; Laird, B. B. Direct Calculation of the Crystal−Melt Interfacial Free Energies for Continuous Potentials: Application to the Lennard-Jones System. J. Chem. Phys. 2003, 118, 7651−7657. (74) Benjamin, R.; Horbach, J. Crystal-Liquid Interfacial Free Energy via Thermodynamic Integration. J. Chem. Phys. 2014, 141, 044715. (75) Ardham, V. R.; Deichmann, G.; van der Vegt, N. F. A.; Leroy, F. Solid-liquid Work of Adhesion of Coarse-Grained Models of n-Hexane on Graphene Layers Derived from the Conditional Reversible Work Method. J. Chem. Phys. 2015, 143, 243135. (76) Shuttleworth, R. The Surface Tension of Solids. Proc. Phys. Soc., London, Sect. A 1950, 63, 444−457. (77) Dingreville, R.; Qu, J. Interfacial Excess Energy, Excess Stress and Excess Strain in Elastic Solids: Planar Interfaces. J. Mech. Phys. Solids 2008, 56, 1944−1954. (78) Sun, Y.; Xia, Y. Large-Scale Synthesis of Uniform Silver Nanowires Through a Soft, Self-Seeding Polyol Process. Adv. Mater. 2002, 14, 833−837. (79) Zeng, J.; Zheng, Y.; Rycenga, M.; Tao, J.; Li, Z.-Y.; Zhang, Q.; Zhu, Y.; Xia, Y. Controlling the Shapes of Silver Nanocrystals with Different Capping Agents. J. Am. Chem. Soc. 2010, 132, 8552−8553. (80) Zhang, Y. C.; Sizemore, J. P.; Doherty, M. F. Shape Evolution of 3-Dimensional Faceted Crystals. AIChE J. 2006, 52, 1906−1915. (81) Torrie, G. M.; Valleau, J. P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187−199. (82) Kästner, J.; Thiel, W. Bridging the Gap between Thermodynamic Integration and Umbrella Sampling Provides a Novel Analysis Method: Umbrella Integration. J. Chem. Phys. 2005, 123, 144104. J


Theoretical Perspectives on the Influence of Solution-Phase Additives

Recommend Documents