Stylianos
Kioumourtzoglou
*a,
Sebastian
Hof
a,
Cécile
Kalk
a,
Viktor
Toth
b,
Mikaela
Görlin
c,
Jaroslava
Nováková
d and
Jacinto
Sá
*ae
aDepartment of Chemistry-Ånsgtröm, Physical-Chemistry Division, Uppsala University, Lägerhyddsvägen 1, Uppsala 751 20, Sweden. E-mail: stylianos.kioumourtzoglou@kemi.uu.se; jacinto.sa@kemi.uu.se
bToptal, LLC, 2810 N. Church St #36879, Wilmington, DE 19802-4447, USA
cDepartment of Chemistry-Ånsgtröm, Structural Chemistry Division, Uppsala University, Lägerhyddsvägen 1, Uppsala 751 20, Sweden
dDepartment of Surface and Plasma Science, Charles University, V holesovickach 2, Prague 8, 18000, Czech Republic
eInstitute of Physical Chemistry, Polish Academy of Sciences, Marcina Kasprzaka 44/52, Warsaw 01-224, Poland
First published on 11th June 2024
Establishing scalable nanomaterials synthesis protocols remains a bottleneck towards their commercialisation and, thus, a topic of intense research and development. Herein, we present an automated machine-learning microfluidic platform capable of synthesising optically active nanomaterials from target spectra originating from prior experience, theorised or published. Implementing unsupervised Bayesian optimisation with Gaussian processes reduces the optimisation time and the need for prior knowledge to initiate the process. Using PTFE tubing and connectors enables facile change in reactor design. Ultimately, the platform substitutes the labour-intensive trial-and-error synthesis and provides a pathway to standardisation and volume synthesis, slowing down the translation and commercialisation of high-quality nanomaterials. As a proof-of-concept, Ag nanoplates and Prussian-blue nanoparticle protocols were optimised and validated for volume production.
New conceptsWe report a new automated machine-learning microfluidic platform capable of synthesising optically active nanomaterials from target spectra originating from prior experience, theorised or published. The differentiation between the published research is twofold: (i) the use of unsupervised Bayesian optimisation with Gaussian processes to reduce the optimisation time and the need for prior knowledge to initiate the process, and (ii) microfluidic reactors made from PTFE tubing that enables customisation of design at low cost. Eventually, the platform substitutes the labour-intensive trial-and-error synthesis and provides a pathway to standardisation and volume synthesis, slowing down the translation and commercialisation of high-quality nanomaterials. The proposed approach's conceptual advantage is that it enables the production of particles with the same quality and morphology as the optimised protocol in volume and anywhere. Ultimately, the innovation driver is the establishment of the Nanomaterials as a Service (NaaS) concept, where end-users request particles with specific characteristics, and the manufacturer can rapidly optimise a protocol and produce in volume the desired nanomaterial or provide the protocol for production at the end-user site, which is yet to be achieved. |
The exponential growth in nanoscience materials synthesis and theoretical capabilities, allied with open access efforts, reduced the challenge in accessing target data necessary to implement machine-learning approaches to reduce labour-intense synthesis protocol optimisation. However, the development of NaaS concepts remains constrained by two significant challenges: lack of low-cost and flexible reactor design and deficiency in reliable and fast unsupervised optimisation algorithms.
Regarding reactor design, microfluidics have emerged as a powerful tool for producing high-quality nanoparticles.3–6 However, microfluidic systems often use pre-established reactor geometries called reactor chips. Their design geometry rigidity demands significant prior knowledge from the user when selecting them or extensive trial-and-error attempts that are unsustainable and expensive. Furthermore, chip integrity, resistance to fouling and ageing are also things to be considered, especially when running unsupervised optimisations. These issues can be mitigated with routine cleaning processes (increasing the optimisation time)7 and multi-point analysis of the reactor's health (increasing cost and complexity).8 Additionally, to be effective, the multi-point analysis must be fast and non-invasive, reducing the range of analytics to in-line optical methods that might not be sensitive to the species, causing reactor fouling or ageing. Therefore, when it comes to reactors, their geometry should be changeable and adaptable, made of durable fouling-resistant materials, easy to integrate with pumping and analytical systems, and cost-effective.
The growing interest in machine learning as a problem-solving tool led to the rapid proliferation and development of optimisation algorithms to a level that has become a deterrent to nano-innovators' use due to the overwhelming amount of choice. The decision of which method to use often comes from weighing the following parameters: how big is the parameter space, how continuous the process is, how much prior knowledge the user has, and how fast it converges. The latter does not relate to the time it takes to optimise, but the number of data points one must acquire to get the desired result.
A large class of nanomaterials manufacturing uses kinetically controlled synthesis methods, where particles are made from precursor salts, reducing or oxidising chemicals and growth-directing agents. Examples of nanomaterials produced via this approach include plasmonic nanoparticles, semiconductors, quantum dots, etc.9 Their synthesis has a characteristically low number of continuous control parameters (less than ten) that are easy to generate data points, making them particularly suitable for Bayesian optimisation with Gaussian processes. Bayesian optimisation with Gaussian processes is an unsupervised algorithm that converges fast and requires relatively low computational power.
This communication reports a development that quells bottlenecks surrounding the discovery, optimisation, standardisation and development of the NaaS concept to fabricate high-quality nanomaterials. The workflow is schematically representation in Fig. 1. The optimised protocols consist of the exact pump flows with specific substrate concentrations and the reactor configuration, which should be transferable to any production facility.
Fig. 1 Schematic workflow representation for optimising a synthesis protocol with the proposed approach. |
As a proof-of-concept, we optimised and standardised protocol for synthesising plasmonic silver nanoplates and Prussian blue nanoparticles, demonstrating the breadth of materials that can be produced with the presented approach. Plasmonic nanoparticles can capture and concentrate distant radiation within subwavelength regions defies diffraction limits,10,11 resulting in powerful near-fields, hot carriers and localised heat,12,13 which spurred the development applications in biosensing,14 cancer therapy,15 photovoltaic,16,17 photodetectors,18 catalysis,19,20 emitting devices,21,22 solar energy converters,23 among others. Prussian blue, which people also call ferric ferrocyanide, is used in catalysis,24 hydrogen storage,25 photochemistry,26 chemical and bio-sensors,27 clinical medicine,28 and most promisingly, high crystalline Prussian blue analogues are used as a cathode on sodium-ion batteries.29,30
Polytetrafluoroethylene (PTFE) is a hydrophobic, non-wetting, high-density, resistant to high temperatures (m.p. 327 °C), versatile and best known for its non-stick properties material. These properties make PTFE an ideal reactor material. Additionally, PTFE tubes with 1/16′′ outer diameter (OD) can be bought with a variable inner diameter (ID) ranging from 0.2–1 mm, ideally suited for microfluidics. The fix OD enables easy connectivity through T and cross-tube fitting connections and offers constant connector dimensions to pumps, the detection system and the reactor. Maintenance costs are dramatically reduced for new tubing and connectors. This versatility and variety permit the creation of tailor-made reactor geometries and continuous experimentation and optimisation of the geometry, as demonstrated by the proof-of-concept examples herein. The precise description of the reactor used for each proof-of-concept synthesis is shown in the Sections 2.3 and 2.4.
Our goal in employing Bayesian optimization is to understand and optimize an unknown function f(x) based on sampled data. We postulate that observations might be noisy, such that an observed value f(x) equals the true function value plus an error term ε, where the expected value of ε is zero and independent across samples, i.e.,
The process involves the following steps: at any given step n, after observing the pairs (x1, f(x1)),…,(xn,f(xn)), we infer the function f at new points not in the existing observations using a Gaussian process. This yields a normal distribution for each untested point, characterized by a calculated mean and variance. We then select the next point xn + 1 to sample by maximizing an acquisition function g(x) which considers both the mean and variance, effectively balancing between exploring new areas and exploiting known ones. An example of such an acquisition function is the upper confidence bound (UCB). After selecting xn + 1, we observe and record the new noisy sample f(xn + 1)).
The choice of Bayesian optimization is justified for nanoparticle synthesis due to the typically small number of controllable parameters, usually less than 10, such as chemical concentrations and temperatures. The continuous nature of f(x), reflecting changes in spectral outputs with variations in synthesis conditions, supports the use of this method. Additionally, the frequency of updating x is deliberately low to manage computational load, though this can be adjusted for faster synthesis methods like hot injections.
In summary to use Baysien optimization for nanomaterials synthesis it must keep on the following aspects:
• Small number k of control parameters, as a rule of thumb k < 20 because a nonconvex problem in is solved to choose the next point xn. This condition is easily fulfilled in nanoparticle synthesis since the variables are chemicals and temperature used, which rarely exceeds 10 controllable parameters.
• f is continuous for the Bayesian optimization to leverage the fact that for x ≈ x′ → f(x) ≈ f(x′). This is also valid in most syntheses since small changes in reactant concentrations rarely produce completely different nanoparticle outputs.
• (x) is defined as the measured spectral fit at control parameters x to the target spectrum. The fit of two spectra may be defined in various ways to establish a continuous function f.
• x is updated at a low frequency because deciding the next point xn is a computationally intensive task. This is once more easily achievable in kinetically controlled synthesis processes but might require the introduction of time stamps for faster synthesis procedures, such as hot injections commonly used in quantum dot synthesis.
• Bayesian optimisation can use prior knowledge, such as historical data and/or knowledge, to speed up the process.
In this context, a Gaussian process helps model the uncertainty about f(x), considering all previous samples to predict the function's behavior at unsampled points. The computation of the mean and variance for each point is automated, using the prior mean mx and the kernel function K(x,x′) which ideally represents the covariance between sampled points.
In the present case, the unknown f(x) is modelled as a Gaussian process. Given previous samples (x1),…, (xn), one can infer f(x) for any via the Gaussian process posterior. As such, at every point f(x) a normal distribution is obtained with a mean and variance predicted. Briefly, given a collection of jointly Gaussian random variables, there is an elegant way to compute the probability of an unknown variables subset from an observable subset through a Gaussian posterior expressed as:
yA = f(x) |
yB = [(x1),…, (xn)] |
Prior mean
Kernel function
K(x,x′). Ideally, K(x,x′) = Cov(f(x),f(x′)) |
Lastly, the integration of Bayesian optimization allows for the inclusion of diverse data sources including historical, published, and theoretical data, enabling robust exploration of the parameter space to quickly assess the proximity of the optimum to the available chemical and temperature conditions. This process is facilitated by the parallel execution of deployment and function inference steps, enhancing the efficiency of the optimization process.
Maggi's tutorial provides a detailed explanation of how this is computed.31 The code was written in Python, which is also used to integrate the hardware and analytics.
The seeds were prepared manually because the process is relatively simple and rapid. Moreover, using NaBH4 (unstable reactant) and the light colour of the seed solution make optimisation using the approach described herein challenging. The protocol for Ag seed production is described in the Methods section. The quality of the seeds was evaluted by dynamic light scattering (DLS) prior to use to confrim they have the correct size.
The reactor layout used for synthesising the Ag nanoplates is shown in Fig. 2a. Pump channel 3 contains the seeds solution diluted 1:10, mixed with the 1 mM ascorbic acid (pump channel 4) via a T-junction connector. The chemicals are combined on a 30 cm 1/16′′ PTFE tube with an ID of 500 μm. After the mixing part, 0.5 mM of AgNO3 is added through pump channel 2. The components are reacted on a 200 cm 1/16′′ PTFE tube with an ID of 500 μm. The formed particles are capped with citrate, added after the reaction via pump channel 1 (trisodium citrate 25 mM). The nanoparticle capping occurs in a 40 cm 1/16′′ PTFE tube with an ID of 500 μm.
The quality of the particles is analysed via online UV-Vis spectrometry through a Z-cell unit and compared to a pre-loaded spectrum of the desired target particles. Several protocols were established and scored using Bayesian optimisation with Gaussian processes. Each protocol has a unique set of pump flows. The highest-scoring protocol (i.e. the best fit) is summarised in Table 1.
Pump channel | Chemicals [concentration] | Flow rate (μL min−1) |
---|---|---|
1 | Trisodium citrate [25 mM] | 49.0 |
2 | AgNO3 [0.5 mM] | 132.0 |
3 | Seeds [1:10] | 101.0 |
4 | Ascorbic acid [1 mM] | 237.0 |
Fig. 2b shows the optical image of the particle prepared with the optimised protocol, consistent with a UV-Vis spectrum with an absorption maximum centred at 645 nm (ESI,† Fig. S1). The DLS measurement revealed two peaks (Fig. 2c), one corresponding to 2.0 nm and another to 43.8 nm. Transmission electron microscopy (TEM) shown in Fig. 2d reveals that the particles are Ag nanoplates with triangular shapes, as expected from Aherne et al.32 studies. The triangle side lengths are between 30–40 nm. The discrepancy with the size estimated from DLS is that DLS only measures the hydrodynamic radius and TEM for the metallic core. The second peak in the DLS is related to the nanoplate thickness, responsible for the UV-Vis absorption peak at around 401 nm (see Fig. S1, ESI†).
The absorption spectrum of the optimized Ag nanoplates protocol is consistent with the spectrum of the nanoparticles produced using the Aherne et al.32 batch protocol (see Fig. S1, ESI†). In both cases, the Ag+ ions are fully consumed based on X-ray fluorescence (XRF) analysis of the supernatant and the observation that unreacted Ag+ ions left in a synthesizing solution will be slowly added to the nanoplates increasing the particle size. However, UV-Vis absorption and DLS analysis after 15 days showed no change in optical absorption and size.
Having established that all the silver ions are completely consumed in both syntheses, the process's yield for metal consumption is 100%. However, as is noticeable in Fig. S1 (ESI†), the optical absorption intensity is different for both methods. The absorption intensity at the maximum in solutions after synthesis reflects the nanoparticle concentration in the solution. In the case of the batch, the peak has a net intensity of 0.74 versus 0.36 for the ones produced in the microfluidic reactor. The batch process generates about 1 mL min−1 of Ag NPs, but it is challenging to scale up production because insufficient stirring creates inhomogeneity in the solution, increasing particle polydispersity. The microfluidic reactor generates nanoparticles at a rate of ca. 0.5 mL min−1 but continuously and reproducibly. This is the proposed approach's conceptual advantage; it enables the production of nanomaterials with the same quality and morphology as the optimised protocol in volume and anywhere.
To demonstrate this, the reactor was configured with the same geometry, and the pump flows were set to the optimised values. The system was allowed to produce particles with these settings for 1 h, and the output was monitored via online UV-Vis. As shown in Fig. 2e, the reactor produced particles with the same optical absorption without noticeable reactor degradation due to fouling or blockage. This is even more evident from the Fig. 2e inset showing the position of the absorption maximum over time, which is very sensitive to changes in particle morphology. The quality of the particles is similar to the ones offered by commercial vendors,33 with the added benefit that the settings can be used in any lab or industrial plant to synthesise this optimised protocol. This is the competitive edge of the proposed approach.
The quality of the particles is analysed via online UV-Vis spectrometry through a Z-cell unit and compared to a pre-loaded spectrum of the desired target particles. Several protocols were established and scored using Bayesian optimisation with Gaussian processes. Each protocol has a unique set of pump flows. The best-fit protocol is summarised in Table 2.
Pump channel | Chemicals [concentration] | Flow rate (μL min−1) |
---|---|---|
1 | Citric acid [25.5 mM] | 220.0 |
2 | FeCl3 [1 mM] | 134.0 |
3 | K4[Fe(CN)6]x 3H2O [1 mM] + citric acid [25.5 mM] | 137.0 |
Fig. S2 (ESI†) shows the UV-Vis spectrum of the optimised nanoparticles with an absorption maximum centred at 697 nm, consistent with the formation of Prussian blue nanomaterial.35 The TEM of the prepared particles is shown in Fig. 4a, revealing small particles that are relatively uniform in size and shape. The histogram in Fig. 4b shows that particles are about 7.2 ± 1.3 nm, a uniformity in size confirmed by DLS analysis (Fig. S3, ESI†). The discrepancy between the average size estimated from DLS and TEM data is related to what is measured with each technique, as highlighted when presenting the Ag nanoplates data. The Prussian blue nanoparticle structure was determined by powder X-ray diffraction (XRD) shown in Fig. S4 (ESI†), which revealed the characteristic cubic framework built from Fe(II)–C–N–Fe(III) sequences, consistent with what has been published.36 As with the Ag nanoplates protocol, once optimized, the protocol was used to synthesize large batches of samples to be tested in applications like batteries and oxygen evolution reactions.
Fig. 4 TEM of the Prussian blue nanoparticles produced with the optimised synthesis protocol. (a) Representative TEM image; (b) histogram showing the distribution of particle sizes. |
The concentration of Prussian blue synthesized can be estimated from the optical absorption. Considering that the produced Prussian blue solution had an optical absorption at 700 nm of ca. 0.7 and the molar absorption coefficient at this wavelength is 3.0 × 104 M−1 cm−1,37 this equates to a concentration of 0.023 mM of Prussian blue. The reactor produces particles at a rate of about 0.5 mL min−1, equating to 8.4 × 10−8 mol of Fe per min in the product considering the formula C18Fe7N18. Considering the amount of iron used per minute, one estimates a yield concerning the metal of about 30%.
The current iteration of the proposed reactor permits reactions up to 3 bar pressure and 150 °C (limited by the heating system). For example, Au nanospheres with a diameter of ca. 8 nm at 85 °C (see ESI,† for reactor scheme Fig. S5 and the characteristic UV-Vis (Fig. S6) and DLS (Fig. S6) data). However, the approach can be adapted to high-pressure and high-temperature synthesis by replacing the tubing material with stainless steel and implementing more advanced heating and pressure regulator systems, which are commercially available.38 In this respect, the limitations are primarily related to available budgets and user willingness to handle more complex and sophisticated hardware. During the discovery and optimization stage, materials with the wrong properties are generated, which can induce reactor fouling and clogging. Thus, when possible, it is advisable to use cheaper and less prone to-fouling reactor tubing materials, like PTFE and PFA.
The unsupervised Bayesian optimization with Gaussian processes reduces the optimization time and the need for prior knowledge to initiate the process. However, the optimization algorithm is limited to twenty optimization parameters, and the math becomes computationally expensive once we go over ten parameters. However, considering the published wet chemical synthesis protocols, this is not a significant impediment because the protocols often use less than ten parameters.
A commercial NaaS concept can take several dimensions, which are difficult to cover in their entirety. However, we would like to mention three examples of how this could be beneficial.
Using the proposed approach to changing the current industrial manufacturing nanomaterials landscape is challenging, even when presented with cheaper, more adaptable, and more sustainable methods. This is because existing producers have established production lines, often involving significant capex expenditure and hired staff with specialization in such methods. Therefore, the use of a new approach will find more substantial penetration in unscaled, underdeveloped, and to-be-discovered processes. However, if one can identify a material producer willing to disrupt, the proposed approach could circumvent potential limitations with already commercialized nanoparticles. For example, establishing alternative flow protocols to batch process enables modularization. Modular concepts that can be produced in series allow for cost-efficient production through economies of scale. As their systems and components can be factory-assembled, they can be transported as modules or whole units to a location, reducing installation costs. However, one needs to examine processes individually and account for the end-user cost limitations and the pain associated with modifying the current production method, which is out of this work scope.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4nh00174e |
This journal is © The Royal Society of Chemistry 2024 |