Sjoerd de Haan, DataDrivenLab · July 2025 · ML for materials and measurement

Nano Particle Synthesis with Machine Learning

Can small machine learning experiments coupled with UV/Vis give insight into the synthesis of nanoparticles?

This paper shows that it can!

Synthesis of Gold nanoparticles for Machine Learning

The authors of “Machine Learning Analysis of Reaction Parameters in UV-Mediated Synthesis of Gold Nanoparticles” @gudaMachineLearningAnalysis2023 employed a UV-mediated synthesis approach where gold nanoparticles were produced by mixing a gold precursor (HAuCl₄·4H₂O), a surfactant (hexadecyltrimethylammonium bromide or CTAB), a reducing agent (L-ascorbic acid), AgNO₃, acetone, and cyclohexane in a quartz vessel. The mixture was irradiated with UV light from a mercury-xenon lamp (386 W) for 60 minutes at a distance of 25 cm.

In this study, the authors systematically varied 3 of the syntisis parameters; the concentrations of

With a improved Latin hypercube sampling, they explored 30 different reaction conditions.

Meanwhile they monitored the reaction progress in real-time using in situ UV-vis spectroscopy. This approach leverages the advantages of UV-Vis over TEM for real-time monitoring.

Reaction stages

The reaction process occurs in the following stages:

1. Initial Complex Formation

2. First Stage - Reduction and Decoloration

3. Second Stage - Nanoparticle Formation

Measurements

The UV/Vis spectra over time Figure 2 show the different stages.

Figure 3 shows the peak at 400 nm from the [AuBr₄]⁻ complex in blue. This peak disappears and a new peak grows at 525 nm, indicating the formation of gold nanoparticles. Figure 4 does not show the initial peak; the first stage goes very fast. It only shows the resultant peak at 525 nm, characteristic for spherical gold particles. Figure 5 shows the initial peak has not disappeared after 1 hour, indicating a slower reaction rate in this case.

Figure 2: UV/Vis spectra measured in situ during the synthesis of gold nanoparticles. The shown spectra are representative for the three observed behavoirs all spectra fall into.

3 groups

Figure 6 shows the UV/Vis spectra of all 30 samples. The spectra can be grouped into three categories based on their behavior:

The peaks at 400 nm are not visible in the theoretical model. The model used here @pashkovQuantitativeAnalysisUV2021 I discussed in a previous post. It deals with spherical and rod-shaped gold nano-particles; not with the formation process.

While some qualitative features are of the theoretical spectra are present, the experimental spectra are clearly different; representing the complex formation process of gold nanoparticles.

Machine learning

If the nanoparticles were formed during the synthesis, what was the shape of the nanoparticles and the time of their formation?

Machine learning was used to answer these questions, roughly:

  1. Did nano particles form during the 60 minutes?
  2. Were the nano particles spherical or elongated?
  3. Did it take more than 1500 s to form nano particles

The spectrum descriptors χ1\chi_1 and χ2\chi_2 were used in constructing labels for the machine learning model. These descriptors are defined in Table 1. So the spectrum forms the output space of the models.

The input of the models were the three synthesis parameters, linear combinations thereof and quadratic combinations; see Table 1.

The classifiers are thus mapping synthesis parameters to spectral properties. The mapping learned.

Table 1: Descriptors of spectra and descriptors of synthesis applied for training the Extra Trees Model
Descriptor comment
Spectra
χ1=A(525)/A(450)\chi_1 = A(525)/A(450) normalized height of the SPR maximum
χ2=A(625)/A(525)\chi_2 = A(625)/A(525) asymmetry of the SPR maximum
Synthesis
Au, AA, CTAB variable components in the synthesis
c1⋅Au+c2⋅AA+c3⋅CTABc_1 \cdot \text{Au} + c_2 \cdot \text{AA} + c_3 \cdot \text{CTAB}; ci=±1c_i = \pm 1 linear combinations of initial parameters
Au⋅\cdotAA, Au⋅\cdotCTAB, AA⋅\cdotCTAB products of the initial parameters
Table 2: Target classes for the machine learning model. The descriptors χ1\chi_1 and χ2\chi_2 are used to classify the spectra into three classes for each target.
# Target Descriptor Class 0 Class 1 Class 2
1 Particle Formation χ1=A(525)A(450)\chi_1 = \frac{A(525)}{A(450)} χ1<1\chi_1 < 1
(none)
1<χ1<1.51 < \chi_1 < 1.5
(intermediate)
χ1>1.5\chi_1 > 1.5
(significant)
2 Particle Shape χ2=A(625)A(525)\chi_2 = \frac{A(625)}{A(525)} χ1<1\chi_1 < 1
(no particles)
χ2<0.7\chi_2 < 0.7
(spherical)
χ2>0.7\chi_2 > 0.7
(elongated)
3 Formation Time τ\tau when χ1>1\chi_1 > 1 τ>4000\tau > 4000 s
(late/none)
1500<τ<40001500 < \tau < 4000 s
(intermediate)
τ<1500\tau < 1500 s
(early)

Feature selection

The extra tree model was used to perform feature selection in combination with leave-one-out cross validation.

Extra trees (Extremely randomized trees) is a tree ensemble method that is more robust than random forest and can be used with fewer data points.

Table 3: Selected Triples of features for synthesis prediction
Target Variable Selected Triple of Synthesis Descriptors Classification Accuracy
χ₁
(particle formation)
CTAB, Au, acid 0.60
CTAB, AA − CTAB·Au − AA 0.69
AA − CTAB, Au·AA·Au − AA − CTAB 0.79
χ₂
(particle shape)
CTAB, Au, Acid 0.57
AA + CTAB, Au, Au·AA 0.75
Au + AA + CTAB, Au − AA, Au − AA − CTAB 0.79
τ
(formation time)
CTAB, Au, acid 0.70
CTAB, Au − AA + CTAB, Au − AA 0.82
Au·AA, Au − AA + CTAB, Au − AA − CTAB 0.83

For χ1\chi_1 and τ\tau classifiers, the feature space could be reduced further to 2 descriptors. In Figure 7, the decision boundary of the 3 classes for χ1\chi_1 are visible.

Leave-One-Out Cross-Validation (LOOCV)

The leave one-out evaluation used here is prone to over-fitting.

Risk of Overfitting

High throughput experiment

The authors experiment with a miniaturized and automated version of the experiment that completed 60 experiments in just 2 hours.

This may be a good way to get started quickly with experimenting and developing machine learning along side. The authors note they there some issues with the miniature version.

For some combinations of flows, we have observed a plasmon peak shift toward lower wavelength, which was not observed in batch experiments. We attribute this shift to the formation of the AuAg alloy

As expected, some new effects are introduced. It is questionable how to extrapolate the results to a large scale production setup.

Lessons

With only 30 data points, the authors managed to distill some knowledge on gold nano particle synthesis. Some crucial parts are

  1. Limit to varying 3 production parameters
  2. Observe 3 groups of synthesis behaviour in the spectra
  3. Characterize the groups in terms of simple spectral quantities
  4. Simplify the spectrum prediction problem to a classification problem
  5. Use of robust learner like extra trees
  6. Feature selection

The paper positions this as a “proof of concept” for using ML with small synthesis datasets

The results should be viewed as demonstrating the feasibility of the ML approach rather than providing a robust, generalizable model. For production use, many more synthesis experiments would be needed.

It is also not clear how to extrapolate results to a production setting, to a different lab or different operator.

Working on a simulation or measurement pipeline that is too slow to iterate on? Book Free Consultation.