Open Access Article
Junxiao Li†
a,
Sizhe Lin†a,
Wei Liua,
Changlong Lib,
Jiacheng Songb,
Zhiwei Liub,
Yue Hu
b and
Shuchen Zhang
*a
aState Key Laboratory of Precision and Intelligent Chemistry, Department of Materials Science and Engineering, School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, Anhui 230026, China. E-mail: zhangshuchen@ustc.edu.cn
bKey Laboratory of Carbon Materials of Zhejiang Province, College of Chemistry and Materials Engineering, Wenzhou University, Wenzhou 325000, P. R. China
First published on 23rd March 2026
Precise control over the diameter and chirality of single-walled carbon nanotubes (SWNTs) is critical for scalable electronic applications, yet it remains a challenge for horizontally aligned arrays. Although low-temperature powder synthesis routinely yields sub-1 nm SWNTs with relatively high chirality selectivity, directly transferring these growth conditions to substrate-guided horizontal arrays has been hindered by poorly understood catalyst evolution and coupled growth factors (gas, temperature, pressure etc). Here, we employ a machine-learning-assisted framework that integrates systematic experiments with interpretable models to optimize the growth parameters, enabling sub-1 nm SWNT growth in horizontal arrays at low temperatures. Our ML model reveals that SWNT formation is governed by kinetic competition between amorphous-carbon deposition and graphitized-carbon incorporation, with catalytic cracking acting as the dominant carbon source under low-temperature conditions. Further, we successfully achieve horizontally aligned SWNT arrays enriched in (6, 5) nanotubes with a chiral selectivity of approximately 70% at 650 °C under an optimized environment predicted by the ML model. This is the first demonstration of sub-1 nm chiral-selective growth in horizontal arrays at such a low temperature. This work bridges the long-standing gap between powder and horizontal-array growth and establishes a rational, data-driven strategy for the deterministic synthesis of SWNTs.
The prevailing understanding attributes the resulting chirality distribution of SWNTs to nucleation thermodynamics, growth kinetics, and the size and crystal structure of catalyst nanoparticles.12–17 Several catalyst-engineering strategies have successfully demonstrated the selective enrichment of specific chiralities. For example, Yang et al. employed W6Co7 catalysts exposing the (0012) facet to produce (12, 6) tubes, whereas Zhang et al. introduced a symmetry-matching design principle that enabled the growth of (12, 6) and (8, 4) tubes using Mo2C and WC catalysts, respectively.18–20 Nonetheless, most selectively synthesized chiralities remain metallic (2n, n) tubes, reflecting the kinetic dominance intrinsic to high-temperature growth conditions, under which only a narrow subset of catalyst species retain structurally stable configurations, thereby suppressing thermodynamically governed assembly pathways due to the limited availability of carbon atoms.
In contrast, low-temperature growth of SWNTs provides a broader design window for chirality control because a wider range of catalyst compositions can maintain stable and well-defined structures.21 Consistent with this rationality, powder samples from chemical vapor deposition (CVD) have long yielded enriched chiralities such as (6, 5),22–25 (7, 5),26–28 (7, 6),29 and (9, 8),30 which are more inclined towards selection by interface thermodynamic between SWNT and catalysts. However, a persistent discrepancy arises when these powder-based strategies are directly applied to the growth of horizontally aligned SWNT arrays on single-crystal substrates. While powder SWNT samples readily produce sub-1 nm nanotubes, substrate-guided horizontal-array growth—essential for wafer-scale device integration—typically results in substantially larger diameters.31 This inconsistency constrains the accessible chiral structure of SWNTs in aligned arrays and obscures the fundamental mechanisms governing catalyst evolution and nanotube nucleation. Resolving this long-standing contradiction requires identifying the true rate-limiting and structure-determining factors in the SWNT growth—a challenge made more difficult by the strong correlations among synthesis variables, many of which cannot be independently tuned through experiments alone.
Recently, machine learning (ML) has presented a powerful route for extracting mechanistic insights from complex and multidimensional synthesis datasets.32,33 By identifying hidden interactions among growth parameters, ML has begun to reshape approaches for catalyst discovery and mechanism elucidation.34 Notably, Zhang et al. developed the CARCO platform, which carried out millions of simulated experiments over 43 days and identified a high-performance Ti–Pt catalyst for dense horizontal-array growth with an experimental prediction accuracy of 56.25%.35 These achievements highlight the transformative potential of data-driven strategies in revealing previously inaccessible governing principles, implying the possibility of addressing the aforementioned issue.
Herein, we establish an ML-assisted framework that integrates systematic experimental datasets with interpretable models to identify the parameters that enable sub-1 nm SWNT growth in horizontal arrays at low temperatures. Our ML model shows that, under low-temperature conditions, SWNT formation is governed by the competition between amorphous-carbon deposition and graphitized-carbon incorporation, with catalytic cracking of carbon sources serving as the dominant carbon supply pathway relative to thermal decomposition. Regulating the carbon-atom flux shifts carbon nanotube formation from kinetically dominated accumulation toward thermodynamically guided assembly, thereby favoring the nucleation of structurally stable SWNT species. Further, we achieved horizontally aligned SWNT arrays enriched in the chiral (6, 5) tubes with a selectivity of 70% at 650 °C for the first time using the best optimized growth factors on catalysts and environment predicted by the ML model. This mechanistic understanding, coupled with data-driven guidance, reconciles the longstanding disparity between powder and horizontal-array growth and establishes a rational pathway toward the controlled synthesis of structurally well-defined SWNT arrays.
:
0, 3
:
1, 2
:
1, 1
:
1 and 0
:
1, respectively. The resulting precursor solutions were subsequently diluted 100-fold with ethanol, yielding a final Co2+ concentration of 0.05 mmol L−1.
Single-crystal ST-cut quartz substrates (single side polished, miscut angle <0.5°, surface roughness <5 Å; purchased from Hefei Kejing Materials Technology Co., Ltd, China) were sequentially ultrasonically cleaned in deionized water, acetone, ethanol and deionized water, each for 15 min. After cleaning, the quartz substrates were annealed in air at 900 °C for 8 h.
For catalyst deposition, 5 μL of the precursor solution was drop-cast onto each quartz substrate. The quartz substrates were then heated in a muffle furnace at 450 °C for 60 min with a ramping rate of 5 °C min−1 to ensure complete oxidization and dispersion of the catalyst precursors.
When ethanol was employed as the carbon source, the procedure was modified as follows. After catalyst reduction, the system was heated to the target temperature under an Ar/H2 atmosphere and held for extra 5 min. The flow rates of Ar and H2 were then each fixed at 200 sccm. Ethanol vapor was introduced by passing an independent Ar stream through an ethanol bubbler at a flow rate of 20, 40, 60, 100, or 150 sccm, and SWNT growth was carried out for 30 min. Following SWNT growth, the ethanol supply was terminated, and the system was cooled to 500 °C under an Ar/H2 atmosphere. The H2 flow was subsequently stopped, and the furnace was cooled to room temperature under a pure Ar flow.
Fitting of growth ability (CO as carbon source), run by Bayesian logistic regression model: Markov chains = 4, cores = 4, iteration = 8000, warmup = 6000, adapt_delta = 0.95, max_treedepth = 15 for running the model.
Fitting of growth ability (ethanol as carbon source), run by Bayesian logistic regression model: Markov chains = 4, cores = 4, iteration = 8000, warmup = 6000, adapt_delta = 0.95, max_treedepth = 15 for running the model.
Fitting of small-diameter selectivity, run by Bayesian nonlinear regression model: Markov chains = 4, cores = 4, iteration = 6000, warmup = 3000, adapt_delta = 0.98, max_treedepth = 15.
The validation process was performed with the BRMS package built-in diagnostic during the posterior analysis of the models. The BRMS package checks detailed posterior distribution information for model parameters, with validation of the stationary and convergent distributed Markov chain and sufficient number of posterior samples.
In addition, the specific content of chiral SWNTs was evaluated by UV-Vis-NIR spectroscopy. For this purpose, an aqueous sodium dodecyl sulfate (SDS, 1 wt%) solution was employed as the dispersion medium. The SWNTs were transferred from 80 quartz substrates into the SDS solution and subsequently dispersed by ultrasonication for 2 h. To ensure accurate background correction, a reference solution was prepared by subjecting an identical number of blank quartz substrates to the same treatment in an SDS solution of the same concentration under identical conditions. The resulting background signal was subtracted from the measured spectra, thereby isolating the intrinsic optical response of the SWNTs for subsequent analysis.
![]() | ||
| Fig. 1 Workflow for model development to guide the synthesis of a small-diameter SWNT horizontal array. | ||
Given their widespread use in powder synthesis and horizontally aligned SWNT growth, cobalt-based catalysts were selected as representative systems, including pure Co and Co–Mo alloy catalysts. Initial exploration experiments were conducted using Co catalysts with carbon monoxide (CO) as the carbon source. As shown in Fig. 2a–d, no SWNT growth was observed at 800 °C under a fixed CO concentration. Reducing the temperature to the range of 750–650 °C resulted in pronounced nanotube growth, whereas further cooling to 600 °C again strongly suppressed growth. These observations clearly indicate that low-temperature SWNT growth cannot be rationalized by temperature alone. Instead, the growth behavior emerges from the coupled effects of temperature, carbon-source partial pressure, and catalyst composition. Consequently, identifying growth windows that explicitly capture these interdependencies is essential for understanding and controlling the synthesis of SWNTs. However, direct experimental mapping of this multi-parameter growth landscape is prohibitively labor-intensive and experimentally inefficient.
To overcome this limitation, we introduce an ML-assisted framework that enables the rapid and quantitative construction of SWNT growth windows from limited experimental data. A structured dataset was generated using an orthogonal experimental design that systematically varied the growth temperature and carbon-source concentration while minimizing experimental redundancy. All samples were characterized by scanning electron microscopy (SEM) (Fig. S1–S4), and the average nanotube length was extracted as a quantitative descriptor of growth behavior (Fig. S5–S7). For each catalyst system, 16 sets of orthogonal experiments were performed, yielding a compact but information-rich database suitable for ML analysis.
Then, a Bayesian logistic regression model was applied to probabilistically classify the growth and non-growth regimes while minimizing uncertainty in the growth boundaries. This model was selected for its robustness in small-sample, high-dimensional parameter spaces, as well as its ability to provide interpretable posterior distributions rather than deterministic predictions. The trained model exhibited excellent convergence, with potential scale reduction factor (È) values below 1.05 for all variables, including catalyst composition, CO partial pressure, and temperature, and sufficiently large effective sample sizes (ESS) for all posterior distributions (Fig. S8a and b). These diagnostics confirm the statistical reliability of the ML-derived growth windows. Detailed descriptions of the ML methodology, including model construction and validation procedures, are provided in the Methods section.
Using Co and Co/Mo (2
:
1) catalysts as representative examples (Fig. 2e and f, respectively), the ML model reveals well-defined decision boundaries separating growth and non-growth regions. Both catalysts exhibit relatively broad growth windows when CO is used as the carbon source. Notably, SWNT growth is consistently suppressed at temperatures above ∼780 °C, whereas within the accessible growth regime, the dependence on temperature is comparatively weak. Comparative analysis across catalyst compositions, CO concentrations, and temperatures reveals a unified physical picture: catalytically active Co sites capable of dissociating CO are essential for enabling low-temperature SWNT growth, whereas excessively high temperatures universally suppress nanotube formation regardless of catalyst composition. This behavior reflects the thermodynamic and kinetic characteristics of CO cracking, in which elevated temperatures inhibit effective carbon incorporation into growing nanotubes. Correlation analysis of the full experimental-ML database further supports these conclusions (Fig. 2g and Fig. S8c). Together, these results demonstrate that low-temperature growth of horizontally aligned SWNT arrays can be efficiently achieved in a CO environment through the utilization of active Co sites.
However, the intrinsic role of temperature in SWNT growth cannot be fully resolved using CO as the sole carbon source, because CO cracking exhibits relatively weak temperature sensitivity within the growth window. To decouple temperature effects from carbon chemistry, additional carbon sources with distinct activation energies and decomposition pathways were introduced. Ethanol, a widely used carbon source for horizontally aligned SWNT arrays, was selected as a contrasting system. Accordingly, 16 sets of orthogonal experiments were conducted for both Co and Co/Mo (2
:
1) catalysts (Fig. S9 and S10, respectively) and analyzed using the same ML-assisted framework. As shown in Fig. 2h and i, the ethanol-based growth system exhibits growth windows that differ markedly from those observed for CO. In particular, the growth window narrows sharply with decreasing temperature, while increasing temperature significantly broadens the growth regime. This behavior reveals a strong and direct temperature dependence of SWNT growth in the ethanol system. Furthermore, pure Mo catalysts can serve as templates for SWNT growth at high temperatures when ethanol is used, but growth is entirely suppressed at low temperatures (Fig. 2j). These observations indicate that ethanol decomposition proceeds through two parallel pathways: thermal cracking and catalytic cracking. At high temperatures, ethanol can thermally decompose without catalytically active sites, generating sufficient carbon species for nanotube growth. In contrast, low-temperature growth relies predominantly on catalytic cracking, consistent with previous studies. Because ethanol has a relatively low activation energy for decomposition, catalytic cracking at low temperatures can generate excessive carbon species. Consequently, a high hydrogen concentration is required to suppress amorphous carbon accumulation and maintain the cleanliness of the catalyst surface. This requirement explains both the pronounced temperature sensitivity and the narrow growth windows observed in ethanol-based systems. Correlation analysis derived from the Bayesian ML model further substantiates these conclusions (Fig. S11). Collectively, these results establish ML-assisted growth-window construction as a powerful and generalizable strategy for decoding complex nanotube growth landscapes and for distinguishing the roles of carbon chemistry and temperature.
Based on the comparative analysis of the CO- and ethanol-based growth systems, we propose a microscopic mechanism governing the growth of horizontally aligned SWNT arrays at low temperatures (<800 °C). Three fundamental processes are involved: carbon atom generation, amorphous carbon accumulation, and SWNT growth. Carbon generation drives the latter two processes, which are intrinsically competitive. At low temperatures, amorphous carbon formation exhibits faster kinetics due to its lower activation barrier compared with SWNT nucleation. As illustrated in Fig. 3a, when the rate of carbon production (RC) via catalytic cracking exceeds the rate of amorphous carbon formation (RAC), rapid accumulation of amorphous carbon occurs on the catalyst surface, ultimately suppressing SWNT growth (Fig. 3b and c). Conversely, when RC is lower than RAC but higher than the SWNT growth rate (RGC), corresponding to a graphitization-dominated regime, SWNT growth is favored (Fig. 3d and e). Thus, low-temperature SWNT growth emerges from a delicate balance among these three competing processes. This framework explains why ethanol-based growth requires high hydrogen concentrations at low temperatures, leading to window confinement, whereas the CO-based system readily enables low-temperature growth with a comparatively broad window.
Having established growth windows for low-temperature SWNT synthesis, we next investigated structural control within these windows. Selective SWNT growth arises from the combined effects of thermodynamic nucleation and kinetic growth on catalyst surfaces. Given the narrow growth window of ethanol, which disfavors thermodynamically prioritized selectivity, subsequent investigations were confined to the CO-based system. Key factors influencing thermodynamic selectivity, including catalyst composition, CO partial pressure, and growth temperature, were systematically explored through structural characterization of as-grown SWNT arrays. SWNT growth was carried out under varying temperatures and CO partial pressures using catalysts of different compositions, and the resulting horizontal arrays were characterized by Raman spectroscopy. As shown in Fig. 4a–f and Fig. S12–S14, increasing the Mo content in the catalyst led to a pronounced narrowing of the radial breathing mode (RBM) peak distribution, whereas modest variations in CO partial pressure produced only minor changes within a given catalyst system. To quantify diameter selectivity, the RBM peak at 220 cm−1, corresponding to SWNTs with a diameter of ∼1 nm, was used as a threshold to distinguish small- (<1 nm) and large-diameter (>1 nm) nanotubes. Statistical analysis (Fig. 4g) reveals that for the pure Co catalyst, the fraction of small-diameter SWNTs decreases significantly with increasing temperature, likely due to catalyst coarsening induced by enhanced migration at elevated temperatures. In contrast, for the Co/Mo (2
:
1) catalyst under a CO/Ar flow ratio of 50
:
100, the relative abundance of small-diameter SWNTs increases from 50.0% to 84.6% as the temperature decreases. Across all the catalyst systems, the diameter selectivity exhibits weak dependence on the CO partial pressure, consistent with the intrinsically broad growth window of CO-based growth systems. These observations suggest that CO cracking proceeds relatively slowly, facilitating the interface-controlled, thermodynamically preferred growth of small-diameter SWNTs. The enhanced small-diameter selectivity observed for the Mo-rich catalysts at low temperatures, particularly at 650 °C, is attributed to the improved catalyst dispersibility and thermal stability of the catalyst structure arising from the high melting point of Mo.
To further quantify the relationships between growth parameters and structural selectivity, the relative abundances of different diameter populations were used as input variables for a Bayesian nonlinear ML model (Fig. 4h). The trained model also demonstrates good predictive accuracy, as confirmed by validation against experimental datasets (Fig. S15). Importantly, the model reveals a clear negative correlation between growth temperature and the abundance of small-diameter SWNTs, together with a positive correlation between small-diameter selectivity and molybdenum content in the catalyst. These statistical trends are fully consistent with the experimental observations and provide a quantitative framework for interpreting structure–growth relationships. Atomic force microscopy (AFM) analysis directly supports this interpretation, showing that catalysts with lower Co/Mo ratios exhibit more uniform nanoparticle dispersion and reduced average particle diameters (Fig. S16 and S17). These smaller and more stable nanoparticles provide energetically favorable nucleation sites for small-diameter SWNTs, thereby explaining the experimentally observed selective enrichment. Collectively, these results demonstrate that the catalyst composition governs the nanotube diameter not only through chemical activity but also by controlling the catalyst morphology and thermal stability during growth.
Guided by insights from the Bayesian nonlinear model, we propose a rational SWNT horizontal array strategy to further enhance the chiral selectivity at low temperatures. This strategy is based on three synergistic principles: minimizing the growth temperature to favor thermodynamic selectivity, maximizing the catalyst dispersion to suppress particle coarsening, and employing CO as a carbon source to provide a slow and controllable carbon flux. To experimentally validate this ML-guided strategy, the Mo content in the Co/Mo alloy catalyst was further decreased to achieve a Co/Mo ratio of 1
:
1. The corresponding AFM characterization reveals a pronounced enhancement in catalyst dispersity, accompanied by a further reduction in the average particle diameter to 1.05 nm at 650 °C (Fig. S18). This behavior is attributed to the high melting point of Mo, which suppresses catalyst migration and aggregation under thermal conditions, as supported by Fig. S19. Transmission electron microscopy (TEM) analysis further confirms the formation of a homogeneous alloy structure (Co2Mo3) (Fig. S20 and 21), indicating intimate mixing of Co and Mo at the nanoscale.36 These alloyed nanoparticles combine the catalytic activity of Co with the structural stabilization provided by Mo, yielding a catalyst architecture well suited for selective SWNT nucleation.
Using this optimized Co/Mo (1
:
1) catalyst, horizontally aligned SWNT arrays were grown on quartz substrates under CO at 650 °C, corresponding to the optimal growth window predicted by the Bayesian logistic regression model (Fig. S22). As shown in Fig. 5a, well-aligned SWNT horizontal arrays were obtained, with a linear density of approximately 1 tube per μm and a high degree of alignment along the substrate direction. Statistical analysis of the nanotube lengths indicates an average length of 9.87 μm, as displayed in Fig. 5b, demonstrating that enhanced selectivity is achieved without sacrificing growth efficiency. Raman spectroscopy measurements reveal pronounced enrichment of (6, 5) SWNTs, manifested as a dominant radial breathing mode (RBM) peak near 270 cm−1 (Fig. 5c). Systematic variation of the CO partial pressure shows that increasing the carbon flux further enhances the (6, 5) selectivity, reflecting the superior interfacial stability of this chirality with this Co/Mo (1
:
1) alloy catalyst (Fig. 5d). However, when the CO partial pressure approaches the upper boundary of the ML-predicted growth window, both alignment and areal density deteriorate, accompanied by the appearance of amorphous carbon. This behavior underscores the critical role of growth-window confinement in balancing the carbon supply, catalyst surface cleanliness, and selective nanotube nucleation. Ultimately, a maximum (6, 5) selectivity of 70% and an overall small-diameter selectivity of 91.2% were achieved (Fig. 5e), highlighting the effectiveness of the ML-guided growth strategy.
Additional spectroscopic characterization provides independent confirmation of the high chirality purity obtained under the optimized conditions. For example, the ultraviolet-visible-near-infrared absorption spectra (Fig. 5f and Fig. S23) of SWNTs dispersed from the substrate into solution exhibit a pronounced S11 transition at 967 nm, a typical characteristic of (6, 5) nanotubes, with minimal contributions from other chiral species. Although a weak feature that may be assigned to (10, 0) nanotubes is discernible in the ensemble absorption spectrum, its markedly lower intensity, together with the absence of statistically significant RBM signals in spatially localized Raman line-scan measurements, consistently indicates that these species, if present, exist only as a negligible minority. This result corroborates the Raman analysis and confirms that the observed enrichment originates from intrinsic growth selectivity rather than post-growth separation effects.
:
1) alloy catalyst was chosen at 650 °C to yield horizontally aligned SWNT arrays enriched in the thermodynamically dominant (6, 5) chirality with 70% selectivity for the first time. More broadly, the ML-assisted growth-window paradigm introduced here offers a general strategy for understanding and controlling complex material growth landscapes in low-dimensional nanomaterials.
Footnote |
| † These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2026 |