Modeling phase separation in polymer-derived carbonitride ceramics through extended machine learning molecular dynamics
Pith reviewed 2026-05-21 07:24 UTC · model grok-4.3
The pith
Machine learning molecular dynamics reveals carbon nucleating into graphene-like sheets within the SiCN ceramic matrix during heating.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central discovery is a phase separation process in which carbon domains nucleate from the amorphous SiCN matrix during thermal treatment, forming graphene-like sheets while the integrity of the ceramic network is preserved. Defective five- and seven-membered carbon rings mediate the conversion to stable six-membered aromatic structures. The simulated atomic configurations reproduce experimental pair distribution functions with high fidelity, supplying atomic-scale explanations for the combination of ceramic thermal stability and graphitic features in these materials.
What carries the argument
A machine learning interatomic potential trained on a diversified database of over 9000 configurations that enables large-scale molecular dynamics simulations of 8000-atom SiCNH systems at high temperature.
If this is right
- The method supplies atomic-level mechanisms for the thermal stability and mixed properties of polymer-derived ceramics.
- Defective carbon rings are identified as the structural intermediates that enable the transition to ordered graphitic domains.
- The same potential framework can be applied to other complex amorphous multicomponent systems at experimentally relevant sizes.
- Microscopic pathways for structural transformation during processing are now accessible for direct comparison with experiment.
Where Pith is reading between the lines
- Processing temperatures or times could be adjusted to control the size and connectivity of the emerging carbon sheets for targeted electrical or mechanical performance.
- The same simulation approach might be extended to predict how dopants or different polymer precursors alter the phase-separation route.
- Longer simulation timescales could reveal whether the graphene-like domains continue to grow or stabilize once formed.
Load-bearing premise
The trained machine learning potential reproduces the energies and dynamics of carbon nucleation and ring changes in large systems without large extrapolation errors or missing interactions.
What would settle it
A simulation run at the same temperature and scale that produces no carbon domain nucleation or that yields atomic pair distribution functions clearly different from experimental measurements would falsify the reported phase-separation mechanism.
Figures
read the original abstract
Polymer-derived ceramics combine the thermal stability of ceramics with the versatile properties of carbon domains, but modeling their atomic-scale evolution during processing remains elusive due to the limitations of traditional computational methods. To address this issue, here we develop and apply a machine learning interatomic potential for silicon carbonitride-based (Si-C-N-H) systems, trained on a diversified database of over 9000 configurations -including amorphous models, high-temperature states, surfaces, and crystal structure predictions - to capture the full complexity of these materials. This potential enables large-scale molecular dynamics simulations of 8000-atom systems revealing the atomic-scale evolution of the polymer-derived ceramic during thermal treatment. A key result of this work is the occurrence of a phase separation where carbon domains progressively nucleate from the amorphous SiCN matrix during thermal processing, forming distinct graphene-like sheets while preserving the integrity of the ceramic network. The resulting models reproduce the experimental atomic pair distribution functions with exceptional fidelity, validating our approach and providing microscopic explanations for the material unique combination of ceramic and graphitic properties. In this process, defective 5- and/or 7-member carbon rings, mediate the transformation to stable 6-member aromatic structures. These findings offer new atomic-scale insights into the thermal stability and structural transformation pathways of polymer-derived ceramics, while our methodology opens avenues for studying complex amorphous systems with first-principles accuracy at experimentally relevant scales.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a machine learning interatomic potential for Si-C-N-H systems trained on a diversified database of over 9000 configurations (including amorphous, high-temperature, surface, and crystal structures). It applies this potential to large-scale MD simulations of 8000-atom polymer-derived SiCN models during thermal processing. The central claim is that carbon domains nucleate from the amorphous matrix to form graphene-like sheets via mediation by defective 5- and 7-membered rings that rearrange into stable 6-membered aromatic structures, while the final atomic configurations reproduce experimental pair distribution functions with high fidelity.
Significance. If the dynamical pathway is faithfully reproduced, the work provides valuable atomic-scale mechanistic insight into phase separation and carbon domain formation in polymer-derived ceramics, explaining their hybrid ceramic-graphitic properties at experimentally relevant scales. The scale-up to 8000-atom systems and reproduction of experimental PDFs are notable strengths that could inform processing strategies.
major comments (1)
- [MD simulation results and validation section] The load-bearing claim is that the ML potential correctly drives carbon nucleation and the 5-/7- to 6-membered ring conversion during high-T annealing. However, validation is reported only via reproduction of experimental PDFs for the final structures (Abstract and results on MD trajectories). PDF agreement constrains equilibrium pair correlations but does not test relative energies or barriers of the defective-ring intermediates identified as mediating the transformation, leaving open the possibility that the observed pathway is an extrapolation artifact rather than a physically faithful mechanism.
minor comments (2)
- [Methods] The training database is described at a high level ('diversified database of over 9000 configurations'); more detail on the distribution of high-temperature states and how extrapolation risk was quantified for the 8000-atom runs would strengthen the methods.
- [Results] Error bars or uncertainty estimates on the PDF comparisons and on the reported ring statistics during annealing are not mentioned; adding these would improve clarity of the validation.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive review of our manuscript. We address the major comment regarding validation of the dynamical pathway below, and we outline revisions that will strengthen the presentation of our results.
read point-by-point responses
-
Referee: [MD simulation results and validation section] The load-bearing claim is that the ML potential correctly drives carbon nucleation and the 5-/7- to 6-membered ring conversion during high-T annealing. However, validation is reported only via reproduction of experimental PDFs for the final structures (Abstract and results on MD trajectories). PDF agreement constrains equilibrium pair correlations but does not test relative energies or barriers of the defective-ring intermediates identified as mediating the transformation, leaving open the possibility that the observed pathway is an extrapolation artifact rather than a physically faithful mechanism.
Authors: We agree that reproduction of experimental PDFs validates the final atomic configurations but does not directly probe the relative energies or barriers associated with the 5-/7- to 6-membered ring rearrangements. Our ML potential was trained on a diversified database of over 9000 configurations that explicitly includes high-temperature amorphous states, surfaces, and crystal structures, which encompass a broad range of carbon ring environments and defective motifs. This training strategy is intended to ensure faithful description of the relevant energetics. Nevertheless, to directly address the concern about possible extrapolation artifacts, the revised manuscript will include new analysis: we will report relative energies (computed with both the ML potential and reference DFT) for representative small-model systems containing 5-, 6-, and 7-membered carbon rings and their interconversions. These additional checks will confirm that the observed mediation pathway is consistent with the underlying potential energy surface learned from the training data. revision: yes
Circularity Check
No significant circularity; derivation relies on external training data and experimental validation
full rationale
The paper trains a machine learning interatomic potential on an external diversified database of over 9000 configurations (including amorphous models and high-temperature states) and then performs large-scale MD simulations on 8000-atom systems. The reported phase separation and ring rearrangements emerge as simulation outcomes, with final structures validated against independent experimental PDFs. No equations, fitted parameters, or self-citations reduce the central claim to the inputs by construction; the workflow remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- ML interatomic potential hyperparameters
axioms (1)
- domain assumption The diversified training database of over 9000 configurations sufficiently samples the relevant configuration space for amorphous, high-temperature, and surface states in Si-C-N-H systems.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
defective 5- and/or 7-member carbon rings mediate the transformation to stable 6-member aromatic structures
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
phase separation where carbon domains progressively nucleate from the amorphous SiCN matrix
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Configurations derived from first-principles molecular dynamics The SiCNH-periodic class contains configurations extracted from FPMD trajectories of periodic models with com- position Si32C25N24H19. This class includes configurations from seven models presented in our previous work [30] in addition to three new models produced with the purpose of increasi...
-
[2]
Initial configurations were obtained from the work of Deringer and Csányi
Amorphous Carbon and Silicon configurations The Amorphous-Carbon-200 and Amorphous-Carbon-250 classes incorporate amorphous carbon structures at den- sities of 2.00 and 2.50 g/cm3, respectively. Initial configurations were obtained from the work of Deringer and Csányi
-
[3]
where melt-quenched models were produced at different densities, resulting in models with carbon sp, sp2, and sp3 hybridizations. Each initial amorphous carbon configuration was subjected to a short thermal cycle of about 20 ps at T = 300 K, from which 249 configurations were extracted from each model. Similarly, for the Amorphous-Silicon class, we collec...
-
[4]
Crystal Structure Predictions The Crystal-Structure-Prediction class contains 1165 configurations generated using the USPEX program [43–45] for various Si-C-N-H compositions and densities. USPEX is a crystal structure prediction (CSP) program employing evolutionary algorithm to generate, generation after generation, a diverse set of candidate structures f...
work page 2000
-
[5]
Partial pair distribution functions The models being validated, we now focus on their structural characterization. Figure 5 presents the partial PDFs for the four 8000-atom models excluding hydrogen-containing pairs, which provide only very small contributions to the total PDF. Before examining the individual partial PDFs, we first point that most major p...
-
[6]
Local atomic environments Table VI shows the dominant Si local environments in the four 8000-atom models and Tables S4, S5, S6 and S7 in the supplementary materials, provide the complete list of local environments for all the models. We find that, in the four models, between 87% and 90% of silicon atoms are 4-fold coordinated, as expected from a Si-based ...
-
[7]
Evolution of C motifs Figure S7 in the supplementary materials shows snapshots of the carbon domains of all models at different stages of the thermals cycles. It reveals that, during the thermal annealing, carbon atoms undergo substantial organisation leading to the formation of dimers, chains and sheets. Specifically, we define five categories of carbon ...
-
[8]
Nucleation and growth mechanisms Figure 8a-f show local atomic configurations in the SiCNH8000-Random model at different times during the T = 1800 K plateau, highlighting the formation of an aromatic C ring. Starting from single carbon atoms and carbon dimers initially dispersed in the amorphous ceramic network (Figure 8a), we observe the progressive form...
-
[9]
P. Colombo, G. Mera, R. Riedel, and G. D. Sorarù, J. Am. Ceram. Soc.93, 1805 (2010)
work page 2010
-
[10]
R. Riedel, Z. Yu, S. Bernard, R. Kumar, and G. Singh, J. Am. Ceram. Soc.108, https://doi.org/10.1111/jace.70187 (2025)
- [11]
-
[12]
F. M. Ainger and J. M. Herbert,Special Ceramics, edited by P. Popper (Academic press, New York, US, 1960) pp. 168–182
work page 1960
-
[13]
P. G. Chantrell and P. Popper,Special Ceramics, edited by P. Popper (Academic press, New York, US, 1965) pp. 87–103
work page 1965
-
[14]
P. Colombo, G. Mera, R. Riedel, and G. D. Sorarù, J. Am. Ceram. Soc.93, 1805 (2010). 20
work page 2010
- [15]
-
[16]
R. P. Chaudhary, C. Parameswaran, M. Idrees, A. S. Rasaki, C. Liu, Z. Chen, and P. Colombo, Prog. Mater. Sci.128, 100969 (2022)
work page 2022
- [17]
-
[18]
Z. C. Eckel, C. Zhou, J. H. Martin, A. J. Jacobsen, W. B. Carter, and T. A. Schaedler, Science351, 58 (2016)
work page 2016
-
[19]
G. Barroso, Q. Li, R. K. Bordia, and G. Motz, J. Mater. Chem. A7, 1936 (2019)
work page 1936
-
[20]
Q. Wen, Z. Yu, and R. Riedel, Prog. Mater. Sci.109, 100623 (2020)
work page 2020
-
[21]
S. Trassl, H.-J. Kleebe, H. Störmer, G. Motz, E. Rössler, and G. Ziegler, J. Am. Ceram. Soc.85, 1268 (2002)
work page 2002
- [22]
-
[23]
Y. Chen, X. Yang, Y. Cao, Z. Gan, and L. An, Acta Mater.72, 22 (2014)
work page 2014
- [24]
-
[25]
O. D. Adigun, E. Ricohermoso, A. A. Daniyan, L. E. Umoru, and E. Ionescu, Ceramics5, 690 (2022)
work page 2022
- [26]
-
[27]
G. Mera, R. Riedel, F. Poli, and K. Müller, J. Eur. Ceram. Soc.29, 2873 (2009)
work page 2009
-
[28]
A. Klausmann, K. Morita, K. E. Johanns, C. Fasel, K. Durst, G. Mera, R. Riedel, and E. Ionescu, J. Eur. Ceram. Soc. 35, 3771 (2015)
work page 2015
-
[29]
J. Ding and W. Li, inAdvanced Ceramic Materials-Emerging Technologies(IntechOpen, 2024)
work page 2024
-
[30]
R. K. Morais Ferreira, M. Ben Miled, R. K. Nishihora, N. Christophe, P. Carles, G. Motz, A. Bouzid, R. Machado, O. Masson, Y. Iwamoto, S. Célérier, A. Habrioux, and S. Bernard, Nanoscale Adv.5, 701 (2023)
work page 2023
-
[31]
M. Ben Miled, M. Fradin, N. Benbakoura, L. Mazière, J. Rousseau, A. Bouzid, P. Carles, Y. Iwamoto, O. Masson, A. Habrioux, and S. Bernard, ChemSusChemn/a, e202400561 (2024)
work page 2024
-
[32]
S. Widgeon, G. Mera, Y. Gao, E. Stoyanov, S. Sen, A. Navrotsky, and R. Riedel, Chem. Mater.24, 1181 (2012)
work page 2012
-
[33]
J. Dürr, P. Lamparter, J. Bill, S. Steeb, and F. Aldinger, J. Non-Cryst. Solids232-234, 155 (1998)
work page 1998
-
[34]
K. Matsunaga, Y. Iwamoto, C. A. J. Fisher, and H. Matsubara, J. Ceram. Soc. Jpn.107, 1025 (1999)
work page 1999
-
[35]
N. Liao, W. Xue, and M. Zhang, J. Eur. Ceram. Soc.32, 1275 (2012)
work page 2012
- [36]
- [37]
-
[38]
F. Mortier, S. Cadars, M. B. Miled., O. Masson, G. Ori, M. Boero, Y. Wang, S. Bernard, and A. Bouzid, Phys. Chem. Chem. Phys. 10.1039/D5CP04954G (2026)
-
[39]
Cadars, Python atomic modeling analyzes (2023)
S. Cadars, Python atomic modeling analyzes (2023)
work page 2023
- [40]
- [41]
-
[42]
I. Batatia, D. P. Kovacs, G. Simm, C. Ortner, and G. Csányi, Advances in neural information processing systems35, 11423 (2022)
work page 2022
-
[43]
I. Batatia, S. Batzner, D. P. Kovács, A. Musaelian, G. N. C. Simm, R. Drautz, C. Ortner, B. Kozinsky, and G. Csányi, Nat. Mach. Intell.7, 56 (2025)
work page 2025
-
[44]
T. D. Kühne, M. Iannuzzi, M. Del Ben, V. V. Rybkin, P. Seewald, F. Stein, T. Laino, R. Z. Khaliullin, O. Schütt, F. Schiffmann, D. Golze, J. Wilhelm, S. Chulkov, M. H. Bani-Hashemian, V. Weber, U. Borštnik, M. Taillefumier, A. S. Jakobovits, A. Lazzaro, H. Pabst, T. Müller, R. Schade, M. Guidon, S. Andermatt, N. Holmberg, G. K. Schenter, A. Hehn, A. Bussy...
work page 2020
- [46]
-
[47]
J. VandeVondele, M. Krack, F. Mohamed, M. Parrinello, T. Chassaing, and J. Hutter, Comput. Phys. Commun.167, 103 (2005)
work page 2005
- [48]
-
[49]
V. L. Deringer and G. Cs’anyi, Phys. Rev. B95, 094203 (2017)
work page 2017
-
[50]
V. L. Deringer, N. Bernstein, A. P. Bartók, M. J. Cliffe, R. N. Kerber, L. E. Marbella, C. P. Grey, S. R. Elliott, and G. Csányi, J. Phys. Chem. Lett.9, 2879 (2018)
work page 2018
-
[51]
A. R. Oganov and C. W. Glass, J. Chem. Phys.124, 244704 (2006)
work page 2006
-
[52]
A. R. Oganov, A. O. Lyakhov, and M. Valle, Acc. Chem. Res.44, 227 (2011)
work page 2011
-
[53]
A. O. Lyakhov, A. R. Oganov, H. T. Stokes, and Q. Zhu, Comput. Phys. Commun.184, 1172 (2013)
work page 2013
- [54]
- [55]
-
[56]
J. B. Kruskal, Psychometrika29, 115 (1964)
work page 1964
-
[57]
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg,et al., J. Mach. Learn. Res.12, 2825 (2011)
work page 2011
- [58]
-
[59]
D. P. Kov’acs, I. Batatia, E. S. Arany, and G. Cs’anyi, J. Chem. Phys.159, https://doi.org/10.1063/5.0155322 (2023)
-
[60]
R. Barrett, C. Ortner, and J. Westermayr, arXiv preprint arXiv:2502.12870 https://doi.org/10.48550/arXiv.2502.12870 (2025)
-
[61]
A. Roy, R. Devanathan, S. Allec, G. Nandipati, A. M. Casella, D. J. Senor, D. D. Johnson, G. Balasubramanian, and A. Soulami, Adv. Intell. Discov. , e202500196 (2026). 21
work page 2026
-
[62]
D. P. Kingma and J. Ba, Adam: A method for stochastic optimization (2017)
work page 2017
-
[63]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Comput. Phys. Commun. 271, 108171 (2022)
work page 2022
-
[64]
H. C. Edwards, C. R. Trott, and D. Sunderland, J. Parallel Distrib. Comput.74, 3202 (2014)
work page 2014
- [65]
-
[66]
S. Nosé, J. Chem. Phys.81, 511 (1984)
work page 1984
-
[67]
W. G. Hoover, Phys. Rev. A31, 1695 (1985)
work page 1985
-
[68]
J. M. Martínez and L. Martínez, J. Comput. Chem.24, 819 (2003)
work page 2003
-
[69]
L. Martínez, R. Andrade, E. G. Birgin, and J. M. Martínez, J. Comput. Chem.30, 2157 (2009)
work page 2009
-
[71]
S. Schempp, J. Dürr, P. Lamparter, J. Bill, and F. Aldinger, Zeitschrift für Naturforschung A53, 127 (1998)
work page 1998
-
[72]
N. Choong Kwet Yive, R. Corriu, D. Leclercq, P. Mutin, and A. Vioux, Chem. Mater.4, 1263 (1992)
work page 1992
- [73]
- [74]
- [75]
-
[76]
G. J. Martyna, M. L. Klein, and M. Tuckerman, J. Chem. Phys.97, 2635 (1992)
work page 1992
-
[77]
T. D. Kühne, M. Krack, F. R. Mohamed, and M. Parrinello, Phys. Rev. Lett.98, 066401 (2007)
work page 2007
-
[78]
T. D. Kühne, WIREs Computational Molecular Science4, 391 (2014)
work page 2014
-
[79]
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett.77, 3865 (1996)
work page 1996
-
[80]
S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys.132, 154104 (2010). 1 Supplementary Materials: Modeling phase separation in polymer-derived carbonitride ceramics through extended machine learning molecular dynamics Fabien Mortier, Sylvian Cadars, Olivier Masson, Mauro Boero, Guido Ori, Yun Wang, Samuel Bernard, and Assil Bouzid. S1. ADDITION...
work page 2010
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.