Christopher B.
Whitehead
and
Richard G.
Finke
*
Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA. E-mail: Richard.Finke@colostate.edu
First published on 27th September 2021
Following a brief description of synchrotron X-ray absorption fine structure (XAFS) spectroscopy and small/wide angle X-ray scattering (SAXS/WAXS), the definition of and four primary criteria for attaining reliable, disproof-based chemical mechanisms of particle formation are given. A total of 74 papers using synchrotron techniques for mechanistic investigation are then analyzed in detail via the construction of four tables provided in the ESI that analyze each of the 74 papers. Six primary case studies are selected out of the 74 total papers for presentation in greater detail in the main text, specifically the illustrative case histories of: (i) palladium nanoparticles studied using SAXS, (ii) rhodium nanocubes examined using XAFS, (iii) iridium nanoparticles studied by XAFS and SAXS, (iv) gold nanoparticles studied by XAFS, SAXS, and XRD, (v) cadmium-selenide nanocrystals studied by XAFS, and (vi) zinc-oxide nanoparticles studied using SAXS/WAXS/UV-vis. Additionally, two shorter case studies are presented that address particle formation shapes: tungstite nanoplatelets and copper nanocrystals. Two summary tables are presented in the main text that present the current state of disproof-based, deliberately minimalistic particle formation mechanisms expressed in terms of generalized, pseudoelementary steps. Finally, a Conclusions section with nine takeaways and an Outlook section are also provided. The goal of the present review is to expedite the use of powerful synchrotron SAXS and XAFS studies to provide reliable, disproof-based chemical mechanisms and their associated quantitative differential equations that can then be employed to predict particle-size distributions via the recent development of Mechanism-Enabled Population Balance Modeling (ME-PBM).
Due to the sub-nanoscale size (<1 nm) of just-formed nuclei, it is very difficult to observe directly in real time the smallest, kinetically first-formed cluster(s), termed the kinetically effective nucleus (KEN).30 Indeed, recent work shows that the KEN consists of just 2–3 atoms in at least some strong-bonding systems.31,32 Hence, there is a need to use synchrotron-based spectroscopic techniques to monitor the nucleation, growth and any agglomeration of particles, notably synchrotron X-ray absorption fine structure (XAFS) spectroscopy (that provides oxidation state, coordination number, and associated structural information) as a function of time, and synchrotron small-angle X-ray scattering (SAXS) (that allows monitoring of particle radius and, significantly, the number of particles) in real-time.33–38 These techniques can be used in tandem with each other or with other methods, and they can be used to monitor directly particle formation reactions in situ.
However, significant gaps remain in the nanoparticle formation literature en route to a deeper understanding of (i) first, what constitutes a reliable, disproof-based chemical reaction mechanism; (ii) the experimental data needed to be able to claim such a mechanism; (iii) how such mechanisms differ from non-disproof-based model(s); (iv) the preferred experimental methods for collecting the needed kinetics data; and then (v) what are the arguably best, prototype mechanistic case studies to date by synchrotron-based methods?
Hence herein we begin by discussing the requirements for a reliable, disproof-based chemical mechanism and then we present six case studies of particle formation monitored by synchrotron-based SAXS, XAFS or tandem SAXS/XAFS methods.39–44 We systemically assess the results and strengths of each case study and, where needed, suggest additional studies en route to a reliable reaction mechanism. The goal is to use these six case studies as pedagogically valuable examples of how to most efficiently achieve reliable, disproof-based reaction mechanisms of particle formation regardless of the exact type or composition of the resultant particles. Two summary tables that present the current state-of-the-art in terms of disproof-based, deliberately minimalistic (i.e., Ockham's razor obeying) particle formation mechanisms are provided as part of the summary and conclusions.
(1) If a complete mass- and charge-balanced reaction stoichiometry is experimentally determined, as that stoichiometry is what the proposed steps of the mechanism must add up to.48 Problems result if even a world class expert mechanistic chemist49 studies a much simpler reaction than nanoparticle formation in the absence of an experimentally established, full, balanced stoichiometry.50
(2) If kinetics data are obtained,51,52 ideally the full rate law, and ideally using multiple, direct physical methods and handles;
(3) If elementary, or if not possible pseudoelementary,53,54 step reactions are written that present the stepwise mechanism that also add up to the experimentally established reaction stoichiometry. Those mechanistic steps also define the rate constants and define the unambiguous words55,56 one can use to discuss the observed chemistry—the latter point being much more important than one might initially imagine for avoiding confusion in the description of the mechanism and associated chemistry;
(4) If, critically, disproof-based57–59 quantitative comparisons of the experimental data to competing, alternative mechanisms are made, specifically how well (or poorly) a given postulated mechanism is able to fit that data. Adherence to Ockham's razor is also important en route to deliberately minimalistic mechanisms that are the initial goal of mechanistic studies.60 Those minimum mechanisms are expected to be replaced—disproved if you like—by more complex mechanisms as more data, from better experiments or improved computational methods, become available over time.
Once one has the minimum, disproof-based chemical mechanism in hand, then one has a good start on a broader and more generalized understanding of the system. That mechanistic understanding of nucleation, growth, and agglomeration that are ubiquitous processes across nature can then be used to further control key properties of the system: average particle size, particle-size distribution, catalytic activity, photochemical properties (that are often size and size-dispersion dependent, e.g., in semiconductor quantum dot nanoparticles), and so on.
Unfortunately, experimental kinetics and rate-law-backed mechanisms for nucleation are rare at the molecular and atomic level, save for a small handful of studies.26,30–32,39,41 The precise mechanism of nucleation is critical as it starts off the particle-formation process; without the precise details of the nucleation process one cannot possibly claim to a complete nor reliable mechanism of particle formation. Growth and agglomeration are in a bit better shape,61–66 but still not well supported kinetically and mechanistically with, again, the required experimental rate law for growth or agglomeration. Hence, advances are needed there as well. Through the basic four requirements outlined above en route to a reliable mechanism illustrated via six case studies, we hope to demystify how one can perform reliable mechanistic investigations that allow elucidation of the most probable particle formation process. We also detail the importance of pseudoelementary steps53,54 in the construction, and attempted disproof, of more complex mechanisms in materials and other chemistries.
With the 2019 advent of Mechanism-Enabled Population Balance Modeling (ME-PBM),67,68 one can track every particle in a proposed particle-formation pathway consisting of even thousands of elementary steps. Disproof-based, deliberately minimum mechanisms are necessary and sufficient to be able to predict both average particle sizes and particle-size distributions (PSDs), including the PSD shape, with ME-PBM. The available evidence to date also illustrates the ability of disproof-based ME-PBM to confirm—or more often to date to refute67,68—a proposed mechanism67,68 en route to a refined mechanism. ME-PBM is also able to extract rate constants for the proposed, pseudoelementary step mechanism from PSD data.67,68 But, to exploit the full generality, predictability and for knowledge transfer to other systems or conditions, a reliable minimum reaction mechanism is required to start, as the name Mechanism-Enabled PBM indicates. Ligand effects are, once again, just now being taken into account as part of ME-PBM,31,32,39vide infra.
The interested reader is encouraged—actually strongly encouraged—to study and analyze for themselves these four tables to understand the variety of systems covered, the status of the approaches, and the overall efforts at nanoparticle reaction development and associated mechanistic investigation summarized in Tables S1–S4 of the ESI.† Indeed, if Tables S1–S4 of the ESI† covered anything smaller than their current 50 pages, we would have placed them upfront in this review to make them even more available and apparent. Self-study of Tables S2–S4 of the ESI† will generate the reader's own assessment and insights regarding the status of the field of particle formation synthesis, kinetics and mechanism using synchrotron methods. Papers that did not present, discuss, or claim a mechanism are not included even if they reported XAFS or SAXS data.
In what follows we have selected 6 case studies for a closer look and critical analysis based, overall, on how they illustrate together the state-of-the-art in the field as well as what else is needed to attain reliable chemical mechanisms for the particle formation reaction under study. Our apologies in advance to authors of the many other interesting studies summarized in ESI,† Tables S2–S4 that either space, or our approach and organization of this review, did not permit us to cover in detail in the main text that follows. Also not covered herein due to scope and space limitations are computational chemistry approaches and contributions that promise to be of increasing importance, perhaps especially to a deeper understanding of nucleation processes.69,70 That said, computations not carefully connected to experimental results can lead to erroneous conclusions, even for chemical systems much simpler than particle formation reactions.71
The four case studies are: (Section 3.1) palladium nanoparticle formation monitored by SAXS;39 (Section 3.2) rhodium nanocube formation monitored by XAFS;40 (Section 3.3) iridium nanoparticle formation monitored by XAFS and SAXS;41 and (Section 3.4) gold nanosphere and nanowire formation monitored by XAFS, SAXS, and XRD.42 Each case study is organized by: (i) a summary of the system and techniques used to study it; (ii) a review of the key kinetics data; (iii) the authors’ proposed formation model or mechanism; and (iv) an analysis of the case's results as compared to the four components required for establishing a minimum, disproof-based mechanism. Overall, our intent is to be supportive of the excellent synchrotron-based studies and encouraging to how those studies can be taken to the next level of mechanistic formulation, analysis, and use.
(1) |
The authors directly monitored the Pd nanoparticle formation by in situ synchrotron SAXS experiments—while making the important point that SAXS counts as “two methods” if one monitors both the particle radius as well as the number of particles vs. time,39,72vide infra. By using a Schultz polydisperse spherical model, they were able to fit the SAXS data and simultaneously extract both the particle sizes and the total number of particles with respect to time.
The data display direct evidence of slow, continuous nucleation with fast (but slower than diffusion-controlled) growth, results in direct contradiction to the classical LaMer model of 1950 that postulates burst nucleation followed by diffusion-controlled growth.39,65,66 Furthermore, the continuous increase in particles suggests that agglomeration, which would reduce the number of particles, is not involved in the Pd(0)n-particle formation process under their conditions that include an excess of strongly coordinating tri-octylphosphine ligand. The authors hypothesized that, at a minimum, continuous nucleation, fast autocatalytic growth, and the incorporation of ligand effects will be needed for any kinetic model and associated minimum mechanism able to fit their SAXS data.
The authors proceeded by constructing a model relying on three assumptions: (i) that the fraction of surface versus core atoms can be ignored, which for particles smaller than 1.5 nm (as seen in Fig. 1) means that the percentage of surface atoms is >75% of the total atoms;39 (ii) that the growth and ligand binding to the particle surface is independent of particle size and ligand coverage, in this first model to be able to account for ligand effects. (“In the absence of experimental or theoretical information, the effects of size, polydispersity and ligand coverage on the rate constants are not included in the model.”39); and (iii) that the ligand–metal precursor binding ratio is assumed to be 1 because isothermal titration calorimetry73 shows, “the equilibrium binding constant for the second TOP binding (A·L + L ⇌ A·L2) is around two orders of magnitude lower than the first binding”,39 where A is the Pd precursor and L is TOP.
Fig. 1 (a) SAXS data at different reaction times after absolute scaling with fitting by Schultz polydisperse spherical model; (b) particle size evolution (blue circles) and number of particles (open orange squares) plotted as a function of time (in seconds). The Pd nanoparticles were prepared from a solution of 10 mM Pd(OAc)2, 20 mM TOP, in a 50:50 solution of toluene and 1-hexanol at 100 °C. Figure reproduced with permission from ref. 39. Copyright 2017 Royal Society of Chemistry. |
By using their proposed mechanism under the above reasonable if not necessary assumptions for an initial, carefully constructed, not over-parameterized model, the authors developed a kinetic model and associated minimum mechanism able to fit quantitatively their experimentally determined number of particles as a function of time data. The proposed kinetics model for nucleation, autocatalytic growth, and accompanying ligand effects used to fit the experimental data in Fig. 2 is given in Scheme 1. The model was constructed using the minimum number of pseudoelementary steps necessary to describe the reaction, that is, by ensuring that Ockham's razor60 is obeyed in their proposed mechanism.
Fig. 2 (a) Experimentally determined concentration of Pd atoms in nanoparticles (molarity, M) as a function of time (seconds, s) is plotted (red open circles) and fit with the kinetic model (black line) accounting for Pd atoms as total nanoparticles (B + B·L). The resultant rate constants are: k1 = 2.45 × 10−5 s−1, k2 = 8.5 × 10−1 s−1 M−1, k3 = 7.9 × 10−3 s−1 M−1, k4 = 2.1 × 10−1 s−1 M−1, Keq,3 = 2.18 × 101 M−1, and Keq,4 = 1.27 × 103 M−1. Based on the kinetic model rate constants, concentrations of A (red line), B (light blue line), A·L (green line), and B·L (dark blue line) were simulated as a function of time. (b) Experimentally determined number of particles as a function of time fit with the kinetic model. The reaction conditions for the experimental data were 10 mM Pd(OAc)2, 50:50 toluene:hexanol solvent, TOP:Pd ratio of 2:1, and 100 °C reaction temperature. Figure reproduced with permission from ref. 39. Copyright 2017 Royal Society of Chemistry. |
Scheme 1 Karim and co-workers39 4-step model for particle formation, autocatalytic surface growth, ligand–precursor interactions (A·L), and ligand–particle interactions (B·L). |
The model contains four pseudoelementary steps that define the corresponding rate constants. The first step is nucleation including reduction (A → B, k1), where Pd(II) is reduced to Pd(0) and Pd(0)n nuclei are formed. Note, as discussed elsewhere,30 under constant initial [A]Total, higher-order nucleation nA → nB (= Bn) is kinetically hidden because (d[A0]/dt)/n = k1[A0]n, where [A]0 is often constant to a ≥99.9% level during the induction period.30 In the present case, unimolecular rate-determining Pd(II) → Pd(0) is arguably likely present given that the slow step is the release of A, that is A − L → A + L. The second step employed is autocatalytic surface growth (A + B → 2B, k2),54 where Pd(II) is reduced and added to the surface of the growing Pd(0)n particle. Next, the equilibrium between the Pd(II) precursor complex and the ligand, TOP, (A + L ⇌ A·L, k3 and Keq,3), and the equilibrium between the Pd(0)n particle and TOP (B + L ⇌ B·L, k4 and Keq,4) are included. The kinetic model's ability to account quantitatively for the experimental data in Fig. 2 demonstrates the value of obtaining simultaneous size and number of particles data by SAXS as Prof. Karim has insightfully emphasized.39,72 The quantitative fits presented Fig. 2 argue strongly that the ligand-based kinetics model, incorporating the species A, B, A·L, B·L, and L, is a satisfactory, quantitative kinetics model for describing this valuable Pd(0)n particle-formation system from the Karim laboratory.
The model can be further visualized using the pictorial representation provided by the authors shown in Fig. 3:
Fig. 3 Schematic representation of the ligand-based model. The figure represents the 4 steps of the model: (1) reduction of Pd(II) to Pd(0) and the formation of the Pd nucleus (A → B); (2) autocatalytic growth of the Pd nucleus by further addition and reduction of Pd(II) (A + B → 2B); (3) reversible ligand binding with the precursor Pd(II) (A + L ⇌ A·L); and (4) reversible ligand binding to the particle surface (B + L ⇌ B·L). Orange circles are the unreduced Pd(II), small black circles with green tails are trioctylphosphine (TOP, ligand), and blue circles are reduced Pd(0)n particles. Figure reproduced with permission from ref. 39. Copyright 2017 Royal Society of Chemistry. |
Noteworthy here is that the Karim and co-authors’ contribution is ahead of most in that it provides a balanced reaction, the detailed kinetic scheme in Scheme 1, quantitative fitting of the SAXS data, and the resultant quantitative rate constants that can be used to quantitatively describe and predict particle size and size-distributions via Population-Balance Modeling. Furthermore, the Karim group tested alternative mechanisms as is needed en route to a reliable, disproof-based mechanism. In the ESI of ref. 39, five models were tested. Included in these models are mechanisms that did not include ligand steps (A → B and A + B → 2B), included only the A·L ligand step (A → B, A + B → 2B, and A + L ⇌ A·L), and included only the B·L ligand step (A → B, A + B → 2B, and B + L ⇌ B·L). These three models, the first of which is the classic 2-step mechanism of just A → B, A + B → 2B,54 were unable to fit the data. Out of the 5 models tested, the only model capable of accounting for the experimental data was the four-step model presented in Scheme 1. Here, the authors disproved four alternative mechanisms en route to the Ockham's Razor60 obeying mechanistic model in Scheme 1 that can quantitatively describe their Pd(0)n nanoparticle system.
Additionally, the authors demonstrated the applicability of their minimal ligand-based kinetic model to other literature systems that exhibit a particle-size dependence based on the ligand concentration employed. They found that their ligand-based kinetic model, that includes reversible ligand binding to the precursor (“A”) and the particle surface (“B”), fit literature data sets for Pd/PVP74 and Au/thiol75 systems as well—a valuable demonstration of the broader generality of their kinetics model.
The authors’ work highlights one of the most important pieces of this mechanistic investigation—the power of pseudoelementary steps.53,54 A pseudoelementary step is a composite of underlying elementary steps,53,54 even if those exact steps are not known to start. Nanoparticle formation consists of at least hundreds (for smaller particles) to often thousands of elementary steps—and millions of elementary steps once one reaches micron-sized particles. In the case of Karim's Pd system, the final particles contain 100s of atoms in a mechanism that must, therefore, correspond minimally to hundreds of elementary steps. Yet, Karim and coworkers were able to describe the Pd nanoparticle formation39 process, and the formation process of others systems,74,75 quantitatively using only 4 pseudoelementary steps as shown in Scheme 1—an instructive demonstration of the pseudoelementary step concept.53,54
Available elsewhere is the most illustrative example presently available where the pseudoelementary step concept has been used to elucidate the intimate, elementary-step catalytic mechanism.48 A review of the pseudoelementary step method for approaching the mechanism of complex reactions is in progress and will be published in due course. There is of course ambiguity in what especially “B” is in Scheme 1 and in general in use of the PEStep concept—B really being a sum of different nanoparticles en route to the final nanoparticles, also grouped under B. But the case illustrated elsewhere48 shows how one can deconvolute B in favorable cases to its precise specie(s)—and indeed in that case, that only by using the PEStep concept and an initial fit to a simple A → B, k1, A + B → 2B, k2, PEStep mechanism were the otherwise stumped authors able to deduce the true catalytic mechanism. That mechanism was then verified by following directly by 1H NMR four the five reaction reagents and products and showing they fit the kinetics predicted by the deduced, elementary step mechanism.48
In summary of this first important case study, the authors used SAXS to simultaneously monitor the total number of particles and the concentration of Pd nanoparticles as a function of time. Their mechanistic model satisfies three of the four criteria for a reliable mechanistic study and partially satisfies the first criteria. For criterion (i)—a complete, balanced reaction stoichiometry—“A” is a complex of the Pd starting material, Pd(OAc)2, the solvent (toluene or pyridine), and/or hexanol,39 so that the precise composition of “A” is not known unequivocally and will be needed for a more detailed understanding of the more intimate nucleation mechanism. The precise compositions of the products, PdnLm(solvent)a, are not known, but proved to not be a major hindrance to the study and its goals. The authors’ TEM studies show that an increase in solvent polarity results in an increase in particle size, yet other studies76,77 find the opposite trend, so that the exact composition of “A” and the resulting effect of the choice of solvent remain topics of interest and potential future study in this interesting system. Second, the authors fully satisfied criterion (ii)—collection of kinetics data—with their direct SAXS kinetics data that yields two observables, particle-size and concentration.39,72 Finally, the authors also satisfied criteria (iii) and (iv)—(pseudo)-elementary steps present the mechanism and disproof of multiple alternative hypotheses, respectively—leading to their development of a minimal, pseudoelementary 4-step kinetic model. The model accurately describes the physical processes, defines the rate constants, and importantly provides the correct words for describing, unambiguously, each pseudoelementary step. The final 4-step kinetics model proposed is the only mechanism out of several tested to account for the experimental data and is an Ockham's Razor-obeying model.
Overall, the excellent work39 by Karim and coworkers is an important case study of nanoparticle formation employing SAXS monitoring and proper mechanistic model building. It is a “must-read” in our opinion for anyone pursuing reliable particle-formation mechanisms and certainly a must read for researchers hoping to understand ligand-effects in their particle formation and stabilization reaction.
(2) |
To begin, the authors directly monitored the EXAFS intensity as a function of time, Fig. 4.
Fig. 4 FT-EXAFS peak intensities as a function of time are shown for Rh–Br (black), Rh–Rh first coordination shell (red), and Rh–Rh higher-order coordination shells (blue). Overlain are four detectable stages for the formation of Rh nanocrystals for the RhCl3–TTAB–PVP–EG system labeled by the authors as: ligand exchange (stage I), Rh2–4 nuclei formation (stage II), nanocrystal growth (stage III), and transformation to high-quality nanocubes (stage IV). Figure reproduced with permission from ref. 40. Copyright 2012 American Chemical Society. |
The authors observed three points where the EXAFS data changed and thereby identified four “Stages” in the reaction. During the first stage (Stage I), they observed a sharp increase in Rh–Br bonds. The second stage (Stage II) was highlighted by a decrease in Rh–Br bonds, an increase in first shell Rh–Rh bonds, but no change in higher shell Rh–Rh interactions. In the third stage (Stage III), they saw a continued decrease in Rh–Br bonds, a sharp increase in first shell Rh–Rh bonds, and a steady increase in higher shell Rh–Rh interactions. The final stage (Stage IV) was characterized by a leveling off for all species. The EXAFS intensity versus time data is consistent with four kinetic stages of the Rh particle formation process. Whether or not it requires three, four (or more) kinetic steps is not clear from just the raw data.
Next, the authors observed the Rh XANES K-edge as it evolved over time (Fig. 5a–c) and separated it for the first three stages observed in Fig. 4. In Fig. 5a, they observed the signal of the reaction solution changes from being consistent with the RhCl3 reference to being consistent with the RhBr3 reference. The change suggests that in Stage I the Rh precursor undergoes a ligand exchange—note that determining if the balanced reaction stoichiometry involves ligands, as it seemingly must, would be useful here. In Fig. 5b, the signal begins to show features consistent with the Rh foil reference, and in Fig. 5c, significant changes are observed to suggest the formation of Rh nanoparticles.
Fig. 5 Rh K-edge XANES evolution for the RhCl3–TTAB–PVP–EG system. (a) Stage I: chloro ligands substituted by bromo ligands, time 0 min to 3 min, with RhCl3 and RhBr3 reference spectra given. (b) Stage II: formation of Rh2–4 nuclei, time 3 min to 15 min, with RhBr3 and Rh foil reference spectra given. (c) Stage III: growth of nuclei to nanoparticles of size 5.2 ± 0.8, time 15 min to 30 min, with RhBr3 and Rh foil reference spectra given. Arrows present in (a)–(c) indicate the direction of spectra change with time. Figure adapted with permission from ref. 40. Copyright 2012 American Chemical Society. |
The EXAFS and XANES data, from Fig. 4 and 5 respectively, are interpreted by the authors as suggesting a four-stage process for formation of their Rh nanoparticles. To investigate their hypothesis, the authors sought to fit their kinetics data, the loss of [Rh3+] as a function of time (Fig. 6). They used a known literature minimum mechanism, the 2-step mechanism of slow, continuous nucleation followed by autocatalytic surface growth,54 to fit their kinetics data (green triangles in Fig. 6).
Fig. 6 Loss of [Rh3+] as a function of time is plotted for RhCl3–PVP–EG with the bromide TTAB (green triangles) and without the bromide source TTAB (red squares). The system without TTAB (red) was fit using a pseudo-first-order rate law. The system with TTAB (green) was fit using the 2-step mechanism of slow continuous nucleation followed by autocatalytic surface growth (A → B with rate constant k1 and A + B → 2B, rate constant k2, where A = Rh3+ and B = Rh0). The fit yielded the following results: k1 = 0.005 min−1 and k2 = 8.77 min−1 M−1. Figure reproduced with permission from ref. 40. Copyright 2012 American Chemical Society. |
The authors demonstrated in Fig. 6 that bromide is critical kinetically to their reaction. The reaction without TTAB (red squares) was fit with an empirical power law79—and could not be fit (not shown) with a 2-step model.54 However, the reaction with TTAB (green triangles) was fit at least roughly by the 2-step mechanistic model of continuous nucleation, A → B with rate constant k1, and autocatalytic growth, A + B → 2B with rate constant k2, where A = Rh3+ and B = Rh(0).54 The data in Fig. 4–6, combined with mass spectrometry data in Fig. S1–S3 (ESI†), led the authors to proposal the pictorial four-stage model for Rh nanocrystal formation shown in Fig. 7.
Fig. 7 Schematic representation of the proposed four-stage formation process of Rh nanocrystals from RhCl3 deduced from XAFS and MALDI-TOF experiments. Stage I is ligand substitution from chloro to bromo ligands. Stage II is the formation of Rh2–4 nuclei. Stage III is proposed to be continued nucleation and growth of nuclei into nanoparticles and cubic nanocrystals. Finally, Stage IV is a final growth stage where imperfect nanoparticles and nanocrystals transform into more crystalline nanocrystals. The authors note, “the Rh species shown in the five boxes represent the predominant but not the exclusive Rh species present, and the chemical transformations illustrated under each arrow represent the major characteristic, but not necessarily the only process involved during each stage.”40 Figure reproduced with permission from ref. 40. Copyright 2012 American Chemical Society. |
Overall, the authors demonstrated the power of directly monitoring the loss of Rh3+ and the formation of Rh0 with DXAFS while also employing ex situ MS and TEM. Their 4-stage pictorial model—while not a chemical mechanism—does satisfy two of the criteria and partially satisfies the other two criteria underlying reliable mechanistic studies. For criterion (i), the authors monitored the transformation of Rh species throughout the reaction from Rh3+–Cl to Rh3+–Br and finally Rh0. They identified, through MS, low molecularity nucleation, where the kinetically effective nucleation (KEN)30 was observed to be Rh2–4—an important result in its own right.
However, the complete, balanced reaction stoichiometry was not reported and would at a minimum suggest the species that could be important in the rate law for particle formation. The true nature of the final product remains unspecified, nominally something like “RhnBra(PVP)b(glycolic acid)c(ethylene glycol)d”, where the unknown coefficients c and d might be small enough to be negligible. Additionally, the Stage I ligand substitution reaction merits additional study—as does the simplification of the system by starting with preformed RhBr3. The fully study of the exchange reaction and the proposed equilibrium constants are of interest as are controls of the effects of added H2O on this reaction, almost surely important as glycolaldehyde is expected to be a kinetically more facile oxidant than ethylene glycol when the H2O needed for the reaction converting glycolaldehyde plus H2O to glycolic acid (and 2e−/2H+) is present—additional reasons why getting the balanced reaction stoichiometry is needed as well as key for this system.
As for the second criterion (ii) of a reliable mechanistic study, the authors collected excellent, direct kinetics data (Fig. 4 and 5). The Rh species were monitored by in situ XANES and EXAFS throughout the entire reaction. (iii) The third criterion of pseudoelementary steps to describe the physical processes has been partially satisfied. Unfortunately, the authors then appear to skip the 2-step mechanism used in Fig. 6 to fit their kinetics data and jump to the proposed a 4-stage representation in Fig. 7 that is ultimately just a pictorial restatement of their hypothesis of a 4-stage process. Helping the community move beyond just such pictorial, qualitative representation of their otherwise powerful synchrotron-radiation-based studies is a primary goal of this review.
Finally, criterion (iv) has been accounted for as the authors disproved alternative hypotheses throughout their investigation. One example was the reaction without TTAB (red squares) in Fig. 6, where the control experiment produced different kinetics and further supported Stage I of the proposed 4-stage model. Additionally, the authors fit their kinetics data (in their ESI†40) with an additional literature mechanism, a 4-step model80 of continuous nucleation, autocatalytic surface growth, bimolecular agglomeration, and autocatalytic agglomerative growth. They found that agglomerative growth is not a significant contributor to the growth of their nanocubes, evidence supportive of their proposed growth pathway (Stage III) and shape correction (Stage IV). Presumably, the 2-step mechanism used to fit their data corresponds to Stages II and III, but without writing out the precise pseudoelementary steps, one does not know what steps are being fit in Fig. 6.
Noteworthy is that the authors are careful and do not claim that they know the mechanism. Rather, they have offered their evidence and conclusions while working towards a disproof-based mechanism, achieving a quite plausible model consistent with and supported by their data. Their study is a good example of acquiring high-quality, direct kinetics data from XAFS and corroborating it with other physical techniques.
Welcome on this valuable system would be additional XAFS as well as SAXS studies (i) starting with RhBr3 to avoid the ligand exchange, (ii) establishing the composition of the “RhnBra(PVP)b(glycolic acid)c(ethylene glycol)d” products, and (iii) establishing the balanced reaction including any role of water and the expected role of glycol as the reductant and the oxidized, glycolaldehyde/glycolic products. Writing plausible balanced reactions teaches that H+ are almost surely produced in the particle-formation process so that the effects of added, non-coordinating bases such as Proton Sponge™ may well significantly accelerate the reaction—and likely lead to better stabilized particles. Then (iv) measuring the full rate law for the nucleation and growth processes ([RhBr3], [Br−], [glycol], [Base], [PVP], Proton Sponge™ and any other dependencies will be needed, along critically with (v) expressing the results in pseudoelementary step mechanisms that can be subjected to further disproof until, ideally, a single proposed mechanism remains, are also needed. The results of those studies promises to be a reliable mechanism for Rh(0)n and by extension to other particle formations starting with the common precursor of the metal-halide and using alcohol as a common reductant in the presence of water.
The experimentally determined, balanced reaction stoichiometry of the HPO42−-stabilized Ir(0)∼150 nanoparticles has been published,32 is given in eqn (3) in its more general form, and makes apparent that the HPO42− nanoparticle-stabilizing ligand is one key in the reaction that, therefore, needs investigation in the kinetics studies:
(3) |
Scheme 2 Formation of a 1:2 HPO42−/IrI(1,5-COD)+ Intermediate, 2, and Then the final 1:1 diphosphate-bridged complex, 3, supported by the 1H NMR titration studies. Reproduced with permission from ref. 32. Copyright 2019 American Chemical Society. |
A convenient, albeit indirect, cyclohexene hydrogenation catalytic reporter reaction (CHCRR)54 was developed54 and used for in-house monitoring of the Ir(0)n particle formation. Another advantage of that method besides its convenience is that it generates thousands of precise, ±0.01 psig kinetics data points. (For a more extensive understanding of the CHCRR and its accompanying assumptions, please refer to ref. 32, 41, and Appendix A of ref. 54 that provides the original derivation for the use of the CHCRR.) In addition, in-house GLC of the cyclooctane product, 1H NMR, and UV-vis were used to monitor the reaction and check on the in-house CHCRR kinetics—the CHCRR proving to be the most indirect, but the quickest and overall quite useful, monitoring method to scout out optimized conditions for subsequent more direct methods of XAFS and SAXS monitoring of the particle formation reaction, vide infra.
The reproducible sigmoidal kinetics curves obtained prior to the XAFS and SAXS studies were fit by the literature 1997 2-step mechanism54 of slow, continuous nucleation (A → B, k1obs) and autocatalytic surface growth (A + B → 2B, k2obs). Reactions were run around a range of [Ir2]Initial from 1.5 mM to 6.0 mM and at five different amounts of HPO42− stabilizer, [HPO42−]Added. A total of 20 different combinations of [Ir2]Initial and [HPO42−]Added were studied and run in triplicate or more—multiple experiments any time they are desired being a significant advantage of in-house methods such as the CHCRR prior to more direct XAFS and SAXS studies. Additionally, the dependence of the nucleation rate constant, k1obs, on both the starting [Ir2]Initial and, importantly, on the amount of [HPO42−]Added (as suggested by the balanced reaction stoichiometry) were examined experimentally. The resultant data are reproduced below in Fig. 8a and b.32
Fig. 8 (a) Plot of k1obs (h−1) vs. the initial {[(1,5-COD)IrI·HPO4][Bu4N]}2 precursor concentration (M). Nucleation rate constants (k1obs) were collected from 0.0015–0.006 M at five ratios of [(Bu4N)2HPO4]:[Ir] (red squares = 1.8:1, orange circles = 2.7:1, green diamonds = 3.6:1, blue triangles = 4.5:1, and purples X's = 5.4:1). “The important result is the slope = 0 in each case within experimental error, indicating that there is no further dependence (i.e., beyond the observed, first-order dependence) of k1obs on the starting concentration of {[(1,5-COD)IrI·HPO4][Bu4N]}2.”32 (b) Plot of k1obs (h−1) vs. the added (Bu4N)2HPO4 beyond a 1:1 ratio of [HPO42−]:[Ir]. The starting iridium concentrations were 3.6 mM (yellow squares), 6.0 mM (blue circles), and 9.6 mM (purple diamonds). Data were fit based on a derivation for KDiss (see Scheme 3 below) that is given in ref. 32. Reproduced with permission from ref. 32. Copyright 2019 American Chemical Society. |
Fig. 8a demonstrates that, for this 2nd-generation iridium nanoparticle formation, nucleation is first-order in the well-characterized, dimeric precursor {[(1,5-COD)IrI·HPO4]2}2−, consistent with a small kinetically effective nucleus (KEN)30 of just two iridium atoms, Ir(0)2. Furthermore, the k1obs nucleation rate constants versus the added [(Bu4N)2(HPO4)] data led to the proposed nucleation mechanism32,41 given in Scheme 3 as the only mechanism of five nucleation mechanisms considered able to fit the data in Fig. 8b. The kinetics analysis also reveal that {[(1,5-COD)IrI·HPO4]2}2− is not the active species that nucleates to Ir(0)2, but instead suggests {[(COD)IrI(solv)]2·HPO4}0 as one top candidates for a kinetically competent intermediate involved in nucleation (eqn (5)) (and possibly also (1,5-COD)IrI(solvent)2+ based on prior precedent31 and studies in progress).Noteworthy here is that multiple attempts over many years to obtain XAFS and SAXS kinetics data for the first-generation {(1,5-COD)IrI·POM}8− system failed due to the large, W-containing POM9− stabilizer, the large ca. 1.5 nm by 1.2 nm “cigar-shaped” POM9− size obscuring SAXS monitoring of early time nucleation events and the W in the POM9− (= P2W15Nb3O629−) interfering with the Ir-XANES. Hence, the {[(1,5-COD)IrI·HPO4]2}2− second generation system was developed specifically with XAFS and SAXS monitoring in mind and in order to check on the convenient, precise-data-generating, but indirect, CHCRR kinetics method.
Scheme 3 Proposed nucleation mechanism32 involving the dissociative equilibrium from {[(1,5-COD)IrI·HPO4][Bu4N]2} in acetone, based on the [HPO42−] dependence of the kinetics and 1H NMR evidence for the neutral, {[(COD)IrI(solv)]2·HPO4}0 Intermediate. Adapted and Reproduced with Permission from ref. 41. Copyright 2021 American Chemical Society. |
Synchrotron XAFS and SAXS experiments were conducted on the second-generation {[(1,5-COD)IrI·HPO4]2}2− precursor to Ir(0)n·(HPO42−)x nanoparticle formation system at initial Ir concentrations from 3.0 mM to 12.0 mM and at five Ir:HPO42− stabilizer ratios of from 1:1.8 to 1:5.4.41 The Ir(0)n nanoparticle formation reactions were run under the standard conditions of ∼40 psig H2 at 22 °C and with 1 molar equivalent of Proton Sponge base per Ir present to absorb the H+ produced in the balanced reaction, eqn (3)—the Proton Sponge also prevents hydrogenation of the acetone solvent to 2-propanol by the highly catalytically active Ir(0)n nanoparticles,84 experimental design that takes advantage of the deep knowledge of the two Ir(0)n nanoparticle formation systems.30–32,54,67,68,80–82,84
Optimized conditions for the XANES ([Ir] = 0.005 M, 2.25 molar equiv. HPO42−) and SAXS ([Ir] = 0.009 M, 3.6 molar equiv. HPO42−) were worked out first with the in-house CHCRR. Those optimized conditions were then employed at the synchrotron along with Tandem CHCRR kinetics data collected simultaneously with (separate) XANES and SAXS studies done at different synchrotrons via separate collaborators.41 Also collected was additional in-house CHCRR data (i.e., under the XANES and SAXS conditions41) as was gas–liquid chromatography quantification of the cyclooctane (COA) product formation vs. time, again under the XANES and SAXS conditions. The composite data in terms of molar concentration of iridium (M) as a function of time (h) are given below in Fig. 9.
Fig. 9 (a) XANES data (hollow purple circles), tandem CHCRR data collected at the synchrotron (hollow black diamonds), in-house CHCRR data collected under otherwise identical conditions (open green squares), and GLC of the cyclooctane product (solid yellow triangles) were fit with the 2-step mechanism54 to yield a global average fit (dashed black line). Resultant rate constants are k1obs(avg,XANES) = (2.2 ± 0.3) × 10−2 h−1 and k2obs(avg,XANES) = (3.7 ± 0.1) × 102 h−1 M−1. (b) SAXS data (hollow blue circles), tandem CHCRR data collected at the synchrotron (hollow black diamonds), in-house CHCRR data collected under otherwise identical conditions (open green squares), and GLC of the cyclooctane produced (solid yellow triangles) were fit with the 2-step mechanism54 to yield a global average fit (dashed black line). Resultant rate constants are k1obs(avg,SAXS) = (1.7 ± 0.2) × 10−2 h−1 and k2obs(avg,SAXS) = (2.0 ± 0.1) × 102 h−1 M−1. Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
The k1obs nucleation rate constants for the XANES, in-house CHCRR, and GLCcyclooctane datasets are equivalent within 1.6-fold. The Tandem CHCRR dataset to the XANES was determined to be an outlier statistically (>5σ) and experimentally (due, apparently, to small amounts of Ir(0) from X-ray radiolysis catalyzing and accelerating the CHCRR).41 The k1obs nucleation rate constants for the SAXS, Tandem CHCRR, in-house CHCRR, and GLCcyclooctane datasets all proved to be within a similar factor of 2. The Tandem CHCRR run simultaneously to the SAXS undercuts the synchrotron data some, but there was evidence of Ir(0) metal fouling found on the SAXS cell window,41 presumably due to actually well-precedented—but little discussed and certainly under recognized—X-ray radiolysis during such synchrotron studies.85–93
There is compelling precedent in the well-studied Ir(0)n nanoparticle system that even trace amounts of adventitious Ir(0) hugely increase the nucleation and subsequent growth rates, sometimes resulting in complete elimination of the induction period of the typically sigmoidal kinetics curve.54,82 A fundamental contribution from as well as these particular synchrotron XANES and SAXS studies is, then, to raise a red flag concerning the involvement in X-ray-based methods that are more often than not assumed to be completely direct without artifacts. Relevant here is that Frenkel and collaborators have shown85 that a photon flux of ∼3 × 1013 photons per second results in an estimated ∼1016 solvated electrons94 even when that radiolysis involves the otherwise improbable, highly energetic removal of an electron from divalent zinc, Zn2+ + hν → Zn3+ + e−. Highly relevant here is that in particle formation systems exhibiting autocatalytic, exponential growth as is commonly seen, any trace nucleation events due to X-ray radiation-induced radiolysis will then be quickly magnified, autocatalytically and, hence, exponentially.
The data from all the monitoring methods (XANES, SAXS, CHCRR, and GLC) were standardized, compiled, and analyzed as a single dataset containing 1178 total data points. These data and the global fit using the 2-step mechanism are reproduced as Fig. 10. The resultant, globally average, fit-determined rate constants yielded the relatively precise rate constants k1obs,globalavg = (1.5 ± 0.1) × 10−2 and k2obs,globalavg = (2.4 ± 0.1) × 102 (±7% and ±4% error, respectively) that are believed to be reasonably accurate41k1obs and k2obs rate constants as well.
Fig. 10 “The seven datasets are given: XANES as hollow purple circles, in-house CHCRR (under XANES conditions) as hollow green squares, COA (under XANES conditions) as hollow yellow triangles, SAXS as solid blue circles, in-house CHCRR (under SAXS conditions) as solid green squares, COA (under SAXS conditions) as solid yellow triangles, and Tandem CHCRR (simultaneous with SAXS measurements) as solid black diamonds. The solid red line represents the global fit to all seven datasets.”41 The resultant, globally average, fit-determined rate constants with fitting error are k1obs,globalavg = (1.5 ± 0.1) × 10−2 and k2obs,globalavg = (2.4 ± 0.1) × 102. All datasets, as mentioned in ref. 41, were standardized to 7.5 mM in IrI. Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
Overall, the XANES, SAXS, in-house CHCRR, and GLC methods all reported equivalent rate constants within either 1.6 or 2.0 orders of magnitude for XANES and SAXS, respectively, a rare comparison of multiple in-house as well as synchrotron XANES plus SAXS particle formation monitoring methods. The composite kinetics data from all available experimental methods (XANES, SAXS, CHCRR, GLC, and 1H NMR), were fit by the 2-step mechanism of slow, continuous nucleation (A → B, k1obs) and autocatalytic surface growth (A + B → 2B, k2obs). Therefore, all experimental methods, including the synchrotron XANES and SAXS, presently support the 2-step mechanism as the minimum mechanism to describe quantitatively Ir(0)n formation.
What is interesting from the well-studied Ir(0)n nanoparticle systems is that, in both the first and second-generation systems now, even with the addition of synchrotron XANES and SAXS and XAFS, the sum of all of the kinetics methods are insufficient to determine the true particle formation mechanism as demonstrated by ME-PBM analysis, discussed next. In addition, even with all the effort that went into this32,41 Ir(0)n system, the exact role of the HPO42− ligand has been elucidated on only the nucleation step, so only partially—so that full incorporation of specific ligand effects into the pseudoelementary step mechanism is of interest and remains to be accomplished. Ligand effects of the POM9− polyoxometalate on the nucleation31 and growth95 steps of the first generation {(1,5-COD)IrI·POM}8− precursor to Ir(0)n·(POM9−)m nanoparticle system are available, however, to help guide the needed studies.
ME-PBM is also able to inform the inverse problem of going from “observations/effects” back to “cause/a mechanism” whereby one tests and thereby refutes—or supports—that minimum mechanism via ME-PBM analysis of the PSD data.67,68 As alluded to above, in both the first67,68 and now second generation Ir(0)n nanoparticle systems, the ME-PBM analysis of the wealth of kinetics data otherwise buried in the PSD has led to disproof of the minimum mechanisms that were previously able to account for all of the experimental data (i.e., other than the PSD). Significantly, ME-PBM analysis of the PSD in the first generation Ir(0)n system also provided a new paradigm for how narrow PSDs can be formed:67,68 smaller particles grow faster than larger particles, thereby catching up in size to them and resulting in near-monodisperse PSDs despite the inherently broadening effects of continuous nucleation.67,68
Hence, initial ME-PBM analysis of the TEM-determined PSDs in the second-generation Ir(0)n nanoparticle formation system was reported in the 2021 paper that reports the XANES and SAXS studies.41 Specifically, ME-PBM was used to analyze the end-time PSD at 10.0 h (shown as the black distribution in Fig. 11), using to start the 2-step mechanism while also including the experimentally determined prior equilibrium nucleation mechanism (given back as eqn (4) and (5)).
Fig. 11 “Particle size distributions for the formation of Ir(0)n nanoparticles at 1.5 h (green), 3.25 h (teal), 5.0 h (purple), and 10.0 h (black) binned in 0.2 nm width bins. At each time point, a new reaction solution was prepared of 5.0 mM [(1,5-COD)IrI(NCCH3)2][BF4] in the presence of 2.25 molar equiv. of (Bu4N)2HPO4 in 3.33 mL acetone and 0.67 mL cyclohexene at 22.0 ± 0.1 °C. Each data point represents >450 measured particles; across the 4 samples >2700 particles were measured.”41 Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
The attempted fit using the 2-step mechanism consistent with the XANES, SAXS and all the other data to this point, along with the simulated [A]t loss curve from the resulting rate constants of that attempted fit, are reproduced as Fig. 12a and b, respectively. The fit is obvious poor—no fit at all, really—disproving the 2-step mechanism that was otherwise consistent with all of the XAFS, SAXS, and other kinetics data.41
Fig. 12 (a) Attempted fit to the end-time, 10.0 h histogram using ME-PBM built off of the 2-step mechanism with experimentally determined nucleation mechanism. The resultant fit-determined rate constants are as follows: k+Diss = 4.0 × 10−1 h−1 M−2, k−Diss = 3.7 × 104 h−1 M−1, k1 = 6.6 × 10−1 h−1, k2 = 9.2 × 103 h−1 M−1, unreliable rate constants given the poor fit. The Best Function Value (BFV)67,68 is 90.0, indicating a poor fit (as lower BFVs indicate67,68 better agreement between the experimental data and the attempted fit). However, the rate constants are provided here because they are used to generate the predicted precursor loss curve discussed next. (b) Calculated precursor loss, [A]t, using the “best-fit” rate constants parameters from the above attempted fit to the PSD, are co-plotted with the experimental global [A]t data from the simulated global fit of data from all the available methods. Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
Given the poor fit in Fig. 12a, the authors considered—that is, attempted to disprove—15 total minimal mechanisms, even if those are but a fraction of the 96 possible mechanisms summarized via Table S5 (ESI†) of the present review. The minimal mechanism the authors report as a current best-fit of the PSD is the 4-step mechanism given in Scheme 4 that contains slow, continuous nucleation (Asolv → B, k1obs), small particle growth (A + B → C, k2obs), bimolecular small particle agglomeration (B + B → C, k3obs), and large particle growth (A + C → 1.5C, k4obs), where A = {[(1,5-COD)IrI·HPO4]2}2−, Asolv = {[(COD)IrI(solv)]2·HPO4}0, B = small Ir(0)m, and C = large Ir(0)n.
Scheme 4 Experimentally determined prior equilibrium and the proposed 4-step mechanism of (i) slow, continuous nucleation, (ii) small particle growth, (iii) bimolecular small particle agglomeration, and (iv) large particle growth. Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
Although the authors demonstrate that the 4-step mechanism given in Scheme 4 can fit the PSD reasonably well, Fig. 13a, the authors note that the fit to the experimental [A]t loss curve in Fig. 12b “is poor” when compared to the simulated [A]t loss curve (i.e., the expected [A]t loss curve simulated using the rate constant parameters extracted from the PDS). This example illustrates both the mechanistic power, as well as the amount of work needed, when using ME-PBM as a now apparently “Gold-Standard Test” of ones proposed particle formation mechanism.41,67,68 As the authors note, they “still have more work to do to uncover the more detailed, even closer to correct particle formation mechanism.”41 Preliminary studies point towards the possibility that Ir(1,5-COD)(solv)2+ may function as a precedented,31 kinetically competent intermediate in the true nucleation mechanism and its associated elementary steps.
Fig. 13 (a) Fit to the end-time, 10.0 h histogram using ME-PBM built off of the new 4-step mechanism with experimentally determined nucleation mechanism. The resultant fit-determined rate constants are as follows: k+Diss = 2.6 × 10−1 h−1 M−2, k−Diss = 2.2 × 104 h−1 M−1, k1 = 2.2 h−1, k2 = 5.4 × 104 h−1 M−1, k3 = 1.6 × 106 h−1 M−1, k4 = 1.0 × 103 h−1 M−1, and the B versus C cutoff value M = 23. The Best Function Value (BFV)67,68 is 25.4, indicating an improved fit over the 2-step shown in Fig. 12—as, again, a lower BFV indicate67,68 better agreement between the experimental data and the attempted fit. (b) Calculated precursor loss, [A]t, using the “best-fit” rate constants parameters from the above attempted fit to the PSD, are co-plotted with the experimental global [A]t data from the simulated global fit of data from all the available methods. Reproduced with permission from ref. 41. Copyright 2021 American Chemical Society. |
Evaluating the second-generation {[(1,5-COD)IrI·HPO4]2}2− precursor to Ir(0)n·(HPO42−)x nanoparticle system32,41 for how well it stacks up to the four criteria for establishing a minimum mechanism, the study provides: (i) an experimentally determined, balanced reaction stoichiometry that predicted the subsequently observed importance of the HPO42− nanoparticle-stabilizing ligand in the particle-formation kinetics; (ii) extensive kinetics data by 4 methods including XANES and SAXS; (iii) (pseudo)-elementary reaction steps for the proposed mechanism that add up to the experimental stoichiometry, define the rate constants, and provide defined, unambiguous words for describing the mechanism. Additionally, the authors provide (iv) disproof of 14 alternative (out of 15 total), possible hypothesized mechanisms.
Somewhat hidden behind the success of the Ir(0)n systems is the amount of synthesis and characterization of the first generation {(1,5-COD)IrI·POM}8− precursor81,82 and then the second-generation {[(1,5-COD)IrI·HPO4]2}2− precursor32,96,97 characterization of their nanoparticle products by as many applicable methods as possible,81,82,97 and the care in establishing the balanced reaction stoichiometry81,82,97 before beginning serious kinetics and mechanistic studies. Strict adherence to a disproof- and Ockham's razor-based approach, as rigorous mechanism demands, is another key underlying the Ir(0)n kinetics and mechanistic studies—as is a philosophical acceptance of the truism that all deliberately minimum mechanisms will eventually be upgraded (“disproved”) as new methods, experiments, or data analyses (such as ME-PBM) become available. At present, a 1-step upgrade, hence 3-step total, PEStep mechanism is what is really a minimum mechanism for particle formation67,68 and despite the (still recommended67,68) use of just the 2-step mechanism for the initial analysis of kinetics data before moving on to the examination and fitting of PSD data (the 2-step rate constants being, for example, useful starting guesses for 3-step mechanism analyses and fits of PSDs).67,68
Needed additional studies for even this relatively well studied, second-generation Ir(O)n nanoparticle formation system include: (i) resolution of the precise nucleation mechanism with the help of the additional ME-PBM; (ii) a better accounting of the reporter reaction kinetics in Fig. 12b (i.e., by an updated nucleation mechanism or some other, needed update of the mechanism); and (iii) elucidation of how ksurface growth varies with nanoparticle size in more detail−that is, elucidation of the more detailed, “growth kernel”, including the likely role of ligand capping effects in that growth kernel, better than the k2 > k4 step-function kernel66,67 used presently (Scheme 4 and Fig. 13).
(6) |
First, the formation of nanospheres from the reactions with 50 mM and 100 mM OY (i.e., and not the nanowires formed under different conditions, vide infra) were investigated by XAFS and SAXS. The XAFS data were interpreted using linear combination analysis (LCA) that allowed for simultaneous determination of precursor concentration (AuIII), intermediate concentration (AuI), and final product concentration (Au0) as a function of time. The XAFS formation data as a function of time and an example of the LCA result are given as Fig. 14a and b, respectively.
Fig. 14 (a) In situ XAS spectra at the Au LIII-edge monitoring the formation of Au nanospheres from the gold precursor solution with [OY] = 50 mM. Spectra are given stepwise from 0–88.6 s and 30 min. The decrease in the intensity of the white line intensity at 11921.2 eV is denoted by “A”, and the increase in the signal corresponding to Au foil at 11948.5 eV is denoted by “B”. (b) The XAFS spectrum collected at 29.2 s, hollow circles, is given as a representative spectrum for linear combination analysis (LCA). The pink line is the LCA fit of the XAFS spectrum, which yields the three components: Au(III) precursor (black line), Au(I) intermediate species (red line), and resultant Au(0) nanospheres (blue line). Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
As CNT did not accurately describe the particle-formation data, a new kinetic model had to be devised in order to fit the data. The authors constructed a multistep kinetic model to account for the induction period, nucleation, and for growth. The proposed nanosphere formation kinetic model is given as Scheme 5 and includes equilibrium expressions between reactive and non-reactive Au(III) and Au(I) species.
Scheme 5 Kinetic model for the formation of Au nanospheres with [OY] = 50 mM. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
The proposed kinetic model contains two reduction steps: Au(III)reac to Au(I)reac with rate constant k1 and Au(I)reac to Au(0) with rate constant kN. That second step with kN is the first step in the two-step nucleation/growth process. The final, autocatalytic process54 involves Au(0) and Au(I)reac to form additional Au(0) with rate constant kAC. The authors report that “a competition between reduction and complexation of the Au(III) and Au(I) by OY is described by two additional equilibria during the induction stage. The two ‘non-reactive’ complexes are not directly involved in the reduction steps but serve as a reservoir for the ‘reactive complexes’.”42
Using the mass-balance expressions of [Au(III)Total] = [Au(III)R] + [Au(III)NR], [Au(I)Total] = [Au(I)R] + [Au(I)NR], and [Au(0)] = 1 − [Au(III)Total] − [Au(I)Total], the authors input their model into MatLab as a series of ordinary differential equations (ODEs), eqn (7)–(10).
(7) |
(8) |
(9) |
(10) |
Fig. 15 The Au(III), Au(I), and Au(0) concentrations determined by LCA analysis of the XAFS spectra (example given in Fig. 14b) are given here from 0–120 s as filled black squares, filled red circles, and filled blue triangles, respectively. The best-fit lines from each Au species are given using the corresponding connected hollow shape. The black (+) represents the reactive Au(III) species, while the black (×) represent the non-reactive Au(III) species. The red (+) represents the reactive Au(I) species, while the red (×) represent the non-reactive Au(I) species. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
The data are well-fit from 0–20 s and then again from 45–120 s. For all species of Au, there is a discrepancy between the data and the fit lines between 20–45 s, somewhat surprising given the considerable flexibility of 7 variable rate-constant parameters, an indication that one or more of the assumed steps in the proposed model is incorrect (e.g., the assumed first-order nucleation?) or that needed additional steps are missing.
Next, the formation of nanowires from the reaction with 400 mM OY were also investigated by XAFS and SAXS. The XAFS data were interpreted using LCA to simultaneous determine precursor concentration (AuIII), intermediate concentration (AuI), and final product concentration (Au0) as a function of time. The XAFS formation data as a function of time is given as Fig. 16.
Fig. 16 (a) In situ XAS spectra at the Au LIII-edge monitoring the formation of Au nanowires from the gold precursor solution with [OY] = 400 mM. Spectra are given from 0 min – 4 days. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
Using the kinetics model presented in Scheme 6 and the mass-balance expressions from before, the authors derived ODEs, reproduced below as eqn (11)–(13), and associate 6 rate-constant parameters.
(11) |
(12) |
(13) |
Scheme 6 Kinetics model for the formation of Au nanowires with [OY] = 400 mM. Reproduced with permission from ref. 42. Copyright Royal Society of Chemistry 2020. |
Fig. 17 The Au(III), Au(I), and Au(0) concentrations determined by LCA analysis of the XAFS spectra are given here from 0–650 min as filled black squares, filled red circles, and filled blue triangles, respectively. The best-fit lines from each Au species are given using the corresponding hollow shape. The black (+) represents the reactive Au(I) species, while the black (×) represent the non-reactive Au(I) species. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
Next, and returning to the case with [OY] = 50 mM that formed nanospheres, the in situ SAXS patterns were closely analyzed during the first 38 s when the authors hypothesize pre-nucleation clusters (PNCs) are formed. Prior to nucleation at 38 s, it was observed that the scattering pattern changed. The scattering pattern shifts again following the onset of nucleation during the primarily growth period from 48 s–30 min. These scattering patterns have been reproduced herein as Fig. 18.
Fig. 18 In situ SAXS patterns for the Au precursor system with [OY] = 50 mM from (a) 4–38 s, where it is believed that pre-nucleation clusters are formed and nucleation begins. Next, from (b) 48 s–30 min, primarily growth is believed to occur. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
These scattering patterns were fit, and the particle diameter versus time plotted. These data are presented in Fig. 19 as diameter in nanometers versus time in seconds, where time is plotted on a log scale.
Fig. 19 Mean particle diameter (nm) is given as a function of log time (s) for the Au precursor system with [OY] = 50 mM, as shown in Fig. 18. In the induction period, the mean diameter was calculated “from the radius of gyration of the Au(III)–Au(I) pre-nucleation clusters”.42 In the growth period, the mean diameter was determined from the Au(0) nanospheres. The shaded region is believed to be the time period when primarily nucleation takes place. Reproduced with permission from ref. 42. Copyright 2020 Royal Society of Chemistry. |
Based on the data presented in Fig. 19, the authors believe during the induction period that the SAXS scattering is reporting “the radius of gyration of the Au(III)–Au(I) pre-nucleation clusters”42 and not of Au(0) nanospheres. Then, following what the authors have labeled as the nucleation period, the implication is that the SAXS is reporting the mean radius of the resultant Au(0) nanospheres.42 Hence, it follows that the true nucleation kinetics are obscured by the much larger PNC, so that only growth data are what is primarily observed in the SAXS measurements. The authors did not attempt to fit the SAXS data with their kinetic model.
Overall, Ramamoorthy and co-workers have reported a very valuable study on the formation of Au(0)n in which they directly monitor the formation of Au(0)n nanospheres or nanowires by in situ XAFS, SAXS, and HE-XRD, and also characterize the resultant particles with ex situ TEM.42 Of perhaps special interest is the observation of Pre Nucleation Clusters (PNCs) and even though those PNCs have not yet been shown by the needed kinetics studies to be on (vs. off) the kinetically dominant pathway to the nanoparticle product. In terms of satisfying the four criteria for a reliable mechanistic investigation, for criterion (i) the authors have proposed a balanced reaction stoichiometry that was reproduced herein as eqn (6). The exact composition of the pre-nucleation clusters would be a useful target for future studies. The authors discussed possible compositions of the Au-based PNCs, but the PNCs were not determined beyond aggregates of {(RNH3+)(AuCl4−)(OY)}n, where R = isopropyl.42 However, after nucleation has been initiated and growth has begun, the authors have been able to determine the Au speciation—a non-trivial task—and plot several Au species as a function of time, another non-trivial achievement. Criterion (ii) was also fully satisfied as the kinetics of Au(0)n were collected in situ by both EXAFS and SAXS over the entire reaction as Fig. 14–18 document for both the nanospheres and the nanowires.
Criterion (iii), namely the task of constructing a pseudoelementary step-based mechanism, has been partially satisfied. For the two cases—nanospheres and nanowires—the authors have presented proposed kinetic models reproduced herein as Schemes 5 and 6, respectively. They have used these kinetic models to write differential equations and then fit their kinetics data with those differential equations. In a future study, once the composition of the PNCs have been determined, then one will be able to write the exact pseudoelementary step reactions that sum to the overall net Au(0)n formation mechanism. Missing at present is an experimentally determined rate law for nucleation. As for criterion (iv), one alternative model was examined. The authors demonstrated that CNT was invalid for their Au(0)n formation system, but they did not report the disproof of any other alternative mechanisms en route to their proposed kinetic models. A greater amount of disproof will be required to provide higher confidence in the resultant, proposed mechanism.
In short, the valuable Ramamoorthy and co-workers Au(0)n system is an illustrative case history demonstrating state-of-the-art, in situ use of three separate synchrotron techniques: XAFS, SAXS, and HE-XRD, a study that also identifies PNCs. Additional studies of this interesting system the community would likely welcome include those hinted at above: (i) determination of the composition(s) of the pre-nucleation clusters; (ii) further work on the experimental rate law especially of nucleation and showing whether or not the observed PNCs are on the kinetically dominant nucleation pathway; and (iii) testing of additional reasonable alternative mechanistic models for formation of the gold nanosphere and nanowires. We look forward to seeing those important, likely exciting studies.
The two case studies are: (Section 4.1) cadmium-selenide nanocrystal formation monitored by XAFS43 and (Section 4.2) zinc-oxide nanoparticle formation monitored by SAXS/WAXS/UV-visible.44 Like before in Section 3, each case study is organized by: (i) a summary of the system and techniques used to study it; (ii) a review of the key kinetics data; (iii) the authors’ proposed formation model or mechanism; and (iv) an analysis of the case's results as compared to the four components required for establishing a minimum, disproof-based mechanism. Following these two main case studies of this section are two additional shorter sections: (Section 4.3) two studies on shape control (tungstite nanoplatelets and copper nanocrystals), and a listing of (Section 4.4) Additional systems of interest (perovskites, quaternary nanocrystals, amorphous nanoparticles, and carbon quantum dots). These latter sections serve to highlight newer systems of interest to the broader nano community even if of those systems are less well investigated at present— our intent being to be supportive and encouraging of those intriguing studies as they strive to achieve the next level of mechanistic analysis.
Fig. 20 “Schematic principle of the in situ XAFS experimental set-up using a microfluidic cell. Fluorescence yield spectra record at positions (x,y) along a microreactor channel provide time-dependent XAFS spectra. In the present reaction, precursor flows (TOP-Se and Cd stock solution) are mixed before introducing to a heated microchannel.”43 Reproduced with permission from ref. 43. Copyright 2011 International Union of Crystallography. |
The proposed reaction stoichiometry, that is, as far as we can determine it from the data provided in ref. 43, is given in eqn (14).
(14) |
The authors analyzed their XANES, EXAFS, UV-vis and photoluminescence (PL) results to determine the early time (0.0–8.1 s) (CdSe)n formation behavior. The XANES results given in Fig. 21a indicate the author's direct observation of (CdSe)n particle formation. Simultaneously, the authors used UV-vis and PL spectroscopies (Fig. 21b and c) to monitor directly both (CdSe)n particle formation and crystallinity.
Fig. 21 (a) Experimental XANES spectra collected from 0.0–8.1 s after “initiation of the reaction with the addition of 5 wt% DDA,” where again DDA is dodecylamine.43 Reference spectra of (CdSe)n were given for (CdSe)n powder crystals of 3 nm. “Open circles indicate simulated XANES spectra.”43 (b) Photoluminescence spectra of (CdSe)n nanoparticles in the presence of 5 wt% DDA.43 (c) UV-vis absorbance spectra of the (CdSe)n nanoparticles for the same solution as (b). Reproduced with permission from ref. 43. Copyright 2011 International Union of Crystallography. |
By XANES, one can see the increase at ca. 12670 eV with time that is indicative of (CdSe)n particle formation, by comparison to authentic (CdSe)n powder. (CdSe)n particle formation is further monitored by the UV-vis absorbance at 350 nm, due to the well-studied cross section99, and the full-wide half-maximum of the PL peaks.43 Notably, prior to 2.4 s, there are no absorption or emission peaks present, indicating a lack of crystalline (CdSe)n. A closer observation of the PL spectra shows a red-shift with time, interpreted by the authors as the loss of deep-trap states due to increased crystallinity.43 The PL spectra also indicate that the particle shape is growing more symmetrical, the narrowing and increased symmetry of the peak supporting the interpretation of increased crystallinity.100 Hence, at the earlier times of 2.4 and 3 s, amorphous (CdSe)n nanoparticles are most likely formed, but by 8.1 s they have increased crystallinity.
The synchrotron EXAFS and XANES results were compiled in Fig. 22, where the Se–Cd coordination number was independently determined three ways as a function of time. The two data sets collected by EXAFS are the same within error, which is consistent with “a two-component model” (i.e., “unreacted Se–TOP does not dissolve leaving Se2− ions”).43 The authors remark that the XANES results are lower in average Se–Cd coordination number due to the reference used in processing the data.43 Furthermore, the authors considered three possible alternative hypotheses for the difference between the EXAFS and XANES data en route to their final Se–Cd coordination numbers.43 The (CdSe)n study43 presented by Maeda and co-workers demonstrates an important approach to observing particle formation at early times. They have designed a direct, tandem UV-vis and XAFS microfluidic technique and used it to monitor the formation of (CdSe)n during the first 10 s of the reaction.
Fig. 22 “Se–Cd coordination number as a function of elapsed time determined by EXAFS independently for the Se–P and Se–Cd pairs (open squares and open circles, respectively) compared with the value determined by the XANES fit (closed circles). The dashed line indicates the Se–Cd coordination limit. The reacted amount of Se and Cd during the reaction is assumed to be equal.”43 Reproduced with permission from ref. 43. Copyright 2011 International Union of Crystallography. |
In future studies the exact reaction stoichiometry would be desirable, thereby satisfying criterion (i) of a reliable mechanistic study. In that endeavor, the considerable literature in the semiconductor/quantum dot field, and the advances the field has made in the last decade, are likely to prove useful.25,101–104 For example, the authors claim43 that ODE was used because it is a non-coordinating solvent. However, a recent 2019 paper has demonstrated that ODE is not innocuous at temperatures greater than 393 K, where ODE will autopolymerize and coat the nanoparticle in a manner that cannot be removed with standard purification methods.105 Hence, at the temperature of the reaction43 (513 K), ODE polymers at least are likely to be coordinating.
Criterion (ii)—namely collecting kinetics data—was well satisfied by the use of XANES, EXAFS, UV-vis, and PL to address the kinetics of (CdSe)n nanoparticle formation. The use of a microfluidic device was effective in providing early time points. That said, future studies are needed to obtain the full rate law for particle formation. Criterion (iii) will also need to be addressed in future studies, that is, the need to write out the complete pseudoelementary steps to describe the (CdSe)n particle-formation mechanism. Finally, regarding criterion (iv), the authors did apply a disproof-based method to a number of their experiments. Importantly, they were careful to not overinterpreted their data, nor have they claimed to know more about the precise (CdSe)n formation mechanism than their data allows.
Overall, the authors provide a noteworthy initial study43 using in situ XAFS to produce Se–Cd bond formation kinetics data and garner insights into the early time of (CdSe)n nanoparticle formation. Future studies of (a) determination of a full rate law for particle formation, including nucleation, and (b) pseudoelementary step construction consistent with that experimental rate law would be most welcome on this valuable, interesting system building off the reported, valuable initial kinetics studies.
(15) |
(16) |
Fig. 23 As noted by the authors,44 “Schematic presentation of the in situ setup employed for real-time SAXS/WAXS/UV-vis measurements during the formation of ZnO nanoparticles. The setup measures SAXS, WAXS, and the UV-vis spectra simultaneously in the same sample volume.”44 Reproduced with permission from ref. 44. Copyright 2019 American Chemical Society. |
Fig. 24 “Time evolution of the average ZnO nanoparticle radii determined from the fitted SAXS curves for two different concentrations and temperature during the first 1000 s of the nucleation and growth reaction. The solid lines are fits to a homogeneous nucleation and growth model.”44 Reproduced with permission from ref. 44. Copyright 2019 American Chemical Society. |
Based on the in situ SAXS data reported in Fig. 24 and eqn (16), a series of differential equations, eqn (17)–(21), were derived to describe the ZnO particle formation process, only eqn (17) of which looks like it is somewhat based on a pseudoelementary step, that is, on eqn (16), albeit with a second-order dependence on [OH−] being implied if eqn (16) were truly elementary.
(17) |
(18) |
(19) |
(20) |
(21) |
The authors report that solving the set of “coupled differential equations yields the concentration [PJ](t) and radii RJ(t) of all particles J that have been formed during the nucleation process.”44 They solved the set of differential equations for each data set in Fig. 24, which resulted in the fits given and values for the precursor rate constant (k1), the growth rate constant (kgr), the interfacial tension (γ), and the saturation concentration ([ZnO]sat). However, even with ≥4 adjustable parameters, the fits appear to match only ca. 50% of the experimental data. This suggests that something important is missing from the model, a point which will be discussed more in a moment.
Of interest is that the authors observe a change in the WAXS signal and UV-vis peaks after 1 h, where higher-order Bragg reflections are identified and a red-shift of the absorption edge is observed (Fig. 25).44 The overall particle-size by SAXS does not change significantly after 1 h, which suggests that a different type of growth process is occurring. The authors conclude that an increase in ZnO crystallinity is taking place during this “second growth phase”, where the ZnO is observed to “develop full crystallinity.”44 Overall, the authors present expertly collected, high-quality kinetics data on the formation of ZnO. They collected it by multiple direct techniques using three direct physical methods, SAXS/WAXS and UV-vis, reasons among others that this study was selected as an illustrative case study.
Fig. 25 “UV-vis spectra as a function of time recorded in situ during the formation of the ZnO nanoparticles at a concentration [Zn2+]0 = 0.133 mol L−1 at temperatures of 40 °C (A) and 50 °C (B).”44 The red-shift of the absorption edge for both spectra is observed after 1 h, after which the primary nucleation and growth processes are hypothesized to have concluded as judged by the authors’ interpretation of their SAXS data. Reproduced with permission from ref. 44. Copyright 2019 American Chemical Society. |
The authors’ kinetics model partially satisfies two of the four minimum criteria for reliable mechanistic study. Unfortunately, the first criterion was not met in that a complete, balanced reaction stoichiometry was not provided.44 We presented a possible mass- and charge-balanced formation reaction back in eqn (15), but further characterization studies are warranted to elucidate the exact values for x and n in eqn (15). Note that even the stoichiometry written eqn (15), suggests a possible higher, perhaps second order in [OH−], rather than the assumed first-order dependence, [OH−]1, used in eqn (17). Next, criterion (ii) was satisfied: the authors collected high-quality size versus time SAXS data and absorbance versus time UV-vis data for the entire ZnO formation. The third criterion was partially satisfied as differential equations for the different processes were derived but, unfortunately, pseudoelementary steps for the reaction were not provided as required for a clearer chemical mechanism. Classical Nucleation Theory (CNT) was assumed over the needed experimental rate law. As a result, the d[ZnON]/dt term is equal to a compilation of CNT-based constants and not the desired concentration terms (i.e., functions of [Zn2+], [OH−], and [C18H33O2−]) and associated rate constants corresponding to a mechanism composed of (pseudo)-elementary step reactions.106 Lastly, criterion (iv) remains to be addressed in that other possible, if not probable, nucleation and growth models need to be tested to determine what model is best accounts for the formation of the observed “ZnO”—really more closer to (ZnO)n·(C18H33O2−)x—nanoparticles.
In short, the authors study zinc-oxide particle formation using state-of-the-art tandem in situ SAXS/WAXS and UV-vis techniques. A single, non-disproof-based physical–chemical model is provided with ≥4 adjustable parameters that doesn’t fully fit the observed kinetics data. Additional studies that would be welcome include: (i) writing out proposed pseudoelementary steps for the reaction; (ii) obtaining an experimental rate law in terms of [ZnII], [OH−], and [C18H33O2−] as a minimum; and then (iii) testing the proposes as well as alternative models to determine the best, disproof-based model for ZnO nanoparticle formation. Importantly, a more reliable, disproof-based mechanism for ZnO nanoparticle formation could result from just one set of additional studies based on this valuable, initial, once again interesting work.
WCl6 + 8.7BnOH → WO3·H2O + 3.3BnCl + 2.7BnOBn + 2.7HCl + 2H2O | (22) |
Fig. 26 EXAFS measurements were collected at the W L3-edge at around 10207 eV and then analyzed using LCA to product the relative amount of each species as a function of time in seconds for the tungstite formation reaction. Reproduced with permission from ref. 110. Copyright 2012 Wiley-VCH Verlag GmbH & Co. |
Next, to investigate the morphological changes during the reaction, the authors used SAXS and TEM. Their reported TEM micrographs are shown in Fig. 27.
Fig. 27 “Time-dependent TEM and HRTEM images of the tungstite nanostructures after different reaction times (a) and (b) 10 min, (c) and (d) 20 min, (e) 60 min and (f) 240 min.”111 Figure reproduced with permission from ref. 111. Copyright 2013 Royal Society of Chemistry. |
From these studies, they observed the formation of spherical particles that assembled into cylinders, and then exfoliated into nanoplatelets. A proposed reaction scheme is provided as Scheme 7.
Scheme 7 Proposed reaction scheme for the formation of tungstite nanoplatelets starting from WCl6 in benzyl alcohol. Reproduced with permission from ref. 111. Copyright 2013 Royal Society of Chemistry. |
Overall, Neiderberger and coworkers have demonstrated the applicability and usefulness of synchrotron XAFS and SAXS measurements to the synthesis of tungstite nanoplatelets. They were able to determine the reaction stoichiometry and then follow the reaction both directly and indirectly. Even with the multiple techniques and depth of analysis, the exact understanding of why, and by what mechanism, the cylinders exfoliate to nanoplatelets is presently unknown. The authors comment that “these two studies exemplarily show how elaborate the investigation of just one system can be and how many ex and in situ techniques have to be applied to monitor both the organic and inorganic processes,”111 sage advice for other workers. Also, “the study of the reaction mechanism over the course of the whole reaction from the dissolution of the precursor to the final nanomaterial is not straightforward and dedicated analytical tools still have to be developed,”110 more truths according to all the present review has uncovered or we know.
Hence, while these two papers uncover a number of insights into formation of non-spherical shapes, mechanistic work remains to be done on these systems. One could conceivably take their proposed reaction scheme, strive to write pseudoelementary step reactions, and see if it is possible to fit their hard-won kinetics data. Furthermore, as the authors state “it is worth the efforts to study as many systems as possible, because every mechanistic study contributes to the overall picture of nanoparticle formation and helps to systematically extend the crystallization models.”111 The value of this system and others like it,113 striving for mechanistic understanding of shape control, is high with a bright future for the needed additional studies.
Fig. 28 (A) “Complexes formed at 80 °C from reaction of CuBr(OLAM)3 with TOPO and TOP, where TOP, TOPO and OLAM ligands are represented by their donor atoms P, O and N, respectively.”112 (B) In situ EXAFS spectra. (C) 1H NMR spectra at room temperature. “OLAM C2N2 resonances are highlighted in red, the broadened C2NH2 resonances are marked by * and shown separately in the inset plot.”112 (D) 31P{1H} NMR spectra. (E) “Mass spectra of reaction mixtures measured by MALDI-TOF for TOPO and TOP complexes and ESI-MS for the OLAM complex; M is the molecular ion, and the fragmentation or adduct formation required to give the simulated isotope patterns (gray) in indicated.”112 Reproduced with permission from ref. 112. Copyright 2019 American Chemical Society. |
With the critical precursor speciation in hand, Buonsanti and coworkers conducted a series of kinetics experiments to try and elucidate the nucleation pathway. The authors propose two potential pathways for Cu(0) formation, but were not able to definitively support one over the other. However, they provide the community with a description of the needed, additional experiments. They state that “EXAFS combined with X-ray scattering and relative PDF analysis can help in catching the nucleation event and in identifying the actual monomer species.”112
The authors also use their understanding of the reaction kinetics and thermodynamics to access different final shapes. Specifically, they correlate “the conversion kinetics of the intermediates to the Cu(0) monomer flux and the final shape of the NC products… to explore the free energy landscape of the reaction.”112 Those experiments were conducted by varying two parameters: the temperature of the injection and the speed of the injection. Reactions were carried out at 260 °C and 335 °C and injections were done either rapid (all at once) or slow (dropwise). The resulting TEMs from these experiments are reproduced herein as Fig. 29.
Fig. 29 Representative TEM images of Cu NCs obtained after (A and B) dropwise and rapid injection of CuBr(OLAM)2(TOPO) at 260 °C, respectively; (C and D) dropwise and rapid injection of CuBr(OLAM)2(TOPO) at 335 °C, respectively; (E and F) dropwise and rapid injection of {CuBr(TOP)2}2; and (G and H) dropwise and rapid injection of {CuBr(TOP)2}2 at 335 °C, respectively.112 Reproduced with permission from ref. 112. Copyright 2019 American Chemical Society. |
Overall, Buonsanti and coworkers have compellingly demonstrated the value and power of elucidating the precursor speciation. By doing this, they are able to probe the kinetic and thermodynamic control of their reaction, which resulted in control over the final shape of the particles. The researchers were able to obtain particles with spherical, cubic, octahedral, truncated-octahedral, and tetrahedra morphologies. Access to these different shapes was due to their in-depth characterization of the reaction speciation using a multitude of techniques (EXAFS, NMR, and MS). They state that “as a general lesson on the synthesis of colloidal metal NCs, we have learned that the choice of the precursor and the identification of the reaction intermediates are crucial to enable a more rational approach to shape control.”112 Well said. Here, the authors have satisfied the first criterion for a mechanistic study of knowing their reaction—which includes knowing the precursor speciation. They also give attention to disproving potential alternative hypotheses during their experiments. In the future, it would be of great interest to try and fit their kinetics data by deriving the equations that correspond to their proposed nucleation pathways. Overall, the study by Buonsanti and coworkers is an exciting and cutting-edge advancement in elucidating shape control, a leading study in an exciting, still wide-open subfield of particle formation synthesis, kinetics and associated mechanism(s).
Amorphous nanoparticles are of interest due to their applications in catalysis, optics, magnetism, bioactivity, electrochemistry, and nanocomposites.128–134 In particular, there are cases within the water oxidation catalysis literature where the amorphous materials have produced superior catalytic activity to their crystalline counterparts.135,136 For the formation of amorphous nanoparticles, there are a limited number of studies that use synchrotron XAFS or SAXS techniques to investigate the local structural changes of what are typically amorphous nanoparticles at least initially.137–140
A good example is a 2013 study on the synthesis of Ni2P, where the authors used to two similar syntheses to produce amorphous and crystalline Ni2P nanoparticles.141 They used EXAFS combined with XRD and TEM results to develop a method to “determine crystalline behavior based on standard deviations in interatomic distance offsets.”141 Then, when combined with XRD, one can determine the level of crystalline versus amorphous character in their sample.141 Overall, there is great potential for synchrotron XAFS in the formation of amorphous nanoparticles, particularly when paired with an addition technique like XRD or PDF.
A final field of particular interest in the past decade, that a referee asked we comment on, is carbon-based quantum dots (CQDs or C dots).142–144 These materials are of interest for their tunable photoluminescent (PL) properties and applications in energy, catalysis, and biomedical engineering.145–148 There are no current studies that we could find that utilize synchrotron XAFS or SAXS to study the formation of CQDs. A number of papers, including an excellent review by Shamsipur and coworkers,149 have detailed the various synthetic strategies to produce CQDs that fluorescence from blue to red.150–152 However, the field of CQDs remains wide open for investigation by synchrotron-based methods. As Shamsipur and coworkers wrote, “it should be noted that most of the reported [CQDs] with unique PL properties were obtained by chance and not by design.”149 In short, a bottom line regarding these additional systems of often high current interest is that each is wide open for kinetics and other studies aimed at a detailed understanding of their formation mechanism.
The LaMer model166 postulating “burst” nucleation and “diffusion-controlled” growth was developed for sulfur sols and has since been used to try to explain hundreds of particle formation systems.65,66,167 Yet, in its 70 years of existence, the LaMer mathematical model and equation has never been used to successfully fit kinetics data save his original case on (S)n sols,166 a remarkable testament to its lack of quantitative utility. Moreover, the LaMer qualitative hypothesis of “burst” nucleation and “diffusion-controlled” growth has been compellingly disputed25,26,168–173 in three recent65,66,167 reviews of 1953 papers citing the LaMer model over the 70+ years since 1950. Relevant here are the wise words of B. Peters that “The classic theories work poorly or not at all when applied beyond the validity of their assumptions.”174
Little known is that Victor LaMer also published a stepwise mechanism for particle formation in 1952,175 but never tested it with nor used it to explain experimental kinetics data for particle formation. Turkevich's 1953 “organizer mechanism”176,177 merits mention in Table 1 as well, even though it lacks specific pseudoelementary steps, because it is the early example of what are now termed Pre-Nucleation Clusters as discussed in a 2021 review.66 Indeed, reading Turkevich's literature reveals he was often far ahead of his time in his intuition and understanding about particle synthesis, physical properties, and probable mechanisms.
Year | Person(s) | Mechanism (A, B, and C formalism) | Mechanism (in words), and comments |
---|---|---|---|
1952 | V. K. LaMer175 | mA1 ⇌ Bm Bm + A1 ⇌ Bm+1 B(i−1) + A1 ⇌ Bi | Stepwise bimolecular addition leading to “burst nucleation” (NB: “diffusion-controlled growth”166 was not included175 as part of the scheme provided) |
[A1 = “kinetically independent unit of phase A”, Bi = “embryo of phase B containing i units”, and “m is usually 2, whereas i in the case of water vapor condensation is about 80—i.e., this reaction is of the 80th order”] | |||
1953 | Turkevich,176,177 Stevenson, Hillier | No pseudoelementary step mechanism was provided in ref. 176 and 177. | Nucleating agents build up a chemical complex, organize macromolecularly, then undergo a molecular rearrangement to produce the metal nucleus, Turkevich's “organizer” mechanism/hypothesis—what can now be seen to be an early version of the prenucleation clusters postulate. |
1997 | Finke, Watzky54 | Slow, kinetically continuous nucleation followed by typically fast, explosive autocatalytic surface growth. Important to note is that higher order nucleations are covered by (i.e., can be hidden in) the pseudoelementary first step, as the subsequent discovery of second-order and termolecular nucleations demonstrates.30,31 | |
[A = metal precursor, B = growing, metal(0) particle] | |||
2004 | Hornstein, et al.184 | Evidence for bimolecular agglomeration as a 3rd step added to the 1997 2-step mechanism54 | |
2005 | Besson, Finney, et al.80,185,186 | Evidence for a 4th step of autocatalytic agglomeration where the smaller particles, B, and larger particles, C, agglomerate with a separate, k4, rate constant | |
[A = metal precursor, B = small particles, C = larger particles] | |||
2005 | Drew, Katsoulakis, Tsapatsis188 | A proposed mechanism for crystal growth by aggregation of meta-stable nanoparticle precursors, one used with Population-Balance Modeling. Three alternative mechanisms were tested; shown is the author's preferred Mechanism II. | |
[A = primary colloidal particles; B = growing crystal; C1 = nucleated primary particles] | |||
[b = nAoKi/k2; c = nAoKij/k2; nAo = initial number of fresh primary particles; ni = particle density of crystal i; y = /nAo; τ = k2t] | |||
2014 | Kent, Mondloch, et al.187 | Evidence for a new 4th step of secondary autocatalytic growth where larger particles, C, growth with a separate, k4′, rate constant. | |
[A = metal precursor, B = smaller particles, C = larger particles] | |||
2017 | Karim, et al.39 | A·L ⇌ A + L | Proposed four step mechanism, including two ligand binding equilibria steps, for the formation of metal nanoparticles, the first A·L step having additional precedent.31,32 |
A set of 5 differential equations are reproduced in the ESI. | |||
B·L ⇌ B + L | |||
[A = metal precursor, B = growing particle, A·L = precursor with ligand, B·L = growing particle with ligand; both ligand steps have associated Keq equilibrium constants] | |||
2019 | Handwerk, Shipman, Whitehead, Özkar, et al.67,68 | A new, 3-step mechanism, discovered by disproof-based Mechanism-Enabled Population-Balance Modeling, in which the pseudoelementary step of the 1997 2-step mechanism was expanded into two growth steps (or, alternatively, and as actually discovered67,68 the 2nd and 3rd steps of the 2014 4-step mechanism are combined into a single, new 2nd step). Eleven alternative mechanisms were disproved. Quantitative fitting of Particle-Size Distributions (PSDs) was demonstrated, including the PSD shape. | |
[A = metal precursor, B = smaller particles, C = larger particles] |
The first, disproof-based minimum mechanism for particle formation mechanism composed of specific pseudoelementary steps (that also add up to an experimentally demonstrated balanced reaction to well-characterized products) did not appear until 1997,54 the 2-step mechanism of continuous nucleation and autocatalytic surface growth listed in Table 1. This simplest version of particle formation has since been used to successfully fit homogeneous nanoparticle formation,178 heterogeneous nanoparticle formation,179 solid state kinetics,180 dye aggregation,181 protein aggregation,182 and other “cooperative” autocatalytic phenomena in nature.183 Six upgraded, themselves deliberately minimum mechanisms for particle formation have since been reported39,54,67,68,80,184–187 as summarized in Table 1. Those specific upgrades to the 2-step mechanism occurred 5 times as additional evidence for more complex pathways became available (entries 2, 3, 5, 6, and 7 – Table 1)—and include disproof57–59 of at least 21 additional alternatively mechanisms en route to those 6 of the 7 mechanisms in Table 1.26,54,67,68,80,184–187 The 2005 crystal-growth model from Tsapatsis and co-workers188 is also included in Table 1 because their model has proposed pseudoelementary steps supported by Population Balance Modeling while disproving one alternative mechanism—and hence is an early example of a mechanistic model partially tested by what is denoted ME-PBM.67,68
An important point growing out of Table 1 is that the 1997 2-step mechanism was upgraded67,68—“disproved” if one likes—in favor of the 2019 3-step mechanism as needed to be able to quantitatively account for PSD data.67,68 Put in other terms, although the 2-step mechanism has been used in several studies in this review—and remains recommended67,68—for initial data analysis of sigmoidal particle-formation kinetics curves, a new 2019 3-step mechanism67,68 is now the minimal PEStep mechanism able to account for PSDs once programed into Mechanism-Enabled Population Balance Modeling.67,68
Of considerable interest is what results if one takes just the 5 basic classes of mechanistic models in Table 1 that have been subjected to extensive disproof and therefore deemed to be more reliable,54,67,68,80,184–187 and then combines them with 3 different precedented nucleation molecularities30–32,54 and up to 3 different types of ligand binding steps based on A. Karim's groundbreaking39 work on ligand-dependent PESteps and A·L and B·L species (plus other precedent for A·L to A + L dissociative steps31,32). Looking at the possible combinations, one now has 5 classes of ≥96 possible mechanisms that one can use to try to use in a disproof-based approach57–59 to fitting one's particle formation data while striving for a disproof-based minimum mechanism for one's particle formation reaction.
• Employing synchrotron X-ray radiation-based experiments provides powerful, largely direct particle formation kinetics data. Indeed, Tables S1–S4 of the ESI† document the wealth of typically expert obtained, in situ, many times creative tandem-methods monitoring of particle formation reactions. One caveat here, that merits closer attention in the future, is that X-ray radiation induced radiolysis can initiate nucleation events that then quickly magnify the loss of precursor by exponential, autocatalytic growth54 in at least some cases.
• The value in completely characterizing ones’ precursors and the overall speciation in one's starting solutions cannot be overemphasized. This follows both from experience as well as the fact that kinetics generally gives one the composition of the activated complex of the rate-determining step, which in turn relies on knowing what species you are actually changing when you vary the concentrations of the starting materials—that is, relies on you being able to confidently write d[Product]/dt in terms of known species.
• Four minimum criteria one needs to satisfy in order to be able to claim a chemical mechanism of particle formation for a given system are the following,65 criteria well-known from physical-organic and physical-inorganic mechanistic chemistry: (i) a complete mass- and charge-balanced, experimentally determined reaction stoichiometry—because the proposed mechanistic steps must add up to that observed reaction, for example the presence of ligand-based steps such as particle capping; (ii) kinetics data, ideally the full rate law for nucleation, growth and any agglomeration that might be present—again in terms of known precursor species; (iii) postulation of pseudoelementary step-based mechanisms, such as those in Table 2, that are also able to add up to the observed reaction stoichiometry, followed by writing down the differential equations corresponding to each mechanism; and critically then also (iv) disproof of ideally all reasonable alternative mechanistic hypotheses by checking each mechanism against the observed kinetics data en route to an Ockham's razor obeying, experimentally supported, disproof-based proposed minimum mechanism of particle formation for the system at hand.
Class | Mechanism | Unimolecular nucleation | Bimolecular nucleation | Termolecular nucleation | Ligand dependence |
---|---|---|---|---|---|
(I) | FW 2-step | A → B | 2A → 2B | 3A → 3B | (i) A·L ⇌ A + L |
A + B → 2B | A + B → 2B | A + B → 2B | (ii) B·L ⇌ B + L | ||
(iii) With A·L and B·L steps | |||||
(II) | Classic 3-step with bimolecular agglomeration | A → B | 2A → 2B | 3A → 3B | (i) A·L ⇌ A + L |
A + B → 2B | A + B → 2B | A + B → 2B | (ii) B·L ⇌ B + L | ||
B + B → C | B + B → C | B + B → C | (iii) C·L ⇌ C + L | ||
(iv) With A·L and B·L steps | |||||
(v) With A·L and C·L steps | |||||
(vi) With B·L and C·L steps | |||||
(III) | New, Handwerk et al. 3-step | A → B | 2A → 2B | 3A → 3B | (i) A·L ⇌ A + L |
A + B → C | A + B → C | A + B → C | (ii) B·L ⇌ B + L | ||
A + C → 1.5C | A + C → 1.5C | A + C → 1.5C | (iii) C·L ⇌ C + L | ||
(iv) With A·L and B·L steps | |||||
(v) With A·L and C·L steps | |||||
(vi) With B·L and C·L steps | |||||
(IV) | 4-Step with B + C agglomerative growth | A → B | 2A → 2B | 2A → 2B | (i) A·L ⇌ A + L |
A + B → 2B | A + B → 2B | A + B → 2B | (ii) B·L ⇌ B + L | ||
B + B → C | B + B → C | B + B → C | (iii) C·L ⇌ C+ L | ||
B + C → 1.5C | B + C → 1.5C | B + C → 1.5C | (iv) With A·L and B·L steps | ||
(v) With A·L and C·L steps | |||||
(vi) With B·L and C·L steps | |||||
(V) | 4-Step with A + C secondary autocatalytic surface growth | A → B | 2A → 2B | 3A → 3B | (i) A·L ⇌ A + L |
A + B → 2B | A + B → 2B | A + B → 2B | (ii) B·L ⇌ B + L | ||
B + B → C | B + B → C | B + B → C | (iii) C·L ⇌ C + L | ||
A + C → 1.5C | A + C → 1.5C | A + C → 1.5C | (iv) With A·L and B·L steps | ||
(v) With A·L and C·L steps | |||||
(vi) With B·L and C·L steps |
• Even a brief survey of Tables S2–S4 of the ESI,† that we hope the reader has done for themselves, teaches several other important points: (i) at least ∼94% of the time to date, the systems examined are not well defined, without for example even a complete, balanced reaction stoichiometry having been determined. In that sense, more in-house, more “synthetic” efforts before “synchrotron spectroscopic” efforts—actually, much, much more—would, in many cases, have provided better defined systems worthier of the sophisticated, expert, often creative synchrotron XAFS and SAXS efforts one finds listed in Tables S2–S4 of the ESI.† Put another way, one can view all 69 cases in Tables S2–S4 (ESI†) as using synchrotron radiation methods to examine particle formation systems directly and see what insights can be obtained. That is fine, but one should not claim that a mechanism has resulted from those largely survey studies. (ii) In 4 of the 69 cases in Tables S2–S4 of the ESI,† all four of the requirements (i)–(iv) for a reliable mechanism are missing; in 39 of the 69 cases (57% of the time) in Tables S2–S4 of the ESI,† three of the four requirements are missing; in 20 cases (29% of the time) two of (i)–(iv) are missing; and in only 4 cases (5.8% of the time) are just one of requirements (i)–(iv) for a reliable mechanism missing.
• Worth emphasizing is that the required but missing first component in ca. 94% of particle formation studies examined by this review is the balanced reaction stoichiometry. Without a balanced reaction stoichiometry, the pseudoelementary steps of the proposed mechanism—that must add up to that net reaction stoichiometry—cannot be reliably nor unequivocally determined. Even for the 6 studies in the main text, only one has (and only 2 more partially have) an experimentally determined, balanced reaction stoichiometry. The studies examined herein, as well as one's intuition, teach that a balanced reaction stoichiometry promises to be especially important for obtaining and understanding ligand effects in particle formation reactions, as A. Karim's valuable, telling example demonstrates.
• Also missing in 66 of 69 examples in Tables S2–S4 (ESI†) and in 3 of the 6 cases in the main text is a full, experimentally determined rate law for nucleation and growth. One cannot possibly have a reliable reaction mechanism without the full rate law. Rate laws for nucleation are especially rare,31 being present in arguably only 1 case32 of the 64 total studies covered in this review—and even then that proposed nucleation mechanism is currently under further scrutiny using ME-PBM.
• Too often presently mechanisms are claimed based on qualitative data or ‘snapshots’ and sometimes cartoons of what size particles are present at what times. These qualitative schematic models or just word-based pictorial models189–197 are of some use, but are not, and should not be confused with, disproof-based, minimum mechanisms expressed in balanced reactions as detailed in Tables 1 and 2.
• Presently, one can write 5 primary classes of pseudoelementary-step based, minimum mechanisms, that include 3 precedented types of nucleation and 3 precedented types of ligand binding, for a total of ≥96 possible mechanisms now available to test out on one's experimental data. Relevant here is that only two studies have elucidated the effects of ligands on the kinetics and mechanism of nucleation and growth.31,39,72,95 Hence, especially important for future studies is to include the specific effects of ligands on particle nucleation, growth, agglomeration and associated ligand-capping (of the particle) processes.
• Mechanism-Enabled Population Balance Modeling has recently emerged as a must-use tool to test one's mechanism, specifically to see if the proposed mechanism can quantitatively account for the proposed particle-size distribution, including its shape.
• One should expect, and indeed hope, that one's initial, minimum mechanism will one day be disproved, replaced by a more complex, experimentally further supported pathway. This is the expectation of minimal, disproof-based mechanisms, where, as time progresses and techniques improve, new insights and new details become possible as new methods, additional experimentation, or applicable computations become available. Relevant here are recommend articles by Platt,57 Chamberlain,58 Scott,59 and Weinfurt198 on employing a proper disproof-based scientific method and avoiding the >400-year-old problem in science of “idols and false notions”.198
Footnote |
† Electronic supplementary information (ESI) available: Summary of supporting information; case study #6: additional figures and information; differential equations for the Karim 4-step mechanism; Table S1 Summary of reviews on synchrotron XAFS and/or SAXS; Table S2 Summary of papers using SAXS to collect kinetics data; Table S3 Summary of papers using XAFS to collect kinetics data; Table S4 Summary of papers using tandem techniques: use of synchrotron XAFS and/or SAXS in combination with another technique. See DOI: 10.1039/d1ma00222h |
This journal is © The Royal Society of Chemistry 2021 |