Nanoparticles of nickel phosphide are finding wide ranging utility as catalysts for hydrodesulfurization, hydrogen evolution reaction, and hydrodeoxygenation of bio-oils. Herein, we present a methodology to tailor monodisperse nickel phosphide nanoparticles in terms of size and phase through the use of a statistical response surface methodology. Colloidal nickel phosphide nanoparticles were synthesized by replacing octadecene (ODE), a commonly used organic solvent, by a more sustainable phosphonium-based ionic liquid (IL). The replacement of ODE with the phosphonium-based IL resulted in faster crystallization at lower temperatures to yield phase-pure, monodisperse Ni2P nanoparticles. Using a first-order design, the PPh3/Ni precursor ratio was identified as the most critical factor influencing the resulting size and phase of the nanoparticles. Optimization using a Doehlert matrix for second-order design yielded a second-degree polynomial equation used to predict the mean diameter of the nanoparticles (over a range of 4–12 nm) as a function of the PPh3/Ni precursor ratio and the temperature used during synthesis. The resulting model was validated by performing reactions using randomly chosen sets of conditions; the experimentally determined nanoparticle sizes were in excellent agreement with the theoretical sizes predicted by our model. This demonstrates the utility of a multivariate experimental design as a powerful tool in the development of synthetic strategies toward the preparation of colloidal nanoparticles with highly controlled size, size distribution, and phase.

The understanding and control of colloidal nanocrystal syntheses are essential for discovery and optimization of desired properties and therefore play a key role in the applications of these materials. Typical one variable at a time (OVAT) methods limit the ability of researchers to achieve such goals by providing one-dimensional insight into a complex, multidimensional experimental domain, wasting precious resources in the process. Design of experiments (DoE) in conjunction with response surface methodology (RSM) offers an accelerated route for multivariate investigation and optimization of nanocrystal syntheses. The method enables systematic analysis and multidimensional modeling of the independent and dependent effects that any number of factors have on chosen responses, resulting in easy optimization of a large synthetic space in a fraction of the experiments. Herein, we will outline the general steps to follow when utilizing DoE and RSM for screening and optimization of nanocrystal syntheses, as well as the background needed to appropriately design an investigation and understand the results.

Thiospinels, such as CoNi2S4, are showing promise for numerous applications, including as catalysts for the hydrogen evolution reaction, hydrodesulfurization, and oxygen evolution and reduction reactions; however, CoNi2S4 has not been synthesized as small, colloidal nanocrystals with high surface-area-to-volume ratios. Traditional optimization methods to control nanocrystal attributes such as size typically rely upon one variable at a time (OVAT) methods that are not only time and labor intensive but also lack the ability to identify higher-order interactions between experimental variables that affect target outcomes. Herein, we demonstrate that a statistical design of experiments (DoE) approach can optimize the synthesis of CoNi2S4 nanocrystals, allowing for control over the responses of nanocrystal size, size distribution, and isolated yield. After implementing a 25–2 fractional factorial design, the statistical screening of five different experimental variables identified temperature, Co:Ni precursor ratio, Co:thiol ratio, and their higher-order interactions as the most critical factors in influencing the aforementioned responses. Second-order design with a Doehlert matrix yielded polynomial functions used to predict the reaction parameters needed to individually optimize all three responses. A multiobjective optimization, allowing for the simultaneous optimization of size, size distribution, and isolated yield, predicted the synthetic conditions needed to achieve a minimum nanocrystal size of 6.1 nm, a minimum polydispersity (σ/d̅) of 10%, and a maximum isolated yield of 99%, with a desirability of 96%. The resulting model was experimentally verified by performing reactions under the specified conditions. Our work illustrates the advantage of multivariate experimental design as a powerful tool for accelerating control and optimization in nanocrystal syntheses.

Transition metal carbides (TMCs) have attracted significant attention because of their applications toward a wide range of catalytic transformations. However, the practicality of their synthesis is still limited because of the harsh conditions in which most TMCs are prepared. Recently, a solutionphase synthesis of phase-pure α-MoC 1−x nanoparticles was presented. While this synthetic route yielded nanoparticles with exceptional catalytic performance, the reaction parameter space was not explored, and catalyst throughput was not optimized for scale-up. Continuous flow platforms coupled with statistical design of experiments (DoE) can provide a powerful method for understanding the reaction parameter space for optimizations. Here, we demonstrate the use of statistical DoE in tandem with response surface methodology for a parametric screening analysis to optimize the throughput of a MoC 1−x nanoparticle synthesis utilizing a millifluidic flow reactor. A full factorial design was implemented to evaluate four input variables (reaction temperature, flow rate, solvent fraction of oleylamine, and precursor concentration) that carry statistically significant effects on three responses (throughput, residence time, and isolated yield). A Doehlert matrix was implemented to investigate each significant variable at a higher number of levels to optimize throughput. Our results give a nonintuitive set of experimental conditions that resulted in an optimized throughput of 2.2 g h −1 . This translates to a 50-fold increase in throughput compared to the previously reported batch method. The catalytic performance of the MoC 1−x nanoparticles produced under optimized throughput was demonstrated in the CO 2 hydrogenation reaction. This DoE screening analysis and throughput optimization of MoC 1−x synthesis open the door to an increased feasibility for scale-up.

Copper selenides are an important family of materials with applications in catalysis, plasmonics, photovoltaics, and thermoelectrics. Despite being a binary material system, the Cu–Se phase diagram is complex and contains multiple crystal structures in addition to several metastable structures that are not found on the thermodynamic phase diagram. Consequently, the ability to synthetically navigate this complex phase space poses a significant challenge. We demonstrate that data-driven learning can successfully map this phase space in a minimal number of experiments. We combine soft chemistry (chimie douce) synthetic methods with multivariate analyses via classification techniques to enable predictive phase determination. A surrogate model was constructed with experimental data derived from a design matrix of four experimental variables: C–Se bond strength of the selenium precursor, time, temperature, and solvent composition. The reactions in the surrogate model resulted in 11 distinct phase combinations of copper selenide. These data were used to train a classification model that predicts the phase with 95.7% accuracy. The resulting decision tree enabled conclusions to be drawn about how the experimental variables affect the phase and provided prescriptive synthetic conditions for specific phase isolation. This guided the accelerated phase targeting in a minimum number of experiments of klockmannite CuSe, which could not be isolated in any of the reactions used to construct the surrogate model. The reaction conditions that the model predicted to synthesize klockmannite CuSe were experimentally validated, highlighting the utility of this approach.